pdf("g04.pdf");options(width=64)
#setwd("C:\\Users\\kolassa\\Class555")
setwd("~/Taught1/960-555/Data")
<pre>
#********************************************************/
#* Yarn strength data from Example Q of Cox & Snell     */
#* (1981).  Variables represent strength of two types of*/
#* yarn collected from six different bobbins.           */
#********************************************************/
yarn<-as.data.frame(scan("yarn.dat",what=
    list(strength=0, bobbin=0,type="")))
#<a name=block1> Block 1</a>
wilcox.test(strength~type,data=yarn[yarn$bobbin==3,],
   exact=FALSE, correct=FALSE)
#<a name=block2> Block 2</a>
wilcox.test(strength~type,data=yarn[yarn$bobbin==3,],
   exact=FALSE)
#<a name=block3> Block 3</a>
wilcox.test(strength~type,data=yarn[yarn$bobbin==3,],
   exact=TRUE)
#<a name=block4> Block 4</a>
arsenic<-as.data.frame(scan('arsenic.dat',
  what=list(age=0,sex=0,drink=0,cook=0,water=0,nails=0)))
wilcox.test(nails~sex,data=arsenic)
#*************************************************************/
# Data from http://lib.stat.cmu.edu/datasets/Arsenic         */
# reformatted into ASCII on the course home page.  Data re-  */
# flect arsenic levels in toenail clippings; covariates in-  */
# clude age, sex (1=M), categorical measures of quantities   */
# used for drinking and cooking, arsenic in the water, and   */
# arsenic in the nails.  To make arsenic.dat from Arsenic, do*/
#antiword Arsenic|awk '((NR>39)&&(NR<61)){print}'>arsenic.dat*/
#*************************************************************/
arsenic<-read.table("arsenic.dat")
names(arsenic)<-c("age","sex","drink","cook","water",
   "nails")
arsenic$status<-(arsenic$water>=.0001)+0
arsenic$water[arsenic$water<.0001]<-.0001
arsenic$lnail<-log(arsenic$nails)
arsenic$lwater<-log(arsenic$water)
#<a name=block5> Block 5</a>
library(exactRankTests)#For savage and vw scores
arsenic$savagenails<-as.numeric(cscores(
   arsenic$nails,type="Savage"))
arsenic$vwnails<-as.numeric(cscores(arsenic$nails,
   type="Normal"))
#<a name=block6> Block 6</a>
library(MultNonParam)#For genscorestat
genscorestat(arsenic$vwnails,arsenic$sex)
genscorestat(arsenic$savagenails,arsenic$sex)
#<a name=block7> Block 7</a>
library(MultNonParam)#For aov.P
aov.P(dattab=arsenic$nails,treatment=arsenic$sex)
#<a name=block8> Block 8</a>
library(NonparametricHeuristic)#For fun.comparepower
fun.comparepower(samp1sz=10,samp2sz=10,altvalue=0)
#<a name=block9> Block 9</a>
fun.comparepower(samp1sz=10,samp2sz=10,altvalue=1)
</pre>