¦ but I am probably missing something]?
Many thanks
Andrew
--
Andrew Beckerman
University of Sheffield
[[alternative HTML version deleted]]
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the p
Dear all -
Should anyone be interested - using
paste(system(paste("awk 'NR==2'", file.genepop), intern = TRUE), collapse="")
pulls together the pieces (list elements) of the returned values.
--
Andrew Beckerman
University of Sheffield
On Tuesday, 17 S
,method="REML") # multiple
# example matrix for predicting against model3
newdat<-expand.grid(c(1,2,3,4,5),c(1,2,3))
predict(model1,newmods=c(1,2,3,4,5))
predict(model2,newmods=c(1,2,3))
predict(model3,newmods=newdat)
--
Andrew Beckerman
Sent with Sparrow (
8
18 0.8397355 8
19 0.8397355 8
20 0.7203888 9
-----
Dr. Andrew Beckerman
Department of Animal and Plant Sciences, University of Sheffield,
Alfred Denny Building, Western Bank, Sheffield S10 2TN, UK
ph +44 (0)114 222 0026; fx +44 (0)114 222 0002
http://www.beckslab.staff.she
rew
-
Dr. Andrew Beckerman
Department of Animal and Plant Sciences, University of Sheffield,
Alfred Denny Building, Western Bank, Sheffield S10 2TN, UK
ph +44 (0)114 222 0026; fx +44 (0)114 222 0002
http://www.beckslab.staff.shef.ac
))
should avoid the error and give a conservative estimate of the
significance of your interaction.
see also:
http://tolstoy.newcastle.edu.au/R/e2/help/06/10/3565.html
and related posts.
A
-----
Dr. Andrew Beckerman
Dear all -
I would like to adjust the line type of specific contours in
contourplot from the lattice package, but it seems like lty does not
take a list in the call.
Here is my call to contourplot:
contourplot(preds~size+trt|Size.Name,
data=pred.dat,layout=c(2,4),
at=c(0.02
at[seq(1, length(at), 2)]
at.e <- at[seq(2, length(at), 2)]
panel.contourplot(at = at.o, lty = 1, col = "blue", ...)
panel.contourplot(at = at.e, lty = 4, col = "red", ...)
},
at = pretty(grid$z, 10))
HTH,
--sundar
evices
methods base
other attached packages:
[1] survival_2.34-1 Hmisc_3.4-3 lattice_0.17-15 MASS_7.2-44
loaded via a namespace (and not attached):
[1] cluster_1.11.11 grid_2.8.0 tools_2.8.0
Andrew
-------
newdata to get just a single estimate? In the above
example I used unique against the model predictions evaluated at every
datapoint. But with no terms on the RHS of the model, how do you (can
you) specify new data?
Best wishes,
Andrew
On 26 Nov 2008, at 16:27, Andrew Beckerman wrot
10 matches
Mail list logo