|
|---|
|
|
|
|
|
|
w*****1
发帖数: 473 | 1 我需要对每个position做一个male和female的t-test,这里只是部分数据。请问如何用R
的LOOP或者apply来实现这个任务呢？谢谢！
pos male male male male female female female female
100156763 84.61538462 90 88.88888889 66.66666667 81.25 81
.25 71.42857143 67.44186047
100156824 66.66666667 46.66666667 50 40 34.54545455 34.
54545455 50 48
100156849 80.95238095 100 95 80 87.27272727 87.27272727
83.33333333 84
100156857 71.42857143 100 95 100 98.18181818 98.18181818
100 96
100258672 100 90.90909091 100 100 90.32258065 90.32258065
100 90
100258694 100 90.90909091 100 100 93.5483871 93.5483871
80 80
100273779 70 28.57142857 83.33333333 100 60 60 0 65
.2173913
10033520 82.14285714 85.71428571 95.23809524 85.71428571 80
.95238095 80.95238095 100 90.90909091
100413619 95.23809524 85 94.11764706 100 81.81818182 81.
81818182 100 86.66666667 | l****s
发帖数: 6 | 2 高手们都不愿意回答你这样的入门问题。
每一行的pos是不是唯一的？如果是，可以就下面loop代码。如果不是，就稍改一下t.
test那行。看你数据结构，不象是配对t-test，但如果是，也稍改一下t.test那一行。
假设你读入的数据对象名为myData。
#++++++++++++++++++++++++++++++++++++++++
output<-vector("list",length(myData$pos))
n=0
for(ii in unique(myData$pos))
{
n=n+1
indx<-which(myData$pos==ii)
output[n]<-t.test(myData[ind,2:5],myData[ind,6:9])
} | a******1
发帖数: 1031 | 3
是不是用mapply替代循环结构更好？
【在 l****s 的大作中提到】
: 高手们都不愿意回答你这样的入门问题。
: 每一行的pos是不是唯一的？如果是，可以就下面loop代码。如果不是，就稍改一下t.
: test那行。看你数据结构，不象是配对t-test，但如果是，也稍改一下t.test那一行。
: 假设你读入的数据对象名为myData。
: #++++++++++++++++++++++++++++++++++++++++
: output<-vector("list",length(myData$pos))
: n=0
: for(ii in unique(myData$pos))
: {
: n=n+1
|
|
|
|
|
|
|
|
