pdf("g07.pdf");options(width=64)
#setwd("C:\\Users\\kolassa\\Class555")
setwd("~/Taught1/960-555/Data")
<pre>
#<a name=block1> Block 1</a>
twinbrain<-as.data.frame(scan("IQ_Brain_Size",
  what=list(CCMIDSA=0,FIQ=0,HC=0,ORDER=0,PAIR=0,SEX=0,
  TOTSA=0, TOTVOL=0,WEIGHT=0),skip=27,nmax=20))
fir<-twinbrain[twinbrain$ORDER==1,]
fir$v1<-fir$TOTVOL
sec<-twinbrain[twinbrain$ORDER==2,]
sec$v2<-sec$TOTVOL
brainpairs<-merge(fir,sec,by="PAIR")[,c("v1","v2")]
brainpairs$diff<-brainpairs$v2-brainpairs$v1
wilcox.test(brainpairs$diff)
#<a name=block2> Block 2</a>
library(BSDA)#For SIGN.test
SIGN.test(brainpairs$diff)
t.test(brainpairs$diff)
#<a name=block3> Block 3</a>
library(MultNonParam)# for symscorestat
cat("Asymptotic test using normal scores\n")
brainpairs$normalscores<-qqnorm(seq(length(
   brainpairs$diff)),plot.it=F)$x
symscorestat(brainpairs$diff,brainpairs$normalscores)
cat("Asymptotic test using savage scores\n")
brainpairs$savagescores<-cumsum(
   1/rev(seq(length(brainpairs$diff))))
symscorestat(brainpairs$diff,brainpairs$savagescores)
#<a name=block4> Block 4</a>
aves<-outer(brainpairs$diff,brainpairs$diff,"+")/2
sort(aves[upper.tri(aves,diag=TRUE)])
#<a name=block5> Block 5</a>
library(NonparametricHeuristic)#For hodgeslehmannexample
hodgeslehmannexample()
#<a name=block6> Block 6</a>
qsignrank(0.025, 10)
#<a name=block7> Block 7</a>
wilcox.test(brainpairs$diff,conf.int=TRUE)
</pre>