similar to: Bug in getVarCov.gls method (PR#9752)

Hello R-users, I am using gls function in R to fit a model with certain correlation structure. The medol as: fit.a<-gls(y~1,data=test.data,correlation=corAR1(form=~1|aa),method="ML") mu<-summary(fit.a)$coefficient With the toy data I made to test, the estimate of mu is exactly equal to the overall mean of y which can not be true. But, if I make a toy data with y more than two

Dear List: I am estimating a gls model and am having to make some rather unconventional modifications to handle a particular problem I have identified. My aim is to fit a GLS with an AR1 structure, obtain the variance-covariance matrix (V), modify it as needed given my research problem, and then reestimate the GLS by brute force using matrix operations. All seems to be working almost perfectly,

Dear List: I am having some trouble reproducing some GLS estimates using matrix operations that I am not having with other R procedures. Here are some sample data to see what I am doing along with all code: mu<-c(100,150,200,250) Sigma<-matrix(c(400,80,16,3.2,80,400,80,16,16,80,400,80,3.2,16,80,400),n c=4) sample.size<-100 temp <-

Hi, all: Two examples are shown below. I want to use the multiple comparison of Dunnett. It succeeded in upper case "example 1". However, the lower case "example 2" went wrong. In "example 2", the function pmvt return NaN, so I cannot show this simtest result. Is there any solution? (I changed the variable "maxpts" to a large number in front of the

Dear R-users I am relatively new to R, i hope my many novice questions are welcome. I have problems accessing some objects (specifically the random effects, correlation structure and variance function) from an object of class gls and lme. I used the following models: yah <- gls (outcome~ -1 + as.factor(Trial):as.factor(endpoint)+

Dr. Andrzej Galecki will present his online course "Mixed Effects Models with Applications" May 12 ? June 9 at statistics.com. This course will explain the basic theory of linear and non-linear mixed effects models, including hierarchical linear models (HLM). It will outline the algorithms used for estimation, primarily for models involving normally distributed errors, and will provide

Rusers, I am trying to build R2.5 on my Kubuntu machine an configure concludes with : configure: WARNING: you cannot build info or HTML versions of the R manuals What do I need to install on my machine to get rid of this warnings? Eryk [[alternative HTML version deleted]]

Dear useR's, I have a problem with the lattice plotting of some symbols: library(lattice) test<-data.frame(x=c(2,3,1,5),u=c(rep(1,2),rep(2,2)),g=c(rep(c(1,2),2))) xyplot(x~u,groups=g, data=test, par.settings=list( superpose.symbol=list(pch=c(22, 23),cex=c(1.7,1.6),col="black") ), key=list( text=list(c("t1","t2")), space =

When I compile Samba, this is on SuSE 7.2, I keep getting a make failure, but it still completes to a prompt. Should I ignore or does anyone have a "fix"? When I configure, I use these options: ./configure --enable-cups --with-smbwrapper --with-automount --with-smbmount --with-pam --with-pam_smbpass --with-winbind --with-winbind-auth-challenge When I run make, this is the error,

>>>>> "VW" == Viechtbauer Wolfgang (STAT) <Wolfgang.Viechtbauer at STAT.unimaas.nl> >>>>> on Fri, 26 Sep 2008 18:00:19 +0200 writes: VW> Hi all, I'm trying to fit a marginal (longitudinal) VW> model with an exponential serial correlation function to VW> the Orange tree data set. However, R crashes frequently VW>

On 9/29/2020 14:37, David Chisnall via llvm-dev wrote: > On 28/09/2020 15:45, Matt P. Dziubinski via llvm-dev wrote: >> Hey, I noticed this talk from the EuroLLVM 2015 >> (https://llvm.org/devmtg/2015-04/slides/pslp_slides_EUROLLVM2015.pdf) >> on the PSLP vectorization algorithm (CGO 2015 paper: >> http://vporpo.me/papers/pslp_cgo2015.pdf). >> >> Is anyone

Dear Xen users! In case you were struggling with recent version of th gplpv drivers, I''ve possibly found an LVM related bug in it, and reported to James Harper. Thats''s the point: I''ve tested the mentioned version of your driver on top of 32bit xen-3.2.1 32bit xen-3.3.0 and 64bit xen xen-3.3.0 hypervisors, dom0 was always an appropriate version of a 32bit PAE kernel,

I wonder if it's the same part as used on the VW Samba! J On 04/01/2019 21:55, Steven Hirsch via samba wrote: > On Fri, 4 Jan 2019, Susan Slayter King via samba wrote: > >> I crashed my car and need the entire turn indicator assembly. Drivers >> side. >> Can anybody help. No local retailers have it. I did get a wiring tube >> but >> still need the seal,

I want to extract the variance-covariance matrix of an nlme model of a dataset. The object is to pass this to mvrnorm to create pseudo- replicates of the original data. I note the nlme package has a getVarCov method available for lme objects but not nlme objects. Is the vcov function in the base stats package suitable? If so, why is the additional getVarCov provided? thank you Rob

Dear all, I noticed the following in the call of lme using msVerbose. fm1 <- lme(distance ~ age, data = Orthodont, control = lmeControl(msVerbose=T)) 9 318.073: -0.567886 0.152479 1.98021 10 318.073: -0.567191 0.152472 1.98009 11 318.073: -0.567208 0.152473 1.98010 fm2 <- lme(distance ~ age, random =~age, data = Orthodont,

Thank you in advance for reading this help request. I am pretty new to R. I am experiencing some issues getting 2.5 installed on my Ubuntu Fiesty system and seek your advice. To the best of my ability I followed the instructions here: http://cran.r-project.org/bin/linux/ubuntu/README Setting this as the last line in my sources.list: deb http://cran.fhcrc.org/bin/linux/ubuntu feisty/ When I

Hello Hello I have a data frame, x, holding 5 persons answering the question which cars they have used: # the data frame x <- as.data.frame( matrix( c('BMW', '', '', 'Mercedes', 'VW', '', 'Skoda', 'VW', 'BMW', '', '', '', 'VW', 'Skoda',

Or at least post your smb.conf, we tend to base our diagnosis on it. On 07.01.2019 11:01, Sven Schwedas via samba wrote: > It might actually, Karmann-Ghia, Beetle and Samba all shared the same > engine blocks, from what I recall. > > …it's probably still better to ask a car mechanic. > > On 04.01.19 23:04, A. James Lewis via samba wrote: >> I wonder if it's the

Debian r2.5 (Potato) rsync version 2.3.2 protocol version 21 ( the latest stable deb version anyway :) RE: absolutely older file keeps getting chosen for update Hi. I have something like: rsync -vxuaz /sendmail/ /mnt/nfs-mount/sendmail In the nfs drive the files are links to other files in a subdirectory of that same directory. Both the links and the files they point

Hi, I'm slowly working through Tsay's "Analysis of Financial Time Series" 3rd ed. ?I'm trying to replicate Table 2.1 on p.47, which gives PACF, AIC, and BIC for the monthly simple returns of the CRSP value-weighted index. The data: http://faculty.chicagobooth.edu/ruey.tsay/teaching/fts3/m-ibm3dx2608.txt > da <-