similar to: nlme installation/ loading problem

I am using the latest version of R on a Windows machine and get the following error when I try to initialize a correlation structure with the function corAR1 in NLME. This example is taken from the book of Pinheiro and Bates, so it should work. What is going wrong? > library(nlme) > data(Orthodont) > cs1AR1 <- corAR1(0.8, form= ~1 | Subject) > cs1AR1 <- initialize(cs1AR1, data =

Could somebody help with the definition of new correlation structure for use with a linear mixed-effects model (package nlme). Specifically, I want to define a Toeplitz type correlation structure, but due to my inexperience with programming in R, I feel a bit overwhelmed with the task at hand. I understand that you can start with a function like corAR1 as template, but I have no idea how I would

Is there an "easy" way to plot a "color key" next to a color image (with image() in graphics package)? The color key should also include a numerical scale, so that the colors can be cross-referenced with image intensity values. I see that levelplot has a facility for color keys, but with image it seems less straightforward, i.e. generating a color key is not an option in

I have data for several rings of a left heart chamber, and which I would like to display in concentric rings, with color-encoding of the values. Each ring corresponds to one slice through the heart, and the rings correspond to positions from the base to the apex of the heart as you move from the outermost ring to the innermost one. The data have a circular pattern. These types of displays are

Let me modify my question about user-defined covariance structures for LME models: Can somebody tell me how I can see the code for the definition of the correlation structures that come with the NLME package. Specifically I like to see the code for the functions coef, corMatrix, and intialize for any of the pre-defined correlation structures, and use this as a template to define a new correlation

When using read.spss (library: 'foreign') I get the following warning message: Warning message: E:/R4win/mesamri.sav: Unrecognized record type 7, subtype 13 encountered in system file. I don't see anything wrong with record #7 in the database I am trying to read in, but I suspect that the warning message does not refer to a specific record, but a "variable" type. Is this

I am using lmList, which takes 'Data' partitioned according to the levels of a grouping factor (gf) and individual 'lm' fits are obtained for each 'data'partition, using a model defined as in "lm". So my call to lmList looks something like: > mg.lis <- lmList(strmbf ~ age1c + gender1 + sysbp.clinic + + diabp.clinic + ldl1 + mets.total + bmi1c |

I have been looking at a draft of Doug Bate's book on mixed effects models (http://lme4.r-forge.r-project.org/book), and in chapter one (section "1.4.2 Matrices and Vectors in the Fitted Model Object", to be exact) he refers to the function env() for checking the environment of a fitted lmer object. This function does not appear under help, and I have not been able to find

I have a problem with nlme on Windows 2000, and I'm having a devil of a time determining whether the problem is with my computer or with something in R. I'm running v1.8.1 on a Dell Pentium III with 512MB of RAM and all of the recommended Windows 2000 updates applied. If I use Rterm, I can run analyses with NLME to my heart's content. But when I run Rgui, I encounter a floating

Hi all, I have a very high dimensional data and apparently there are several columns that contain similar information (some columns are equal). I want to form a matrix/data frame consisting of unique columns. Does anyone have an efficient way of getting out these columns. A small section of the data frame is given below. Thanks for helping. Stephen. > newdata [,1] [,2] [,3] [,4] [,5]

Hi and thanks to all R developers and helpers, I wanrt to take the 1 and 3 dimensions of this object to create a new one without the years and with the next functions to do some specific calculations: > str(PMpa) num [1:499105, 1:60, 1:12] 29.8 55.8 29.7 25.1 25 ... - attr(*, "dimnames")=List of 3 ..$ punt: NULL ..$ any : chr [1:60] "1950" "1951"

I see that feature freeze for a better than ever R has been announced. May I ask if there will be or could be "merge" function for data frames? I believe it would be very nice if this were "centrally" implemented (with regard to my humble steps to writing such a function and the related error and other precautions). I know from VR2R.pdf (courtesy Profs. Venables and Ripley)

Hello, could You please help: I am looking for a way to formulate test accuracy measures such as test sensitivity, specificity, predictive values, and correct classification rate using p/d/qnorm. The tests' primary values follow a bimodal distribution, which is modelled by a mixture of two normal distributions: p * dnorm ((x - u1) / s1) / s1 + (1 - p) * dnorm ((x - u2) / s2) / s2)

Diese Nachricht wurde automatisch von einem Virenschutzprogramm erzeugt. Es wurde ein Virus in Ihrer Mail mit dem erw?hnten Subject vom Thu, 26 Jun 2003 13:03:32 --0400 an herold at med.uni-marburg.de gefunden. Ihre Nachricht wurde gel?scht und nicht weitergeleitet. Bitte l?sen Sie dieses Problem und versuchen Sie es erneut. This messages was automatically created by a virus-scanning program.

R-help list, I ported some R code from MS Windows PC to sun-solaris. Both do queries on the same MySQL database. PC version of R is 1.6.1 and solaris version is 1.6.0. I use RMySQL_0.5-0 (on solaris) and RODBC 1.0-1 on PC to connect to MySQL. The PC version took about 65 mins. to run wheras the Solaris (SunOS 5.7) version took 375 mins. (> 6 hrs). R and RMySQL are resident on the same host

Hi all, is the Matrix package no longer available for download via install.packages? When I try to install it (from any mirror), I get the following error message: install.packages("Matrix",lib="/home/kn52/R-2.5.0/library") Warning message: package 'Matrix' is not available in install.packages("Matrix",lib="/home/kn52/R-2.5.0/library") Other

Dear R-Users, Is it possible to get the covariance matrix from an lme model that did not converge ? I am doing a simulation which entails fitting linear mixed models, using a "for loop". Within each loop, i generate a new data set and analyze it using a mixed model. The loop stops When the "lme function" does not converge for a simulated dataset. I want to

this is from a fresh svn checkout: [stvjc at stvjc1 R-devel]$ ls -tl Makefile -rw-rw-r-- 1 stvjc stvjc 12833 Jun 26 07:24 Makefile [stvjc at stvjc1 R-devel]$ svn up At revision 45988. [stvjc at stvjc1 R-devel]$ bin/R R version 2.8.0 Under development (unstable) (--) Copyright (C) The R Foundation for Statistical Computing ISBN 3-900051-07-0 .. > sessionInfo() R version 2.8.0 Under

http://bugzilla.mindrot.org/show_bug.cgi?id=917 Summary: ssh_exchange_identification: Connection closed by remote host Product: Portable OpenSSH Version: 3.6.1p2 Platform: ix86 OS/Version: Linux Status: NEW Severity: normal Priority: P3 Component: sshd AssignedTo: openssh-bugs

Short form: How do I transform the output of grid.locator() (or grid.locator(unit='npc')) to the native (or npc) coordinates of a viewport other than the top-level viewport? Thanks in advance. -Ben Long form: I would like the user to be able to click anywhere on the R graphics window. I would then like to figure out which leaf viewport the user clicked in and compute the native