# # plot the COWS data from OSWALD # # ---------------------------------------------------------------------- # data <- matrix( scan("cows.data"), ncol=6, byrow=T ) # cows<-data.frame( id = data[,3], log.weight = data[,5], tx.group = data[,2], day = data[,4], week = data[,4]/7, year = data[,4]/365 ) # attach( cows ) # ### postscript(file="cows.plot1.ps",horiz=F) ### par( pty="s" ) plot( jitter(day), log.weight, pch="o" ) n<-length( table( id ) ) name.id <- as.integer( names( table( id ) ) ) for( j in 1:n ){ lines( day[ id==name.id[j] ], log.weight[ id==name.id[j] ], lty=j ) } title("Weight vs Day") # graphics.off() ### postscript(file="cows.plot2.ps",horiz=F) ### # # by TX group # par( pty="s", mfrow=c(2,2) ) # plot( day, log.weight, pch="o", type="n" ) group <- 1 for( j in 1:n ){ groupj <- mean( tx.group[ id==name.id[j] ] ) if( groupj == group ){ lines( day[ id==name.id[j] ], log.weight[ id==name.id[j] ], lty=j ) } } title( paste( "TX Group", group ) ) # plot( day, log.weight, pch="o", type="n" ) group <- 2 for( j in 1:n ){ groupj <- mean( tx.group[ id==name.id[j] ] ) if( groupj == group ){ lines( day[ id==name.id[j] ], log.weight[ id==name.id[j] ], lty=j ) } } title( paste( "TX Group", group ) ) # plot( day, log.weight, pch="o", type="n" ) group <- 3 for( j in 1:n ){ groupj <- mean( tx.group[ id==name.id[j] ] ) if( groupj == group ){ lines( day[ id==name.id[j] ], log.weight[ id==name.id[j] ], lty=j ) } } title( paste( "TX Group", group ) ) # plot( day, log.weight, pch="o", type="n" ) group <- 4 for( j in 1:n ){ groupj <- mean( tx.group[ id==name.id[j] ] ) if( groupj == group ){ lines( day[ id==name.id[j] ], log.weight[ id==name.id[j] ], lty=j ) } } title( paste( "TX Group", group ) ) # graphics.off() # # end-of-file...