SourceCode_HW3

SourceCode_HW3 - aspect = 1.0) qqmath(~ li | session,...

Info iconThis preview shows pages 1–2. Sign up to view the full content.

View Full Document Right Arrow Icon
Sheet1 Page 1 library(lattice) load("C:\\lattice.RData") attach(keytime) kt14 <- keytime[keytime$subject==14,] kt14 <- keytime[keytime$subject==14,-1] dim(kt14) detach(keytime) attach(kt14) tapply(delay, list(session,run),sum) array<-tapply(delay, list(session,run),sum) dim(array) list(run, session) m<-matrix(tapply(delay,list(run,session), sum),400 , 21) m li<-as.vector(t(m)) li trellis.device(pdf, file="C:\\Users\\phuppo\\Desktop\\Plot.pdf") xyplot(li~run|session, data = kt14, prepanel=function(x,y) prepanel.loess(x,y, degree = 1, family = "s", span = 3, evaluation = 50), panel = function(x,y){ panel.xyplot(x,y) panel.loess(x,y, degree = 1, family = "s", span = 3, evaluation = 50) }, layout = c(4,2), scales = list(y = "free"), pch = ".", cex = .2, sub = list("XY PLOT WITH LOESS FITTING",cex=.8), xlab = "RUN", ylab = "log2 ofSUMMED DELAY",
Background image of page 1

Info iconThis preview has intentionally blurred sections. Sign up to view the full version.

View Full DocumentRight Arrow Icon
Background image of page 2
This is the end of the preview. Sign up to access the rest of the document.

Unformatted text preview: aspect = 1.0) qqmath(~ li | session, data=kt14, prepanel = prepanel.qqmathline, panel = function(x, . ..) { panel.grid() panel.qqmathline(x, distribution = qnorm) panel.qqmath(x, . ..) }, scales = "free", layout = c(2,2), pch=".", cex=2, Sheet1 Page 2 sub = list("Normal QQ Plot",cex=.8), xlab = "Unit Normal Quantile", ylab="log2 of SummedDelay") rDelay<-resid(loess(li~run) qqmath(~ rDelay | session, data=kt14, prepanel = prepanel.qqmathline, panel = function(x, . ..) { panel.grid() panel.qqmathline(x, distribution = qnorm) panel.qqmath(x, . ..) }, scales = "free", pch=".", cex=2, layout = c(2,2), sub = list("Normal QQ Plot of residuals of loess fit",cex=.8), xlab = "Unit Normal Quantile", ylab="residules of loess fit") graphics.off()...
View Full Document

This note was uploaded on 07/10/2011 for the course STAT 695 taught by Professor Staff during the Spring '08 term at Purdue.

Page1 / 2

SourceCode_HW3 - aspect = 1.0) qqmath(~ li | session,...

This preview shows document pages 1 - 2. Sign up to view the full document.

View Full Document Right Arrow Icon
Ask a homework question - tutors are online