similar to: How to get p-values, seperate vectors of regression coefficients and their s.e. from the "yags" output?

Dear R-Users: I am trying to find the robust (or sandwich) estimates of the standard error of fixed effects parameter estimates using the package "lmer2". In model-1, I used "robust=TRUE" on the other, in model-2, I used "robust=FALSE". Both models giving me the same estimates. So my question is, does the robust option works in lmer2 to get the robust estimates of

Hello R-users: Could you please tell me how can I extract the "df (degree of freedom)" and "estfun (estimating functions)" for the following lmer (or lmer2) model? wtd.mixed<-lmer(ddimer~race+steroid+psi+sofa+apache + (1|subject), method="ML", data=final, cluster="id", weights=w) I tried the following codes: - for the degree of freedom (erorr

Hi A quick question. Standard errors reported by gee/yags differs from the ones in geeglm (geepack). require(gee) require(geepack) require(yags) mm <- gee(breaks ~ tension, id=wool, data=warpbreaks, corstr="exchangeable") mm2 <- geeglm(breaks ~ tension, id=wool, data=warpbreaks, corstr="exchangeable", std.err = "san.se") mm3 <- yags(breaks ~

Here is an example given from ?yags library(methods) data(stackloss) Y1 <- yags(stack.loss~Air.Flow,id=1:21, data=stackloss) How can I access parts of the output. I tried: > str(Y1) Formal class 'yagsResult' [package "yags"] with 25 slots ..@ coefficients : num [1:2] -44.13 1.02 ..@ coefnames : chr(0) > Y1$coefnames Error in Y1$coefnames : $

Dear R users, I search the R archives and noted that the same problem has been posted but without solution. I know there IS instructions by the author of yags, but I just couldn't figured out. I know gee and geepack can also perform generalized estimating equation, but the reason why I need yags is I want to perform a Small Sample Adjustments by Michael P. Fay and Barry I. Graubard in

Dear R-ers, I am just a simple 'end-user' of R and am trying to analyse data with a binary response variable (dead or alive) in relation to weight and sex (of young birds). As some of the birds have the same biological mother, I am using mixed models with the identity of the mother as a random factor. (please, Mick Crawley, when are you going to write a chapter on mixed models with binary

Dear R-ers, I am just a simple 'end-user' of R and am trying to analyse data with a binary response variable (dead or alive) in relation to weight and sex (of young birds). As some of the birds have the same biological mother, I am using mixed models with the identity of the mother as a random factor. (please, Mick Crawley, when are you going to write a chapter on mixed models with binary

I have included data at the bottom of this email. It can be read in by highlighting the data and then using this command: dat <- read.table("clipboard", header = TRUE,sep="\t") I can obtain solutions with both of these: library(gee) fit.gee<-gee(score ~ chem + time, id=id, family=gaussian,corstr="exchangeable",data=dat) and library(yags) fit.yags <-

i'm trying to learn R. i have no programing experience. i tried to add on a package and just picked "yags". i can't get it to work or get any help. --- > install.packages("yags") Warning in install.packages("yags") : argument 'lib' is missing: using 'C:\Documents and Settings\j\My Documents/R/win-library/2.6' trying URL

El 15/11/17 a las 3:11, david escribi?: > Folks > > I have a Centos7 system (SOFA) and want to install a Samba share named > "STUFF" for the machines inside my home.? All users in my home have > read access to the share, but only one user "me" has write > permission.? The configuration below worked just fine when the Samba > system was on Centos 6, but

> -----Original Message----- > From: CentOS [mailto:centos-bounces at centos.org] On Behalf Of david > Sent: Tuesday, November 14, 2017 9:11 PM > To: CentOS mailing list <centos at centos.org> > Subject: [CentOS] Samba help > > Folks > > I have a Centos7 system (SOFA) and want to install a Samba share > named "STUFF" for the machines inside my home.

Folks I have a Centos7 system (SOFA) and want to install a Samba share named "STUFF" for the machines inside my home. All users in my home have read access to the share, but only one user "me" has write permission. The configuration below worked just fine when the Samba system was on Centos 6, but did not work under Centos 7. The client machine is Windows 10. I have

I have read in a file (call it myData). The actual file is about 3000x30,000 columns and object.size() says myData takes: > 737910472/(1024^2) [1] 703.7263 Unfortunately, my program ends up using 40GB as indicated by maxvmem on Unix, which causes my department's cluster to stop working. Perhaps, I have some copying going on that I cannot find. I have created an example below that mimics

Full_Name: Markus Jantti Version: 1.5.1 OS: GNU/Linux Unstable (Sid) Submission from: (NULL) (193.166.59.44) I was trying to get the package yags 3.20 by Vince Carey <stvjc@channing.harvard.edu> to run, but - after successfull installation -- loading the library failed with the message > library(yags) Error in stop(paste("Duplicated slot names: ", paste(slots[duplicated(slo

markus.jantti@iki.fi wrote: > > Full_Name: Markus Jantti > Version: 1.5.1 > OS: GNU/Linux Unstable (Sid) > Submission from: (NULL) (193.166.59.44) > > I was trying to get the package yags 3.20 by Vince Carey > <stvjc@channing.harvard.edu> > to run, but - after successfull installation -- loading the library failed with > the > message > > >

En el correo anterior se me olvidó mencionar que trabajo con Rstudio El 28 de mayo de 2015, 15:18, MªLuz Morales <mlzmrls en gmail.com> escribió: > Hola, > soy nueva en esta lista y también en R. Yo he realizado un script en R > que carga dos archivos csv, uno de ellos con casi 2 millones de filas. El > programa carga esos archivos a data frame, y se trata simplemente de >

Hola, quiero compartir con vosotros mi problema y la solución que me han planteado. Mi programa carga Outcomes.csv y Set-A.csv (descargados de http://garrickadenbuie.com/blog/2013/04/11/visualize-physionet-data-with-r/, apartado Getting Started --> the code and the data set) de unos 50MB entre los dos. Mi código era: # Transforma csv a data frame seta <- read.csv('Set-A.csv');

Hola, Si no tienes inconveniente en compartir tu conjunto de datos (puedes dejarlo en un Dropbox y compartir enlace) o incluir una salida de la variables: "seta" y "outcomes" (función "save.image()") con eso podemos darte alguna solución mucho más rápida que la que planteas. En tu código con un bucle estás tratando de rellenar una lista que son los diferentes

I need capabilities, for my data analysis, like the Pinheiro & Bates S-Plus/R package nlme() but with binomial family and logit link. I need multiple crossed, possibly interacting fixed effects (age cohort of twin when entered study, sex of twin, sampling method used to acquire twin pair, and twin zygosity), a couple of random effects other than the cluster variable, and the ability to

I just installed R 2.12.1, and when I went to run a few old programs with it, nothing worked. I got a ton of error messages saying such and such package was built before R 2.10.0 and needed to be reinstalled. These were not just warning messages, but error messages that prevent the programs from running when they were running just fine with R 2.10.1 For some of those packages, such as deSolve, I