#~New session: Version: S Wed Mar 29 17:35:08 PST 1995, Time: Wed Jan 17 09:33:55 1996
#~get "/usr/statlocal/lib/Splus3.3/splus/.Datasets/..min.script.id" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Datasets/.Copyright" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Datasets/version" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/.laenv" -1 ...
#~
reading <- read.table("reading", head = T)
#~error: Can't open file "reading" for reading
#~put ".Data.hp/last.dump" 821889274 ...
#~
reading <- read.table("reading.dat", head = T)
#~put ".Data.hp/reading" 821889282 ...
#~
reading
#~get ".Data.hp/reading" 821889282 ...
#~
attach(reading, pos = 1)
#~get ".Data.hp/reading" 821889282 ...
#~put "reading$.Last.value" 0 ...
#~
ls()
#~put "reading$.Last.value" 0 ...
#~
?box
#~put "reading$.Last.value" 0 ...
#~
motif()
#~get "/usr/statlocal/lib/Splus3.3/splus/.Datasets/.PostScript.Options" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Datasets/ps.preamble.std" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Datasets/ps.setcolor.rgb" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Datasets/ps.setfont.std" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Datasets/ps.paper.regions" -1 ...
#~put "reading$.Last.value" 0 ...
#~
?boxplot
#~put "reading$.Last.value" 0 ...
#~
plot(reading)
#~get ".Data.hp/reading" 821889282 ...
#~put "reading$.Last.value" 0 ...
#~
boxplot(Pre1, group = Group)
#~put "reading$.Last.value" 0 ...
#~
boxplot(list(Pre1, Group))
#~put "reading$.Last.value" 0 ...
#~
?boxplot
#~put "reading$.Last.value" 0 ...
#~
boxplot(Pre1, names = Group)
#~put "reading$last.dump" 0 ...
#~
boxplot(split(Pre1, Group))
#~put "reading$.Last.value" 0 ...
#~
pars(mfrow(2, 2))
#~error: couldn't find function "pars"
#~
par(mfrow(2, 2))
#~error: couldn't find function "mfrow"
#~put "reading$last.dump" 0 ...
#~
par(mfrow = c(2, 2))
#~put "reading$.Last.value" 0 ...
#~
plot(Pre1)
#~
plot(Split(Pre1, Group))
#~error: couldn't find function "Split"
#~put "reading$last.dump" 0 ...
#~
plot(split(Pre1, Group))
#~error: Cannot find x and y in list
#~put "reading$last.dump" 0 ...
#~
qqnorm(split(Pre1, Group))
#~error: ordering not defined for mode "list"
#~put "reading$last.dump" 0 ...
#~
split(Pre1, Group)
#~put "reading$.Last.value" 0 ...
#~
qqnorm(split(Pre1, Group)$Basal)
#~put "reading$.Last.value" 0 ...
#~
apropos(frame)
#~error: couldn't find function "apropos"
#~
clearframe()
#~error: couldn't find function "clearframe"
#~
newframe()
#~error: couldn't find function "newframe"
#~
clearplot
#~
clearplot
#~
clearplot()
#~put "reading$.Last.value" 0 ...
#~
qqnorm(split(Pre1, Group)$Basal)
#~put "reading$.Last.value" 0 ...
#~
qqnorm(split(Pre1, Group)$DRTA)
#~put "reading$.Last.value" 0 ...
#~
qqnorm(split(Pre1, Group)$Strat)
#~put "reading$.Last.value" 0 ...
#~
clearplot()
#~put "reading$.Last.value" 0 ...
#~
qqnorm(split(Pre1, Group)$Basal)
#~put "reading$.Last.value" 0 ...
#~
qqnorm(Pre1[Group == Basal])
#~error: Object "Basal" not found
#~put "reading$last.dump" 0 ...
#~
qqnorm(Pre1[Group == "Basal"])
#~put "reading$.Last.value" 0 ...
#~
clearplot()
#~put "reading$.Last.value" 0 ...
#~
qqnorm(split(Pre1, Group)$Basal)
#~put "reading$.Last.value" 0 ...
#~
qqnorm(Pre1[Group == "DRTA"])
#~put "reading$.Last.value" 0 ...
#~
qqnorm(Pre1[Group == "Strat"])
#~put "reading$.Last.value" 0 ...
#~
for(x in split(Pre1, Group)) {
	print(summary(x))
	print(var(x))
}
#~
reading.aov <- aov(Pre1 ~ Group)
#~get "/usr/statlocal/lib/Splus3.3/s/.Datasets/.Machine" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/stat/.Datasets/store.co.helmert" -1 ...
#~put "reading$reading.aov" 0 ...
#~
summary(reading.aov)
#~put "reading$.Last.value" 0 ...
#~
anova(reading.aov)
#~put "reading$.Last.value" 0 ...
#~
ls()
#~put "reading$.Last.value" 0 ...
#~
post.aov <- aov(Post1 ~ Group)
#~get "/usr/statlocal/lib/Splus3.3/s/.Datasets/.Machine" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/stat/.Datasets/store.co.helmert" -1 ...
#~put "reading$post.aov" 0 ...
#~
summary(post.aov)
#~put "reading$.Last.value" 0 ...
#~
anova(post.aov)
#~put "reading$.Last.value" 0 ...
#~
clearplot
#~
clearplot()
#~put "reading$.Last.value" 0 ...
#~
boxplot(split(residuals(post.aov), Group))
#~put "reading$.Last.value" 0 ...
#~
title("Pre Residuals")
#~
boxplot(split(Post1, Group))
#~put "reading$.Last.value" 0 ...
#~
title("Post Data")
#~
boxplot(split(residuals(post.aov), Group))
#~put "reading$.Last.value" 0 ...
#~
title("Post Residuals")
#~
clearplot()
#~put "reading$.Last.value" 0 ...
#~
boxplot(split(residuals(reading.aov), Group))
#~put "reading$.Last.value" 0 ...
#~
title("Pre Residuals")
#~
boxplot(split(Post1, Group))
#~put "reading$.Last.value" 0 ...
#~
title("Post Data")
#~
boxplot(split(residuals(post.aov), Group))
#~put "reading$.Last.value" 0 ...
#~
title("Post Residuals")
#~
split(reading)
#~error: Argument "group" is missing, with no default
#~put "reading$last.dump" 0 ...
#~
split(reading, Group)
#~get ".Data.hp/reading" 821889282 ...
#~put "reading$.Last.value" 0 ...
#~
reading
#~get ".Data.hp/reading" 821889282 ...
#~
Pre1[Group == "Basal"]
#~put "reading$.Last.value" 0 ...
#~
c(Pre1[Group == "Basal"], Pre1[Group == "DRTA"], Pre1[Group == "Strat"])
#~put "reading$.Last.value" 0 ...
#~
table(Pre1[Group == "Basal"], Pre1[Group == "DRTA"], Pre1[Group == "Strat"])
#~put "reading$.Last.value" 0 ...
#~
data.frame(Pre1[Group == "Basal"], Pre1[Group == "DRTA"], Pre1[Group == "Strat"
	])
#~put "reading$.Last.value" 0 ...
#~
data.frame(Post1[Group == "Basal"], Post1[Group == "DRTA"], Post1[Group == 
	"Strat"])
#~put "reading$.Last.value" 0 ...
#~
mymod <- lm(Post1 ~ Group)
#~get "/usr/statlocal/lib/Splus3.3/s/.Datasets/.Machine" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/stat/.Datasets/store.co.helmert" -1 ...
#~put "reading$mymod" 0 ...
#~
summary(mymod)
#~put "reading$.Last.value" 0 ...
#~
anova(mymod)
#~put "reading$.Last.value" 0 ...
#~
Group
#~
split(Group)
#~error: Argument "group" is missing, with no default
#~put "reading$last.dump" 0 ...
#~
x1 <- Group == Basal
#~error: Object "Basal" not found
#~
x1 <- Group == "Basal"
#~put "reading$x1" 0 ...
#~
x1
#~
x2 <- Group == "DRTA"
#~put "reading$x2" 0 ...
#~
x3 <- Group == "Strat"
#~put "reading$x3" 0 ...
#~
x1 <- as.integer(x1)
#~put "reading$x1" 0 ...
#~
x2 <- as.integer(x2)
#~put "reading$x2" 0 ...
#~
x3 <- as.integer(x3)
#~put "reading$x3" 0 ...
#~
data.frame(Group, Post1, x1, x2, x3)
#~put "reading$.Last.value" 0 ...
#~
mymod2 <- lm(Post1 ~ x1 + x2)
#~get "/usr/statlocal/lib/Splus3.3/s/.Datasets/.Machine" -1 ...
#~put "reading$mymod2" 0 ...
#~
summary(mymod2)
#~put "reading$.Last.value" 0 ...
#~
anova(mymod2)
#~put "reading$.Last.value" 0 ...
#~
ls()
#~put "reading$.Last.value" 0 ...
#~
q()
#~End session: Time: 821893510; Process: 12549
#~New session: Version: S Wed Mar 29 17:35:08 PST 1995, Time: Wed Jan 24 09:55:41 1996
#~get "/usr/statlocal/lib/Splus3.3/splus/.Datasets/..min.script.id" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Datasets/.Copyright" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Datasets/version" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/.laenv" -1 ...
#~
ls()
#~put ".Data.hp/.Last.value" 822495345 ...
#~
blood <- read.table("blood.dat", h = T)
#~put ".Data.hp/blood" 822495357 ...
#~
motif()
#~get "/usr/statlocal/lib/Splus3.3/splus/.Datasets/.PostScript.Options" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Datasets/ps.preamble.std" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Datasets/ps.setcolor.rgb" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Datasets/ps.setfont.std" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Datasets/ps.paper.regions" -1 ...
#~put ".Data.hp/.Last.value" 822495369 ...
#~
plot.factor(blood)
#~get ".Data.hp/blood" 822495357 ...
#~put ".Data.hp/.Last.value" 822495386 ...
#~
mymod <- aov(Time ~ C(Type, c(-1, 1, 1, -1), 1) + Type)
#~error: Object "Time" not found
#~put ".Data.hp/last.dump" 822495571 ...
#~
attach(blood, pos = 1)
#~get ".Data.hp/blood" 822495357 ...
#~put "blood$.Last.value" 0 ...
#~
mymod <- aov(Time ~ C(Type, c(-1, 1, 1, -1), 1) + Type)
#~get "/usr/statlocal/lib/Splus3.3/s/.Datasets/.Machine" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/stat/.Datasets/store.co.helmert" -1 ...
#~put "blood$mymod" 0 ...
#~
plot(mymod)
#~get "/usr/statlocal/lib/Splus3.3/s/.Datasets/gr.pars" -1 ...
#~put "blood$.Last.value" 0 ...
#~
par(mfrow = c(3, 3))
#~put "blood$.Last.value" 0 ...
#~
plot(mymod)
#~get "/usr/statlocal/lib/Splus3.3/s/.Datasets/gr.pars" -1 ...
#~put "blood$.Last.value" 0 ...
#~
par(mfrow = c(2, 3))
#~put "blood$.Last.value" 0 ...
#~
plot(mymod)
#~get "/usr/statlocal/lib/Splus3.3/s/.Datasets/gr.pars" -1 ...
#~put "blood$.Last.value" 0 ...
#~
plot.factor(data.frame(resid(mymod), Type))
#~put "blood$.Last.value" 0 ...
#~
par(mfrow = c(1, 1))
#~put "blood$.Last.value" 0 ...
#~
plot.factor(data.frame(resid(mymod), Type))
#~put "blood$.Last.value" 0 ...
#~
ls()
#~put "blood$.Last.value" 0 ...
#~
detach(1, save = F)
#~put ".Data.hp/.Last.value" 822496388 ...
#~
ls()
#~put ".Data.hp/.Last.value" 822496389 ...
#~
q()
#~End session: Time: 822496392; Process: 25986
#~interrupt
#~New session: Version: S Wed Mar 29 17:35:08 PST 1995, Time: Tue Jan 30 14:17:00 1996
#~get "/usr/statlocal/lib/Splus3.3/splus/.Datasets/..min.script.id" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Datasets/.Copyright" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Datasets/version" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/.laenv" -1 ...
#~
corn <- read.table("bhh-corn.dat", header = T)
#~put ".Data.hp/corn" 823029442 ...
#~
motif()
#~get "/usr/statlocal/lib/Splus3.3/splus/.Datasets/.PostScript.Options" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Datasets/ps.preamble.std" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Datasets/ps.setcolor.rgb" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Datasets/ps.setfont.std" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Datasets/ps.paper.regions" -1 ...
#~put ".Data.hp/.Last.value" 823029451 ...
#~
plot.factor(corn)
#~get ".Data.hp/corn" 823029442 ...
#~put ".Data.hp/.Last.value" 823029466 ...
#~
par(mfrow = c(2, 2))
#~put ".Data.hp/.Last.value" 823029476 ...
#~
plot.factor(corn)
#~get ".Data.hp/corn" 823029442 ...
#~put ".Data.hp/.Last.value" 823029483 ...
#~
?interaction.plot
#~put ".Data.hp/.Last.value" 823029577 ...
#~
attach(corn, pos = 1)
#~get ".Data.hp/corn" 823029442 ...
#~put "corn$.Last.value" 0 ...
#~
interaction.plot(Blend, Tmp, Yld)
#~error: Object "Tmp" not found
#~put "corn$last.dump" 0 ...
#~
interaction.plot(Blend, Tmt, Yld)
#~get "/usr/statlocal/lib/Splus3.3/s/.Datasets/letters" -1 ...
#~
interaction.plot(Tmt, Blend, Yld)
#~get "/usr/statlocal/lib/Splus3.3/s/.Datasets/letters" -1 ...
#~
blood2 <- read.table("blood2.dat", header = T)
#~put "corn$blood2" 0 ...
#~
search()
#~put "corn$.Last.value" 0 ...
#~
detach(1, save = F)
#~put ".Data.hp/.Last.value" 823030607 ...
#~
blood2
#~error: Object "blood2" not found
#~
blood2 <- read.table("blood2.dat", header = T)
#~put ".Data.hp/blood2" 823030613 ...
#~
attach(blood1, pos = 1)
#~error: Object "blood1" not found
#~put ".Data.hp/last.dump" 823030620 ...
#~
attach(blood2, pos = 1)
#~get ".Data.hp/blood2" 823030613 ...
#~put "blood2$.Last.value" 0 ...
#~
plot.factor(blood2)
#~get ".Data.hp/blood2" 823030613 ...
#~put "blood2$.Last.value" 0 ...
#~
plot.factor(blood2)
#~get ".Data.hp/blood2" 823030613 ...
#~put "blood2$.Last.value" 0 ...
#~
interaction.plot(Drug, Biofeedback, Resp)
#~get "/usr/statlocal/lib/Splus3.3/s/.Datasets/letters" -1 ...
#~
interaction.plot(Biofeedback, Drug, Resp)
#~get "/usr/statlocal/lib/Splus3.3/s/.Datasets/letters" -1 ...
#~
search()
#~put "blood2$.Last.value" 0 ...
#~
q()
#~End session: Time: 823033100; Process: 537
#~New session: Version: S Wed Mar 29 17:35:08 PST 1995, Time: Tue Jan 30 15:48:54 1996
#~get "/usr/statlocal/lib/Splus3.3/splus/.Datasets/..min.script.id" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Datasets/.Copyright" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Datasets/version" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/.laenv" -1 ...
#~
delme <- scan("coxdata.dat")
#~put ".Data.hp/delme" 823034951 ...
#~
delme
#~get ".Data.hp/delme" 823034951 ...
#~
#~interrupt
#~
rep(c("T1", "T2", "T3", "T4", "T5"), 3)
#~put ".Data.hp/.Last.value" 823034998 ...
#~
delmeT <- rep(c("T1", "T2", "T3", "T4", "T5"), 3)
#~put ".Data.hp/delmeT" 823035012 ...
#~
?rep
#~put ".Data.hp/.Last.value" 823035049 ...
#~
rep(c("B1", "B2", "B3"), c(5, 5, 5))
#~put ".Data.hp/.Last.value" 823035075 ...
#~
delmeB <- rep(c("B1", "B2", "B3"), c(5, 5, 5))
#~put ".Data.hp/delmeB" 823035082 ...
#~
data.frame(delme, delmeT, delmeB)
#~get ".Data.hp/delme" 823034951 ...
#~get ".Data.hp/delmeT" 823035012 ...
#~get ".Data.hp/delmeB" 823035082 ...
#~put ".Data.hp/.Last.value" 823035092 ...
#~
objects(p = "del*")
#~put ".Data.hp/.Last.value" 823035163 ...
#~
remove(objects(p = "del*"))
#~remove ".Data.hp/delme" 0 ...
#~remove ".Data.hp/delmeB" 0 ...
#~remove ".Data.hp/delmeT" 0 ...
#~
objects(p = "del*")
#~put ".Data.hp/.Last.value" 823035169 ...
#~
ls()
#~put ".Data.hp/.Last.value" 823035171 ...
#~
coxdat <- read.table("coxdata.dat")
#~put ".Data.hp/coxdat" 823035191 ...
#~
coxdat
#~get ".Data.hp/coxdat" 823035191 ...
#~
attach(coxdat, pos = 1)
#~get ".Data.hp/coxdat" 823035191 ...
#~put "coxdat$.Last.value" 0 ...
#~
mymod <- aov(Yield ~ Block + Tmt)
#~error: Object "Yield" not found
#~put "coxdat$last.dump" 0 ...
#~
ls()
#~put "coxdat$.Last.value" 0 ...
#~
detach(1, save = F)
#~put ".Data.hp/.Last.value" 823035317 ...
#~
coxdat <- read.table("coxdata.dat", h = T)
#~put ".Data.hp/coxdat" 823035324 ...
#~
attach(coxdat, pos = 1)
#~get ".Data.hp/coxdat" 823035324 ...
#~put "coxdat$.Last.value" 0 ...
#~
mymod <- aov(Yield ~ Block + Tmt)
#~get "/usr/statlocal/lib/Splus3.3/s/.Datasets/.Machine" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/stat/.Datasets/store.co.helmert" -1 ...
#~put "coxdat$mymod" 0 ...
#~
coxdat
#~get ".Data.hp/coxdat" 823035324 ...
#~
anova(mymod)
#~put "coxdat$.Last.value" 0 ...
#~
motif()
#~get "/usr/statlocal/lib/Splus3.3/splus/.Datasets/.PostScript.Options" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Datasets/ps.preamble.std" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Datasets/ps.setcolor.rgb" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Datasets/ps.setfont.std" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Datasets/ps.paper.regions" -1 ...
#~put "coxdat$.Last.value" 0 ...
#~
plot.factor(coxdat)
#~get ".Data.hp/coxdat" 823035324 ...
#~put "coxdat$.Last.value" 0 ...
#~
interaction.plot(coxdat)
#~error: Argument "trace.factor" is missing, with no default
#~put "coxdat$last.dump" 0 ...
#~
plot.factor(coxdat)
#~get ".Data.hp/coxdat" 823035324 ...
#~put "coxdat$.Last.value" 0 ...
#~
fitted(mymod)
#~put "coxdat$.Last.value" 0 ...
#~
residual(mymod)
#~error: couldn't find function "residual"
#~
residuals(mymod)
#~put "coxdat$.Last.value" 0 ...
#~
effects(mymod)
#~put "coxdat$.Last.value" 0 ...
#~
ls()
#~put "coxdat$.Last.value" 0 ...
#~
?mean
#~put "coxdat$.Last.value" 0 ...
#~
?by
#~put "coxdat$.Last.value" 0 ...
#~
by(Yeild, Block, mean)
#~error: Object "Yeild" not found
#~put "coxdat$last.dump" 0 ...
#~
by(Yield, Block, mean)
#~error: Cannot coerce mode list to double
#~put "coxdat$last.dump" 0 ...
#~
by(Yeild, Block, mean)
#~error: Object "Yeild" not found
#~put "coxdat$last.dump" 0 ...
#~
?by
#~put "coxdat$.Last.value" 0 ...
#~
ls()
#~put "coxdat$.Last.value" 0 ...
#~
by(Yield, Block, mean)
#~error: Cannot coerce mode list to double
#~put "coxdat$last.dump" 0 ...
#~
?by
#~put "coxdat$.Last.value" 0 ...
#~
by(coxdat, Block, mean)
#~error: Cannot coerce mode list to double
#~put "coxdat$last.dump" 0 ...
#~
by(data.frame(Yield), Block, summary)
#~put "coxdat$.Last.value" 0 ...
#~
mean(data.frame(Yield))
#~error: Cannot coerce mode list to double
#~put "coxdat$last.dump" 0 ...
#~
by(Yield, Tmt, summary)
#~put "coxdat$.Last.value" 0 ...
#~
delme <- by(Yield, Tmt, summary)
#~put "coxdat$delme" 0 ...
#~
delme
#~
attributes(delme)
#~put "coxdat$.Last.value" 0 ...
#~
delme$T1
#~put "coxdat$.Last.value" 0 ...
#~
attributes(delme$T1)
#~put "coxdat$.Last.value" 0 ...
#~
delme$T1[3]
#~error: No data to interpret as numerical value
#~
?tapply
#~put "coxdat$.Last.value" 0 ...
#~
tapply(Yield, Block, mean)
#~put "coxdat$.Last.value" 0 ...
#~
tapply(Yield, Tmt, mean)
#~put "coxdat$.Last.value" 0 ...
#~
mymeans <- tapply(Yield, Tmt, mean)
#~put "coxdat$mymeans" 0 ...
#~
mymeans
#~
?apply
#~put "coxdat$.Last.value" 0 ...
#~
?tapply
#~put "coxdat$.Last.value" 0 ...
#~
mymeans$T2 - mymeans$T5
#~put "coxdat$.Last.value" 0 ...
#~
mymeans
#~
mymeans[2] - mymeans[5]
#~put "coxdat$.Last.value" 0 ...
#~
diff(mymeans)
#~put "coxdat$.Last.value" 0 ...
#~
anova(mymod)
#~put "coxdat$.Last.value" 0 ...
#~
(0.043685 * 2)/3
#~put "coxdat$.Last.value" 0 ...
#~
sqrt((0.043685 * 2)/3)
#~put "coxdat$.Last.value" 0 ...
#~
4.17 * 0.17
#~put "coxdat$.Last.value" 0 ...
#~
4.17/sqrt(2)
#~put "coxdat$.Last.value" 0 ...
#~
(4.17 * 0.17)/sqrt(2)
#~put "coxdat$.Last.value" 0 ...
#~
qt(0.95, 8)
#~put "coxdat$.Last.value" 0 ...
#~
qt(0.95, 8) * 0.17
#~put "coxdat$.Last.value" 0 ...
#~
?qf
#~put "coxdat$.Last.value" 0 ...
#~
qf(0.9, 4, 8)
#~put "coxdat$.Last.value" 0 ...
#~
sqrt(4 * qf(0.9, 4, 8))
#~put "coxdat$.Last.value" 0 ...
#~
sqrt(4 * qf(0.9, 4, 8)) * 0.17
#~put "coxdat$.Last.value" 0 ...
#~
2.95 * 0.17
#~put "coxdat$.Last.value" 0 ...
#~
#~End session: Time: 823097609; Process: 821
#~New session: Version: S Wed Mar 29 17:35:08 PST 1995, Time: Mon Feb  5 09:10:42 1996
#~get "/usr/statlocal/lib/Splus3.3/splus/.Datasets/..min.script.id" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Datasets/.Copyright" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Datasets/version" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/.laenv" -1 ...
#~
caster <- read.table("castle.dat", h = T)
#~put ".Data.hp/caster" 823529491 ...
#~
rm(caster)
#~remove ".Data.hp/caster" 0 ...
#~
castle <- read.table("castle.dat", h = T)
#~put ".Data.hp/castle" 823529503 ...
#~
castle <- data.frame(Resp = castle$Resp, Height = as.factor(castle$Height), 
	Width = as.factor(castle$Width))
#~get ".Data.hp/castle" 823529503 ...
#~put ".Data.hp/castle" 823529546 ...
#~
castle
#~get ".Data.hp/castle" 823529546 ...
#~
attach(castle, pos = 1)
#~get ".Data.hp/castle" 823529546 ...
#~put "castle$.Last.value" 0 ...
#~
par(mfrow = c(2, 2))
#~error: Device not active
#~put "castle$last.dump" 0 ...
#~
motif()
#~get "/usr/statlocal/lib/Splus3.3/splus/.Datasets/.PostScript.Options" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Datasets/ps.preamble.std" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Datasets/ps.setcolor.rgb" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Datasets/ps.setfont.std" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Datasets/ps.paper.regions" -1 ...
#~put "castle$.Last.value" 0 ...
#~
par(mfrow = c(2, 2))
#~put "castle$.Last.value" 0 ...
#~
plot.factor(castle)
#~get ".Data.hp/castle" 823529546 ...
#~put "castle$.Last.value" 0 ...
#~
interaction.plot(Width, Weight, Resp)
#~error: Object "Weight" not found
#~put "castle$last.dump" 0 ...
#~
ls()
#~put "castle$.Last.value" 0 ...
#~
interaction.plot(Width, Width, Resp)
#~get "/usr/statlocal/lib/Splus3.3/s/.Datasets/letters" -1 ...
#~
interaction.plot(Width, Width, Resp)
#~get "/usr/statlocal/lib/Splus3.3/s/.Datasets/letters" -1 ...
#~
plot.factor(castle)
#~get ".Data.hp/castle" 823529546 ...
#~put "castle$.Last.value" 0 ...
#~
interaction.plot(Width, Height, Resp)
#~get "/usr/statlocal/lib/Splus3.3/s/.Datasets/letters" -1 ...
#~
interaction.plot(Height, Width, Resp)
#~get "/usr/statlocal/lib/Splus3.3/s/.Datasets/letters" -1 ...
#~
objects()
#~put "castle$.Last.value" 0 ...
#~
detach(1, save = F)
#~put ".Data.hp/.Last.value" 823531755 ...
#~
objects()
#~put ".Data.hp/.Last.value" 823531757 ...
#~
#~End session: Time: 823531760; Process: 1740
#~New session: Version: S Wed Mar 29 17:35:08 PST 1995, Time: Mon Feb  5 09:49:24 1996
#~get "/usr/statlocal/lib/Splus3.3/splus/.Datasets/..min.script.id" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Datasets/.Copyright" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Datasets/version" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/.laenv" -1 ...
#~
objects()
#~put ".Data.hp/.Last.value" 823531767 ...
#~
motif()
#~get "/usr/statlocal/lib/Splus3.3/splus/.Datasets/.PostScript.Options" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Datasets/ps.preamble.std" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Datasets/ps.setcolor.rgb" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Datasets/ps.setfont.std" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Datasets/ps.paper.regions" -1 ...
#~put ".Data.hp/.Last.value" 823531776 ...
#~
pars(mfrow = c(2, 2))
#~error: couldn't find function "pars"
#~
par(mfrow = c(2, 2))
#~put ".Data.hp/.Last.value" 823531810 ...
#~
attach(poisson.frame, pos = 1)
#~error: Object "poisson.frame" not found
#~put ".Data.hp/last.dump" 823531820 ...
#~
attach(poison.frame, pos = 1)
#~get ".Data.hp/poison.frame" 823531736 ...
#~put "poison.frame$.Last.value" 0 ...
#~
plot.factor(poison.frame)
#~get ".Data.hp/poison.frame" 823531736 ...
#~put "poison.frame$.Last.value" 0 ...
#~
interaction.plot(Poison, Tmt, poison)
#~get "/usr/statlocal/lib/Splus3.3/s/.Datasets/letters" -1 ...
#~
interaction.plot(Tmt, Poison, poison)
#~get "/usr/statlocal/lib/Splus3.3/s/.Datasets/letters" -1 ...
#~
anova(mymod <- aov(poison ~ Poison + Tmt + Poison * Tmt))
#~get "/usr/statlocal/lib/Splus3.3/s/.Datasets/.Machine" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/stat/.Datasets/store.co.helmert" -1 ...
#~put "poison.frame$mymod" 0 ...
#~put "poison.frame$.Last.value" 0 ...
#~
plot.factor(residuals(mymod))
#~error: can't interpret as a factor or data fram
#~put "poison.frame$last.dump" 0 ...
#~
plot(residuals(mymod))
#~
?boxplot
#~put "poison.frame$.Last.value" 0 ...
#~
boxplot(split(residuals(mymod), Poison))
#~put "poison.frame$.Last.value" 0 ...
#~
boxplot(split(residuals(mymod), Tmt))
#~put "poison.frame$.Last.value" 0 ...
#~
anova(mymod <- aov((poison)^(-1) ~ Poison + Tmt + Poison * Tmt))
#~get "/usr/statlocal/lib/Splus3.3/s/.Datasets/.Machine" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/stat/.Datasets/store.co.helmert" -1 ...
#~put "poison.frame$mymod" 0 ...
#~put "poison.frame$.Last.value" 0 ...
#~
newframe()
#~error: couldn't find function "newframe"
#~
clearplot
#~
clearplot()
#~put "poison.frame$.Last.value" 0 ...
#~
false -  - nostringval -  -  -  - nostringval -  - false
#~error: Object "false" not found
#~
plot.factor(data.frame(resid(mymod), Poison, Tmt))
#~put "poison.frame$.Last.value" 0 ...
#~
plot(fitted(mymod), resid(mymod))
#~
anova(mymod <- aov((poison)^(0.5) ~ Poison + Tmt + Poison * Tmt))
#~get "/usr/statlocal/lib/Splus3.3/s/.Datasets/.Machine" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/stat/.Datasets/store.co.helmert" -1 ...
#~put "poison.frame$mymod" 0 ...
#~put "poison.frame$.Last.value" 0 ...
#~
clearplot()
#~put "poison.frame$.Last.value" 0 ...
#~
plot.factor(data.frame(resid(mymod), Poison, Tmt))
#~put "poison.frame$.Last.value" 0 ...
#~
plot(fitted(mymod), resid(mymod))
#~
clearplot()
#~put "poison.frame$.Last.value" 0 ...
#~
plot.factor(poison.frame)
#~error: couldn't open file ".Data.hp/poison.frame" to read S object "poison.frame"
#~put "poison.frame$last.dump" 0 ...
#~
objects()
#~put "poison.frame$.Last.value" 0 ...
#~
objects(pos = 2)
#~put "poison.frame$last.dump" 0 ...
#~
objects(w = 2)
#~put "poison.frame$.Last.value" 0 ...
#~
poison.frame
#~error: couldn't open file ".Data.hp/poison.frame" to read S object "poison.frame"
#~
poison.frame
#~get ".Data.hp/poison.frame" 823535045 ...
#~
plot.factor(poison.frame)
#~get ".Data.hp/poison.frame" 823535045 ...
#~put "poison.frame$.Last.value" 0 ...
#~
interaction.plot(Poison, Tmt, poison)
#~get "/usr/statlocal/lib/Splus3.3/s/.Datasets/letters" -1 ...
#~
interaction.plot(Tmt, Poison, poison)
#~get "/usr/statlocal/lib/Splus3.3/s/.Datasets/letters" -1 ...
#~
anova(mymod <- aov(poison ~ Poison + Tmt + Poison * Tmt))
#~get "/usr/statlocal/lib/Splus3.3/s/.Datasets/.Machine" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/stat/.Datasets/store.co.helmert" -1 ...
#~put "poison.frame$mymod" 0 ...
#~put "poison.frame$.Last.value" 0 ...
#~
plot.factor(data.frame(resid(mymod), Poison, Tmt))
#~put "poison.frame$.Last.value" 0 ...
#~
plot(fitted(mymod), resid(mymod))
#~
object()
#~error: couldn't find function "object"
#~
gas.frame
#~error: Object "gas.frame" not found
#~
ls()
#~put "poison.frame$.Last.value" 0 ...
#~
objects(w = 2)
#~put "poison.frame$.Last.value" 0 ...
#~
detach(1, save = F)
#~put ".Data.hp/.Last.value" 823537172 ...
#~
ls()
#~put ".Data.hp/.Last.value" 823537174 ...
#~
#~End session: Time: 823537177; Process: 1925
#~New session: Version: S Wed Mar 29 17:35:08 PST 1995, Time: Mon Feb  5 11:19:39 1996
#~get "/usr/statlocal/lib/Splus3.3/splus/.Datasets/..min.script.id" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Datasets/.Copyright" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Datasets/version" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/.laenv" -1 ...
#~
attach(gas.frame, pos = 1)
#~get ".Data.hp/gas.frame" 823537122 ...
#~put "gas.frame$.Last.value" 0 ...
#~
ls()
#~put "gas.frame$.Last.value" 0 ...
#~
motif()
#~get "/usr/statlocal/lib/Splus3.3/splus/.Datasets/.PostScript.Options" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Datasets/ps.preamble.std" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Datasets/ps.setcolor.rgb" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Datasets/ps.setfont.std" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Datasets/ps.paper.regions" -1 ...
#~put "gas.frame$.Last.value" 0 ...
#~
par(mfrow = c(2, 2))
#~put "gas.frame$.Last.value" 0 ...
#~
plot.factor(gas.frame)
#~get ".Data.hp/gas.frame" 823537122 ...
#~put "gas.frame$.Last.value" 0 ...
#~
anova(mymod <- aov(resp ~ gas + car + driver))
#~get "/usr/statlocal/lib/Splus3.3/s/.Datasets/.Machine" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/stat/.Datasets/store.co.helmert" -1 ...
#~put "gas.frame$mymod" 0 ...
#~put "gas.frame$.Last.value" 0 ...
#~
plot(fitted(mymod), residual(mymod))
#~error: couldn't find function "residual"
#~put "gas.frame$last.dump" 0 ...
#~
plot(fitted(mymod), residuals(mymod))
#~
options()
#~put "gas.frame$.Last.value" 0 ...
#~
?plot
#~put "gas.frame$.Last.value" 0 ...
#~
par
#~
par()
#~get "/usr/statlocal/lib/Splus3.3/s/.Datasets/gr.pars" -1 ...
#~put "gas.frame$.Last.value" 0 ...
#~
options$ask
#~put "gas.frame$.Last.value" 0 ...
#~
options()
#~put "gas.frame$.Last.value" 0 ...
#~
options()$ask
#~put "gas.frame$.Last.value" 0 ...
#~
par(mfrow = c(1, 1))
#~put "gas.frame$.Last.value" 0 ...
#~
plot(fitted(mymod), residuals(mymod))
#~
plot.factor(gas.frame)
#~get ".Data.hp/gas.frame" 823537122 ...
#~put "gas.frame$.Last.value" 0 ...
#~
plot.factor
#~
?plot.factor
#~put "gas.frame$.Last.value" 0 ...
#~
q()
#~End session: Time: 823537566; Process: 2280
#~New session: Version: S Wed Mar 29 17:35:08 PST 1995, Time: Wed Feb 14 12:52:02 1996
#~get "/usr/statlocal/lib/Splus3.3/splus/.Datasets/..min.script.id" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Datasets/.Copyright" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Datasets/version" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/.laenv" -1 ...
#~
chem <- read.table("chemical.dat", header = T)
#~put ".Data.hp/chem" 824320330 ...
#~
mychem <- data.frame(chem$Run, Temp = as.factor(chem$Temp), Conc = as.factor(
	chem$Conc), Caty = chem$Caty, Yield = chem$Yield)
#~get ".Data.hp/chem" 824320330 ...
#~put ".Data.hp/mychem" 824320405 ...
#~
mychem <- data.frame(Run = chem$Run, Temp = as.factor(chem$Temp), Conc = 
	as.factor(chem$Conc), Caty = chem$Caty, Yield = chem$Yield)
#~get ".Data.hp/chem" 824320330 ...
#~put ".Data.hp/mychem" 824320416 ...
#~
mychem
#~get ".Data.hp/mychem" 824320416 ...
#~
chem <- mychem
#~get ".Data.hp/mychem" 824320416 ...
#~put ".Data.hp/chem" 824320483 ...
#~
rm(mychem)
#~remove ".Data.hp/mychem" 0 ...
#~
chem
#~get ".Data.hp/chem" 824320483 ...
#~
motif()
#~get "/usr/statlocal/lib/Splus3.3/splus/.Datasets/.PostScript.Options" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Datasets/ps.preamble.std" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Datasets/ps.setcolor.rgb" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Datasets/ps.setfont.std" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Datasets/ps.paper.regions" -1 ...
#~put ".Data.hp/.Last.value" 824320503 ...
#~
plot.factor(chem, y = Yield)
#~error: Object "Yield" not found
#~put ".Data.hp/last.dump" 824320515 ...
#~
attach(chem, pos = 1)
#~get ".Data.hp/chem" 824320483 ...
#~put "chem$.Last.value" 0 ...
#~
plot.factor(chem, y = Yield)
#~get ".Data.hp/chem" 824320483 ...
#~put "chem$.Last.value" 0 ...
#~
pars(mfrow = c(2, 2))
#~error: couldn't find function "pars"
#~
par(mfrow = c(2, 2))
#~put "chem$.Last.value" 0 ...
#~
plot.factor(chem, y = Yield)
#~get ".Data.hp/chem" 824320483 ...
#~put "chem$.Last.value" 0 ...
#~
newframe()
#~error: couldn't find function "newframe"
#~
search()
#~put "chem$.Last.value" 0 ...
#~
objects(w = 3)
#~put "chem$.Last.value" 0 ...
#~
clearplot
#~
clearplot()
#~put "chem$.Last.value" 0 ...
#~
interaction.plot(Temp, Conc, Yield)
#~get "/usr/statlocal/lib/Splus3.3/s/.Datasets/letters" -1 ...
#~
interaction.plot(Temp, Caty, Yield)
#~get "/usr/statlocal/lib/Splus3.3/s/.Datasets/letters" -1 ...
#~
interaction.plot(Conc, Caty, Yield)
#~get "/usr/statlocal/lib/Splus3.3/s/.Datasets/letters" -1 ...
#~
mymod <- aov(Yield ~ Temp + Conc + Caty + Temp * Caty)
#~get "/usr/statlocal/lib/Splus3.3/s/.Datasets/.Machine" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/stat/.Datasets/store.co.helmert" -1 ...
#~put "chem$mymod" 0 ...
#~
plot(mymod)
#~get "/usr/statlocal/lib/Splus3.3/s/.Datasets/gr.pars" -1 ...
#~put "chem$.Last.value" 0 ...
#~
plot(mymod)
#~get "/usr/statlocal/lib/Splus3.3/s/.Datasets/gr.pars" -1 ...
#~put "chem$.Last.value" 0 ...
#~
clearplot
#~
clearplot()
#~put "chem$.Last.value" 0 ...
#~
plot(mymod)
#~get "/usr/statlocal/lib/Splus3.3/s/.Datasets/gr.pars" -1 ...
#~put "chem$.Last.value" 0 ...
#~
par(mfrow = c(2, 3))
#~put "chem$.Last.value" 0 ...
#~
clearplot()
#~put "chem$.Last.value" 0 ...
#~
plot(mymod)
#~get "/usr/statlocal/lib/Splus3.3/s/.Datasets/gr.pars" -1 ...
#~put "chem$.Last.value" 0 ...
#~
mymod <- aov(Yield ~ Temp + Conc + Caty + Temp * Caty + Temp * Conc + Caty * 
	Conc, data = chem)
#~get ".Data.hp/chem" 824320483 ...
#~get "/usr/statlocal/lib/Splus3.3/s/.Datasets/.Machine" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/stat/.Datasets/store.co.helmert" -1 ...
#~put "chem$mymod" 0 ...
#~
mymod <- aov(Yield ~ Temp + Conc + Caty + Temp * Caty)
#~get "/usr/statlocal/lib/Splus3.3/s/.Datasets/.Machine" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/stat/.Datasets/store.co.helmert" -1 ...
#~put "chem$mymod" 0 ...
#~
clearplot()
#~put "chem$.Last.value" 0 ...
#~
qqnorm(mymod)
#~put "chem$.Last.value" 0 ...
#~
par(mfrow = c(1, 1))
#~put "chem$.Last.value" 0 ...
#~
qqnorm(mymod)
#~put "chem$.Last.value" 0 ...
#~
q()
#~End session: Time: 824321316; Process: 11484
#~New session: Version: S Wed Mar 29 17:35:08 PST 1995, Time: Wed Feb 14 13:40:51 1996
#~get "/usr/statlocal/lib/Splus3.3/splus/.Datasets/..min.script.id" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Datasets/.Copyright" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Datasets/version" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/.laenv" -1 ...
#~
sum(c(5, 6, 4, 6, 8, 8, 1.5, 1, 1, 3, 5, 3, 12, 8, 4, 1))
#~put ".Data.hp/.Last.value" 824323278 ...
#~
sum(c(5, 6, 4, 6, 8, 8, 1.5, 1, 1, 3, 5, 3, 12, 8, 4, 1)) * 250
#~put ".Data.hp/.Last.value" 824323281 ...
#~
#~End session: Time: 824323549; Process: 11621
#~New session: Version: S Wed Mar 29 17:35:08 PST 1995, Time: Wed Feb 21 08:54:00 1996
#~get "/usr/statlocal/lib/Splus3.3/splus/.Datasets/..min.script.id" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Datasets/.Copyright" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Datasets/version" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/.laenv" -1 ...
#~
hay <- scan("hayfever.dat")
#~put ".Data.hp/hay" 824910852 ...
#~
hay
#~get ".Data.hp/hay" 824910852 ...
#~
fac.design(c(3, 3, 3))
#~get "/usr/statlocal/lib/Splus3.3/s/.Datasets/LETTERS" -1 ...
#~put ".Data.hp/.Last.value" 824910893 ...
#~
fac.design(c(3, 3, 4))
#~get "/usr/statlocal/lib/Splus3.3/s/.Datasets/LETTERS" -1 ...
#~put ".Data.hp/.Last.value" 824910908 ...
#~
fac.design(c(4, 3, 3))
#~get "/usr/statlocal/lib/Splus3.3/s/.Datasets/LETTERS" -1 ...
#~put ".Data.hp/.Last.value" 824910921 ...
#~
fac.design(c(3, 3, 4))
#~get "/usr/statlocal/lib/Splus3.3/s/.Datasets/LETTERS" -1 ...
#~put ".Data.hp/.Last.value" 824910999 ...
#~
fac.design(c(3, 4, 4))
#~get "/usr/statlocal/lib/Splus3.3/s/.Datasets/LETTERS" -1 ...
#~put ".Data.hp/.Last.value" 824911110 ...
#~
fac.design(c(3, 4, 3))
#~get "/usr/statlocal/lib/Splus3.3/s/.Datasets/LETTERS" -1 ...
#~put ".Data.hp/.Last.value" 824911127 ...
#~
?fac.design
#~put ".Data.hp/.Last.value" 824911171 ...
#~
fac.design(c(3, 4, 3), fac.names = c("I1", "Rep", "I2"))
#~put ".Data.hp/last.dump" 824911196 ...
#~
fac.design(c(3, 4, 3), factor.names = c("I1", "Rep", "I2"))
#~put ".Data.hp/.Last.value" 824911206 ...
#~
fac.design(c(3, 4, 3), factor.names = c("I-one", "Rep", "I-two"))
#~put ".Data.hp/.Last.value" 824911230 ...
#~
hayfever <- fac.design(c(3, 4, 3), factor.names = c("I-one", "Rep", "I-two"))
#~put ".Data.hp/hayfever" 824911245 ...
#~
hayfever$Resp <- hay
#~get ".Data.hp/hay" 824910852 ...
#~get ".Data.hp/hayfever" 824911245 ...
#~put ".Data.hp/hayfever" 824911260 ...
#~
hayfever
#~get ".Data.hp/hayfever" 824911260 ...
#~
attach(hayfever, pos = 1)
#~get ".Data.hp/hayfever" 824911260 ...
#~put "hayfever$.Last.value" 0 ...
#~
tapply(hayfever, I - one, mean)
#~error: Object "one" not found
#~put "hayfever$last.dump" 0 ...
#~
search()
#~put "hayfever$.Last.value" 0 ...
#~
tapply(hayfever, I - one, mean)
#~error: Object "one" not found
#~put "hayfever$last.dump" 0 ...
#~
ls()
#~put "hayfever$.Last.value" 0 ...
#~
detach(1, save = F)
#~put ".Data.hp/.Last.value" 824911351 ...
#~
hayfever <- fac.design(c(3, 4, 3), factor.names = c("IO", "Rep", "IT"))
#~put ".Data.hp/hayfever" 824911363 ...
#~
hayfever$Resp <- hay
#~get ".Data.hp/hay" 824910852 ...
#~get ".Data.hp/hayfever" 824911363 ...
#~put ".Data.hp/hayfever" 824911368 ...
#~
attach(hayfever, pos = 1)
#~get ".Data.hp/hayfever" 824911368 ...
#~put "hayfever$.Last.value" 0 ...
#~
tapply(hayfever, I - one, mean)
#~error: Object "one" not found
#~put "hayfever$last.dump" 0 ...
#~
tapply(hayfever, IO, mean)
#~error: Data and all indices must have same length
#~put "hayfever$last.dump" 0 ...
#~
ls()
#~put "hayfever$.Last.value" 0 ...
#~
Resp
#~
tapply(Resp, IO, mean)
#~put "hayfever$.Last.value" 0 ...
#~
tapply(Resp, IT, mean)
#~put "hayfever$.Last.value" 0 ...
#~
anova(Resp ~ IO + IT)
#~error: No specific or default method found
#~put "hayfever$last.dump" 0 ...
#~
anova(aov(Resp ~ IO + IT))
#~get "/usr/statlocal/lib/Splus3.3/s/.Datasets/.Machine" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/stat/.Datasets/store.co.helmert" -1 ...
#~put "hayfever$.Last.value" 0 ...
#~
anova(aov(Resp ~ IO * IT))
#~get "/usr/statlocal/lib/Splus3.3/s/.Datasets/.Machine" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/stat/.Datasets/store.co.helmert" -1 ...
#~put "hayfever$.Last.value" 0 ...
#~
#~End session: Time: 824911744; Process: 2930
#~New session: Version: S Wed Mar 29 17:35:08 PST 1995, Time: Thu Feb 22 10:36:51 1996
#~get "/usr/statlocal/lib/Splus3.3/splus/.Datasets/..min.script.id" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Datasets/.Copyright" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Datasets/version" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/.laenv" -1 ...
#~
?plot.factor
#~put ".Data.hp/.Last.value" 825003432 ...
#~
#~End session: Time: 825003491; Process: 8282
#~New session: Version: S Wed Mar 29 17:35:08 PST 1995, Time: Mon Feb 26 13:48:31 1996
#~get "/usr/statlocal/lib/Splus3.3/splus/.Datasets/..min.script.id" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Datasets/.Copyright" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Datasets/version" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/.laenv" -1 ...
#~
Reiss <- fac.design(c(8, 8))
#~get "/usr/statlocal/lib/Splus3.3/s/.Datasets/LETTERS" -1 ...
#~put ".Data.hp/Reiss" 825360542 ...
#~
Reiss
#~get ".Data.hp/Reiss" 825360542 ...
#~
?fac.design
#~put ".Data.hp/.Last.value" 825360568 ...
#~
Reiss <- fac.design(c(8, 8), factor.names = c("First", "Second"))
#~put ".Data.hp/Reiss" 825360584 ...
#~
Reiss
#~get ".Data.hp/Reiss" 825360584 ...
#~
Reiss <- fac.design(c(8, 8), factor.names = c("Second", "First"))
#~put ".Data.hp/Reiss" 825360627 ...
#~
?fac.design
#~put ".Data.hp/.Last.value" 825360635 ...
#~
Reiss <- fac.design(c(8, 8), factor.names = c("Second", "First"))
#~put ".Data.hp/Reiss" 825360636 ...
#~
Reiss$Resp <- scan("reiss.dat")
#~get ".Data.hp/Reiss" 825360636 ...
#~put ".Data.hp/Reiss" 825360681 ...
#~
Reiss[1:10,  ]
#~get ".Data.hp/Reiss" 825360681 ...
#~put ".Data.hp/.Last.value" 825360688 ...
#~
mymod <- glm(Resp ~ First + Second, error = poisson, data = Reiss)
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.variances" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.links" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/stat/.Functions/glm.weights" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/stat/.Datasets/store.co.helmert" -1 ...
#~get ".Data.hp/Reiss" 825360681 ...
#~put ".Data.hp/mymod" 825360747 ...
#~
summary(mymod)
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.variances" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.links" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/stat/.Functions/glm.weights" -1 ...
#~get ".Data.hp/mymod" 825360747 ...
#~put ".Data.hp/.Last.value" 825360754 ...
#~
anova(mymod)
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.variances" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.links" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/stat/.Functions/glm.weights" -1 ...
#~get ".Data.hp/mymod" 825360747 ...
#~get ".Data.hp/Reiss" 825360681 ...
#~put ".Data.hp/.Last.value" 825360764 ...
#~
mymod <- glm(Resp ~ First + Second, error = "poisson", data = Reiss)
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.variances" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.links" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/stat/.Functions/glm.weights" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/stat/.Datasets/store.co.helmert" -1 ...
#~get ".Data.hp/Reiss" 825360681 ...
#~put ".Data.hp/mymod" 825360791 ...
#~
anova(mymod)
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.variances" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.links" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/stat/.Functions/glm.weights" -1 ...
#~get ".Data.hp/mymod" 825360791 ...
#~get ".Data.hp/Reiss" 825360681 ...
#~put ".Data.hp/.Last.value" 825360793 ...
#~
fitted(mymod)
#~get ".Data.hp/mymod" 825360791 ...
#~put ".Data.hp/.Last.value" 825360801 ...
#~
?glm
#~put ".Data.hp/.Last.value" 825360832 ...
#~
mymod <- glm(Resp ~ First + Second, family = poisson, data = Reiss)
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.variances" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.links" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/stat/.Functions/glm.weights" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/s/.Datasets/.Machine" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/stat/.Datasets/store.co.helmert" -1 ...
#~get ".Data.hp/Reiss" 825360681 ...
#~put ".Data.hp/mymod" 825360847 ...
#~
anova(mymod)
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.variances" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.links" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/stat/.Functions/glm.weights" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/s/.Datasets/.Machine" -1 ...
#~get ".Data.hp/mymod" 825360847 ...
#~get ".Data.hp/Reiss" 825360681 ...
#~put ".Data.hp/.Last.value" 825360851 ...
#~
?table
#~put ".Data.hp/.Last.value" 825360921 ...
#~
table(8, 8, fitted(mymod))
#~error: All arguments must have the same length
#~put ".Data.hp/last.dump" 825360935 ...
#~
table(8, 8)
#~put ".Data.hp/.Last.value" 825360966 ...
#~
table(1:8, 1:8)
#~put ".Data.hp/.Last.value" 825360986 ...
#~
?table
#~put ".Data.hp/.Last.value" 825361030 ...
#~
?matrix
#~put ".Data.hp/.Last.value" 825361046 ...
#~
matrix(fitted(mymod), 8, 8)
#~get ".Data.hp/mymod" 825360847 ...
#~put ".Data.hp/.Last.value" 825361058 ...
#~
matrix(round(fitted(mymod)), 8, 8)
#~get ".Data.hp/mymod" 825360847 ...
#~put ".Data.hp/.Last.value" 825361085 ...
#~
?round
#~put ".Data.hp/.Last.value" 825361090 ...
#~
matrix(round(fitted(mymod), digits = 2), 8, 8)
#~get ".Data.hp/mymod" 825360847 ...
#~put ".Data.hp/.Last.value" 825361099 ...
#~
matrix(round(fitted(mymod), digits = 3), 8, 8)
#~get ".Data.hp/mymod" 825360847 ...
#~put ".Data.hp/.Last.value" 825361103 ...
#~
deviance(mymod)
#~get ".Data.hp/mymod" 825360847 ...
#~put ".Data.hp/.Last.value" 825361379 ...
#~
anova(mymod)
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.variances" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.links" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/stat/.Functions/glm.weights" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/s/.Datasets/.Machine" -1 ...
#~get ".Data.hp/mymod" 825360847 ...
#~get ".Data.hp/Reiss" 825360681 ...
#~put ".Data.hp/.Last.value" 825361385 ...
#~
df(mymod)
#~error: no df1 arg
#~put ".Data.hp/last.dump" 825361396 ...
#~
?deviance
#~put ".Data.hp/.Last.value" 825361411 ...
#~
mymod
#~get ".Data.hp/mymod" 825360847 ...
#~
anova(mymod)
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.variances" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.links" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/stat/.Functions/glm.weights" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/s/.Datasets/.Machine" -1 ...
#~get ".Data.hp/mymod" 825360847 ...
#~get ".Data.hp/Reiss" 825360681 ...
#~put ".Data.hp/.Last.value" 825361433 ...
#~
Pearson <- function(mymod)
{
 + sum((mymod$y - fitted(mymod))^2/fitted(mymod))
}
#~put ".Data.hp/Pearson" 825361560 ...
#~
Pearson
#~
Pearson
#~
names(mymod)
#~get ".Data.hp/mymod" 825360847 ...
#~put ".Data.hp/.Last.value" 825361847 ...
#~
mymod$rank
#~get ".Data.hp/mymod" 825360847 ...
#~put ".Data.hp/.Last.value" 825361884 ...
#~
mymod
#~get ".Data.hp/mymod" 825360847 ...
#~
glm.print
#~error: Object "glm.print" not found
#~
print.glm
#~
mymod$nobs - mymod$rank
#~get ".Data.hp/mymod" 825360847 ...
#~put ".Data.hp/.Last.value" 825361949 ...
#~
Pearson
#~
mymod$df.resid
#~get ".Data.hp/mymod" 825360847 ...
#~put ".Data.hp/.Last.value" 825362735 ...
#~
Pearson <- function(mymod)
{
# This is not completely bullet proof.
	nobs <- length(mymod$residuals)
	rdf <- mymod$df.resid
	if(is.null(rdf))
		rdf <- nobs - rank
	cat("\nDegrees of Freedom:", nobs, "Total;", rdf, "Residual\n")
	cat("Residual Pearson X^2", format(sum((mymod$y - fitted(mymod))^2/
		fitted(mymod))), "\n")
	invisible(mymod)
}
#~put ".Data.hp/Pearson" 825362841 ...
#~
Pearson(mymod)
#~get ".Data.hp/mymod" 825360847 ...
#~put ".Data.hp/.Last.value" 825362847 ...
#~
fix(Pearson)
#~put ".Data.hp/Pearson" 825362873 ...
#~put ".Data.hp/.Last.fixed" 825362873 ...
#~put ".Data.hp/.Last.value" 825362874 ...
#~
?fix
#~put ".Data.hp/.Last.value" 825362882 ...
#~
fix
#~
?fix
#~put ".Data.hp/.Last.value" 825362897 ...
#~
search()
#~put ".Data.hp/.Last.value" 825362901 ...
#~
objects(w = 2)
#~put ".Data.hp/.Last.value" 825362907 ...
#~
emacs
#~
fix
#~
fix(Pearson)
#~put ".Data.hp/Pearson" 825362947 ...
#~put ".Data.hp/.Last.fixed" 825362947 ...
#~put ".Data.hp/.Last.value" 825362947 ...
#~
Pearson(mymod)
#~get ".Data.hp/mymod" 825360847 ...
#~put ".Data.hp/.Last.value" 825362954 ...
#~
Pearson
#~
resid(mymod)
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.variances" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.links" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/stat/.Functions/glm.weights" -1 ...
#~get ".Data.hp/mymod" 825360847 ...
#~put ".Data.hp/.Last.value" 825363038 ...
#~
names(mymod)
#~get ".Data.hp/mymod" 825360847 ...
#~put ".Data.hp/.Last.value" 825363050 ...
#~
mymod$residuals
#~get ".Data.hp/mymod" 825360847 ...
#~put ".Data.hp/.Last.value" 825363061 ...
#~
?glm
#~put ".Data.hp/.Last.value" 825363107 ...
#~
help.start()
#~put ".Data.hp/.Last.value" 825363112 ...
#~
residuals(mymod, type = "pearson")
#~get ".Data.hp/mymod" 825360847 ...
#~put ".Data.hp/.Last.value" 825363156 ...
#~
matrix(round(residuals(mymod, type = "pearson"), dig = 3), 8, 8)
#~get ".Data.hp/mymod" 825360847 ...
#~put ".Data.hp/.Last.value" 825363335 ...
#~
trans(matrix(round(residuals(mymod, type = "pearson"), dig = 3), 8, 8))
#~error: couldn't find function "trans"
#~
tran(matrix(round(residuals(mymod, type = "pearson"), dig = 3), 8, 8))
#~error: couldn't find function "tran"
#~
t(matrix(round(residuals(mymod, type = "pearson"), dig = 3), 8, 8))
#~get ".Data.hp/mymod" 825360847 ...
#~put ".Data.hp/.Last.value" 825363482 ...
#~
t(matrix(round(fitted(mymod), digits = 3), 8, 8))
#~get ".Data.hp/mymod" 825360847 ...
#~put ".Data.hp/.Last.value" 825363528 ...
#~
residuals(mymod, type = "pearson")[1]^2
#~get ".Data.hp/mymod" 825360847 ...
#~put ".Data.hp/.Last.value" 825363601 ...
#~
#~End session: Time: 825363987; Process: 19258
#~New session: Version: S Wed Mar 29 17:35:08 PST 1995, Time: Tue Feb 27 08:37:07 1996
#~get "/usr/statlocal/lib/Splus3.3/splus/.Datasets/..min.script.id" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Datasets/.Copyright" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Datasets/version" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/.laenv" -1 ...
#~
help.start()
#~put ".Data.hp/.Last.value" 825428235 ...
#~
#~End session: Time: 825428290; Process: 20278
#~New session: Version: S Wed Mar 29 17:35:08 PST 1995, Time: Mon Mar  4 10:14:42 1996
#~get "/usr/statlocal/lib/Splus3.3/splus/.Datasets/..min.script.id" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Datasets/.Copyright" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Datasets/version" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/.laenv" -1 ...
#~
fac.design(2, 2, 2)
#~put ".Data.hp/.Last.value" 825952492 ...
#~
fac.design(c(2, 2, 2))
#~get "/usr/statlocal/lib/Splus3.3/s/.Datasets/LETTERS" -1 ...
#~put ".Data.hp/.Last.value" 825952496 ...
#~
fac.design(c(2, 2, 2), c("Dose", "Age", "Cancer"))
#~put ".Data.hp/.Last.value" 825952531 ...
#~
reddye <- fac.design(c(2, 2, 2), c("Dose", "Age", "Cancer"))
#~put ".Data.hp/reddye" 825952562 ...
#~
reddye$Resp <- c(4, 7, 0, 7, 26, 16, 14, 14)
#~get ".Data.hp/reddye" 825952562 ...
#~put ".Data.hp/reddye" 825952583 ...
#~
reddye
#~get ".Data.hp/reddye" 825952583 ...
#~
anova(mymod <- glm(Resp ~ Dose, family = Poisson, data = reddye))
#~error: Object "Poisson" not found
#~put ".Data.hp/last.dump" 825952759 ...
#~
anova(mymod <- glm(Resp ~ Dose, family = poisson, data = reddye))
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.variances" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.links" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/stat/.Functions/glm.weights" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/s/.Datasets/.Machine" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/stat/.Datasets/store.co.helmert" -1 ...
#~get ".Data.hp/reddye" 825952583 ...
#~put ".Data.hp/mymod" 825952767 ...
#~put ".Data.hp/.Last.value" 825952767 ...
#~
Pearson(mymod)
#~get ".Data.hp/mymod" 825952767 ...
#~put ".Data.hp/.Last.value" 825952793 ...
#~
anova(mymod <- glm(Resp ~ Dose + Age, family = poisson, data = reddye))
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.variances" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.links" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/stat/.Functions/glm.weights" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/s/.Datasets/.Machine" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/stat/.Datasets/store.co.helmert" -1 ...
#~get ".Data.hp/reddye" 825952583 ...
#~put ".Data.hp/mymod" 825953235 ...
#~put ".Data.hp/.Last.value" 825953235 ...
#~
anova(mymod <- glm(Resp ~ Dose * Age, family = poisson, data = reddye))
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.variances" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.links" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/stat/.Functions/glm.weights" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/s/.Datasets/.Machine" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/stat/.Datasets/store.co.helmert" -1 ...
#~get ".Data.hp/reddye" 825952583 ...
#~put ".Data.hp/mymod" 825953265 ...
#~put ".Data.hp/.Last.value" 825953265 ...
#~
anova(mymod <- glm(Resp ~ Dose * Age + Age * Cancer + Dose * Cancer, family = 
	poisson, data = reddye))
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.variances" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.links" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/stat/.Functions/glm.weights" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/s/.Datasets/.Machine" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/stat/.Datasets/store.co.helmert" -1 ...
#~get ".Data.hp/reddye" 825952583 ...
#~put ".Data.hp/mymod" 825953367 ...
#~put ".Data.hp/.Last.value" 825953367 ...
#~
anova(update(mymod, . ~ . - Age:Cancer))
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.variances" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.links" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/stat/.Functions/glm.weights" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/s/.Datasets/.Machine" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/stat/.Datasets/store.co.helmert" -1 ...
#~get ".Data.hp/mymod" 825953367 ...
#~get ".Data.hp/reddye" 825952583 ...
#~put ".Data.hp/.Last.value" 825953459 ...
#~
qchi(1, 0.05)
#~error: couldn't find function "qchi"
#~
?chi
#~put ".Data.hp/.Last.value" 825953555 ...
#~
help.start()
#~put ".Data.hp/.Last.value" 825953575 ...
#~
qchisq(0.05, 1)
#~put ".Data.hp/.Last.value" 825953615 ...
#~
qchisq(0.95, 1)
#~put ".Data.hp/.Last.value" 825953622 ...
#~
anova(update(mymod, . ~ . - Age:Dose))
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.variances" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.links" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/stat/.Functions/glm.weights" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/s/.Datasets/.Machine" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/stat/.Datasets/store.co.helmert" -1 ...
#~get ".Data.hp/mymod" 825953367 ...
#~get ".Data.hp/reddye" 825952583 ...
#~put ".Data.hp/.Last.value" 825953712 ...
#~
anova(mymod <- update(mymod, . ~ . - Age:Cancer))
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.variances" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.links" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/stat/.Functions/glm.weights" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/s/.Datasets/.Machine" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/stat/.Datasets/store.co.helmert" -1 ...
#~get ".Data.hp/mymod" 825953367 ...
#~put ".Data.hp/mymod" 825953728 ...
#~get ".Data.hp/reddye" 825952583 ...
#~put ".Data.hp/.Last.value" 825953728 ...
#~
anova(mymod <- update(mymod, . ~ . - Age:Dose))
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.variances" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.links" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/stat/.Functions/glm.weights" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/s/.Datasets/.Machine" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/stat/.Datasets/store.co.helmert" -1 ...
#~get ".Data.hp/mymod" 825953728 ...
#~put ".Data.hp/mymod" 825953736 ...
#~get ".Data.hp/reddye" 825952583 ...
#~put ".Data.hp/.Last.value" 825953736 ...
#~
Pearson(mymod)
#~get ".Data.hp/mymod" 825953736 ...
#~put ".Data.hp/.Last.value" 825953832 ...
#~
qchisq(0.95, 3)
#~put ".Data.hp/.Last.value" 825953902 ...
#~
motif()
#~get "/usr/statlocal/lib/Splus3.3/splus/.Datasets/.PostScript.Options" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Datasets/ps.preamble.std" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Datasets/ps.setcolor.rgb" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Datasets/ps.setfont.std" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Datasets/ps.paper.regions" -1 ...
#~put ".Data.hp/.Last.value" 825957124 ...
#~
plot(mymod)
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.variances" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.links" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/stat/.Functions/glm.weights" -1 ...
#~get ".Data.hp/mymod" 825953736 ...
#~put ".Data.hp/.Last.value" 825957129 ...
#~
par(mfrow = c(3, 3))
#~put ".Data.hp/.Last.value" 825957142 ...
#~
plot(mymod)
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.variances" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.links" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/stat/.Functions/glm.weights" -1 ...
#~get ".Data.hp/mymod" 825953736 ...
#~put ".Data.hp/.Last.value" 825957144 ...
#~
par(mfrow = c(2, 3))
#~put ".Data.hp/.Last.value" 825957150 ...
#~
plot(mymod)
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.variances" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.links" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/stat/.Functions/glm.weights" -1 ...
#~get ".Data.hp/mymod" 825953736 ...
#~put ".Data.hp/.Last.value" 825957152 ...
#~
fac.design(c(6, 2, 2))
#~get "/usr/statlocal/lib/Splus3.3/s/.Datasets/LETTERS" -1 ...
#~put ".Data.hp/.Last.value" 825959539 ...
#~
fac.design(c(2, 2, 6))
#~get "/usr/statlocal/lib/Splus3.3/s/.Datasets/LETTERS" -1 ...
#~put ".Data.hp/.Last.value" 825959552 ...
#~
fac.design(c(2, 2, 6), c("Admit", "Sex", "Major"))
#~put ".Data.hp/.Last.value" 825959621 ...
#~
berk <- fac.design(c(2, 2, 6), c("Admit", "Sex", "Major"))
#~put ".Data.hp/berk" 825959650 ...
#~
berk$Resp <- scan("/tmp/berk")
#~get ".Data.hp/berk" 825959650 ...
#~put ".Data.hp/berk" 825959707 ...
#~
berk
#~get ".Data.hp/berk" 825959707 ...
#~
?fac.design
#~put ".Data.hp/.Last.value" 825959760 ...
#~
berk$Resp <- scan("/tmp/berk")
#~get ".Data.hp/berk" 825959707 ...
#~put ".Data.hp/berk" 825959834 ...
#~
berk <- fac.design(c(2,2,6),fnames=list(Admit=c("Yes","No"),
Sex=c("Men","Women"), Major=c("A","B","C","D",E","F"))
#~error: statement with syntax error
#~
berk <- fac.design(c(2, 2, 6), fnames = list(Admit = c("Yes", "No"), Sex = c(
	"Men", "Women"), Major = c("A", "B", "C", "D", "E", "F")))
#~put ".Data.hp/last.dump" 825959950 ...
#~
berk <- fac.design(c(2, 2, 6), fnames = list(Admit = c("Yes", "No"), Sex = c(
	"Men", "Women"), Major = c("A", "B", "C", "D", "E", "F")))
#~put ".Data.hp/last.dump" 825960022 ...
#~
berk <- fac.design(c(2, 2, 6), fac = list(Admit = c("Yes", "No"), Sex = c("Men",
	"Women"), Major = c("A", "B", "C", "D", "E", "F")))
#~put ".Data.hp/berk" 825960048 ...
#~
berk$Resp <- scan("/tmp/berk")
#~get ".Data.hp/berk" 825960048 ...
#~put ".Data.hp/berk" 825960057 ...
#~
berk
#~get ".Data.hp/berk" 825960057 ...
#~
attach(berk, pos = 1)
#~get ".Data.hp/berk" 825960057 ...
#~put "berk$.Last.value" 0 ...
#~
ls()
#~put "berk$.Last.value" 0 ...
#~
Resp[Admit == "Yes"] + Resp[Admit == "No"]
#~put "berk$.Last.value" 0 ...
#~
applied <- Resp[Admit == "Yes"] + Resp[Admit == "No"]
#~put "berk$applied" 0 ...
#~
applied
#~
Resp[Admit == "Yes"]/applied
#~put "berk$.Last.value" 0 ...
#~
Resp[Sex = "Men"]
#~put "berk$.Last.value" 0 ...
#~
Resp[Sex == "Men"]
#~put "berk$.Last.value" 0 ...
#~
Resp[Sex == "Men" && Admit == "Yes"]
#~put "berk$.Last.value" 0 ...
#~
Resp[Sex == "Men" & Admit == "Yes"]
#~put "berk$.Last.value" 0 ...
#~
sum(Resp[Sex == "Men"])
#~put "berk$.Last.value" 0 ...
#~
sum(Resp[Sex == "Men" & Admit == "Yes"])
#~put "berk$.Last.value" 0 ...
#~
sum(Resp[Sex == "Women" & Admit == "Yes"])
#~put "berk$.Last.value" 0 ...
#~
sum(Resp[Sex == "Women"])
#~put "berk$.Last.value" 0 ...
#~
1198/2691
#~put "berk$.Last.value" 0 ...
#~
557/1835
#~put "berk$.Last.value" 0 ...
#~
round(Resp[Admit == "Yes"]/applied, dig = 3)
#~put "berk$.Last.value" 0 ...
#~
dodo <- round(Resp[Admit == "Yes"]/applied, dig = 3)
#~put "berk$dodo" 0 ...
#~
odd
#~error: Object "odd" not found
#~
even
#~error: Object "even" not found
#~
dodo[(0:5) * 2 + 1]
#~put "berk$.Last.value" 0 ...
#~
dodo[(0:5) * 2 + 2]
#~put "berk$.Last.value" 0 ...
#~
dodo[Sex = "Men"]
#~put "berk$.Last.value" 0 ...
#~
berk
#~get ".Data.hp/berk" 825960057 ...
#~
ls()
#~put "berk$.Last.value" 0 ...
#~
mymod <- glm(Resp ~ Sex * Major, family = poisson)
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.variances" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.links" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/stat/.Functions/glm.weights" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/s/.Datasets/.Machine" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/stat/.Datasets/store.co.helmert" -1 ...
#~put "berk$mymod" 0 ...
#~
poisson
#~
anova(mymod)
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.variances" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.links" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/stat/.Functions/glm.weights" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/s/.Datasets/.Machine" -1 ...
#~put "berk$.Last.value" 0 ...
#~
anova(mymod <- glm(Resp ~ Sex * Major, family = poisson))
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.variances" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.links" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/stat/.Functions/glm.weights" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/s/.Datasets/.Machine" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/stat/.Datasets/store.co.helmert" -1 ...
#~put "berk$mymod" 0 ...
#~put "berk$.Last.value" 0 ...
#~
anova(mymod <- glm(Resp ~ Sex * Major + Admit, family = poisson))
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.variances" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.links" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/stat/.Functions/glm.weights" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/s/.Datasets/.Machine" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/stat/.Datasets/store.co.helmert" -1 ...
#~put "berk$mymod" 0 ...
#~put "berk$.Last.value" 0 ...
#~
anova(mymod <- glm(Resp ~ Sex * Major + Admit * Sex, family = poisson))
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.variances" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.links" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/stat/.Functions/glm.weights" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/s/.Datasets/.Machine" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/stat/.Datasets/store.co.helmert" -1 ...
#~put "berk$mymod" 0 ...
#~put "berk$.Last.value" 0 ...
#~
anova(mymod <- glm(Resp ~ Sex * Major + Admit * Major, family = poisson))
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.variances" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.links" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/stat/.Functions/glm.weights" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/s/.Datasets/.Machine" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/stat/.Datasets/store.co.helmert" -1 ...
#~put "berk$mymod" 0 ...
#~put "berk$.Last.value" 0 ...
#~
anova(mymod <- glm(Resp ~ Sex * Major + Admit * Major + Admit * Sex, family = 
	poisson))
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.variances" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.links" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/stat/.Functions/glm.weights" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/s/.Datasets/.Machine" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/stat/.Datasets/store.co.helmert" -1 ...
#~put "berk$mymod" 0 ...
#~put "berk$.Last.value" 0 ...
#~
anova(mymod <- glm(Resp ~ Sex * Major + Admit * Major, family = poisson))
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.variances" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.links" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/stat/.Functions/glm.weights" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/s/.Datasets/.Machine" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/stat/.Datasets/store.co.helmert" -1 ...
#~put "berk$mymod" 0 ...
#~put "berk$.Last.value" 0 ...
#~
qchisq(0.95, 6)
#~put "berk$.Last.value" 0 ...
#~
pchisq(21, 6)
#~put "berk$.Last.value" 0 ...
#~
anova(mymod <- glm(Resp ~ Sex * Major + Admit * Major + Admit * Sex, family = 
	poisson))
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.variances" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.links" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/stat/.Functions/glm.weights" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/s/.Datasets/.Machine" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/stat/.Datasets/store.co.helmert" -1 ...
#~put "berk$mymod" 0 ...
#~put "berk$.Last.value" 0 ...
#~
anova(mymod <- glm(Resp ~ Sex * Major + Admit * Major + Admit * Sex, family = 
	poisson))
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.variances" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.links" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/stat/.Functions/glm.weights" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/s/.Datasets/.Machine" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/stat/.Datasets/store.co.helmert" -1 ...
#~put "berk$mymod" 0 ...
#~put "berk$.Last.value" 0 ...
#~
pchisq(1.531, 1)
#~put "berk$.Last.value" 0 ...
#~
plot(mymod)
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.variances" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.links" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/stat/.Functions/glm.weights" -1 ...
#~put "berk$.Last.value" 0 ...
#~
plot(mymod)
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.variances" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.links" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/stat/.Functions/glm.weights" -1 ...
#~put "berk$.Last.value" 0 ...
#~
newplot()
#~error: couldn't find function "newplot"
#~
clearplot()
#~put "berk$.Last.value" 0 ...
#~
plot(mymod)
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.variances" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.links" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/stat/.Functions/glm.weights" -1 ...
#~put "berk$.Last.value" 0 ...
#~
residuals(mymod)
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.variances" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.links" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/stat/.Functions/glm.weights" -1 ...
#~put "berk$.Last.value" 0 ...
#~
min(residuals(mymod))
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.variances" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.links" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/stat/.Functions/glm.weights" -1 ...
#~put "berk$.Last.value" 0 ...
#~
fitted(mymod)
#~put "berk$.Last.value" 0 ...
#~
Resp
#~
search()
#~put "berk$.Last.value" 0 ...
#~
detach(1, save = F)
#~put ".Data.hp/.Last.value" 825962947 ...
#~
#~End session: Time: 825962949; Process: 9451
#~New session: Version: S Wed Mar 29 17:35:08 PST 1995, Time: Mon Mar  4 14:17:46 1996
#~get "/usr/statlocal/lib/Splus3.3/splus/.Datasets/..min.script.id" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Datasets/.Copyright" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Datasets/version" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/.laenv" -1 ...
#~
GG <- read.table("midterm.grades", sep = ",")
#~put ".Data.hp/GG" 825967085 ...
#~
GG
#~get ".Data.hp/GG" 825967085 ...
#~
?edit
#~put ".Data.hp/.Last.value" 825967151 ...
#~
?edit.data
#~put ".Data.hp/.Last.value" 825967156 ...
#~
help.start()
#~put ".Data.hp/.Last.value" 825967161 ...
#~
data.ed(GG)
#~get ".Data.hp/GG" 825967085 ...
#~put ".Data.hp/.Last.value" 825967210 ...
#~
?data.ed
#~put ".Data.hp/.Last.value" 825967213 ...
#~
data.ed(GG)
#~get ".Data.hp/GG" 825967085 ...
#~put ".Data.hp/.Last.value" 825967237 ...
#~
data.ed(GG, width = 15, heigh = 30)
#~get ".Data.hp/GG" 825967085 ...
#~put ".Data.hp/.Last.value" 825967295 ...
#~
rm(GG)
#~remove ".Data.hp/GG" 0 ...
#~
#~End session: Time: 825967302; Process: 9894
#~New session: Version: S Wed Mar 29 17:35:08 PST 1995, Time: Mon Mar 11 12:11:33 1996
#~get "/usr/statlocal/lib/Splus3.3/splus/.Datasets/..min.script.id" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Datasets/.Copyright" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Datasets/version" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/.laenv" -1 ...
#~
options()$prompt
#~put ".Data.hp/.Last.value" 826564294 ...
#~
search()
#~put ".Data.hp/.Last.value" 826564294 ...
#~
fac.design(c(2, 2, 2, 3), list(Temp = c("H", "L"), Pre = c("M", "X"), Brand = c(
	"X", "M"), Soft = c("Soft", "Med", "Hard")))
#~put ".Data.hp/.Last.value" 826564430 ...
#~
ries = fac.design(c(2,2,2,3),list(Temp=c("H","L"),Pre=c("M","X"),Brand=c("X","M"), So#~error: statement with syntax error
#~
ries <- fac.design(c(2, 2, 2, 3), list(Temp = c("H", "L"), Pre = c("M", "X"), 
	Brand = c("X", "M"), Soft = c("Soft", "Med", "Hard")))
#~put ".Data.hp/ries" 826564451 ...
#~
ries$Resp <- scan("ries-smith.dat")
#~get ".Data.hp/ries" 826564451 ...
#~put ".Data.hp/ries" 826564470 ...
#~
ries
#~get ".Data.hp/ries" 826564470 ...
#~
deviance(glm(Resp ~ Temp + Pre + Brand + Soft, family = poisson, data = ries))
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.variances" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.links" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/stat/.Functions/glm.weights" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/s/.Datasets/.Machine" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/stat/.Datasets/store.co.helmert" -1 ...
#~get ".Data.hp/ries" 826564470 ...
#~put ".Data.hp/.Last.value" 826564587 ...
#~
summary(glm(Resp ~ Temp + Pre + Brand + Soft, family = poisson, data = ries))
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.variances" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.links" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/stat/.Functions/glm.weights" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/s/.Datasets/.Machine" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/stat/.Datasets/store.co.helmert" -1 ...
#~get ".Data.hp/ries" 826564470 ...
#~put ".Data.hp/.Last.value" 826564610 ...
#~
(glm(Resp ~ Temp + Pre + Brand + Soft, family = poisson, data = ries))
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.variances" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.links" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/stat/.Functions/glm.weights" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/s/.Datasets/.Machine" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/stat/.Datasets/store.co.helmert" -1 ...
#~get ".Data.hp/ries" 826564470 ...
#~put ".Data.hp/.Last.value" 826564622 ...
#~
names(glm(Resp ~ Temp + Pre + Brand + Soft, family = poisson, data = ries))
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.variances" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.links" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/stat/.Functions/glm.weights" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/s/.Datasets/.Machine" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/stat/.Datasets/store.co.helmert" -1 ...
#~get ".Data.hp/ries" 826564470 ...
#~put ".Data.hp/.Last.value" 826564635 ...
#~
glm(Resp ~ Temp + Pre + Brand + Soft, family = poisson, data = ries)$
	df.residual
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.variances" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.links" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/stat/.Functions/glm.weights" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/s/.Datasets/.Machine" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/stat/.Datasets/store.co.helmert" -1 ...
#~get ".Data.hp/ries" 826564470 ...
#~put ".Data.hp/.Last.value" 826564656 ...
#~
summary.glm
#~
(glm(Resp ~ Temp + Pre + Brand + Soft, family = poisson, data = ries))
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.variances" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.links" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/stat/.Functions/glm.weights" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/s/.Datasets/.Machine" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/stat/.Datasets/store.co.helmert" -1 ...
#~get ".Data.hp/ries" 826564470 ...
#~put ".Data.hp/.Last.value" 826564693 ...
#~
summary(glm(Resp ~ Temp + Pre + Brand + Soft, family = poisson, data = ries))
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.variances" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.links" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/stat/.Functions/glm.weights" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/s/.Datasets/.Machine" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/stat/.Datasets/store.co.helmert" -1 ...
#~get ".Data.hp/ries" 826564470 ...
#~put ".Data.hp/.Last.value" 826564731 ...
#~
glm(Resp ~ Temp + Pre + Brand + Soft, family = poisson, data = ries)
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.variances" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.links" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/stat/.Functions/glm.weights" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/s/.Datasets/.Machine" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/stat/.Datasets/store.co.helmert" -1 ...
#~get ".Data.hp/ries" 826564470 ...
#~put ".Data.hp/.Last.value" 826564758 ...
#~
glm(Resp ~ Temp + Brand + Pre * Soft, family = poisson, data = ries)
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.variances" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.links" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/stat/.Functions/glm.weights" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/s/.Datasets/.Machine" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/stat/.Datasets/store.co.helmert" -1 ...
#~get ".Data.hp/ries" 826564470 ...
#~put ".Data.hp/.Last.value" 826564845 ...
#~
glm(Resp ~ Temp + Brand * Pre + Soft, family = poisson, data = ries)
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.variances" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.links" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/stat/.Functions/glm.weights" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/s/.Datasets/.Machine" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/stat/.Datasets/store.co.helmert" -1 ...
#~get ".Data.hp/ries" 826564470 ...
#~put ".Data.hp/.Last.value" 826564918 ...
#~
glm(Resp ~ Soft + Pre * Brand + Temp * Brand, family = poisson, data = ries)
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.variances" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.links" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/stat/.Functions/glm.weights" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/s/.Datasets/.Machine" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/stat/.Datasets/store.co.helmert" -1 ...
#~get ".Data.hp/ries" 826564470 ...
#~put ".Data.hp/.Last.value" 826565020 ...
#~
glm(Resp ~ Soft * Temp + Pre * Brand + Temp * Brand, family = poisson, data = 
	ries)
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.variances" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.links" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/stat/.Functions/glm.weights" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/s/.Datasets/.Machine" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/stat/.Datasets/store.co.helmert" -1 ...
#~get ".Data.hp/ries" 826564470 ...
#~put ".Data.hp/.Last.value" 826565098 ...
#~
glm(Resp ~ Soft * Temp + Pre * Brand * Temp, family = poisson, data = ries)
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.variances" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.links" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/stat/.Functions/glm.weights" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/s/.Datasets/.Machine" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/stat/.Datasets/store.co.helmert" -1 ...
#~get ".Data.hp/ries" 826564470 ...
#~put ".Data.hp/.Last.value" 826565218 ...
#~
glm(Resp ~ Pre * Soft * Temp + Pre * Brand * Temp, family = poisson, data = 
	ries)
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.variances" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.links" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/stat/.Functions/glm.weights" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/s/.Datasets/.Machine" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/stat/.Datasets/store.co.helmert" -1 ...
#~get ".Data.hp/ries" 826564470 ...
#~put ".Data.hp/.Last.value" 826565348 ...
#~
help("stepwise", pager = "cat", window = F)
#~put ".Data.hp/.Last.value" 826565564 ...
#~
help("step", pager = "cat", window = F)
#~put ".Data.hp/.Last.value" 826565597 ...
#~
mymod <- glm(Resp ~ Pre + Soft + Temp + Brand, family = poisson, data = ries)
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.variances" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.links" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/stat/.Functions/glm.weights" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/s/.Datasets/.Machine" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/stat/.Datasets/store.co.helmert" -1 ...
#~get ".Data.hp/ries" 826564470 ...
#~put ".Data.hp/mymod" 826565666 ...
#~
step(mymod)
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.variances" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.links" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/stat/.Functions/glm.weights" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/s/.Datasets/.Machine" -1 ...
#~get ".Data.hp/mymod" 826565666 ...
#~get ".Data.hp/ries" 826564470 ...
#~put ".Data.hp/.Last.value" 826565714 ...
#~
step(mymod, scope = list(upper =  ~ .^4))
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.variances" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.links" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/stat/.Functions/glm.weights" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/s/.Datasets/.Machine" -1 ...
#~get ".Data.hp/mymod" 826565666 ...
#~get ".Data.hp/ries" 826564470 ...
#~put ".Data.hp/.Last.value" 826565798 ...
#~
step(mymod, scope = list(upper =  ~ .^4, lower =  ~ .^1))
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.variances" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.links" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/stat/.Functions/glm.weights" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/s/.Datasets/.Machine" -1 ...
#~get ".Data.hp/mymod" 826565666 ...
#~get ".Data.hp/ries" 826564470 ...
#~put ".Data.hp/.Last.value" 826565878 ...
#~
mymod <- glm(Resp ~ ^4, family=poisson, data=ries)
#~error: statement with syntax error
#~
mymod <- glm(Resp ~ Soft * Pre * Temp * Brand, family = poisson, data = ries)
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.variances" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.links" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/stat/.Functions/glm.weights" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/s/.Datasets/.Machine" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/stat/.Datasets/store.co.helmert" -1 ...
#~get ".Data.hp/ries" 826564470 ...
#~put ".Data.hp/mymod" 826566149 ...
#~
estimates(mymod)
#~error: couldn't find function "estimates"
#~
mymod
#~get ".Data.hp/mymod" 826566149 ...
#~
coef(mymod)
#~get ".Data.hp/mymod" 826566149 ...
#~put ".Data.hp/.Last.value" 826566170 ...
#~
summary(mymod)
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.variances" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.links" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/stat/.Functions/glm.weights" -1 ...
#~get ".Data.hp/mymod" 826566149 ...
#~put ".Data.hp/.Last.value" 826566175 ...
#~
ries
#~get ".Data.hp/ries" 826564470 ...
#~
summary(mymod)
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.variances" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.links" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/stat/.Functions/glm.weights" -1 ...
#~get ".Data.hp/mymod" 826566149 ...
#~put ".Data.hp/.Last.value" 826566203 ...
#~
mymod <- glm(Resp ~ (Soft + Pre + Temp + Brand)^2, family = poisson, data = 
	ries)
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.variances" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.links" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/stat/.Functions/glm.weights" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/s/.Datasets/.Machine" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/stat/.Datasets/store.co.helmert" -1 ...
#~get ".Data.hp/ries" 826564470 ...
#~put ".Data.hp/mymod" 826566362 ...
#~
summary(mymod)
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.variances" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.links" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/stat/.Functions/glm.weights" -1 ...
#~get ".Data.hp/mymod" 826566362 ...
#~put ".Data.hp/.Last.value" 826566369 ...
#~
summary(mymod <- update(mymod, . ~ . - Pre:Temp - Soft:Brand))
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.variances" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.links" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/stat/.Functions/glm.weights" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/s/.Datasets/.Machine" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/stat/.Datasets/store.co.helmert" -1 ...
#~get ".Data.hp/mymod" 826566362 ...
#~put ".Data.hp/mymod" 826566485 ...
#~get ".Data.hp/ries" 826564470 ...
#~put ".Data.hp/.Last.value" 826566485 ...
#~
summary(mymod <- update(mymod, . ~ . - Pre:Soft))
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.variances" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.links" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/stat/.Functions/glm.weights" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/s/.Datasets/.Machine" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/stat/.Datasets/store.co.helmert" -1 ...
#~get ".Data.hp/mymod" 826566485 ...
#~put ".Data.hp/mymod" 826566576 ...
#~get ".Data.hp/ries" 826564470 ...
#~put ".Data.hp/.Last.value" 826566576 ...
#~
residuals(mymod)
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.variances" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.links" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/stat/.Functions/glm.weights" -1 ...
#~get ".Data.hp/mymod" 826566576 ...
#~put ".Data.hp/.Last.value" 826569569 ...
#~
cooksd
#~error: Object "cooksd" not found
#~
residuals.glm
#~
residuals(mymod, "deviance")
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.variances" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.links" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/stat/.Functions/glm.weights" -1 ...
#~get ".Data.hp/mymod" 826566576 ...
#~put ".Data.hp/.Last.value" 826569626 ...
#~
residuals(mymod, "pearson")
#~get ".Data.hp/mymod" 826566576 ...
#~put ".Data.hp/.Last.value" 826569634 ...
#~
motif()
#~get "/usr/statlocal/lib/Splus3.3/splus/.Datasets/.PostScript.Options" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Datasets/ps.preamble.std" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Datasets/ps.setcolor.rgb" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Datasets/ps.setfont.std" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Datasets/ps.paper.regions" -1 ...
#~put ".Data.hp/.Last.value" 826569657 ...
#~
pars(mfrow = c(2, 3))
#~error: couldn't find function "pars"
#~
plot(mymod)
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.variances" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.links" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/stat/.Functions/glm.weights" -1 ...
#~get ".Data.hp/mymod" 826566576 ...
#~put ".Data.hp/.Last.value" 826569672 ...
#~
par(mfrow = c(2, 3))
#~put ".Data.hp/.Last.value" 826569677 ...
#~
plot(mymod)
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.variances" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.links" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/stat/.Functions/glm.weights" -1 ...
#~get ".Data.hp/mymod" 826566576 ...
#~put ".Data.hp/.Last.value" 826569678 ...
#~
#~End session: Time: 826572109; Process: 5541
#~New session: Version: S Wed Mar 29 17:35:08 PST 1995, Time: Wed Mar 13 12:32:18 1996
#~get "/usr/statlocal/lib/Splus3.3/splus/.Datasets/..min.script.id" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Datasets/.Copyright" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Datasets/version" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/.laenv" -1 ...
#~
?glm
#~New session: Version: S Wed Mar 29 17:35:08 PST 1995, Time: Wed Mar 13 12:47:24 1996
#~get "/usr/statlocal/lib/Splus3.3/splus/.Datasets/..min.script.id" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Datasets/.Copyright" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Datasets/version" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/.laenv" -1 ...
#~
options()$prompt
#~put ".Data.hp/.Last.value" 826739245 ...
#~
search()
#~put ".Data.hp/.Last.value" 826739245 ...
#~
monkey <- fac.design(c(6, 6), list(Watch = c("R", "S", "T", "U", "V", "W"), 
	Display = c("R", "S", "T", "U", "V", "W")))
#~put ".Data.hp/monkey" 826739322 ...
#~
monkey
#~get ".Data.hp/monkey" 826739322 ...
#~
monkey$Resp <- scan("monkey.dat")
#~get ".Data.hp/monkey" 826739322 ...
#~put ".Data.hp/monkey" 826739376 ...
#~
monkey
#~get ".Data.hp/monkey" 826739376 ...
#~
monkey(monkey$Watch == "R")
#~error: couldn't find function "monkey"
#~
monkey(monkey$Watch == "R",  )
#~error: couldn't find function "monkey"
#~
monkey[monkey$Watch == "R",  ]
#~get ".Data.hp/monkey" 826739376 ...
#~put ".Data.hp/.Last.value" 826739487 ...
#~
monkey[monkey$Display == "R",  ]
#~get ".Data.hp/monkey" 826739376 ...
#~put ".Data.hp/.Last.value" 826739500 ...
#~
monkey[monkey$Display == "T",  ]
#~get ".Data.hp/monkey" 826739376 ...
#~put ".Data.hp/.Last.value" 826739518 ...
#~
attach(monkey, pos = 1)
#~get ".Data.hp/monkey" 826739376 ...
#~put "monkey$.Last.value" 0 ...
#~
search()
#~put "monkey$.Last.value" 0 ...
#~
names(monkey)
#~get ".Data.hp/monkey" 826739376 ...
#~put "monkey$.Last.value" 0 ...
#~
www <
)
#~error: statement with syntax error
#~
wei <- sign(Resp)
#~put "monkey$wei" 0 ...
#~
wei
#~
wei <- sign(Resp + 1)
#~put "monkey$wei" 0 ...
#~
monkey[Display == "T",  ]
#~get ".Data.hp/monkey" 826739376 ...
#~put "monkey$.Last.value" 0 ...
#~
monkey[Display == Watch,  ]
#~get ".Data.hp/monkey" 826739376 ...
#~put "monkey$.Last.value" 0 ...
#~
wei[Display == "T",  ] <- 0
#~error: No dim attribute for array subset
#~
wei[Display == "T"] <- 0
#~put "monkey$wei" 0 ...
#~
wei[Display == Watch] <- 0
#~put "monkey$wei" 0 ...
#~
wei
#~
mymod <- glm(Resp ~ Watch + Display, family = poisson, weight = wei)
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.variances" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.links" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/stat/.Functions/glm.weights" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/s/.Datasets/.Machine" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/stat/.Datasets/store.co.helmert" -1 ...
#~put "monkey$mymod" 0 ...
#~
summary(mymod)
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.variances" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.links" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/stat/.Functions/glm.weights" -1 ...
#~put "monkey$.Last.value" 0 ...
#~
anova(mymod)
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.variances" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.links" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/stat/.Functions/glm.weights" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/s/.Datasets/.Machine" -1 ...
#~put "monkey$.Last.value" 0 ...
#~
Pearson(mymod)
#~put "monkey$.Last.value" 0 ...
#~
Pearson
#~
?residuals.glm
#~New session: Version: S Wed Mar 29 17:35:08 PST 1995, Time: Wed Mar 13 13:00:31 1996
#~get "/usr/statlocal/lib/Splus3.3/splus/.Datasets/..min.script.id" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Datasets/.Copyright" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Datasets/version" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/.laenv" -1 ...
#~
options()$prompt
#~put ".Data.hp/.Last.value" 826740031 ...
#~
search()
#~put ".Data.hp/.Last.value" 826740032 ...
#~
attach(monkey, pos = 1)
#~get ".Data.hp/monkey" 826739376 ...
#~put "monkey$.Last.value" 0 ...
#~
search()
#~put "monkey$.Last.value" 0 ...
#~
names(monkey)
#~get ".Data.hp/monkey" 826739376 ...
#~put "monkey$.Last.value" 0 ...
#~
ls()
#~put "monkey$.Last.value" 0 ...
#~
wei <- sign(Resp + 1)
#~put "monkey$wei" 0 ...
#~
wei[Display == "T"] <- 0
#~put "monkey$wei" 0 ...
#~
wei[Display == Watch] <- 0
#~put "monkey$wei" 0 ...
#~
mymod <- glm(Resp ~ Watch + Display, family = poisson, weight = wei)
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.variances" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.links" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/stat/.Functions/glm.weights" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/s/.Datasets/.Machine" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/stat/.Datasets/store.co.helmert" -1 ...
#~put "monkey$mymod" 0 ...
#~
help(residuals.glm)
#~New session: Version: S Wed Mar 29 17:35:08 PST 1995, Time: Wed Mar 13 13:02:17 1996
#~get "/usr/statlocal/lib/Splus3.3/splus/.Datasets/..min.script.id" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Datasets/.Copyright" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Datasets/version" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/.laenv" -1 ...
#~
options()$prompt
#~put ".Data.hp/.Last.value" 826740138 ...
#~
search()
#~put ".Data.hp/.Last.value" 826740138 ...
#~
help("residuals", pager = "cat", window = F)
#~put ".Data.hp/.Last.value" 826740169 ...
#~
residuals
#~
glm.residuals
#~error: Object "glm.residuals" not found
#~
residuals.glm
#~
#~put ".Data.hp/.Last.value" 826740213 ...
#~
#~End session: Time: 826740214; Process: 12150
attach(monkey, pos = 1)
#~get ".Data.hp/monkey" 826739376 ...
#~put "monkey$.Last.value" 0 ...
#~
search()
#~put "monkey$.Last.value" 0 ...
#~
names(monkey)
#~get ".Data.hp/monkey" 826739376 ...
#~put "monkey$.Last.value" 0 ...
#~
wei=1
#~error: statement with syntax error
#~
wei <- rep(1, length(Resp))
#~put "monkey$wei" 0 ...
#~
wei[Display == "T"] <- 0
#~put "monkey$wei" 0 ...
#~
wei[Display == Watch] <- 0
#~put "monkey$wei" 0 ...
#~
mymod <- glm(Resp ~ Watch + Display, family = poisson, weight = wei)
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.variances" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.links" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/stat/.Functions/glm.weights" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/s/.Datasets/.Machine" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/stat/.Datasets/store.co.helmert" -1 ...
#~put "monkey$mymod" 0 ...
#~
anova(mymod)
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.variances" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.links" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/stat/.Functions/glm.weights" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/s/.Datasets/.Machine" -1 ...
#~put "monkey$.Last.value" 0 ...
#~
sum(residuals(mymod, type = "pearson")^2)
#~put "monkey$.Last.value" 0 ...
#~
residuals(mymod)
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.variances" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.links" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/stat/.Functions/glm.weights" -1 ...
#~put "monkey$.Last.value" 0 ...
#~
data.frame(Resp, fitted(mymod), residuals(mymod))
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.variances" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.links" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/stat/.Functions/glm.weights" -1 ...
#~put "monkey$.Last.value" 0 ...
#~
data.frame(Display, Watch, Resp, fitted(mymod), residuals(mymod))
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.variances" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.links" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/stat/.Functions/glm.weights" -1 ...
#~put "monkey$.Last.value" 0 ...
#~
#~End session: Time: 826741015; Process: 12242
#~New session: Version: S Wed Mar 29 17:35:08 PST 1995, Time: Mon Mar 18 12:59:00 1996
#~get "/usr/statlocal/lib/Splus3.3/splus/.Datasets/..min.script.id" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Datasets/.Copyright" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Datasets/version" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/.laenv" -1 ...
#~
help.start()
#~put ".Data.hp/.Last.value" 827171951 ...
#~
#~End session: Time: 827172076; Process: 23913
#~New session: Version: S Wed Mar 29 17:35:08 PST 1995, Time: Mon Mar 18 14:48:04 1996
#~get "/usr/statlocal/lib/Splus3.3/splus/.Datasets/..min.script.id" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Datasets/.Copyright" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Datasets/version" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/.laenv" -1 ...
#~
start.help()
#~error: couldn't find function "start.help"
#~
help.start()
#~put ".Data.hp/.Last.value" 827178498 ...
#~
dgamma((1:100)/100, 2)
#~put ".Data.hp/.Last.value" 827178597 ...
#~
motif()
#~get "/usr/statlocal/lib/Splus3.3/splus/.Datasets/.PostScript.Options" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Datasets/ps.preamble.std" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Datasets/ps.setcolor.rgb" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Datasets/ps.setfont.std" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Datasets/ps.paper.regions" -1 ...
#~put ".Data.hp/.Last.value" 827178604 ...
#~
xx <- (1:200)/200
#~put ".Data.hp/xx" 827178619 ...
#~
plot(xx, dgamma(xx, 2))
#~get ".Data.hp/xx" 827178619 ...
#~
xx <- (1:200)/20
#~put ".Data.hp/xx" 827178646 ...
#~
plot(xx, dgamma(xx, 2))
#~get ".Data.hp/xx" 827178646 ...
#~
points(xx, dgamma(xx, 3))
#~get ".Data.hp/xx" 827178646 ...
#~
points(xx, dgamma(xx, 4))
#~get ".Data.hp/xx" 827178646 ...
#~
points(xx, dgamma(xx, 1))
#~get ".Data.hp/xx" 827178646 ...
#~
points(xx, dgamma(xx, 0.5))
#~get ".Data.hp/xx" 827178646 ...
#~
plot(xx, dgamma(xx, 0.5))
#~get ".Data.hp/xx" 827178646 ...
#~
plot(xx, dgamma(xx, 0.5), lty = 1)
#~get ".Data.hp/xx" 827178646 ...
#~
plot(xx, dgamma(xx, 0.5), lty = 2)
#~get ".Data.hp/xx" 827178646 ...
#~
plot(xx, dgamma(xx, 0.5), lty = 3)
#~get ".Data.hp/xx" 827178646 ...
#~
plot(xx, dgamma(xx, 0.5))
#~get ".Data.hp/xx" 827178646 ...
#~
plot(xx, dgamma(xx, 0.5), lty = 3)
#~get ".Data.hp/xx" 827178646 ...
#~
plot(xx, dgamma(xx, 0.5), type = "l")
#~get ".Data.hp/xx" 827178646 ...
#~
lines(xx, dgamma(xx, 1))
#~get ".Data.hp/xx" 827178646 ...
#~
lines(xx, dgamma(xx, 1.1))
#~get ".Data.hp/xx" 827178646 ...
#~
plot(xx, dgamma(xx, 0.5), type = "l")
#~get ".Data.hp/xx" 827178646 ...
#~
lines(xx, dgamma(xx, 1))
#~get ".Data.hp/xx" 827178646 ...
#~
lines(xx, dgamma(xx, 2))
#~get ".Data.hp/xx" 827178646 ...
#~
lines(xx, dgamma(xx, 3))
#~get ".Data.hp/xx" 827178646 ...
#~
xx <- xx + 0.1
#~get ".Data.hp/xx" 827178646 ...
#~put ".Data.hp/xx" 827178957 ...
#~
plot(xx, dgamma(xx, 0.5), type = "l")
#~get ".Data.hp/xx" 827178957 ...
#~
lines(xx, dgamma(xx, 1))
#~get ".Data.hp/xx" 827178957 ...
#~
lines(xx, dgamma(xx, 2))
#~get ".Data.hp/xx" 827178957 ...
#~
lines(xx, dgamma(xx, 3))
#~get ".Data.hp/xx" 827178957 ...
#~
help.start()
#~put ".Data.hp/.Last.value" 827179045 ...
#~
#~End session: Time: 827179123; Process: 24114
#~New session: Version: S Wed Mar 29 17:35:08 PST 1995, Time: Mon Mar 18 15:01:33 1996
#~get "/usr/statlocal/lib/Splus3.3/splus/.Datasets/..min.script.id" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Datasets/.Copyright" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Datasets/version" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/.laenv" -1 ...
#~
motif()
#~get "/usr/statlocal/lib/Splus3.3/splus/.Datasets/.PostScript.Options" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Datasets/ps.preamble.std" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Datasets/ps.setcolor.rgb" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Datasets/ps.setfont.std" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Datasets/ps.paper.regions" -1 ...
#~put ".Data.hp/.Last.value" 827179296 ...
#~
#~End session: Time: 827179300; Process: 24324
#~New session: Version: S Wed Mar 29 17:35:08 PST 1995, Time: Mon Mar 18 15:02:00 1996
#~get "/usr/statlocal/lib/Splus3.3/splus/.Datasets/..min.script.id" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Datasets/.Copyright" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Datasets/version" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/.laenv" -1 ...
#~
motif()
#~get "/usr/statlocal/lib/Splus3.3/splus/.Datasets/.PostScript.Options" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Datasets/ps.preamble.std" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Datasets/ps.setcolor.rgb" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Datasets/ps.setfont.std" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Datasets/ps.paper.regions" -1 ...
#~put ".Data.hp/.Last.value" 827179325 ...
#~
#~End session: Time: 827179328; Process: 24349
#~New session: Version: S Wed Mar 29 17:35:08 PST 1995, Time: Mon Mar 18 15:03:08 1996
#~get "/usr/statlocal/lib/Splus3.3/splus/.Datasets/..min.script.id" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Datasets/.Copyright" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Datasets/version" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/.laenv" -1 ...
#~
motif()
#~get "/usr/statlocal/lib/Splus3.3/splus/.Datasets/.PostScript.Options" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Datasets/ps.preamble.std" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Datasets/ps.setcolor.rgb" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Datasets/ps.setfont.std" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Datasets/ps.paper.regions" -1 ...
#~put ".Data.hp/.Last.value" 827179394 ...
#~
#~End session: Time: 827179403; Process: 24373
#~New session: Version: S Wed Mar 29 17:35:08 PST 1995, Time: Mon Mar 18 15:04:11 1996
#~get "/usr/statlocal/lib/Splus3.3/splus/.Datasets/..min.script.id" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Datasets/.Copyright" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Datasets/version" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/.laenv" -1 ...
#~
motif()
#~get "/usr/statlocal/lib/Splus3.3/splus/.Datasets/.PostScript.Options" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Datasets/ps.preamble.std" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Datasets/ps.setcolor.rgb" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Datasets/ps.setfont.std" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Datasets/ps.paper.regions" -1 ...
#~put ".Data.hp/.Last.value" 827179457 ...
#~
#~End session: Time: 827179458; Process: 24400
#~New session: Version: S Wed Mar 29 17:35:08 PST 1995, Time: Mon Mar 18 15:25:47 1996
#~get "/usr/statlocal/lib/Splus3.3/splus/.Datasets/..min.script.id" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Datasets/.Copyright" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Datasets/version" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/.laenv" -1 ...
#~
motif()
#~get "/usr/statlocal/lib/Splus3.3/splus/.Datasets/.PostScript.Options" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Datasets/ps.preamble.std" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Datasets/ps.setcolor.rgb" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Datasets/ps.setfont.std" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Datasets/ps.paper.regions" -1 ...
#~put ".Data.hp/.Last.value" 827180755 ...
#~
plot(xx, dgamma(xx, 0.5), type = "l")
#~get ".Data.hp/xx" 827178957 ...
#~
lines(xx, dgamma(xx, 1))
#~get ".Data.hp/xx" 827178957 ...
#~
lines(xx, dgamma(xx, 2))
#~get ".Data.hp/xx" 827178957 ...
#~
lines(xx, dgamma(xx, 3))
#~get ".Data.hp/xx" 827178957 ...
#~
dexp
#~
plot(xx, dexp(xx, 0.5), type = "l")
#~get ".Data.hp/xx" 827178957 ...
#~
#~End session: Time: 827182163; Process: 24616
#~New session: Version: S Wed Mar 29 17:35:08 PST 1995, Time: Wed Mar 20 13:04:21 1996
#~get "/usr/statlocal/lib/Splus3.3/splus/.Datasets/..min.script.id" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Datasets/.Copyright" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Datasets/version" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/.laenv" -1 ...
#~
help.start()
#~put ".Data.hp/.Last.value" 827345070 ...
#~
#~End session: Time: 827345643; Process: 26298
#~New session: Version: S Wed Mar 29 17:35:08 PST 1995, Time: Wed Mar 20 13:23:53 1996
#~get "/usr/statlocal/lib/Splus3.3/splus/.Datasets/..min.script.id" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Datasets/.Copyright" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Datasets/version" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/.laenv" -1 ...
#~
help(wilcox.test)
#~put ".Data.hp/.Last.value" 827346242 ...
#~
#~End session: Time: 827346245; Process: 26412
#~New session: Version: S Wed Mar 29 17:35:08 PST 1995, Time: Mon Apr  1 21:29:53 1996
#~get "/usr/statlocal/lib/Splus3.3/splus/.Datasets/..min.script.id" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Datasets/.Copyright" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Datasets/version" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/.laenv" -1 ...
#~
leuk
#~get ".Data.hp/leuk" 828412180 ...
#~
attach(leuk)
#~get ".Data.hp/leuk" 828412180 ...
#~put ".Data.hp/.Last.value" 828412214 ...
#~
motif()
#~get "/usr/statlocal/lib/Splus3.3/splus/.Datasets/.PostScript.Options" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Datasets/ps.preamble.std" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Datasets/ps.setcolor.rgb" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Datasets/ps.setfont.std" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Datasets/ps.paper.regions" -1 ...
#~put ".Data.hp/.Last.value" 828412220 ...
#~
plot(survfit(Surv(time) ~ ag), lty = c(2, 3))
#~put ".Data.hp/.Last.value" 828412246 ...
#~
#~End session: Time: 828412304; Process: 384
#~New session: Version: S Wed Mar 29 17:35:08 PST 1995, Time: Mon Apr  1 21:34:49 1996
#~get "/usr/statlocal/lib/Splus3.3/splus/.Datasets/..min.script.id" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Datasets/.Copyright" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Datasets/version" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/.laenv" -1 ...
#~
leuk
#~get ".Data.hp/leuk" 828412180 ...
#~
motif()
#~get "/usr/statlocal/lib/Splus3.3/splus/.Datasets/.PostScript.Options" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Datasets/ps.preamble.std" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Datasets/ps.setcolor.rgb" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Datasets/ps.setfont.std" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Datasets/ps.paper.regions" -1 ...
#~put ".Data.hp/.Last.value" 828412939 ...
#~
attach(leuk)
#~get ".Data.hp/leuk" 828412180 ...
#~put ".Data.hp/.Last.value" 828412957 ...
#~
plot(survfit(Surv(time) ~ ag), lty = c(2, 3))
#~put ".Data.hp/.Last.value" 828412971 ...
#~
help(survfit)
#~put ".Data.hp/.Last.value" 828412986 ...
#~
legend(80, 0.8, c("ag absent", "ag present"), lty = c(2, 3))
#~
lwbc <- log(wbc)
#~put ".Data.hp/lwbc" 828413370 ...
#~
leuk.glm <- glm(time ~ ag * lwbc, Gamma(log))
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.variances" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.links" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/stat/.Functions/glm.weights" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/s/.Datasets/.Machine" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/stat/.Datasets/store.co.helmert" -1 ...
#~get ".Data.hp/lwbc" 828413370 ...
#~put ".Data.hp/leuk.glm" 828413380 ...
#~
summary(leuk.glm, dispersion = 1)	# force exponential data summary
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.variances" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.links" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/stat/.Functions/glm.weights" -1 ...
#~get ".Data.hp/leuk.glm" 828413380 ...
#~put ".Data.hp/.Last.value" 828413395 ...
#~
anova(leuk.glm)
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.variances" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.links" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/stat/.Functions/glm.weights" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/s/.Datasets/.Machine" -1 ...
#~get ".Data.hp/leuk.glm" 828413380 ...
#~get ".Data.hp/lwbc" 828413370 ...
#~put ".Data.hp/.Last.value" 828413405 ...
#~
leuk.glm <- update(leuk.glm,  ~ . - ag:lwbc)
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.variances" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.links" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/stat/.Functions/glm.weights" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/s/.Datasets/.Machine" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/stat/.Datasets/store.co.helmert" -1 ...
#~get ".Data.hp/leuk.glm" 828413380 ...
#~put ".Data.hp/leuk.glm" 828413418 ...
#~get ".Data.hp/lwbc" 828413370 ...
#~
summary(leuk.glm, dispersion = 1)
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.variances" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.links" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/stat/.Functions/glm.weights" -1 ...
#~get ".Data.hp/leuk.glm" 828413418 ...
#~put ".Data.hp/.Last.value" 828413426 ...
#~
leuk.glmi <- glm(time ~ ag * wbc, Gamma(inverse))
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.variances" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.links" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/stat/.Functions/glm.weights" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/stat/.Datasets/store.co.helmert" -1 ...
#~put ".Data.hp/leuk.glmi" 828413437 ...
#~
summary(leuk.glmi, dispersion = 1)
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.variances" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.links" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/stat/.Functions/glm.weights" -1 ...
#~get ".Data.hp/leuk.glmi" 828413437 ...
#~put ".Data.hp/.Last.value" 828413444 ...
#~
anova(leuk.glmi)
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.variances" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.links" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/stat/.Functions/glm.weights" -1 ...
#~get ".Data.hp/leuk.glmi" 828413437 ...
#~put ".Data.hp/.Last.value" 828413452 ...
#~
help(glm)
#~put ".Data.hp/.Last.value" 828413478 ...
#~
help(survreg)
#~put ".Data.hp/.Last.value" 828413500 ...
#~
plot(log(time) ~ ag + log(wbc))
#~put ".Data.hp/.Last.value" 828413953 ...
#~
pars(mfrow = c(2, 1))
#~error: couldn't find function "pars"
#~
par(mfrow = c(2, 1))
#~put ".Data.hp/.Last.value" 828413974 ...
#~
plot(log(time) ~ ag + log(wbc))
#~put ".Data.hp/.Last.value" 828413977 ...
#~
help(plot)
#~put ".Data.hp/.Last.value" 828418317 ...
#~
help(plot.formula)
#~put ".Data.hp/.Last.value" 828418347 ...
#~
plot.formula
#~
help(plot.data.frame)
#~put ".Data.hp/.Last.value" 828418386 ...
#~
help(plot.default)
#~put ".Data.hp/.Last.value" 828418415 ...
#~
par(mfrow = c(1, 1))
#~put ".Data.hp/.Last.value" 828418643 ...
#~
plot(log(time) ~ ag + log(wbc), mfrow = c(2, 1))
#~put ".Data.hp/last.dump" 828418652 ...
#~
help(plot)
#~put ".Data.hp/.Last.value" 828418677 ...
#~
man(par)
#~error: couldn't find function "man"
#~
help(par)
#~put ".Data.hp/.Last.value" 828418734 ...
#~
plot(log(time) ~ ag + log(wbc), mfrow = c(2, 1))
#~put ".Data.hp/last.dump" 828418744 ...
#~
plot(log(time) ~ ag + log(wbc), mfrow = c(2, 1))
#~put ".Data.hp/last.dump" 828418756 ...
#~
par(mfrow = c(2, 1))
#~put ".Data.hp/.Last.value" 828418765 ...
#~
plot(log(time) ~ ag + log(wbc))
#~put ".Data.hp/.Last.value" 828418769 ...
#~
source("/tmp/dodo")
lwbc <- log(wbc)
       par(mfrow=c(2,1))
       plot(log(time) ~ ag + log(wbc))
       options(contrasts=c("contr.treatment","contr.poly")
       leuk.glm <- glm(time ~ ag*lwbc, Gamma(log))
#~error: statement with syntax error
#~put ".Data.hp/last.dump" 828419004 ...
#~
source("/tmp/dodo")
lwbc <- log(wbc)
       par(mfrow=c(2,1))
       plot(log(time) ~ ag + log(wbc))
       options(contrasts=c("contr.treatment","contr.poly")
       leuk.glm <- glm(time ~ ag*lwbc, Gamma(log))
#~error: statement with syntax error
#~put ".Data.hp/last.dump" 828419028 ...
#~
ls()
#~put ".Data.hp/.Last.value" 828419048 ...
#~
rm(leuk.glm)
#~remove ".Data.hp/leuk.glm" 0 ...
#~
rm(leuk.glmi)
#~remove ".Data.hp/leuk.glmi" 0 ...
#~
source("/tmp/dodo")
lwbc <- log(wbc)
       par(mfrow=c(2,1))
       plot(log(time) ~ ag + log(wbc))
       options(contrasts=c("contr.treatment","contr.poly")
       leuk.glm <- glm(time ~ ag*lwbc, Gamma(log))
#~error: statement with syntax error
#~put ".Data.hp/last.dump" 828419108 ...
#~
ls()
#~put ".Data.hp/.Last.value" 828419111 ...
#~
lwbc <- log(wbc)
#~put ".Data.hp/lwbc" 828419120 ...
#~
plot(log(time) ~ ag + log(wbc))
#~put ".Data.hp/.Last.value" 828419131 ...
#~
options(contrasts = c("contr.treatment", "contr.poly"))
#~put ".Data.hp/.Last.value" 828419149 ...
#~
source("/tmp/dodo")
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.variances" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.links" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/stat/.Functions/glm.weights" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/s/.Datasets/.Machine" -1 ...
#~put ".Data.hp/lwbc" 828419162 ...
#~put ".Data.hp/leuk.glm" 828419162 ...
#~put ".Data.hp/leuk.glmi" 828419162 ...
#~put ".Data.hp/.Last.value" 828419162 ...
#~
source("/tmp/dodo")
lwbc <- log(wbc)
       par(mfrow=c(2,1))
       plot(log(time) ~ ag + log(wbc))
       options(contrasts=c("contr.treatment","contr.poly"))
       leuk.glm <- glm(time ~ ag*lwbc, Gamma(log))
       print(summary(leuk.glm, dispersion=1) # force exponential data summary)
       print(anova(leuk.glm))
#~error: statement with syntax error
#~put ".Data.hp/last.dump" 828419194 ...
#~
source("/tmp/dodo")
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.variances" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.links" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/stat/.Functions/glm.weights" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/s/.Datasets/.Machine" -1 ...
#~put ".Data.hp/lwbc" 828419209 ...
#~put ".Data.hp/leuk.glm" 828419209 ...
#~put ".Data.hp/leuk.glmi" 828419210 ...
#~put ".Data.hp/.Last.value" 828419210 ...
#~
source("/tmp/dodo")
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.variances" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.links" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/stat/.Functions/glm.weights" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/s/.Datasets/.Machine" -1 ...
#~put ".Data.hp/lwbc" 828419367 ...
#~put ".Data.hp/leuk.glm" 828419368 ...
#~put ".Data.hp/leuk.glmi" 828419368 ...
#~put ".Data.hp/.Last.value" 828419368 ...
#~
#~End session: Time: 828419567; Process: 421
#~New session: Version: S Wed Mar 29 17:35:08 PST 1995, Time: Wed Apr  3 09:22:38 1996
#~get "/usr/statlocal/lib/Splus3.3/splus/.Datasets/..min.script.id" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Datasets/.Copyright" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Datasets/version" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/.laenv" -1 ...
#~
chal <- read.table("chal.dat")
#~put ".Data.hp/chal" 828541370 ...
#~
attach(chal, pos = 1)
#~get ".Data.hp/chal" 828541370 ...
#~put "chal$.Last.value" 0 ...
#~
ls
#~
ls()
#~put "chal$.Last.value" 0 ...
#~
detach(1, save = F)
#~put ".Data.hp/.Last.value" 828541391 ...
#~
chal <- read.table("chal.dat", head = T)
#~put ".Data.hp/chal" 828541397 ...
#~
attach(chal, pos = 1)
#~get ".Data.hp/chal" 828541397 ...
#~put "chal$.Last.value" 0 ...
#~
ls()
#~put "chal$.Last.value" 0 ...
#~
rep(1, 6)
#~put "chal$.Last.value" 0 ...
#~
n <- rep(6, length(Temp))
#~put "chal$n" 0 ...
#~
mymod <- glm(cbind(Fail, Fail - n) ~ Temp, family = binomial)
#~error: subroutine glmfit: 23 missing value(s) in argument 5
#~put "chal$last.dump" 0 ...
#~
warnings()
#~
Temp
#~
Fail
#~
n
#~
mymod <- glm(cbind(Fail, n - Fail) ~ Temp, family = binomial)
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.variances" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.links" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/stat/.Functions/glm.weights" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/s/.Datasets/.Machine" -1 ...
#~put "chal$mymod" 0 ...
#~
summary(mymod)
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.variances" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.links" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/stat/.Functions/glm.weights" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/s/.Datasets/.Machine" -1 ...
#~put "chal$.Last.value" 0 ...
#~
pred
#~error: Object "pred" not found
#~
predicted
#~error: Object "predicted" not found
#~
predict
#~
help(predict)
#~put "chal$.Last.value" 0 ...
#~
predict(mymod, 31)
#~error: Invalid frame number, 31
#~put "chal$last.dump" 0 ...
#~
predict(mymod)
#~put "chal$.Last.value" 0 ...
#~
help(predict.glm)
#~put "chal$.Last.value" 0 ...
#~
predict(mymod, se.fit = T)
#~put "chal$.Last.value" 0 ...
#~
fitted(mymod)
#~put "chal$.Last.value" 0 ...
#~
new <- data.frame(Fail = c(30, 31, 32, 33, 34))
#~put "chal$new" 0 ...
#~
predict(mymod, new)
#~error: Length of variable 1 is 23 != length of row names (5)
#~put "chal$last.dump" 0 ...
#~
new
#~
ls()
#~put "chal$.Last.value" 0 ...
#~
help.start()
#~put "chal$.Last.value" 0 ...
#~
new
#~
Fail
#~
new <- data.frame(Temp = c(30, 31, 32, 33, 34))
#~put "chal$new" 0 ...
#~
predict(mymod, new)
#~put "chal$.Last.value" 0 ...
#~
fitted(mymod, new)
#~put "chal$last.dump" 0 ...
#~
unlogit <- function(x){ exp(x) / (1 + exp(x)}
#~error: statement with syntax error
#~
exp(1)
#~put "chal$.Last.value" 0 ...
#~
unlogit <- function(x)
{
	exp(x)/(1 + exp(x))
}
#~put "chal$unlogit" 0 ...
#~
unlogit(predict(mymod, new))
#~put "chal$.Last.value" 0 ...
#~
fitted(mymod)
#~put "chal$.Last.value" 0 ...
#~
mymod <- glm(cbind(Fail, n - Fail) ~ Temp + Temp^2, family = binomial)
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.variances" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.links" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/stat/.Functions/glm.weights" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/s/.Datasets/.Machine" -1 ...
#~put "chal$mymod" 0 ...
#~
summary(mymod)
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.variances" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.links" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/stat/.Functions/glm.weights" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/s/.Datasets/.Machine" -1 ...
#~put "chal$.Last.value" 0 ...
#~
unlogit(predict(mymod, new))
#~put "chal$.Last.value" 0 ...
#~
mymod <- glm(cbind(Fail, n - Fail) ~ Temp, family = binomial)
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.variances" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.links" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/stat/.Functions/glm.weights" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/s/.Datasets/.Machine" -1 ...
#~put "chal$mymod" 0 ...
#~
summary(mymod)
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.variances" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.links" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/stat/.Functions/glm.weights" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/s/.Datasets/.Machine" -1 ...
#~put "chal$.Last.value" 0 ...
#~
ls()
#~put "chal$.Last.value" 0 ...
#~
motif()
#~get "/usr/statlocal/lib/Splus3.3/splus/.Datasets/.PostScript.Options" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Datasets/ps.preamble.std" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Datasets/ps.setcolor.rgb" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Datasets/ps.setfont.std" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Datasets/ps.paper.regions" -1 ...
#~put "chal$.Last.value" 0 ...
#~
plot(Temp[Fail > 0], Fail[Fail > 0])
#~
plot(Temp, Fail)
#~
detach(1, save = F)
#~put ".Data.hp/.Last.value" 828544075 ...
#~
chal <- read.table("chal.dat", head = T)
#~put ".Data.hp/chal" 828544091 ...
#~
attach(chal, pos = 1)
#~get ".Data.hp/chal" 828544091 ...
#~put "chal$.Last.value" 0 ...
#~
attach(chal, pos=1)ls()
#~error: statement with syntax error
#~
ls()
#~put "chal$.Last.value" 0 ...
#~
n <- rep(6, length(Fail))
#~put "chal$n" 0 ...
#~
plot(Temp[Fail > 0], Fail[Fail > 0])
#~
plot(Temp, Fail)
#~
ls()
#~put "chal$.Last.value" 0 ...
#~
FN <- cbind(Fail, n - Fail)
#~put "chal$FN" 0 ...
#~
mymod <- glm(FN ~ Temp, family = binomial)
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.variances" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.links" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/stat/.Functions/glm.weights" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/s/.Datasets/.Machine" -1 ...
#~put "chal$mymod" 0 ...
#~
mymod <- glm(cbind(Fail, n - Fail) ~ Temp, family = binomial)
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.variances" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.links" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/stat/.Functions/glm.weights" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/s/.Datasets/.Machine" -1 ...
#~put "chal$mymod" 0 ...
#~
mymod <- glm(FN ~ Temp, family = binomial)
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.variances" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.links" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/stat/.Functions/glm.weights" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/s/.Datasets/.Machine" -1 ...
#~put "chal$mymod" 0 ...
#~
summary(mymod)
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.variances" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.links" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/stat/.Functions/glm.weights" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/s/.Datasets/.Machine" -1 ...
#~put "chal$.Last.value" 0 ...
#~
anova(mymod)
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.variances" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.links" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/stat/.Functions/glm.weights" -1 ...
#~put "chal$.Last.value" 0 ...
#~
fitted(mymod)
#~put "chal$.Last.value" 0 ...
#~
predicted(mymod)
#~error: couldn't find function "predicted"
#~
predict(mymod)
#~put "chal$.Last.value" 0 ...
#~
unlogit(predict(mymod))
#~error: couldn't find function "unlogit"
#~
search()
#~put "chal$.Last.value" 0 ...
#~
ls()
#~put "chal$.Last.value" 0 ...
#~
help(assign)
#~put "chal$.Last.value" 0 ...
#~
dodo <- function(x)
exp(x)/(1 + exp(x))
#~put "chal$dodo" 0 ...
#~
assign(unlogit, dodo, where = 2)
#~error: Object "unlogit" not found
#~put "chal$last.dump" 0 ...
#~
assign(dodo, "unlogit", where = 2)
#~error: object of mode function and length 2 used as string value
#~put "chal$last.dump" 0 ...
#~
help(assign)
#~put "chal$.Last.value" 0 ...
#~
assign("unlogit", dodo, where = 2)
#~put ".Data.hp/unlogit" 828554289 ...
#~
rm(dodo)
#~remove "chal$dodo" 0 ...
#~
unlogit(predict(mymod))
#~put "chal$.Last.value" 0 ...
#~
fitted(mymod)
#~put "chal$.Last.value" 0 ...
#~
predict(mymod)
#~put "chal$.Last.value" 0 ...
#~
fitted(mymod) - unlogit(predict(mymod))
#~put "chal$.Last.value" 0 ...
#~
unlogit
#~
predict(mymod, se = T)
#~put "chal$.Last.value" 0 ...
#~
unlogit(-1.041893)
#~put "chal$.Last.value" 0 ...
#~
unlogit(-1.041893 + 2 * 0.7682786)
#~put "chal$.Last.value" 0 ...
#~
unlogit(-1.041893 - 2 * 0.7682786)
#~put "chal$.Last.value" 0 ...
#~
predict(mymod, data.frame(Temp = c(30, 31, 32, 33), se = T))
#~error: arguments imply differing number of rows: 4, 1
#~put "chal$last.dump" 0 ...
#~
predict(mymod, data.frame(Temp = c(30, 31, 32, 33)))
#~put "chal$.Last.value" 0 ...
#~
predict(mymod, data.frame(Temp = c(30, 31, 32, 33)), se = T)
#~put "chal$.Last.value" 0 ...
#~
unlogit(1.616887)
#~put "chal$.Last.value" 0 ...
#~
unlogit(1.616887 + 2 * 1.972255)
#~put "chal$.Last.value" 0 ...
#~
unlogit(1.616887 - 2 * 1.972255)
#~put "chal$.Last.value" 0 ...
#~
round(unlogit(1.616887 - 2 * 1.972255), 2)
#~put "chal$.Last.value" 0 ...
#~
round(unlogit(1.616887 + 2 * 1.972255), 2)
#~put "chal$.Last.value" 0 ...
#~
round(unlogit(1.616887 + 2 * 1.972255), 3)
#~put "chal$.Last.value" 0 ...
#~
round(unlogit(1.616887 - 2 * 1.972255), 3)
#~put "chal$.Last.value" 0 ...
#~
plot*mymod)
#~error: statement with syntax error
#~
plot(mymod)
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.variances" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.links" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/stat/.Functions/glm.weights" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/s/.Datasets/.Machine" -1 ...
#~put "chal$.Last.value" 0 ...
#~
plot(mymod, mfrow = c(2, 2))
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.variances" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.links" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/stat/.Functions/glm.weights" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/s/.Datasets/.Machine" -1 ...
#~put "chal$.Last.value" 0 ...
#~
plo
#~error: Object "plo" not found
#~
plot
#~
plot.glm
#~
par(mfrow = c(2, 3))
#~put "chal$.Last.value" 0 ...
#~
plot(mymod)
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.variances" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.links" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/stat/.Functions/glm.weights" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/s/.Datasets/.Machine" -1 ...
#~put "chal$.Last.value" 0 ...
#~
par(mfrow = c(2, 2))
#~put "chal$.Last.value" 0 ...
#~
plot(mymod)
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.variances" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.links" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/stat/.Functions/glm.weights" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/s/.Datasets/.Machine" -1 ...
#~put "chal$.Last.value" 0 ...
#~
mymod <- update(mymod, . ~ . + Temp^2)
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.variances" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.links" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/stat/.Functions/glm.weights" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/s/.Datasets/.Machine" -1 ...
#~put "chal$mymod" 0 ...
#~
plot(mymod)
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.variances" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.links" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/stat/.Functions/glm.weights" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/s/.Datasets/.Machine" -1 ...
#~put "chal$.Last.value" 0 ...
#~
mymod <- update(mymod, . ~ . - Temp^2)
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.variances" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.links" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/stat/.Functions/glm.weights" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/s/.Datasets/.Machine" -1 ...
#~put "chal$mymod" 0 ...
#~
plot(update(mymod, . ~ . + sqrt(Temp)))
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.variances" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.links" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/stat/.Functions/glm.weights" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/s/.Datasets/.Machine" -1 ...
#~put "chal$.Last.value" 0 ...
#~
summary(mymod)
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.variances" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.links" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/stat/.Functions/glm.weights" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/s/.Datasets/.Machine" -1 ...
#~put "chal$.Last.value" 0 ...
#~
mymod <- update(mymod, . ~ . - Temp^2)
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.variances" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.links" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/stat/.Functions/glm.weights" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/s/.Datasets/.Machine" -1 ...
#~put "chal$mymod" 0 ...
#~
summary(mymod)
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.variances" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.links" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/stat/.Functions/glm.weights" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/s/.Datasets/.Machine" -1 ...
#~put "chal$.Last.value" 0 ...
#~
mymod <- glm(FN ~ Temp, family = binomial)
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.variances" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.links" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/stat/.Functions/glm.weights" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/s/.Datasets/.Machine" -1 ...
#~put "chal$mymod" 0 ...
#~
plot(update(mymod, . ~ . + sqrt(Temp)))
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.variances" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.links" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/stat/.Functions/glm.weights" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/s/.Datasets/.Machine" -1 ...
#~put "chal$.Last.value" 0 ...
#~
plot(update(mymod, . ~ . + log(Temp)))
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.variances" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.links" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/stat/.Functions/glm.weights" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/s/.Datasets/.Machine" -1 ...
#~put "chal$.Last.value" 0 ...
#~
plot(update(mymod, . ~ . + exp(Temp)))
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.variances" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.links" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/stat/.Functions/glm.weights" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/s/.Datasets/.Machine" -1 ...
#~put "chal$.Last.value" 0 ...
#~
FN
#~
sign(Fail)
#~put "chal$.Last.value" 0 ...
#~
Bn <- cbind(sign(Fail), 1 - sign(Fail))
#~put "chal$Bn" 0 ...
#~
Bn
#~
mymod <- glm(Bn ~ Temp, family = binomial)
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.variances" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.links" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/stat/.Functions/glm.weights" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/s/.Datasets/.Machine" -1 ...
#~put "chal$mymod" 0 ...
#~
summary(mymod)
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.variances" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.links" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/stat/.Functions/glm.weights" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/s/.Datasets/.Machine" -1 ...
#~put "chal$.Last.value" 0 ...
#~
anova(mymod)
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.variances" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.links" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/stat/.Functions/glm.weights" -1 ...
#~put "chal$.Last.value" 0 ...
#~
plot(mymod)
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.variances" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.links" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/stat/.Functions/glm.weights" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/s/.Datasets/.Machine" -1 ...
#~put "chal$.Last.value" 0 ...
#~
detach(1, save = F)
#~put ".Data.hp/.Last.value" 828557707 ...
#~
stroke <- read.table("stroke.data", head = T)
#~error: Can't open file "stroke.data" for reading
#~put ".Data.hp/last.dump" 828557724 ...
#~
stroke <- read.table("stroke.dat", head = T)
#~put ".Data.hp/stroke" 828557729 ...
#~
stroke[1:10,  ]
#~get ".Data.hp/stroke" 828557729 ...
#~put ".Data.hp/.Last.value" 828557737 ...
#~
mymod <- glm(Status ~ Side + Gender + Age + LOSD + WAB1, family = binom, data
	 = stroke)
#~error: Object "binom" not found
#~put ".Data.hp/last.dump" 828557828 ...
#~
mymod <- glm(Status ~ Side + Gender + Age + LOSD + WAB1, family = binomial, 
	data = stroke)
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.variances" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.links" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/stat/.Functions/glm.weights" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/s/.Datasets/.Machine" -1 ...
#~get ".Data.hp/stroke" 828557729 ...
#~put ".Data.hp/mymod" 828557838 ...
#~
summary(mymod)
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.variances" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.links" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/stat/.Functions/glm.weights" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/s/.Datasets/.Machine" -1 ...
#~get ".Data.hp/mymod" 828557838 ...
#~put ".Data.hp/.Last.value" 828557860 ...
#~
anova(mymod)
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.variances" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.links" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/stat/.Functions/glm.weights" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/s/.Datasets/.Machine" -1 ...
#~get ".Data.hp/mymod" 828557838 ...
#~get ".Data.hp/stroke" 828557729 ...
#~put ".Data.hp/.Last.value" 828557873 ...
#~
plot(mymod)
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.variances" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.links" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/stat/.Functions/glm.weights" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/s/.Datasets/.Machine" -1 ...
#~get ".Data.hp/mymod" 828557838 ...
#~put ".Data.hp/.Last.value" 828558205 ...
#~
q()
#~End session: Time: 828558702; Process: 813
#~New session: Version: S Wed Mar 29 17:35:08 PST 1995, Time: Mon Apr  8 13:13:27 1996
#~get "/usr/statlocal/lib/Splus3.3/splus/.Datasets/..min.script.id" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Datasets/.Copyright" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Datasets/version" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/.laenv" -1 ...
#~
ls()
#~put ".Data.hp/.Last.value" 828983612 ...
#~
stroke
#~get ".Data.hp/stroke" 828557729 ...
#~
attach(stroke, pos = 1)
#~get ".Data.hp/stroke" 828557729 ...
#~put "stroke$.Last.value" 0 ...
#~
mymod <- glm(Status ~ Side, family = binomial)
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.variances" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.links" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/stat/.Functions/glm.weights" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/s/.Datasets/.Machine" -1 ...
#~put "stroke$mymod" 0 ...
#~
summary(mymod)
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.variances" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.links" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/stat/.Functions/glm.weights" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/s/.Datasets/.Machine" -1 ...
#~put "stroke$.Last.value" 0 ...
#~
motif()
#~get "/usr/statlocal/lib/Splus3.3/splus/.Datasets/.PostScript.Options" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Datasets/ps.preamble.std" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Datasets/ps.setcolor.rgb" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Datasets/ps.setfont.std" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Datasets/ps.paper.regions" -1 ...
#~put "stroke$.Last.value" 0 ...
#~
plot(stroke)
#~get ".Data.hp/stroke" 828557729 ...
#~put "stroke$.Last.value" 0 ...
#~
scatmat(stroke)
#~error: couldn't find function "scatmat"
#~
scatter(stroke)
#~error: couldn't find function "scatter"
#~
pairs(stroke)
#~get ".Data.hp/stroke" 828557729 ...
#~put "stroke$.Last.value" 0 ...
#~
plot(WAB1, Status)
#~
help.start()
#~put "stroke$.Last.value" 0 ...
#~
scatter.smootht(WAB1, Status)
#~error: couldn't find function "scatter.smootht"
#~
scatter.smooth(WAB1, Status)
#~put "stroke$.Last.value" 0 ...
#~
scatter.smooth(WAB1, Status, span = 0.3)
#~put "stroke$.Last.value" 0 ...
#~
scatter.smooth(WAB1, Status, span = 0.5)
#~put "stroke$.Last.value" 0 ...
#~
scatter.smooth(WAB1, Status, span = 0.4)
#~error: all data on boundary of neighborhood. make span bigger
#~put "stroke$last.dump" 0 ...
#~
scatter.smooth(WAB1, Status, span = 0.3)
#~put "stroke$.Last.value" 0 ...
#~
scatter.smooth(WAB1, Status, span = 0.2)
#~put "stroke$.Last.value" 0 ...
#~
scatter.smooth(WAB1, Status, span = 0.1)
#~error: subroutine loess_raw: 61 missing value(s) in argument 1
#~put "stroke$last.dump" 0 ...
#~
scatter.smooth(WAB1, Status, span = 0.2)
#~put "stroke$.Last.value" 0 ...
#~
scatter.smooth(WAB1, Status, span = 0.4)
#~error: all data on boundary of neighborhood. make span bigger
#~put "stroke$last.dump" 0 ...
#~
scatter.smooth(WAB1, Status, span = 0.35)
#~put "stroke$.Last.value" 0 ...
#~
scatter.smooth(WAB1, Status, span = 0.5)
#~put "stroke$.Last.value" 0 ...
#~
scatter.smooth(WAB1, Status, span = 0.65)
#~put "stroke$.Last.value" 0 ...
#~
par(mfrow = c(2, 2))
#~put "stroke$.Last.value" 0 ...
#~
scatter.smooth(WAB1, Status, span = 0.2)
#~put "stroke$.Last.value" 0 ...
#~
scatter.smooth(WAB1, Status, span = 0.35)
#~put "stroke$.Last.value" 0 ...
#~
scatter.smooth(WAB1, Status, span = 0.5)
#~put "stroke$.Last.value" 0 ...
#~
scatter.smooth(WAB1, Status, span = 0.65)
#~put "stroke$.Last.value" 0 ...
#~
scatter.smooth(WAB1, Status, span = 0.2)
#~put "stroke$.Last.value" 0 ...
#~
title("Span = 0.2")
#~
scatter.smooth(WAB1, Status, span = 0.35)
#~put "stroke$.Last.value" 0 ...
#~
title("Span = 0.35")
#~
scatter.smooth(WAB1, Status, span = 0.5)
#~put "stroke$.Last.value" 0 ...
#~
title("Span = 0.5")
#~
scatter.smooth(WAB1, Status, span = 0.65)
#~put "stroke$.Last.value" 0 ...
#~
title("Span = 0.65")
#~
par(mfrow = c(1, 1))
#~put "stroke$.Last.value" 0 ...
#~
plot(WAB1, Status)
#~
plot(WAB1, Status + 0.1 * urand(length(Status)))
#~error: couldn't find function "urand"
#~put "stroke$last.dump" 0 ...
#~
randu(10)
#~error: couldn't find function "randu"
#~
rand(10)
#~error: couldn't find function "rand"
#~
plot(WAB1, Status + 0.1 * runif(length(Status)))
#~get "stroke$.Random.seed" 0 ...
#~
.Random.seed <- c(1, 7, 54, 58, 19, 2, 55, 11, 42, 62, 41, 2)
#~put "stroke$.Random.seed" 0 ...
#~
ls()
#~put "stroke$.Last.value" 0 ...
#~
LOSD
#~
plot(LOSD, Status)
#~
scatter.smooth(LOSD, Status, span = 0.65)
#~put "stroke$.Last.value" 0 ...
#~
scatter.smooth(Age, Status, span = 0.65)
#~put "stroke$.Last.value" 0 ...
#~
par(mfrow = c(2, 2))
#~put "stroke$.Last.value" 0 ...
#~
foreach i (c(0.2, 0.35, .5, 0.65)){
#~error: statement with syntax error
#~
for i (c(0.2, 0.35, .5, 0.65)){
#~error: statement with syntax error
#~
#~interrupt
#~
for



#~error: statement with syntax error
#~
for (i in c(0.2, 0.35, .5, 0.65)){
scatter.smooth(Age, Status, span=$i)
#~error: statement with syntax error
#~
for(i in c(0.2, 0.35, 0.5, 0.65)) {
	scatter.smooth(Age, Status, span = i)
	title(paste("Span = ", i))
}
#~
mymod <- glm(Status ~ Side + Gender + Age + Age^2 + LOSD + WAB1, family = 
	binomial)
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.variances" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.links" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/stat/.Functions/glm.weights" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/s/.Datasets/.Machine" -1 ...
#~put "stroke$mymod" 0 ...
#~
summary(mymod)
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.variances" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.links" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/stat/.Functions/glm.weights" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/s/.Datasets/.Machine" -1 ...
#~put "stroke$.Last.value" 0 ...
#~
mymod <- glm(Status ~ Side + Gender + Age^2 + LOSD + WAB1, family = binomial)
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.variances" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.links" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/stat/.Functions/glm.weights" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/s/.Datasets/.Machine" -1 ...
#~put "stroke$mymod" 0 ...
#~
summary(mymod)
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.variances" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.links" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/stat/.Functions/glm.weights" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/s/.Datasets/.Machine" -1 ...
#~put "stroke$.Last.value" 0 ...
#~
mymod <- glm(Status ~ Side + Gender + Age + Age^2 + LOSD + WAB1, family = 
	binomial)
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.variances" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.links" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/stat/.Functions/glm.weights" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/s/.Datasets/.Machine" -1 ...
#~put "stroke$mymod" 0 ...
#~
summary(mymod)
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.variances" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.links" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/stat/.Functions/glm.weights" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/s/.Datasets/.Machine" -1 ...
#~put "stroke$.Last.value" 0 ...
#~
anova(mymod)
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.variances" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.links" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/stat/.Functions/glm.weights" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/s/.Datasets/.Machine" -1 ...
#~put "stroke$.Last.value" 0 ...
#~
mymod <- glm(Status ~ Side + Gender + LOSD + WAB1 + Age + Age^2, family = 
	binomial)
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.variances" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.links" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/stat/.Functions/glm.weights" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/s/.Datasets/.Machine" -1 ...
#~put "stroke$mymod" 0 ...
#~
anova(mymod)
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.variances" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.links" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/stat/.Functions/glm.weights" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/s/.Datasets/.Machine" -1 ...
#~put "stroke$.Last.value" 0 ...
#~
summary(mymod)
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.variances" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.links" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/stat/.Functions/glm.weights" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/s/.Datasets/.Machine" -1 ...
#~put "stroke$.Last.value" 0 ...
#~
mymod <- glm(Status ~ Side + Gender + LOSD + WAB1 + Age, family = binomial)
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.variances" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.links" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/stat/.Functions/glm.weights" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/s/.Datasets/.Machine" -1 ...
#~put "stroke$mymod" 0 ...
#~
residuals(mymod)
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.variances" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.links" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/stat/.Functions/glm.weights" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/s/.Datasets/.Machine" -1 ...
#~put "stroke$.Last.value" 0 ...
#~
residuals.glm
#~
residuals(mymod, type = "deviance")
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.variances" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.links" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/stat/.Functions/glm.weights" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/s/.Datasets/.Machine" -1 ...
#~put "stroke$.Last.value" 0 ...
#~
par(mfrow = c(1, 1))
#~put "stroke$.Last.value" 0 ...
#~
plot(Age, residuals(mymod))
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.variances" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.links" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/stat/.Functions/glm.weights" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/s/.Datasets/.Machine" -1 ...
#~
scatter.smooth(Age, residuals(mymod))
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.variances" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.links" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/stat/.Functions/glm.weights" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/s/.Datasets/.Machine" -1 ...
#~put "stroke$.Last.value" 0 ...
#~
scatter.smooth(Age, residuals(mymod, type = "pearson"))
#~put "stroke$.Last.value" 0 ...
#~
scatter.smooth(Age, residuals(mymod))
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.variances" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.links" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/stat/.Functions/glm.weights" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/s/.Datasets/.Machine" -1 ...
#~put "stroke$.Last.value" 0 ...
#~
identify(Age, residuals(mymod))
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.variances" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.links" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/stat/.Functions/glm.weights" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/s/.Datasets/.Machine" -1 ...
#~put "stroke$.Last.value" 0 ...
#~
Age[59]
#~put "stroke$.Last.value" 0 ...
#~
mymod <- glm(Status ~ Side + Gender + LOSD + WAB1 + Age, family = binomial, 
	data = stroke[-59,  ])
#~get ".Data.hp/stroke" 828557729 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.variances" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.links" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/stat/.Functions/glm.weights" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/s/.Datasets/.Machine" -1 ...
#~put "stroke$mymod" 0 ...
#~
stroke
#~get ".Data.hp/stroke" 828557729 ...
#~
stroke[-59,  ]
#~get ".Data.hp/stroke" 828557729 ...
#~put "stroke$.Last.value" 0 ...
#~
stroke[59,  ]
#~get ".Data.hp/stroke" 828557729 ...
#~put "stroke$.Last.value" 0 ...
#~
mymod <- glm(Status ~ Side + Gender + LOSD + WAB1, family = binomial, data = 
	stroke[-59,  ])
#~get ".Data.hp/stroke" 828557729 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.variances" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.links" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/stat/.Functions/glm.weights" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/s/.Datasets/.Machine" -1 ...
#~put "stroke$mymod" 0 ...
#~
mymod <- glm(Status ~ Side + Gender + LOSD, family = binomial, data = stroke[
	-59,  ])
#~get ".Data.hp/stroke" 828557729 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.variances" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.links" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/stat/.Functions/glm.weights" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/s/.Datasets/.Machine" -1 ...
#~put "stroke$mymod" 0 ...
#~
summary(mymod)
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.variances" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.links" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/stat/.Functions/glm.weights" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/s/.Datasets/.Machine" -1 ...
#~put "stroke$.Last.value" 0 ...
#~
cooksd(mymod)
#~error: couldn't find function "cooksd"
#~
scatter.smooth(Age[-59], residuals(mymod)[-59])
#~error: arguments imply differing number of rows: 59, 60
#~put "stroke$last.dump" 0 ...
#~
scatter.smooth(Age[-59], (residuals(mymod))[-59])
#~error: arguments imply differing number of rows: 59, 60
#~put "stroke$last.dump" 0 ...
#~
length(Age[-59])
#~put "stroke$.Last.value" 0 ...
#~
Age
#~
length(Age)
#~put "stroke$.Last.value" 0 ...
#~
mymod <- glm(Status ~ Side + Gender + LOSD + WAB1 + Age, family = binomial)
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.variances" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.links" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/stat/.Functions/glm.weights" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/s/.Datasets/.Machine" -1 ...
#~put "stroke$mymod" 0 ...
#~
scatter.smooth(Age[-59], (residuals(mymod))[-59])
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.variances" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.links" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/stat/.Functions/glm.weights" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/s/.Datasets/.Machine" -1 ...
#~put "stroke$.Last.value" 0 ...
#~
mymod <- glm(Status ~ Side + Gender + LOSD + WAB1 + Age, family = binomial)
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.variances" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.links" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/stat/.Functions/glm.weights" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/s/.Datasets/.Machine" -1 ...
#~put "stroke$mymod" 0 ...
#~
scatter.smooth(Age, (residuals(mymod)))
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.variances" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.links" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/stat/.Functions/glm.weights" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/s/.Datasets/.Machine" -1 ...
#~put "stroke$.Last.value" 0 ...
#~
identify(Age, (residuals(mymod)))
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.variances" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.links" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/stat/.Functions/glm.weights" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/s/.Datasets/.Machine" -1 ...
#~put "stroke$.Last.value" 0 ...
#~
scatter.smooth(Age[-59], (residuals(mymod))[-59])
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.variances" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.links" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/stat/.Functions/glm.weights" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/s/.Datasets/.Machine" -1 ...
#~put "stroke$.Last.value" 0 ...
#~
par(mfrow = c(2, 1))
#~put "stroke$.Last.value" 0 ...
#~
mymod <- glm(Status ~ Side + Gender + LOSD + WAB1 + Age, family = binomial)
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.variances" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.links" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/stat/.Functions/glm.weights" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/s/.Datasets/.Machine" -1 ...
#~put "stroke$mymod" 0 ...
#~
scatter.smooth(Age, (residuals(mymod)))
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.variances" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.links" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/stat/.Functions/glm.weights" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/s/.Datasets/.Machine" -1 ...
#~put "stroke$.Last.value" 0 ...
#~
identify(Age, (residuals(mymod)))
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.variances" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.links" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/stat/.Functions/glm.weights" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/s/.Datasets/.Machine" -1 ...
#~put "stroke$.Last.value" 0 ...
#~
scatter.smooth(Age[-59], (residuals(mymod))[-59])
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.variances" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.links" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/stat/.Functions/glm.weights" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/s/.Datasets/.Machine" -1 ...
#~put "stroke$.Last.value" 0 ...
#~
par(mfrow = c(4, 4))
#~put "stroke$.Last.value" 0 ...
#~
plot(mymod)
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.variances" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.links" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/stat/.Functions/glm.weights" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/s/.Datasets/.Machine" -1 ...
#~put "stroke$.Last.value" 0 ...
#~
par(mfrow = c(2, 2))
#~put "stroke$.Last.value" 0 ...
#~
plot(mymod)
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.variances" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.links" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/stat/.Functions/glm.weights" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/s/.Datasets/.Machine" -1 ...
#~put "stroke$.Last.value" 0 ...
#~
scatter.smooth(Age[-59], Status[-59])
#~put "stroke$.Last.value" 0 ...
#~
scatter.smooth(Age[-59], Status[-59], span = 0.5)
#~put "stroke$.Last.value" 0 ...
#~
ls
#~
ls()
#~put "stroke$.Last.value" 0 ...
#~
WW <- rep(1, Side)
#~error: rep() only defined for length(times)==1 or length(x)
#~put "stroke$last.dump" 0 ...
#~
WW <- rep(1, length(Side))
#~put "stroke$WW" 0 ...
#~
WW[59] <- 0
#~put "stroke$WW" 0 ...
#~
mymod <- glm(Status ~ Side + Gender + LOSD + WAB1 + Age, weight = WW, family = 
	binomial)
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.variances" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.links" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/stat/.Functions/glm.weights" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/s/.Datasets/.Machine" -1 ...
#~put "stroke$mymod" 0 ...
#~
options()
#~put "stroke$.Last.value" 0 ...
#~
mymod <- glm(Status ~ Side + Gender + LOSD + WAB1 + Age, weight = WW, iter = 30,
	family = binomial)
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.variances" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.links" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/stat/.Functions/glm.weights" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/s/.Datasets/.Machine" -1 ...
#~put "stroke$mymod" 0 ...
#~
help(glm)
#~put "stroke$.Last.value" 0 ...
#~
?control
#~put "stroke$.Last.value" 0 ...
#~
?glm
#~put "stroke$.Last.value" 0 ...
#~
?glm.control
#~put "stroke$.Last.value" 0 ...
#~
mymod <- glm(Status ~ Side + Gender + LOSD + WAB1 + Age, weight = WW, maxit = 
	30, family = binomial)
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.variances" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.links" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/stat/.Functions/glm.weights" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/s/.Datasets/.Machine" -1 ...
#~put "stroke$mymod" 0 ...
#~
summary(mymod)
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.variances" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.links" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/stat/.Functions/glm.weights" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/s/.Datasets/.Machine" -1 ...
#~put "stroke$.Last.value" 0 ...
#~
fitted(mymod)
#~put "stroke$.Last.value" 0 ...
#~
predict(mymod)
#~put "stroke$.Last.value" 0 ...
#~
round(fitted(mymod))
#~put "stroke$.Last.value" 0 ...
#~
fitted(mymod <- glm(Status ~ Side + Gender + LOSD + WAB1 + Age, maxit = 30, 
	family = binomial))
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.variances" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.links" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/stat/.Functions/glm.weights" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/s/.Datasets/.Machine" -1 ...
#~put "stroke$mymod" 0 ...
#~put "stroke$.Last.value" 0 ...
#~
round(fitted(mymod <- glm(Status ~ Side + Gender + LOSD + WAB1 + Age, maxit = 
	30, family = binomial)))
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.variances" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.links" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/stat/.Functions/glm.weights" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/s/.Datasets/.Machine" -1 ...
#~put "stroke$mymod" 0 ...
#~put "stroke$.Last.value" 0 ...
#~
round(fitted(mymod <- glm(Status ~ Side + Gender + LOSD + WAB1 + Age, maxit = 
	30, family = binomial)), 3)
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.variances" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.links" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/stat/.Functions/glm.weights" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/s/.Datasets/.Machine" -1 ...
#~put "stroke$mymod" 0 ...
#~put "stroke$.Last.value" 0 ...
#~
round(fitted(mymod <- glm(Status ~ Side + Gender + LOSD + WAB1 + Age, weight = 
	WW, maxit = 30, family = binomial)), 3)
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.variances" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.links" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/stat/.Functions/glm.weights" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/s/.Datasets/.Machine" -1 ...
#~put "stroke$mymod" 0 ...
#~put "stroke$.Last.value" 0 ...
#~
mymod <- glm(Status ~ Side + Gender + LOSD + WAB1 + Age, weight = WW, family = 
	binomial)
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.variances" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.links" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/stat/.Functions/glm.weights" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/s/.Datasets/.Machine" -1 ...
#~put "stroke$mymod" 0 ...
#~
WW <- rep(1, length(Side))
#~put "stroke$WW" 0 ...
#~
WW[59] <- 0
#~put "stroke$WW" 0 ...
#~
mymod <- glm(Status ~ Side + Gender + LOSD + WAB1 + Age, weight = WW, family = 
	binomial)
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.variances" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.links" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/stat/.Functions/glm.weights" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/s/.Datasets/.Machine" -1 ...
#~put "stroke$mymod" 0 ...
#~
mymod <- glm(Status ~ Side + Gender + LOSD + WAB1 + Age, weight = WW, maxit = 
	30, family = binomial)
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.variances" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.links" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/stat/.Functions/glm.weights" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/s/.Datasets/.Machine" -1 ...
#~put "stroke$mymod" 0 ...
#~
summary(mymod)
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.variances" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.links" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/stat/.Functions/glm.weights" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/s/.Datasets/.Machine" -1 ...
#~put "stroke$.Last.value" 0 ...
#~
anova(mymod)
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.variances" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.links" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/stat/.Functions/glm.weights" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/s/.Datasets/.Machine" -1 ...
#~put "stroke$.Last.value" 0 ...
#~
anova(mymod, maxit = 30)
#~error: invalid formula
#~put "stroke$last.dump" 0 ...
#~
mymod
#~
mymod <- glm(Status ~ Side + Gender + LOSD + WAB1 + Age, weight = WW, maxit = 
	50, family = binomial)
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.variances" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.links" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/stat/.Functions/glm.weights" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/s/.Datasets/.Machine" -1 ...
#~put "stroke$mymod" 0 ...
#~
summary(mymod)
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.variances" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.links" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/stat/.Functions/glm.weights" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/s/.Datasets/.Machine" -1 ...
#~put "stroke$.Last.value" 0 ...
#~
anova(mymod)
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.variances" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.links" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/stat/.Functions/glm.weights" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/s/.Datasets/.Machine" -1 ...
#~put "stroke$.Last.value" 0 ...
#~
mymod
#~
mymod <- glm(Status ~ Side + Gender + LOSD + WAB1 + Age + Age^2, weight = WW, 
	maxit = 50, family = binomial)
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.variances" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.links" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/stat/.Functions/glm.weights" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/s/.Datasets/.Machine" -1 ...
#~put "stroke$mymod" 0 ...
#~
warnings()
#~
summary(mymod)
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.variances" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.links" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/stat/.Functions/glm.weights" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/s/.Datasets/.Machine" -1 ...
#~put "stroke$.Last.value" 0 ...
#~
fitted(mymod)
#~put "stroke$.Last.value" 0 ...
#~
anova(mymod)
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.variances" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.links" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/stat/.Functions/glm.weights" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/s/.Datasets/.Machine" -1 ...
#~put "stroke$.Last.value" 0 ...
#~
summary(mymod)
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.variances" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.links" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/stat/.Functions/glm.weights" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/s/.Datasets/.Machine" -1 ...
#~put "stroke$.Last.value" 0 ...
#~
predict(mymod)
#~put "stroke$.Last.value" 0 ...
#~
fitted(mymod)
#~put "stroke$.Last.value" 0 ...
#~
round(fitted(mymod), 4)
#~put "stroke$.Last.value" 0 ...
#~
Status
#~
model(mymod)
#~error: couldn't find function "model"
#~
formula(mymod)
#~put "stroke$.Last.value" 0 ...
#~
m2 <- update(mymod, ~ . ~ - Side)
#~error: statement with syntax error
#~
m2 <- update(mymod, . ~ . - Side)
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.variances" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.links" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/stat/.Functions/glm.weights" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/s/.Datasets/.Machine" -1 ...
#~put "stroke$m2" 0 ...
#~
anova(m2)
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.variances" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.links" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/stat/.Functions/glm.weights" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/s/.Datasets/.Machine" -1 ...
#~put "stroke$.Last.value" 0 ...
#~
anova(m2 <- update(mymod, . ~ . - Gender))
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.variances" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.links" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/stat/.Functions/glm.weights" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/s/.Datasets/.Machine" -1 ...
#~put "stroke$m2" 0 ...
#~put "stroke$.Last.value" 0 ...
#~
anova(m2 <- update(mymod, . ~ . - LOSD))
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.variances" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.links" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/stat/.Functions/glm.weights" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/s/.Datasets/.Machine" -1 ...
#~put "stroke$m2" 0 ...
#~put "stroke$.Last.value" 0 ...
#~
anova(m2 <- update(mymod, . ~ . - WAB1))
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.variances" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.links" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/stat/.Functions/glm.weights" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/s/.Datasets/.Machine" -1 ...
#~put "stroke$m2" 0 ...
#~put "stroke$.Last.value" 0 ...
#~
anova(mymod)
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.variances" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.links" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/stat/.Functions/glm.weights" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/s/.Datasets/.Machine" -1 ...
#~put "stroke$.Last.value" 0 ...
#~
formula(mymod)
#~put "stroke$.Last.value" 0 ...
#~
anova(mymod)
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.variances" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.links" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/stat/.Functions/glm.weights" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/s/.Datasets/.Machine" -1 ...
#~put "stroke$.Last.value" 0 ...
#~
m <- update(mymod, . ~ . - Age^2 + WAB1^2)
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.variances" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.links" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/stat/.Functions/glm.weights" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/s/.Datasets/.Machine" -1 ...
#~put "stroke$m" 0 ...
#~
anova(m)
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.variances" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.links" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/stat/.Functions/glm.weights" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/s/.Datasets/.Machine" -1 ...
#~put "stroke$.Last.value" 0 ...
#~
summary(mymod)
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.variances" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.links" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/stat/.Functions/glm.weights" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/s/.Datasets/.Machine" -1 ...
#~put "stroke$.Last.value" 0 ...
#~
formula(mymod)
#~put "stroke$.Last.value" 0 ...
#~
mm <- glm(Status ~ Side + Gender + LOSD + WAB1 + Age + WAB1^2, weight = WW, 
	family = binomial)
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.variances" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.links" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/stat/.Functions/glm.weights" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/s/.Datasets/.Machine" -1 ...
#~put "stroke$mm" 0 ...
#~
mm <- glm(Status ~ Side + Gender + LOSD + WAB1 + Age + WAB1^2, maxit = 50, 
	weight = WW, family = binomial)
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.variances" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.links" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/stat/.Functions/glm.weights" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/s/.Datasets/.Machine" -1 ...
#~put "stroke$mm" 0 ...
#~
anova(mymod)
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.variances" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.links" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/stat/.Functions/glm.weights" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/s/.Datasets/.Machine" -1 ...
#~put "stroke$.Last.value" 0 ...
#~
anova(mm)
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.variances" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.links" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/stat/.Functions/glm.weights" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/s/.Datasets/.Machine" -1 ...
#~put "stroke$.Last.value" 0 ...
#~
sum(Status)
#~put "stroke$.Last.value" 0 ...
#~
par(mfrow = c(1, 1))
#~put "stroke$.Last.value" 0 ...
#~
plot(Status, WAB1)
#~
plot(WAB1, Status)
#~
WAB1[Status > 1]
#~put "stroke$.Last.value" 0 ...
#~
WAB1[Status = 1]
#~put "stroke$.Last.value" 0 ...
#~
WAB1[Status == 1]
#~put "stroke$.Last.value" 0 ...
#~
Status
#~
WAB1[Status == 1]
#~put "stroke$.Last.value" 0 ...
#~
sort(WAB1[Status == 1])
#~put "stroke$.Last.value" 0 ...
#~
length(WAB1[WAB1 >= 81.6])
#~put "stroke$.Last.value" 0 ...
#~
length(sort(WAB1[Status == 1]))
#~put "stroke$.Last.value" 0 ...
#~
sort(WAB1[Status == 0])
#~put "stroke$.Last.value" 0 ...
#~
WAB1[WAB1 >= 81.6 && WAB1 <= 91.3]
#~put "stroke$.Last.value" 0 ...
#~
WAB1[(WAB1 >= 81.6) && (WAB1 <= 91.3)]
#~put "stroke$.Last.value" 0 ...
#~
WAB1[(WAB1 >= 81.6) & (WAB1 <= 91.3)]
#~put "stroke$.Last.value" 0 ...
#~
length(WAB1[(WAB1 >= 81.6) & (WAB1 <= 91.3)])
#~put "stroke$.Last.value" 0 ...
#~
#~End session: Time: 828990294; Process: 4993
#~New session: Version: S Wed Mar 29 17:35:08 PST 1995, Time: Wed Apr 10 13:00:55 1996
#~get "/usr/statlocal/lib/Splus3.3/splus/.Datasets/..min.script.id" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Datasets/.Copyright" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Datasets/version" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/.laenv" -1 ...
#~
stroke
#~get ".Data.hp/stroke" 828557729 ...
#~
attach(stroke, pos = 1)
#~get ".Data.hp/stroke" 828557729 ...
#~put "stroke$.Last.value" 0 ...
#~
ls()
#~put "stroke$.Last.value" 0 ...
#~
help.start()
#~put "stroke$.Last.value" 0 ...
#~
glm.influence()
#~error: couldn't find function "glm.influence"
#~
glm.influence
#~error: Object "glm.influence" not found
#~
lm.influence
#~
influence
#~error: Object "influence" not found
#~
objects()
#~put "stroke$.Last.value" 0 ...
#~
mymod <- glm(Status ~ Gender + LOSD + Side + Age + WAB1, family = binomial)
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.variances" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.links" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/stat/.Functions/glm.weights" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/s/.Datasets/.Machine" -1 ...
#~put "stroke$mymod" 0 ...
#~
fitted(mymod)
#~put "stroke$.Last.value" 0 ...
#~
summary(mymod)
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.variances" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.links" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/stat/.Functions/glm.weights" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/s/.Datasets/.Machine" -1 ...
#~put "stroke$.Last.value" 0 ...
#~
mymod
#~
attributes(mymod)
#~put "stroke$.Last.value" 0 ...
#~
coefficients(mymod)
#~put "stroke$.Last.value" 0 ...
#~
WW <- rep(1, length(Side))
#~put "stroke$WW" 0 ...
#~
mymod <- glm(Status ~ Gender + LOSD + Side + Age + WAB1, weight = WW, family = 
	binomial)
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.variances" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.links" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/stat/.Functions/glm.weights" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/s/.Datasets/.Machine" -1 ...
#~put "stroke$mymod" 0 ...
#~
W <- WW
#~put "stroke$W" 0 ...
#~
for (i in 1:length(Side)){
WW<-W; WW[i] <- 0;
print(coefficients(update(mymod, ~ . ~)))
#~error: statement with syntax error
#~
for(i in 1:length(Side)) {
	WW <- W
	WW[i] <- 0
	print(coefficients(update(mymod, . ~ .)))
}
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.variances" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.links" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/stat/.Functions/glm.weights" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/s/.Datasets/.Machine" -1 ...
#~put "stroke$WW" 0 ...
#~
length(Side)
#~put "stroke$.Last.value" 0 ...
#~
matrix(3, 10)
#~put "stroke$.Last.value" 0 ...
#~
?matrix
#~put "stroke$.Last.value" 0 ...
#~
matrix(0, 3, 10)
#~put "stroke$.Last.value" 0 ...
#~
attributes(mymod)
#~put "stroke$.Last.value" 0 ...
#~
summary.glm
#~
source("poorinf")
#~error: can't open source file "poorinf"
#~put "stroke$last.dump" 0 ...
#~
source("poorinf.S")
WW <- rep(1,length(Side))
mymod <- glm(Status ~ Gender + LOSD + Side + Age + WAB1, weight=WW,family=>
#~error: statement with syntax error
#~put "stroke$last.dump" 0 ...
#~
source("poorinf.S")
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.variances" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.links" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/stat/.Functions/glm.weights" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/s/.Datasets/.Machine" -1 ...
#~put "stroke$WW" 0 ...
#~put "stroke$mymod" 0 ...
#~put "stroke$W" 0 ...
#~put "stroke$infs" 0 ...
#~put "stroke$.Last.value" 0 ...
#~
infs
#~
summary.glm
#~
infs
#~
args(apply)
#~put "stroke$.Last.value" 0 ...
#~
apply(infs, 2, var)
#~put "stroke$.Last.value" 0 ...
#~
sqrt(apply(infs, 2, var))
#~put "stroke$.Last.value" 0 ...
#~
infs/sqrt(apply(infs, 2, var))
#~put "stroke$.Last.value" 0 ...
#~
sqrt(apply(infs, 2, var))
#~put "stroke$.Last.value" 0 ...
#~
sqrt(var(infs[, 1]))
#~put "stroke$.Last.value" 0 ...
#~
infs
#~
source("poorinf.S")
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.variances" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.links" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/stat/.Functions/glm.weights" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/s/.Datasets/.Machine" -1 ...
#~put "stroke$WW" 0 ...
#~put "stroke$mymod" 0 ...
#~put "stroke$base.coef" 0 ...
#~put "stroke$W" 0 ...
#~put "stroke$infs" 0 ...
#~put "stroke$.Last.value" 0 ...
#~
matrix(base.coef, 6, 51)
#~put "stroke$.Last.value" 0 ...
#~
matrix(base.coef, 51, 6)
#~put "stroke$.Last.value" 0 ...
#~
t(matrix(base.coef, 6, 51))
#~put "stroke$.Last.value" 0 ...
#~
infs - t(matrix(base.coef, 6, 51))
#~error: Dimension attributes do not match
#~put "stroke$last.dump" 0 ...
#~
infs - t(matrix(base.coef, 6, 61))
#~put "stroke$.Last.value" 0 ...
#~
motif()
#~get "/usr/statlocal/lib/Splus3.3/splus/.Datasets/.PostScript.Options" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Datasets/ps.preamble.std" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Datasets/ps.setcolor.rgb" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Datasets/ps.setfont.std" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Datasets/ps.paper.regions" -1 ...
#~put "stroke$.Last.value" 0 ...
#~
apply(inf, 2, plot, 1:length(Side))
#~error: Object "inf" not found
#~put "stroke$last.dump" 0 ...
#~
apply(infs, 2, plot, 1:length(Side))
#~error: Invalid value for dimension 1
#~put "stroke$last.dump" 0 ...
#~
par(mfrow = c(2, 3))
#~put "stroke$.Last.value" 0 ...
#~
apply(infs, 2, plot, 1:length(Side))
#~error: Invalid value for dimension 1
#~put "stroke$last.dump" 0 ...
#~
for(i in 1:6) {
	plot(1:length(Side), infs[, i])
}
#~
?plot
#~put "stroke$.Last.value" 0 ...
#~
?symbols
#~put "stroke$.Last.value" 0 ...
#~
?par
#~put "stroke$.Last.value" 0 ...
#~
for(i in 1:6) {
	plot(1:length(Side), infs[, i])
	identify(1:length(Side), infs[, i])
}
#~
attributes(mymod)
#~put "stroke$.Last.value" 0 ...
#~
mymod$formula
#~put "stroke$.Last.value" 0 ...
#~
summary.glm
#~
names(coefficients(mymod))
#~put "stroke$.Last.value" 0 ...
#~
for(i in 1:6) {
	plot(1:length(Side), infs[, i])
	title(names(coefficients(mymod))[i])
	identify(1:length(Side), infs[, i])
}
#~
source("poorinf.S")
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.variances" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.links" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/stat/.Functions/glm.weights" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/s/.Datasets/.Machine" -1 ...
#~put "stroke$WW" 0 ...
#~put "stroke$mymod" 0 ...
#~put "stroke$base.coef" 0 ...
#~put "stroke$W" 0 ...
#~put "stroke$infs" 0 ...
#~put "stroke$.Last.value" 0 ...
#~
par(mfrow = c(1, 1))
#~put "stroke$.Last.value" 0 ...
#~
plot(WAB1, Status)
#~
identify(WAB1, Status)
#~put "stroke$.Last.value" 0 ...
#~
plot(Age, Status)
#~
identify(Age, Status)
#~put "stroke$.Last.value" 0 ...
#~
WW <- W
#~put "stroke$WW" 0 ...
#~
WW[59] <- 0
#~put "stroke$WW" 0 ...
#~
summary(update(mymod))
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.variances" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.links" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/stat/.Functions/glm.weights" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/s/.Datasets/.Machine" -1 ...
#~put "stroke$.Last.value" 0 ...
#~
library(MASS)
#~put "stroke$.Last.value" 0 ...
#~
collett
#~error: Object "collett" not found
#~
search()
#~put "stroke$.Last.value" 0 ...
#~
ls(pos = 10)
#~put "stroke$.Last.value" 0 ...
#~
budworm
#~error: Object "budworm" not found
#~
search()
#~put "stroke$.Last.value" 0 ...
#~
detach(1, save = F)
#~put ".Data.hp/.Last.value" 829164544 ...
#~
bud <- read.table("budworm.dat", head = T)
#~put ".Data.hp/bud" 829164553 ...
#~
bud
#~get ".Data.hp/bud" 829164553 ...
#~
attach(bud, pos = 1)
#~get ".Data.hp/bud" 829164553 ...
#~put "bud$.Last.value" 0 ...
#~
N <- rep(20, length(Response))
#~put "bud$N" 0 ...
#~
plot(Response/N, Dose)
#~
plot((Response/N)[Gender == M], Dose[Gender == M], pch = "M")
#~error: Object "M" not found
#~put "bud$last.dump" 0 ...
#~
Gender == M
#~error: Object "M" not found
#~
Gender == "M"
#~put "bud$.Last.value" 0 ...
#~
points((Response/N)[Gender == "M"], Dose[Gender == "M"], pch = "M")
#~
points((Response/N)[Gender == "F"], Dose[Gender == "F"], pch = "F")
#~
plot(Response/N, Dose, log = "x")
#~
Dose
#~
plot(Dose, Response/N, log = "x")
#~
plot(Dose, Response/N)
#~
points(Dose[Gender == "F"], (Response/N)[Gender == "F"], pch = "F")
#~
points(Dose[Gender == "M"], (Response/N)[Gender == "M"], pch = "M")
#~
plot(Dose, Response/N, log = "x")
#~
points(Dose[Gender == "F"], (Response/N)[Gender == "F"], pch = "F")
#~
points(Dose[Gender == "M"], (Response/N)[Gender == "M"], pch = "M")
#~
Response/N
#~put "bud$.Last.value" 0 ...
#~
logit
#~error: Object "logit" not found
#~
logit <- function(x)
(log(x/(1 - x)))
#~put "bud$logit" 0 ...
#~
plot(Dose, logit(Response/N), log = "x")
#~
plot(Dose, logit(Response/N))
#~
par(mfrow = c(2, 1))
#~put "bud$.Last.value" 0 ...
#~
plot(Dose, logit(Response/N), log = "x")
#~
title("Logit vs Log(Dose)")
#~
plot(Dose, logit(Response/N))
#~
title("Logit vs Dose")
#~
opitions(contrasts = c("cont.treatment", "contr.poly"))
#~error: couldn't find function "opitions"
#~
options(contrasts = c("cont.treatment", "contr.poly"))
#~put "bud$.Last.value" 0 ...
#~
log(Dose)
#~put "bud$.Last.value" 0 ...
#~
Dose
#~
ln(Dose)
#~error: couldn't find function "ln"
#~
Dose
#~
ldose <- log(Dose)
#~put "bud$ldose" 0 ...
#~
SF <- cbind(Response, N - Response)
#~put "bud$SF" 0 ...
#~
mymod <- glm(SF ~ Gender * ldose, family = binomial)
#~error: couldn't find function "cont.treatment"
#~put "bud$last.dump" 0 ...
#~
options(contrasts = c("contr.treatment", "contr.poly"))
#~put "bud$.Last.value" 0 ...
#~
mymod <- glm(SF ~ Gender * ldose, family = binomial)
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.variances" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.links" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/stat/.Functions/glm.weights" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/s/.Datasets/.Machine" -1 ...
#~put "bud$mymod" 0 ...
#~
summary(mymod)
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.variances" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.links" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/stat/.Functions/glm.weights" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/s/.Datasets/.Machine" -1 ...
#~put "bud$.Last.value" 0 ...
#~
anova(mymod)
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.variances" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.links" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/stat/.Functions/glm.weights" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/s/.Datasets/.Machine" -1 ...
#~put "bud$.Last.value" 0 ...
#~
levels(Gender)
#~put "bud$.Last.value" 0 ...
#~
anova(mymod, test = "Chisq")
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.variances" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.links" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/stat/.Functions/glm.weights" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/s/.Datasets/.Machine" -1 ...
#~put "bud$.Last.value" 0 ...
#~
anova(mymod <- glm(SF ~ Gender, family = binomial))
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.variances" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.links" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/stat/.Functions/glm.weights" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/s/.Datasets/.Machine" -1 ...
#~put "bud$mymod" 0 ...
#~put "bud$.Last.value" 0 ...
#~
anova(mymod <- glm(SF ~ ldose, family = binomial))
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.variances" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.links" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/stat/.Functions/glm.weights" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/s/.Datasets/.Machine" -1 ...
#~put "bud$mymod" 0 ...
#~put "bud$.Last.value" 0 ...
#~
anova(mymod <- glm(SF ~ ldose + Gender, family = binomial))
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.variances" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.links" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/stat/.Functions/glm.weights" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/s/.Datasets/.Machine" -1 ...
#~put "bud$mymod" 0 ...
#~put "bud$.Last.value" 0 ...
#~
summary(mymod <- glm(SF ~ ldose + Gender, family = binomial))
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.variances" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.links" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/stat/.Functions/glm.weights" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/s/.Datasets/.Machine" -1 ...
#~put "bud$mymod" 0 ...
#~put "bud$.Last.value" 0 ...
#~
summary(mymod <- glm(SF ~ ldose + Gender + ldose * Gender, family = binomial))
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.variances" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.links" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/stat/.Functions/glm.weights" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/s/.Datasets/.Machine" -1 ...
#~put "bud$mymod" 0 ...
#~put "bud$.Last.value" 0 ...
#~
anova(mymod <- glm(SF ~ ldose + Gender + ldose * Gender, family = binomial))
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.variances" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.links" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/stat/.Functions/glm.weights" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/s/.Datasets/.Machine" -1 ...
#~put "bud$mymod" 0 ...
#~put "bud$.Last.value" 0 ...
#~
anova(mymod <- glm(SF ~ ldose + Gender, family = binomial))
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.variances" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.links" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/stat/.Functions/glm.weights" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/s/.Datasets/.Machine" -1 ...
#~put "bud$mymod" 0 ...
#~put "bud$.Last.value" 0 ...
#~
anova(mymod <- glm(SF ~ Gender + ldose, family = binomial))
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.variances" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.links" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/stat/.Functions/glm.weights" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/s/.Datasets/.Machine" -1 ...
#~put "bud$mymod" 0 ...
#~put "bud$.Last.value" 0 ...
#~
summary(mymod)
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.variances" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.links" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/stat/.Functions/glm.weights" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/s/.Datasets/.Machine" -1 ...
#~put "bud$.Last.value" 0 ...
#~
plot(mymod)
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.variances" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.links" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/stat/.Functions/glm.weights" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/s/.Datasets/.Machine" -1 ...
#~put "bud$.Last.value" 0 ...
#~
par(mfrow = c(2, 2))
#~put "bud$.Last.value" 0 ...
#~
plot(mymod)
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.variances" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.links" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/stat/.Functions/glm.weights" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/s/.Datasets/.Machine" -1 ...
#~put "bud$.Last.value" 0 ...
#~
par(mrfow = c(1, 1))
#~put "bud$last.dump" 0 ...
#~
par(mfrow = c(1, 1))
#~put "bud$.Last.value" 0 ...
#~
plot(1:32, 0:1, type = "n", xlab = "Dose", ylab = "Prob", log = "x")
#~put "bud$last.dump" 0 ...
#~
plot(c(1, 32), 0:1, type = "n", xlab = "Dose", ylab = "Prob", log = "x")
#~
predict(mymod)
#~put "bud$.Last.value" 0 ...
#~
points(Dose[Gender == "F"], (Response/N)[Gender == "F"], pch = "F")
#~
points(Dose[Gender == "M"], (Response/N)[Gender == "M"], pch = "M")
#~
lines(Dose[Gender == "M"], predict(mymod)[Gender == "M"])
#~
lines(Dose[Gender == "M"], fitted(mymod)[Gender == "M"])
#~
lines(Dose[Gender == "F"], fitted(mymod)[Gender == "F"])
#~
text(Dose, Response/N, as.character(Gender))
#~
seq(0, 5, 0.1)
#~put "bud$.Last.value" 0 ...
#~
ld <- seq(0, 5, 0.1)
#~put "bud$ld" 0 ...
#~
lines(2^ld, predict(mymod, data.frame(ldose = ld, sex = factor(rep("M", length(
	ld)), levels = levels(sex))), type = "response"))
#~error: Object "sex" not found
#~put "bud$last.dump" 0 ...
#~
lines(2^ld, predict(mymod, data.frame(ldose = ld, Gender = factor(rep("M", 
	length(ld)), levels = levels(Gender))), type = "response"))
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.variances" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.links" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/stat/.Functions/glm.weights" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/s/.Datasets/.Machine" -1 ...
#~
lines(2^ld, predict(mymod, data.frame(ldose = ld, Gender = factor(rep("F", 
	length(ld)), levels = levels(Gender))), type = "response"))
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.variances" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.links" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/stat/.Functions/glm.weights" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/s/.Datasets/.Machine" -1 ...
#~
plot(c(1, 32), 0:1, type = "n", xlab = "Dose", ylab = "Prob", log = "x")
#~
text(Dose, Response/N, as.character(Gender))
#~
ls()
#~put "bud$.Last.value" 0 ...
#~
ldose
#~
lines( 2^ld , predict(mymod, data.frame(ldose=ld, Gender=factor(rep("F", length(ld))))), type="response"))
#~error: statement with syntax error
#~
lines(2^ld, predict(mymod, data.frame(ldose = ld, Gender = factor(rep("F", 
	length(ld)))), type = "response"))
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.variances" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.links" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/stat/.Functions/glm.weights" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/s/.Datasets/.Machine" -1 ...
#~
lines(2^ld, predict(mymod, data.frame(ldose = ld, Gender = factor(rep("M", 
	length(ld)))), type = "response"))
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.variances" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.links" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/stat/.Functions/glm.weights" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/s/.Datasets/.Machine" -1 ...
#~
lines(2^ld, predict(mymod, data.frame(ldose = ld, Gender = factor(rep("F", 
	length(ld)))), type = "response"))
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.variances" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.links" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/stat/.Functions/glm.weights" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/s/.Datasets/.Machine" -1 ...
#~
data.frame(ldose = ld, Gender = factor(rep("F", length(ld))))
#~put "bud$.Last.value" 0 ...
#~
predict(mymod, data.frame(ldose = ld, Gender = factor(rep("F", length(ld)))))
#~put "bud$.Last.value" 0 ...
#~
predict(mymod)
#~put "bud$.Last.value" 0 ...
#~
ld
#~
predict(mymod, data.frame(ldose = log(2^ld), Gender = factor(rep("F", length(ld
	)))))
#~put "bud$.Last.value" 0 ...
#~
predict(mymod, data.frame(ldose = log(2^ld), Gender = factor(rep("F", length(ld
	)))))
#~put "bud$.Last.value" 0 ...
#~
lines(2^ld, predict(mymod, data.frame(ldose = log(2^ld), Gender = factor(rep(
	"F", length(ld))))))
#~
lines(2^ld, predict(mymod, data.frame(ldose = log(2^ld), Gender = factor(rep(
	"F", length(ld)))), type = "response"))
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.variances" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.links" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/stat/.Functions/glm.weights" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/s/.Datasets/.Machine" -1 ...
#~
lines(2^ld, predict(mymod, data.frame(ldose = log(2^ld), Gender = factor(rep(
	"M", length(ld)))), type = "response"))
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.variances" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.links" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/stat/.Functions/glm.weights" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/s/.Datasets/.Machine" -1 ...
#~
plot(c(1, 32), 0:1, type = "n", xlab = "Dose", ylab = "Prob", log = "x")
#~
text(Dose, Response/N, as.character(Gender))
#~
lines(2^ld, predict(mymod, data.frame(ldose = log(2^ld), Gender = factor(rep(
	"M", length(ld)))), type = "response"))
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.variances" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.links" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/stat/.Functions/glm.weights" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/s/.Datasets/.Machine" -1 ...
#~
lines(2^ld, predict(mymod, data.frame(ldose = log(2^ld), Gender = factor(rep(
	"F", length(ld)))), type = "response"))
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.variances" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.links" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/stat/.Functions/glm.weights" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/s/.Datasets/.Machine" -1 ...
#~
family(mymod)
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.variances" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.links" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/stat/.Functions/glm.weights" -1 ...
#~put "bud$.Last.value" 0 ...
#~
family(mymod)$link
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.variances" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.links" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/stat/.Functions/glm.weights" -1 ...
#~put "bud$.Last.value" 0 ...
#~
family(mymod)$link(p)
#~error: Object "p" not found
#~put "bud$last.dump" 0 ...
#~
family(mymod)$link(0.5)
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.variances" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.links" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/stat/.Functions/glm.weights" -1 ...
#~put "bud$.Last.value" 0 ...
#~
family(mymod)$link(0.1)
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.variances" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.links" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/stat/.Functions/glm.weights" -1 ...
#~put "bud$.Last.value" 0 ...
#~
family(mymod)$link
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.variances" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.links" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/stat/.Functions/glm.weights" -1 ...
#~put "bud$.Last.value" 0 ...
#~
birthwt
#~get "/usr/statlocal/lib/Splus3.3/library/MASS/.Data/birthwt" 814043417 ...
#~
birthwt
#~get "/usr/statlocal/lib/Splus3.3/library/MASS/.Data/birthwt" 814043417 ...
#~
ptl
#~error: Object "ptl" not found
#~
attach(birthwt, pos = 1)
#~get "/usr/statlocal/lib/Splus3.3/library/MASS/.Data/birthwt" 814043417 ...
#~put "birthwt$.Last.value" 0 ...
#~
ptl
#~
summary(ptl)
#~put "birthwt$.Last.value" 0 ...
#~
factor(ptl > 0)
#~put "birthwt$.Last.value" 0 ...
#~
ftv
#~
ptd <- factor(ptl)
#~put "birthwt$ptd" 0 ...
#~
ftv <- factor(ftv)
#~put "birthwt$ftv" 0 ...
#~
ftv
#~
levels(ftv)[ - (1:2)] <- "2+"
#~put "birthwt$ftv" 0 ...
#~
ftv
#~
low
#~
low <- factor(low)
#~put "birthwt$low" 0 ...
#~
smoke
#~
smoke > 0
#~put "birthwt$.Last.value" 0 ...
#~
ht
#~
mymod <- glm(low ~ ., family = binomial)
#~error: Object "." not found
#~put "birthwt$last.dump" 0 ...
#~
bwt <- data.frame(low, age, lwt, race, smoke, ptd, ht, ui, ftv)
#~put "birthwt$bwt" 0 ...
#~
assign("btw", btw, pos = 0)
#~put "birthwt$last.dump" 0 ...
#~
?assign
#~put "birthwt$.Last.value" 0 ...
#~
assign("btw", btw, frame = 0)
#~error: Object "btw" not found
#~put "birthwt$last.dump" 0 ...
#~
assign("bwt", bwt, frame = 0)
#~put "birthwt$.Last.value" 0 ...
#~
bwt
#~
detach(1, save = F)
#~put "bud$.Last.value" 0 ...
#~
bwt
#~
mymod <- glm(low ~ ., family = binomial, data = bwt)
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.variances" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.links" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/stat/.Functions/glm.weights" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/s/.Datasets/.Machine" -1 ...
#~put "bud$mymod" 0 ...
#~
summary(mymod)
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.variances" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/glm.links" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/stat/.Functions/glm.weights" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/s/.Datasets/.Machine" -1 ...
#~put "bud$.Last.value" 0 ...
#~
#~End session: Time: 829187356; Process: 7125
#~New session: Version: S Wed Mar 29 17:35:08 PST 1995, Time: Mon Apr 15 13:26:15 1996
#~get "/usr/statlocal/lib/Splus3.3/splus/.Datasets/..min.script.id" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Datasets/.Copyright" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Datasets/version" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/.laenv" -1 ...
#~
motif()
#~get "/usr/statlocal/lib/Splus3.3/splus/.Datasets/.PostScript.Options" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Datasets/ps.preamble.std" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Datasets/ps.setcolor.rgb" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Datasets/ps.setfont.std" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Datasets/ps.paper.regions" -1 ...
#~put ".Data.hp/.Last.value" 829589184 ...
#~
#~End session: Time: 829589186; Process: 562
#~New session: Version: S Wed Mar 29 17:35:08 PST 1995, Time: Mon Apr 22 11:41:23 1996
#~get "/usr/statlocal/lib/Splus3.3/splus/.Datasets/..min.script.id" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Datasets/.Copyright" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Datasets/version" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/.laenv" -1 ...
#~
motif()
#~get "/usr/statlocal/lib/Splus3.3/splus/.Datasets/.PostScript.Options" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Datasets/ps.preamble.std" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Datasets/ps.setcolor.rgb" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Datasets/ps.setfont.std" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Datasets/ps.paper.regions" -1 ...
#~put ".Data.hp/.Last.value" 830187698 ...
#~
par(mrfow = c(4, 1))
#~put ".Data.hp/last.dump" 830187708 ...
#~
par(mfrow = c(4, 1))
#~put ".Data.hp/.Last.value" 830187713 ...
#~
tsplot(dd <- arima.sim(200, model = list(ar = c(0.7, 0.2))))
#~get ".Data.hp/.Random.seed" 0 ...
#~put ".Data.hp/dd" 830187726 ...
#~
.Random.seed <- c(61, 6, 20, 23, 59, 3, 42, 12, 40, 52, 45, 3)
#~put ".Data.hp/.Random.seed" 0 ...
#~
acf(dd)
#~get "/usr/statlocal/lib/Splus3.3/s/.Datasets/gr.pars" -1 ...
#~get ".Data.hp/dd" 830187726 ...
#~put ".Data.hp/.Last.value" 830187732 ...
#~
tsplot(dd <- arima.sim(200, model = list(ar = c(0.7))))
#~get ".Data.hp/.Random.seed" 0 ...
#~put ".Data.hp/dd" 830187743 ...
#~
.Random.seed <- c(53, 36, 49, 20, 61, 3, 48, 46, 60, 19, 41, 2)
#~put ".Data.hp/.Random.seed" 0 ...
#~
acf(dd)
#~get "/usr/statlocal/lib/Splus3.3/s/.Datasets/gr.pars" -1 ...
#~get ".Data.hp/dd" 830187743 ...
#~put ".Data.hp/.Last.value" 830187750 ...
#~
acf(arima.sim(200, model = list(ar = c(0.7, 0.2))), type = "p")
#~get "/usr/statlocal/lib/Splus3.3/s/.Datasets/gr.pars" -1 ...
#~get ".Data.hp/.Random.seed" 0 ...
#~put ".Data.hp/.Last.value" 830187799 ...
#~
.Random.seed <- c(29, 17, 34, 38, 18, 2, 43, 59, 4, 60, 33, 3)
#~put ".Data.hp/.Random.seed" 0 ...
#~
acf(arima.sim(200, model = list(ar = c(0.7, 0.2))), type = "p")
#~get "/usr/statlocal/lib/Splus3.3/s/.Datasets/gr.pars" -1 ...
#~get ".Data.hp/.Random.seed" 0 ...
#~put ".Data.hp/.Last.value" 830187802 ...
#~
.Random.seed <- c(37, 45, 15, 7, 33, 3, 2, 36, 27, 60, 37, 1)
#~put ".Data.hp/.Random.seed" 0 ...
#~
acf(arima.sim(200, model = list(ar = c(0.7))), type = "p")
#~get "/usr/statlocal/lib/Splus3.3/s/.Datasets/gr.pars" -1 ...
#~get ".Data.hp/.Random.seed" 0 ...
#~put ".Data.hp/.Last.value" 830187809 ...
#~
.Random.seed <- c(37, 13, 18, 42, 5, 3, 5, 43, 44, 10, 18, 3)
#~put ".Data.hp/.Random.seed" 0 ...
#~
acf(arima.sim(200, model = list(ar = c(0.7))), type = "p")
#~get "/usr/statlocal/lib/Splus3.3/s/.Datasets/gr.pars" -1 ...
#~get ".Data.hp/.Random.seed" 0 ...
#~put ".Data.hp/.Last.value" 830187811 ...
#~
.Random.seed <- c(45, 4, 26, 16, 37, 2, 60, 38, 2, 46, 38, 0)
#~put ".Data.hp/.Random.seed" 0 ...
#~
acf(arima.sim(200, model = list(ma = c(-0.9))), type = "p")
#~get "/usr/statlocal/lib/Splus3.3/s/.Datasets/gr.pars" -1 ...
#~get ".Data.hp/.Random.seed" 0 ...
#~put ".Data.hp/.Last.value" 830187865 ...
#~
.Random.seed <- c(45, 62, 28, 36, 52, 0, 2, 33, 43, 5, 55, 2)
#~put ".Data.hp/.Random.seed" 0 ...
#~
acf(arima.sim(200, model = list(ma = c(0.9))), type = "p")
#~get "/usr/statlocal/lib/Splus3.3/s/.Datasets/gr.pars" -1 ...
#~get ".Data.hp/.Random.seed" 0 ...
#~put ".Data.hp/.Last.value" 830187876 ...
#~
.Random.seed <- c(21, 31, 52, 19, 55, 3, 11, 59, 48, 29, 4, 1)
#~put ".Data.hp/.Random.seed" 0 ...
#~
acf(arima.sim(200, model = list(ma = c(0.9, 0.6))), type = "p")
#~get "/usr/statlocal/lib/Splus3.3/s/.Datasets/gr.pars" -1 ...
#~get ".Data.hp/.Random.seed" 0 ...
#~put ".Data.hp/.Last.value" 830187887 ...
#~
.Random.seed <- c(5, 17, 41, 63, 5, 3, 36, 32, 36, 58, 60, 1)
#~put ".Data.hp/.Random.seed" 0 ...
#~
acf(arima.sim(200, model = list(ma = c(0.9, 0.6))), type = "p")
#~get "/usr/statlocal/lib/Splus3.3/s/.Datasets/gr.pars" -1 ...
#~get ".Data.hp/.Random.seed" 0 ...
#~put ".Data.hp/.Last.value" 830187892 ...
#~
.Random.seed <- c(53, 63, 47, 0, 52, 2, 51, 15, 22, 44, 57, 3)
#~put ".Data.hp/.Random.seed" 0 ...
#~
tsplot(dd <- arima.sim(200, model = list(ma = c(0.6))))
#~get ".Data.hp/.Random.seed" 0 ...
#~put ".Data.hp/dd" 830188026 ...
#~
.Random.seed <- c(45, 7, 51, 61, 27, 0, 22, 36, 22, 25, 53, 0)
#~put ".Data.hp/.Random.seed" 0 ...
#~
acf(dd)
#~get "/usr/statlocal/lib/Splus3.3/s/.Datasets/gr.pars" -1 ...
#~get ".Data.hp/dd" 830188026 ...
#~put ".Data.hp/.Last.value" 830188032 ...
#~
tsplot(dd <- arima.sim(200, model = list(ma = c(0.6, 0.4))))
#~get ".Data.hp/.Random.seed" 0 ...
#~put ".Data.hp/dd" 830188042 ...
#~
.Random.seed <- c(37, 54, 17, 12, 18, 3, 42, 6, 15, 42, 6, 1)
#~put ".Data.hp/.Random.seed" 0 ...
#~
acf(dd)
#~get "/usr/statlocal/lib/Splus3.3/s/.Datasets/gr.pars" -1 ...
#~get ".Data.hp/dd" 830188042 ...
#~put ".Data.hp/.Last.value" 830188044 ...
#~
#~End session: Time: 830189254; Process: 3500
#~New session: Version: S Wed Mar 29 17:35:08 PST 1995, Time: Tue May 14 11:08:22 1996
#~get "/usr/statlocal/lib/Splus3.3/splus/.Datasets/..min.script.id" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Datasets/.Copyright" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Datasets/version" -1 ...
#~get "/usr/statlocal/lib/Splus3.3/splus/.Functions/.laenv" -1 ...
#~
fff <- read.table("dodo", head = T)
#~put ".Data.hp/fff" 832086521 ...
#~
fff
#~get ".Data.hp/fff" 832086521 ...
#~
fff <- read.table("dodo", head = T)
#~put ".Data.hp/fff" 832087052 ...
#~
ff
#~error: Object "ff" not found
#~
fff
#~get ".Data.hp/fff" 832087052 ...
#~
?apply
#~put ".Data.hp/.Last.value" 832087078 ...
#~
apply(fff, 1, mean)
#~get ".Data.hp/fff" 832087052 ...
#~put ".Data.hp/.Last.value" 832087086 ...
#~
data.frame(apply(fff, 1, mean))
#~get ".Data.hp/fff" 832087052 ...
#~put ".Data.hp/.Last.value" 832087096 ...
#~
mean.fun
#~error: Object "mean.fun" not found
#~
mean.fun <- function(x)
mean(x[!is.na(x)])
#~put ".Data.hp/mean.fun" 832087136 ...
#~
data.frame(apply(fff, 1, mean.fun))
#~get ".Data.hp/fff" 832087052 ...
#~put ".Data.hp/.Last.value" 832087141 ...
#~
mean.fun2 <- function(x)
{
	x[is.na(x)] <- 0
	mean(x)
}
#~put ".Data.hp/mean.fun2" 832087208 ...
#~
data.frame(apply(fff, 1, mean.fun2))
#~get ".Data.hp/fff" 832087052 ...
#~put ".Data.hp/.Last.value" 832087215 ...
#~
data.frame(S = apply(fff, 1, mean.fun), T = apply(fff, 1, mean.fun2))
#~get ".Data.hp/fff" 832087052 ...
#~put ".Data.hp/.Last.value" 832087249 ...
#~
fix(mean.fun)
#~put ".Data.hp/mean.fun" 832087283 ...
#~put ".Data.hp/.Last.fixed" 832087284 ...
#~put ".Data.hp/.Last.value" 832087284 ...
#~
?fix
#~put ".Data.hp/.Last.value" 832087337 ...
#~
options$editor
#~put ".Data.hp/.Last.value" 832087345 ...
#~
options
#~
options()
#~put ".Data.hp/.Last.value" 832087351 ...
#~
options(editor = "emacs")
#~put ".Data.hp/.Last.value" 832087363 ...
#~
fix(mean.fun)
#~put ".Data.hp/mean.fun" 832087391 ...
#~put ".Data.hp/.Last.fixed" 832087391 ...
#~put ".Data.hp/.Last.value" 832087391 ...
#~
fix(mean.fun2)
#~put ".Data.hp/mean.fun2" 832087441 ...
#~put ".Data.hp/.Last.fixed" 832087441 ...
#~put ".Data.hp/.Last.value" 832087442 ...
#~
data.frame(S = apply(fff, 1, mean.fun), T = apply(fff, 1, mean.fun2))
#~get ".Data.hp/fff" 832087052 ...
#~put ".Data.hp/.Last.value" 832087446 ...
#~
round(apply(fff, 1, mean.fun2),  )
#~error: No data to interpret as numerical value
#~put ".Data.hp/last.dump" 832087883 ...
#~
round(apply(fff, 1, mean.fun2), 2)
#~get ".Data.hp/fff" 832087052 ...
#~put ".Data.hp/.Last.value" 832087900 ...
#~
data.frame(T = round(apply(fff, 1, mean.fun2), 2))
#~get ".Data.hp/fff" 832087052 ...
#~put ".Data.hp/.Last.value" 832087912 ...
#~
#~End session: Time: 832088072; Process: 4801
