Thanks much Dr. Ripley. Looks like the setMethod
statment was unneeded. If anyone can help, Im still
looking for some insight into the error message I
receive using nlme:
Warning in conLin$Xy * varWeights(object) : longer
object length is not a multiple of shorter object
length Error in
Hello, I am trying to design new variance structures
by using fixed effects variables in combination with
the VarPower function. That is, I would like to
create and evaluate my own variance function in the
data frame and then incorporate it into the model
using varPower, with value=.5.
As a
Thanks much, Dr. Bates. That helps. Would you have
any comment on the other part of my other question
regarding the manual recreation of varConstPower
function using new fixed effects? Im guessing I
obtain a singularity error because the additional
fixed effects may not alter the return value
Ah, I think Ive found a solution. A variance
function of the form:
(d1 + (covariant)^d2)^2 + U,
can be made by altering the line in varWeights.varComb
to read:
apply(as.data.frame(lapply(object, varWeights)), 1,
sum)
Then the function varComb can be used as in:
varComb(varPower(form =
Hello:
I would like to flag certain warnings to make them
into errors. I thought I could do this by giving
warning.expression a function to evaluate, in which I
could check the warning message and call stop if I
want.
I tried this:
options(warning.expression=expression(myfunction()))
but it
Hello:
I recently asked the list a question regarding
warning.expression. This is a different statement of
the problem.
I am doing a large simulation experiment using nlme
and most of the data realizations run fine. The
simulations stall, however, on a few particularly
noisy data sets. What
Excellent Gabor. Thanks. I had seen the
withCallingHandlers function, but I did not
understand how it was to be used. Thanks for the
example.
John
__
[EMAIL PROTECTED] mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the
Hello. Im hoping someone can help with a grouping
question related to the random= statement within the
nlme function. How do you specify that some grouping
levels are inner to others? I tried several things,
given below.
Lets say I have a data frame with five variables,
resp, cov1, ran1, ran2,
: The specified module could not be found.
Error in fun(...) : Could not load the rggobi library - please ensure
GGobi is on the library path
Error: .onLoad failed in 'loadNamespace' for 'rggobi'
Error in eval(expr, envir, enclos) : could not find function ggobi
Need your help
JJ
Hi,
I have a question about the function lme() in R.
I have a 2*2*3 layout with some missing data (labelled as *). These 3
factors are labelled as A,B,C, the response is Score. The layout is as
follows:-
A B CScore
1 1
I run R 1.9.0 on windows 2000, and have the following libraries installed:
affydata_1.3.1
affy_1.4.23
Biobase_1.4.10
DynDoc_1.3.14
gcrma_1.0.6
hgu133acdf_1.4.3
hgu95av2cdf_1.4.3
hgu95av2probe_1.0
matchprobes_1.0.7
moe430acdf_1.4.3
multcomp_0.4-6
mvtnorm_0.6-6
rae230acdf_1.4.3
reposTools_1.3.29
11 matches
Mail list logo