/picb/share/huangxin/ph525
Download FileZilla,
Unix/linux SCP
Mans SCP
Open FZ
Host: sftp://gate2.picb.ac.cn
Enter PICB user and password
Port: 22
Remote site:/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)
Not finished
~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
Sample Variance
Sample Standard deviations
Confidence intervals
T-statics:
Start editing
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 (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)
R language Drawing