求助R语言如何运算permutation test

Python039

求助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语言

说起统计学非常厉害的学校就不得不说中国人民大学了,他们的统计学专业在中国所有的大学中是最好的。

特色是最优情况能拿两个博士,在中国的两年学分可以转,资格考tamu那边也认可,免GRE成绩。直博按导师组招生,经过一年级学习后再选导师。

优势是学院老师学术水平都很高,全都发过顶刊。课程设置难,要求高,经过几年调整教学质量有保障。19级开始分为数据科学和人工智能两个方向,有理论同时也有应用。本身人大统计招牌就不错,学位含金量不会低。

图为中国人民大学统计学院院长赵彦云做客人民网。

还有一个很好的选择是东北师范大学,东师的统计学是东师6个双一流学科之一(数学也是),2017年教学评估A,次于北大、人大。如果看就业率的话,数学与统计学院应该是常年100%的。

清华大学的统计学也是非常棒了,老师基本上都十分nice,备课也十分认真,用的基本都是各方向的经典教材,认真学的话基本都是干货满满。

学统计必然还是有一些压力的,毕竟是数学课,期中、期末和每周or每两周的作业基本都是逃不掉的,因此还是得根据自己的时间和精力量力而行。

不过要注意的是有些选修课可能两年甚至三年才开一次(比如随机过程就一直没有开出来),因此极有可能出现过了这村就没这店的情况,因此看到心仪的选修课,可能需要咬咬牙也要跟着上。

q1, dexp是指数分布的概率密度;

r*是生成*分布的随机数,p是*分布的概率分布,d是概率密度,还有q是d的反函数,貌似

R可向量操作,数乘向量是用数乘向量内的每个元素

q2, pi在R里是一个常量,可直接用,3.141593,没有你要求的精度高,你可pi<-3.1415926535...

q3, 不会

q4, integrate 函数能算数值积分;蒙特卡洛积分我不会,甚至都没听过,不是数学专业的,正常正常...

以 ∫ e^(-x) dx 上限: 无限, 下限:0为例啊;

f<-function(x) {exp(-x)}

integrate(f, lower=0, upper=Inf)

若是 ∫ e^(-x^3) dx的话,

f<-function(x) {exp(-x^3)}

就行了