首页 >> 知识 >> R语言 by()用法

R语言 by()用法

今天在学习GEO数据分析时,看到以下代码,看了几遍解读,还是不懂,究其原因是不会by()的用法,在网页上搜了一下,在github上看到一个比较详细的解释。果然这种还是要有例子

tmp = by(exp_id, ids$symbol, function(x) rownames(x)[which.max(rowMeans(x))])

使用的是R自带的数据warpbreaks,大家如果对这个by函数不理解的话,也最好看着下面的例子操作一下

示例1: 对warpbreaks这个数据集的第1,2列按“tension"这一列分成若干data.frame,然后进行summary

by(warpbreaks[,1:2],warpbreaks[,"tension"],summary)

结果如下: 示例2: 对warpbreaks这个数据集按“tension"这一列分成若干data.frame,然后分别计算每个data.frame中breaks与wool的相关。

by(warpbreaks,warpbreaks[,"tension"],function(x)lm(breaks~wool,data=x))

将data.frame拆分验证了一下,结果是对的

a=subset(warpbreaks,tension=="L")lm(a$breaks~a$wool)

验证时也又温习了一下subset()的用法,选字符型的时候要用“==”。

PS:always on the way 感谢各种学习帖!

网站地图