similar to: R-help: gls with correlation=corARMA

Hi, I'm a new member here in the list. I am a graduate from University of Georgia. Recently in doing analysis using lme on a dataset, I found several questions: 1. How to express the equation when the correlation structure is very complicated. For exmaple, if the fixed is y(t)=0.03x1(t)+1.5x2(t)(I omitted "hat" and others). And the model with corARMA(p=2,q=3) is proper. What will be

I am currently trying to get a "lme" analyses running to correct for the non-independence of residuals (using e.g. corAR1, corARMA) for a larger data set (>10000 obs) for an independent (lgeodisE) and dependent variable (gendis). Previous attempts using SAS failed. In addition we were told by SAS that our data set was too large to be handled by this procedure anyway (!!). SAS script

Full_Name: Benjamin Tyner Version: 2.6.0 RC 2007-10-01 r43043 OS: WinXP Submission from: (NULL) (171.161.224.10) platform i386-pc-mingw32 arch i386 os mingw32 system i386, mingw32 status RC major

Trying following example from Pinheiro and Bates in order to fit an ARMA(1,1) model: library(nlme) fm1Ovary.lme<-lme(follicles~sin(2*pi*Time)+cos(*pi*Time),data=Ovary,random=p dDiag(~sin(2*pi*Time))) fm5Ovary.lme<-update(fm1Ovary.lme,corr=corARMA(p=1,q=1)) I get follwing error message: Error in "coef<-.corARMA"(`*tmp*`, value = c(62.3428455941166, 62.3428517930051 :

Dear all I am tryiing to fit the following lme with an ARMA correlation structure: test <- lme(fixed=fev1f~year, random=~1|id2, data=pheno2, correlation=corARMA(value=0.2, form=~year|id2), na.action=na.omit) But I get the following error message: Error in getGroupsFormula.default(correlation, asList = TRUE) : "Form" argument must be a formula I have used this same form

R-helpers, n <- 100 arcoefs <- c(0.8) macoefs <- c(-0.6) p <- length(arcoefs) q <- length(macoefs) require(nlme) tmp <- corARMA(value=c(arcoefs,macoefs), form=~1, p=p, q=q) Sigma <- corMatrix(tmp, covariate = 1:n) # results in segfault Have I used these commands in an improper way? Thanks Ben

Hi R-sters, Just wondering what I might be doing wrong. I'm trying to fit a multiple linear regression model, and being ever mindful about the possibilities of autocorrelation in the errors (it's a time series), the errors appear to follow an AR1 process (ar(ts(glsfit$residuals)) selected order 1). So, when I go back and try to do the simultaneous regression and error fit with gls,

Dear All: I am trying to fit a special case of a 2-banded Toeplitz correlation structure. A 2-banded Toeplitz has ones on the diagonal, a correlation, RHO1, on the first off-diagonal, and a correlation, RHO2, on the second off-diagonal, with zeros on all subsequent off-diagonals. After reading relevant sections in Mixed-Effects Models in S and S-PLUS (Pinheiro & Bates, 2000) and searching

I have a data analysis job for which lme may be used. Prof. Spencer Graves had helped me much on that. I'm really appreciated for that. Could anybody else in the list give me some hints from other perspectives? I hope I can learn as much as possible for this complicated real data. Thanks in advance. Hanhan To briefly describe my data: My data is health effect measurements (y) and personal

Dear R-help, I don't understand why Gls gives me an error when trying to fit a model with AR(2) errors, while gls (from nlme) does not. For example: library(nlme) library(rms) set.seed(1) d <- data.frame(x = rnorm(50), y = rnorm(50)) gls(y ~ x, data=d, correlation = corARMA(p=2)) #This works Gls(y ~ x, data=d, correlation = corARMA(p=2)) # Gives error # Error in

Hello, samba v3.0.7 on suse linux v9.1 os/2 v4.x I have started the process of replacing our os/2 hosts with linux, starting with the server(s). Most of the problems have been firewall related, a topic only briefly discussed in the docs (nudge). Some of the hosts can now interact with the linux samba shares; not all though. I have run ip traces and have gotten to a point that is a

I tried different versions of sma**.tar.gz file either with "R CMD INSTALL -l packagepath sma**.tar.gz", "R INSTALL -l path sma**.tar.gz" and under R windows by "install.packages("sma_0.5.15.tar.gz",repos=NULL)" . I got same error ====== * installing *source* package 'sma' ... ** Creating default NAMESPACE file Warning in .write_description(db,

hi I am going to crazy.I want to use beadarray package in R, but I cant. when I want install it, I encounter to this message "Loading required package: hwriter Loading required package: sma Error in library(pkg, character.only = TRUE, logical.return = TRUE, lib.loc = lib.loc) : 'sma' is not a valid installed package" I cant find sma and hwriter packages, it seems that they are

rsync v3.1.0 linux v4.4.104-39-default x86_64 Found in the system log: 2018-02-22T05:02:00-0700 sma-server3 python3[31371]: backintime (sma-user3x/3): WARNING: Command "rsync -rtDHh --links --no-p --no-g --no-o --info=progress2 --no-i-r --delete --delete-excluded -i --dry-run --out-format="BACKINTIME: %i %n%L" --chmod=Du+wx --exclude="/bkp/cgate-backintime"

HELLO ALL, I AM GETTING AN ERROR MESSAGE WHEN TRYING TO RUN A GAMM MODEL LIKE THE ONE BELOW. I AM USING R VERSION 2.14.1 (2011-12-22) AND MGCV 1.7-12. M1 <-gamm(DepVar ~ Treatment + s(Year, by =Treatment), random=list(Block=~1), na.action=na.omit, data = mydata, correlation = corARMA(form =~ Year|Treatment, p = 1, q = 0)) THIS IS THE ERROR MESSAGE Error in `*tmp*`[[k]] : attempt to

-----BEGIN PGP SIGNED MESSAGE----- Hash: SHA1 Hello, I have an iMac v10.6. I'd like to get it to connect to an old os/2 server, v4.5 fp3. Using smbutils on the iMac: $ smbutil -v view //sma-server1.sma.com smbutil: server connection failed: RPC struct is bad Can anyone suggest what changes may be made that will correct this error? Is it even possible? - -- James Moe moe dot james at

Hi all I have just downloaded the sma package from CRAN. On installing on my linux machine, I get the message > library(sma) > data(MouseArray) Warning message: file 'MouseArray.RData' has magic number 'RDX1' Use of save versions prior to 2 is deprecated Hereafter, MouseArray is not found: > MouseArray Error: object "MouseArray" not found We were

-----BEGIN PGP SIGNED MESSAGE----- Hash: SHA1 Hello, linux 2.6.34.7-0.7-desktop x86_64 smbclient 3.5.4-5.1.2-2426-SUSE-SL11.3 I issue this command: smbclient -L SMA-STN14L -U jmoe I get: Connection to SMA-STN14L failed (Error NT_STATUS_BAD_NETWORK_NAME) In <smb.conf> [globals] section: workgroup = SOHNEN-MOE netbios name = SMA-STN14L Adding the -I option made no difference.

openSUSE v13.1 linux 3.11.10-11-desktop x86_64 samba 4.1.6 I recently upgraded an openSUSE server from v12.3 to v13.1. In the older version there was no problem about mounting a remote volume offered by an ancient OS/2 system. Now there is. Trying to mount it by command line: root:/home/sma-user3x> mount /t2 Retrying with upper case share name mount error(6): No such device or address

Probably a very simple query: When I try to plot a curve from a fitted polynomial, it comes out rather jagged, not smooth like fitted curves in other stats software. Is there a way of getting a smooth curve in R? What I'm doing at the moment (for the sake of example) is: > x <- c(1,2,3,4,5,6,7,8,9,10) > y <- c(10,9,8,7,6,6.5,7,8,9,10) > b <- data.frame(cbind(x,y)) >