Re: [R] Re: creating new varFunc classes in nlme .. error: Don't know how to get coefficients for .. object

2004-10-07 Thread JJ
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

[R] manual recreation of varConstPower using new fixed effects variables in nlme

2004-10-17 Thread JJ
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

Re: [R] manual recreation of varConstPower using new fixed effects variables in nlme

2004-10-18 Thread JJ
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

Re: [R] manual recreation of varConstPower using new fixed effects variables in nlme

2004-10-18 Thread JJ
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 =

[R] question on using warning.expression

2004-12-20 Thread JJ
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

[R] hang-up during nlme call

2004-12-20 Thread JJ
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

Re: [R] hang-up during nlme call

2004-12-21 Thread JJ
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

[R] How to conctruct an inner grouping for nlme random statement?

2005-02-22 Thread JJ
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,

[R] rattle()-RData-Explore-GGobi-libggobi-0.dll--Error

2007-03-08 Thread JJ
: 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

[R] How to specify arguments in lme() ?

2007-01-18 Thread w 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

[R] rma and gcrma do not work in R 1.9.0

2004-05-14 Thread Li, James [PRDUS Non-JJ]
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