I have no real idea what you are trying to do, but if a table is
what you want, you can probably get it using the table() function.
Or, more likely, the xtabs() function.
Using your example from an earlier post (adjusted to make it
comprehensible to the human mind):
set.seed(1000)
time <- factor(rep(c("Pre","Post"),each=200))
treatment <- factor(rep(rep(c("Control","Treatment"),each=100),2))
mu <- rep(rep(1:2,each=100),2)
response <- rnorm(400,mean=mu)
xmpldata <- data.frame(time=time,treatment=treatment,response)
mod <- lm(response~time*treatment,data=xmpldata)
yhat <- fitted(mod)
xtb <- with(xmpldata,xtabs(yhat ~ time+treatment))
print(xtb)
> treatment
> time Control Treatment
> Post 94.10501 201.99112
> Pre 101.63792 210.04248
Is that (something like) what you want?
cheers,
Rolf Turner
P.S. You said: "I usually use GUI software". Now *there* lies your
problem. A GUI is a black box that removes all control over what is
going on, from your hands.
R. T.
--
Honorary Research Fellow
Department of Statistics
University of Auckland
Stats. Dep't. (secretaries) phone:
+64-9-373-7599 ext. 89622
Home phone: +64-9-480-4619
______________________________________________
R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.
