求助R语言如何运算permutation test

Python011

求助R语言如何运算permutation test,第1张

原始数据

HAPLOTYPE CASEFREQUENCE CONTROLFREQUENCE

10.412027 0.421346

20.206045 0.233228

30.349419 0.296320

40.032509 0.049106

R语言运算Permutation test

dat=as.matrix(read.table(textConnection('X1 X2)

1 0.412027 0.421346

2 0.206045 0.233228

3 0.349419 0.296320

4 0.032509 0.049106

')))

d=rep(c(1,-1),each=1)

obs=drop(dat%*%d)

perms=replicate(99,dat[,sample(2)]%*%d)

>(nominal.p.value=rowMeans(abs(cbind(obs,perms))>abs(obs)))

本文来自: 人大经济 R语言

用R作图，比用EXCEL要灵活的多。 散点图，直接用plot()即可多类别，在R中就是多变量，用pionts() 加类别 拟合曲线用 fit<- lm()lines(fit)添加文字用 text()