R语言绘图

/picb/share/huangxin/PH525

下载filezilla，

Unix/Linux scp

man scp

打开FZ

主机：sftp://gate2.picb.ac.cn

输入picb用户和密码

端口：22

远程站点：/picb/share/huangxin/PH525

P-Value

Central Lmit Theorem(CLT)

mean(null>diff)

hist(null)

qqnorm(null)

qqline(null)

pops<-read.cssv("mice_pheno.csv")

hed(pops)

hf<- pops[popsSDiet=="hf"&popsSSex=="F",3]

chow<-pops[popsSDiet=="chow"&popsSSex=="F",3]

mean(hf)-mean(chow)

x<- sample(hf,12)

y<-sample(chow,12)

mean(x)_mean(y)

Ns<-c(3,5,10,25)

B<-10000

res<-sapply(Ns,funtion(n){sapply(1:8,function(j){mean(sample(hf,n))})})

lbrary(rafalib)

mypar2(2,2)

未完

~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~

sample variance

sample standard deviations

confidence intervals

t-statics:

开始编辑

dot<-read.csv("femaleMiceWeights.csv")

dot

contrlo <- dot[1:12,2]

treatment<-dot[12+1:12,2]

diff <-mean(treatment)-meant(control)

print(diff)

t.test(treatment,control)

sd(control)

sd(contrlo)/sqrt(length(control))

s <- sqrt(var(treatmen)/length(treatment)+var(control)/length(control))

tstat <-diff/se

1-pnorm(tstat)+pnorm(-tstat)

qqnorm(treatment)

qqline(treatment)

t.test(treatment,control)