David,
It would be helpful to give an example of what you would like to extract.
I guess you know how to extract elements from vectors and lists.
However, sometimes the objects returned by functions can be rather complex
(output of coxph() is...)
A general method to capture printed output is vi
What about the other way round?
x <- 1:10
y <- rnorm(10,10,1)
x2 <- 3*x + 2
plot(y ~ x2, xaxt = "n")
axis(side=1,at = x2, labels = x)
Florent Bresson
x <- 1:10
y <- rnorm(10,10,1)
x2 <- 3*x + 2
plot(y ~ x, xaxt = "n")
axis(side=1,at = x, labels = x2)
Joris
Florent Bresson
<[EMAIL PROTECTE
Martin,
Have you checked ?intervals.gls
This intervals are approximate, but this would be the obvious starting
point to me.
Joris
"Martin Henry H.
The AMMI senso strictu part starts from the corrected data for genotype and
environment. In most cases where AMMI is applied (maybe also in the
agricolae package, I haven't checked), starts from the interaction effects
obtained through a general linear model anova.
It should be possible to replac
Mustapha,
You should tell us what fm1 actually is (which class), before we can give
you a correct answer.
It is helpfull to use str() in this case.
str(fm1) will give you the structure of fm1, and with that information you
extract any value you want.
I suppose that this is what you need:
summ
or
Cat <- c('a','a','a','b','b','b','a','a','b')
C1 <- (Cat=='a')*1
"ONKELINX,
Thierry"
gls() from the package nlme is similar to lm but is meant for models
without random effects.
Joris
Spencer Graves
<[EMAIL PROTECTED]
Guilio,
Have a look at Rnew volume 5/1 (http://cran.r-project.org/doc/Rnews/) where
Doug Bates explains this nicely. Condider using lme4 for your purpose.
Joris
"Giulio Di
to backtransform 'estimate':
if (lambda == 0 ) {
log(estimate)
} else {
estimate^(1/lambda)
}
"Des Callaghan"
<[EMAIL PROTECT
mod <- lmer(resp3~b$age+b$size+b$pcfat+(1|sex), data=b)
coef(mod)[1]$Subject[1,1]
[EMAIL PROTECTED] wrote on 27/02/2007 16:52:54:
> Hi,
>
> I have been dabbling with str() to extract values from outputs such as
> lmer etc and have found it very helpful sometimes.
>
> but only seem to manage to
Romain,
Look for info on mixed models.
In R you do this either with the library nlme or lme4.
A good starting point is an article by Doug Bates in Rnews
http://cran.r-project.org/doc/Rnews/Rnews_2005-1.pdf
Joris
[EMAIL PROTECTED] wrote on 13/02/2007 17:16:41:
>
> Hello,
>
> Does somebody cou
"R --vanilla" ?
>
> Regards,
> Martin Maechler, ETH Zurich
>
> >>>>> "joris" == joris dewolf <[EMAIL PROTECTED]>
> >>>>> on Sun, 14 Jan 2007 19:15:43 +0100 writes:
>
> joris> Dear all,
>
> joris>
Dear all,
I try to run the example of lmer and get the following error message.
> library(lme4)
> example(lmer)
lmer> (fm1 <- lmer(Reaction ~ Days + (Days | Subject), sleepstudy))
[[1]]
Error in get(x, envir, mode, inherits) : variable "as.dpoMatrix" was not
found
This error message is similar
I advice you strongly to use VarCorr(), but if you insist
tmp <- as.matrix(m1$modelStruct$reStruct$MI)
c(sqrt(diag(tmp)), Residual = 1) * m1$sigma
Joris
[EMAIL PROTECTED] wrote on 02/01/2007 12:50:05:
> Here is a piece of code fitting a model to a (part) of a dataset, just
> for
> illustration
Hello,
Has somebody an idea of how to fit contrast involving random terms and
obtain their standard errors with lme?
I am referring to
Welham, S.; Cullis, B.; Gogel, B.; Gilmour, A. and Thompson, R. (2004).
Prediction in linear mixed models.
Australian and New Zealand Journal of Statistics 46, 32
Hi,
I have a lattice xyplot that contains panels according to FactorA, and
curves for the 2 levels of Factor B within a panel.
I try to add text in the panels of a lattice graph. I suppose I have to
write a custom function (panel.txt).
What I really would like is to adapt the text in the panel ac
Could anybody explain this difference in the function order() between R
under windows and R under Linux ?
It seems that under Linux the order in of character dsitinguishes
between caps and normal letters and sorts them starting with the
capitals (see below).
How can I avoid this (I mean: get the
This is not a PS problem. The barplot you made has its horizontl axis
suppressed by default. See argument axis.lty in ?barplot
Joris
Frank Gerrit Zoellner wrote:
Hi all!
I have a little problem with saving plots to file.
I use the command postscript() followed by the plotting command and a dev.
check rbind
Joris
John Sweval wrote:
I have a series of data frames that are identical structurally, i.e. -
made with the same code, but I need to add them together so that they
become one, longer, data frame, i.e. - each of the slot vectors are
increased in length by the length of the added data
20 matches
Mail list logo
