similar to: plot groupedData in nlme

Dear all, I've got a problem in including square variables in lme function. I've tried to work on Dialyzer data of Pinheiro and Bates'book. We fit the heteroscedastic model with: > data(Dialyzer) > fm2Dial.lme<-lme(rate~(pressure+pressure^2+pressure^3+pressure^4)*QB, + Dialyzer,~pressure+pressure^2,weights=varPower(form=~pressure)) We Obtain > fm2Dial.lme Linear

I met today a computer crash and our maintenance officer had to reinstall some components of the OS (MS Windows XP Pro) as well as the Internet browser (among other things). Now, I cannot install packages from a distant repository: > utils:::menuInstallPkgs() Error in .readRDS(pfile) : unknown input format > traceback() 5: .readRDS(pfile) 4: .packages(all.available = TRUE) 3:

  Hola,   Por favor, ¿podríais indicarme qué recursos (librerías o ideas) pueden resultar de utilidad para crear un gráfico del estilo del de la figura 3.8 del siguiente link?   http://www.tsc.uvigo.es/BIO/Bioing/ChrLDoc3.html#3.5   Actualmente estoy utilizando funciones muy básicas y la verdad es que no me encuentro muy satisfecha con el resultado.   Muchas gracias.   Eva [[alternative HTML

Hi, I have not been able to calculate residuals in the lme4 package. I've been trying the resid() function after I ran a GLMM with the lmer() function, but I get an error message that says "residuals are not inserted yet". I looked it up in the "help" history and I realized that several people have had this problem in the past, related to some bug in this function and

Hi, I am very intrigued by the idea of integrating statistical analysis directly with a paper as Sweave does it. But as I am collaborating with several people and they won't have set up R and also they're very unlikely interested in learning it, I am concerned about how such an Latex/Sweave document will work on their side. Can they Latex-compile the document without having R / Sweave

How to add mean,sd, number of observation in each panel for lattice histogram? Aimin

Dear all, I've just downloaded R.12.0 from CRAN (http://cran.r-project.org/bin/windows/base/release.htm). I am unable to find RGUI.exe as well as Rterm.exe. Conversely to what is explained in readme.R.2.12.0: There are two versions of the R executable in R-2.12.0\bin\i386 (32-bit) or R-2.12.0\bin\x64 (64-bit). I can't find directories i386 or x64 in R-2.12.0\bin -- Renaud Lancelot EDEN

Hi, The standard error from logistic regression is slightly different from the naive SE from GEE under independence working correlation structure. Shouldn't they be identical? Anyone has insight about this? Thanks, Qiong a<-rbinom(1000,1) b<-rbinom(1000,2,0.1) c<-rbinom(1000,10,0.5) summary(gee(a~b, id=c,family="binomial",corstr="independence"))$coef

Dear all, I have installed the latest version of R 2.12.0 available on CRAN (http://cran.r-project.org). When I try to load the recommended package lattice: > library(lattice) Error: package 'lattice' is not installed for 'arch=i386' I am running Rgui using C:\R\R-2.12.0\bin\i386\Rgui.exe --vanilla > sessionInfo() R version 2.12.0 (2010-10-15) Platform:

windows XP R 2.8.1 I would like to learn to use Sweave with R. I have no experience using Latex. I have looked at the sweave manual (http://www.statistik.lmu.de/~leisch/Sweave/Sweave-manual.pdf) and have found it difficult to understand, probably because I don't know Latex. Can someone recommend a document that will get me started with Sweave? Thanks John Confidentiality Statement: This

Hi, I could not use 'exlcude=' option in factor() to exclude a level from a existing factor. x is a factor: > x [1] a b c Levels: a b c > factor(x,exclude="c") [1] a b c Levels: a b c Warning message: NAs introduced by coercion However, "c" is not coded as NA. The following does not work either: >

I've been using par() to check the graphics parameters associated with both plot(<fitted linear model>) and plot(<grouped data object>). AFAIK the only differences are in the $cxy, $usr, $xaxp, and $yaxp parameters but the background color for the grouped data plot is grey while the linear model plot has a white background. When I've tried par(bg = "white") prior to

Is it possible to read jpeg files into R? If yes please guide, Thanks.. I tried to search many time but failed to do. Thankis in advance.. with Best Regards, Malik Shahzad Visiting Researcher National Institute of Informatics (NII) Tokyo, Japan Doctoral Student Asian Institute of Technology (AIT) Bangkok, Thailand +66-8-7676-5616 [[alternative HTML version deleted]]

Dear all, I am new to R. I wish to use R's multiple imputation to deal with missing data. I have a data set with the size around 300 observations and 150 variables. I checked the help function in R and could not locate how to write the codes for this. can anyone give a hand? Do appreciate your time and kindness! Q. [[alternative HTML version deleted]]

Hi, I'm running CentOS 6.5 as NFS server (v3 and v4) and exporting Ext4 and XFS filesystem. After many months that all works fine today the server crash: Jan 30 09:46:13 qb-storage kernel: ------------[ cut here ]------------ Jan 30 09:46:13 qb-storage kernel: WARNING: at lib/list_debug.c:26 __list_add+0x6d/0xa0() (Not tainted) Jan 30 09:46:13 qb-storage kernel: Hardware name: PowerEdge

I did testing on generating this report from a 130 MB Quickbooks file. Locally, 15 seconds. Samba Server, 1:35. Shared off of my win98 Workstation, 2+ mins. I do not understand the reason it is SOO slow. I know it is a fairly large file. The only explanation I can come up with is the way QuickBooks is accessing the files. Maybe jumping all around the file, requiring many operations on the

This is set up with Win2K boxes sharing a QB file on a SAMBA share. It may be quite different than with 95-98x-ME(Miserable Edition). There are two or more problems - one thing I identified is that it is possible for two computers to connect to the QB(QuickBooks) share using the same quick books logon - something that will cause data corruption. (Something that could have been easily

I have spent much of the day today researching performance tuning with samba. I have tried everything that I can find out about how to make performance faster. I checked disk performance with Bonnie, installed FTP and tested a transfer that way, achieving 6-7 MB / second. about 30 seconds for 150 MB file. I was finally able to achieve those speeds on a file transfer to the server through

help.search("groupedData") says that it's part of the lme4 package, but it appears not to be there (details below). Is this because lme4 is new and (perhaps) still under development? > update.packages() trying URL `http://cran.r-project.org/bin/windows/contrib/1.9/PACKAGES' Content type `text/plain; charset=iso-8859-1' length 19113 bytes opened URL downloaded 18Kb >

I have a directory shared out via Samba for Quickbooks and seem to have some issues with permissions. The directory being shared is a subdirectory in an ext3 partition being mounted with the acl option. It has been setup as follows: chown root:DOMAIN\AD_Group /mnt/Intuit_Data/ chmod 2770 /mnt/Intuit_Data/ And the Samba share config is has: create mask = 0660 directory mask =