similar to: corCompSymm in nlme package

Dear R Friends, I need help with an error associated with corCompSymm in an lme fit. I am using a mixed effects model to analyze a split-plot with repeated measures and would like to fit with the compound symmetry correlation structure. This problem doesn't occur when using corAR1 or any of the other structures. I would greatly appreciate help on how to solve this issue. Here's my

Dear R-devel members, While writing a new correlation structure similar to 'corCompSymm' and intended to be used with 'gnls', I got puzzled with the 'Initialize' method. Using 'Initialize' before 'gnls' may be regarded as a mean to set an initial value for the corStruct parameter. However 'gnls' does not work properly with a

Dear Friends, I am seeking for any help on an error message in lme functions. I use mixed model to analyze a data with compound symmetric correlation structure. But I get an error message: "Error in corMatrix.corCompSymm(object) : NA/NaN/Inf in foreign function call (arg 1)". If I change the correlation structure to corAR1, then no error. I have no clue how to solve this problem.

Hi, I have confirmed temporal correlation problems in my data. Is there a possibility to use corCompSymm for a gamm()? I am an R-beginner. I have very short time series. There are three years and within each year, there are 10 weeks. he 10 weeks are the same every year and have not unique values, I seem not to be able to use AR-1 (I assume that I have too little data for autoregression

Dear R-list, Realising the importance of updating the R base, I tried to compile the source file of R2.2.1 in my machine, but failed. According to "R for Mac OS X FAQ", i've got the basic dependencies which are gcc 3.3 and g77 3.4. ./configure is fine Here are the errors I got from make: >make ... gcc -I../../src/extra/zlib -I../../src/extra/bzip2 -I../../src/extra/pcre

Dear all, I just found that some of the packages are not able to load any more, after I installed R2.0.1 in my Mac, it even affects my old R1.8 installs. It gives me errors when I load packages that contains "myEnv" settings. such as: RMySQL, DBI, Rggobi, etc. But others that does not require "myENV" is all right, like tcltk that only calls the c functions. The errors

Hello and thanks in advance I am running a glm in R the code is as follows with residual diagnostic code below model4<-glm(Biomass~(Treatment+Time+Site)^2, data=bobB, family=quasi(link="log", variance="mu")) par(mfrow=c(2,2)) plot(model2) to test the effect of grazing exclusion of feral horses for a Phd with following factors: Treatment - 3 levels which are grazed

Sorry, I just posted the email below but realised I did not give a name or details, would it be possible to adjust before posting and send what is below, sorry again, first time user... From: Joanne Lenehan [mailto:jlenehan@une.edu.au] Sent: Tuesday, 1 December 2009 3:51 PM To: 'r-help@r-project.org' Subject: GLM Repeated measures test of assumptions: e.g. test for sphericity e.g.

Dear All -- I am trying to use within a little table producing code an anova comparison of two gls fitted objects, contained in a list of such object, obtained using nlme function gls. The anova procedure fails to locate the second of the objects. The following code, borrowed from the help page of anova.gls, exemplifies: --------------- start example code --------------- library(nlme) ##

Dear R developers, first of all let me join the chorus of congratulations for the release of R 1.0.0. Well, done! Unfortunately, I find it necessary to e-mail in a bug report regarding the `nlme' package. On my office machine I experience the following trouble: bossiaea:/opt/R$ R CMD check -c nlme Checking package `nlme' ... Massaging examples into `nlme-Ex.R' ... Running

I am using the gls function of the nlme package to analyze data sets of soil respiration which have the following design: 3 complete blocks x 5 sampling dates (time from fertilization) x 3 fertilization levels. The fertilization dates are equal for all subjects (blocks) but not periodical (-46, 10, 24, 53, 123 days from the event). The code that I've been using is: fit.csnC<- gls(dno.C

Hi, Well after trying to fix the assign commands in version 3.1.7 of your nlme package I still cannot get the plot functions to work. Basically I have no other packages installed (apart from those provided with R) so nlme is the first. my R version info is; platform sparc-sun-solaris2.7 arch sparc os solaris2.7 system sparc, solaris2.7 status

Hi there, I try to compare the mixed model package "lme" by Splus and R. I used the dataset "Ovary" and the following code assuming AR(1) model for the error term: lme(follicles ~ sin(2*pi*Time) + cos(2*pi*Time), data=Ovary, random = pdDiag(~sin(2*pi*Time) ) , correlation=corAR1() ) But I got different results! And then I used a simpler model: lme(follicles ~

I'm encountering excruciatingly slow load times for character vectors in R 2.6.0-- up to 30sec for a 15K file that contains a no-attributes character vector of length ~1e4 and object size ~0.5MB. In R 2.5.1, repeated loads of the same set of files are near-instantaneous. The problem is proving tricky to reproduce consistently from scratch, so I have attached the 3 files used in the examples

I am trying to use the Windows COM interface under R 2.0, and have encountered the following difficulties: - the package RDCOMClient installs, loads and works under R 1.9.1, installs under R2.0, but does not load or work under R2.0 - the package SWinTypeLibs does not install or load under either R 1.9.1 or under R2.0 For the moment I am concentrating on the RDCOMClient package as it seems to

Dear all, I am trying to define my own corStruct which is different from the classical one available in nlme. The structure of this correlation is given below. I am wondering to know how to continue with this structure by using specific functions (corMatrix, getCovariate, Initialize,...) in order to get a structure like corAR1, corSymm which will be working for my data. Thanks in advance.

Hello, I ran my simulations on the Unix verson of R1.9.1 and the Windows version of R2.0.1 on XP. I kept getting different values for the mean of the same column of the same matrix, and am perplexed. I would appreciate if anyone could help explain the difference? Here is a sample code: set.seed(7293) z1v <- rnorm(1000, mean=68, sd=13) z1v <- (z1v-mean(z1v))/sd(z1v) Using R1.9.1 on Unix,

Hello - I am a fairly new user of R. I currently have R1.7 installed. I tried to install the KernSmooth package but a message tells me I need R1.8. Is there a way to update to R1.8 while in R1.7 in the same manner that I can use "update.packages", or do I have to install R1.8 in a separate directory and re-install the many packages I have already installed in R1.7? Thanks in advance for

Has anyone managed to get this working? Here's what I did: I got the binary build for R2.0.1 from the Omegahat download page, and made a small change to the registerClassID function ( to make it use the right path to RDCOMServer.dll). Then I tried to replicate the simple TTest example from the same web site. The COM class definition and registration worked fine for me. To test this, I

Full_Name: Vidar Hjellvik Version: 2.0.1 OS: windows Submission from: (NULL) (82.134.28.194) I have an tcltk application that runs without problems in R1.9.1, but when I press the "run"-button (and other buttons as well) in R2.0.1, I get the message: "Error in function () : can't change value of a locked binding". I have another tcltk application that runs fine on