similar to: Problem with Trellis graphics in nlme

Hi, As part of the analysis of a longitudinal study, I have a list, called wg, on R 1.3.1 for windows created as follows - wg <- by(xx, xx$mrn, function(xx) lm(xx$maternal.wt ~ xx$gestation, data=xx)) xx is a data.frame, and so wg contains the results of 293 linear regressions, as lm objects. Now what I want is to extract the 293 intercepts and regression coefficients as a vector, and

Dear alltogether, how can a delay like possible with par(ask=TRUE) be attained while using trellis-plots within a loop or something like that? the following draws each plot without waiting for a signal (mouse-klick), so par() does not work for that: library(nlme) for(i in 1:3) { fitlme <- lme(Orthodont) par(ask=TRUE) # does not work with trellis.... print(

Dear colleagues, I have a problem which has bitten me occasionally. I often need to prepare graphs for many variables in a data set, but seldom for all. or for any large number of sequential or sequentially named variables. Often I need several graphs for different subsets of the dataset for a given variable. I run into similar problems with other needs besides graphing. What I would like to

Hello, I am trying to graph a regression line using the followings: Age <- c(39, 47, 45, 47, 65, 46, 67, 42, 67, 56, 64, 56, 59, 34, 42, 48, 45, 17, 20, 19, 36, 50, 39, 21, 44, 53, 63, 29, 25, 69) BloodPressure <- c(144, 220, 138, 145, 162, 142, 170, 124, 158, 154, 162, 150, 140, 110, 128, 130, 135, 114, 116, 124, 136, 142, 120, 120, 160, 158, 144, 130, 125, 175) SimpleLinearReg1=lm(Age ~

Hi folks, Pls help me to understand follow; An Introduction to R 2.4 Logical vectors http://cran.r-project.org/doc/manuals/R-intro.html#R-and-statistics 1) > x [1] 1 2 3 4 5 > temp <- x != 1 > temp [1] FALSE TRUE TRUE TRUE TRUE > 2) > x [1] 1 2 3 4 5 > temp <- x > 1 > temp [1] FALSE TRUE TRUE TRUE TRUE Why NOT > temp [1] TRUE FALSE FALSE FALSE

Hi folks, Pls help me to understand follow; An Introduction to R 2.4 Logical vectors http://cran.r-project.org/doc/manuals/R-intro.html#R-and-statistics 1) > x [1] 1 2 3 4 5 > temp <- x != 1 > temp [1] FALSE TRUE TRUE TRUE TRUE > 2) > x [1] 1 2 3 4 5 > temp <- x > 1 > temp [1] FALSE TRUE TRUE TRUE TRUE Why NOT > temp [1] TRUE FALSE FALSE FALSE

Hi, I think I'm missing something very obvious, but I am missing it, so I would be very grateful for help. I'm using rpart to analyse data on skull base morphology, essentially predicting sex from one or several skull base measurements. The sex of the people whose skulls are being studied is known, and lives as a factor (M,F) in the data. I want to get back predictions of gender, and

Hi, Anne, assign Age and Bloodpressure in the correct order to the axes in your call to plot as in: plot(y = Age, x = BloodPressure) abline(SimpleLinearReg1) Hth -- Gerrit --------------------------------------------------------------------- Dr. Gerrit Eichner Mathematical Institute, Room 212 gerrit.eichner at math.uni-giessen.de Justus-Liebig-University Giessen Tel:

Dear Gerrit, This issue is discussed in a vignette in the car package (both for functions in the car and effects packages): vignette("embedding", package="car") . The solution suggested there is the essentially the one that you used. I hope this helps, John ----------------------------- John Fox, Professor Emeritus McMaster University Hamilton, Ontario, Canada Web:

Dear All, I am analysing a dataset on levels of herbivory in seedlings in an experimental setup in a rainforest. I have seven classes/categories of seedling damage/herbivory that I want to analyse, modelling each separately. There are twenty maternal trees, with eight groups of seedlings around each. Each tree has a TreeID, which I use as the random effect (blocking factor). There are two

Hi, again, I have to modify my query since my first (too simple) example doesn't reflect my actual problem. Second try: When asking Effect() inside a function to compute an effect of an lmer-fit which uses a data frame local to the body of the function, as in the following example (simplifying my actual application), I get the "Error in is.data.frame(data) : object 'X' not

Hello, everyody, when asking, e.g., Effect() to compute the effects of a fitted, e.g., linear model after having deleted the data frame from the workspace for which the model was obtained an error is reported: > myair <- airquality > fm <- lm(Ozone ~ Temp, data = myair) > rm(myair) > Effect("Temp", fm) Error in eval(model$call$data, envir) : object 'myair'

Dear all, I work with R Version 2.0.0 on Machine hardware: sun4u OS version: 5.9 Processor type: sparc Hardware: SUNW,Sun-Blade-1000 and I have a very simple data frame (called OR) with the following variables: > sapply( OR, class) X ci FTyp "factor" "numeric" "factor" (In OR$ci there are some Inf-values. OR's

Dear R-help, I am trying to express myself as best as I can here. If you also use Latex to edit math reports or other languages with similar editing method, you'll see what I'm talking about. My sincere appologies if my question is not clear enough to some extend, as also I'm not able to provide my code here because I don`t know which one I can use... When editing the title in R

Dear colleagues, This may be a question with a really obvious answer, but I can't find it. I have access to a large file with real medical record identifiers (mixed strings of characters and numbers) in it. These represent medical events for many thousands of people. It's important to be able to link events for the same people. It's much more important that the real record numbers

Hi everybody! I'm sure that I overlook something and feel quite stupid to ask, but I have not found an easy solution to the following problem: Take e.g. the Orthodont data from the nlme package: > head(Orthodont) Grouped Data: distance ~ age | Subject distance age Subject Sex 1 26.0 8 M01 Male 2 25.0 10 M01 Male 3 29.0 12 M01 Male 4 31.0 14 M01 Male

Consider the following plot that shows separate regression lines ~ age for each subject in the Pothoff-Roy Orthodont data, with separate panels by Sex: library(nlme) #plot(Orthodont) xyplot(distance ~ age|Sex, data=Orthodont, type='r', groups=Subject, col=gray(.50), main="Individual linear regressions ~ age") I'd like to also show in each panel the pooled OLS

>From owner-r-help at stat.math.ethz.ch Wed Jan 9 10:08:44 2002 > >Does anyone know how to access/list the value labels stored in a SPSS .sav >file? >At the moment I can read in variables using read.spss() (foreign package), >but have to power up SPSS to check value labels. > >Stuart There may be a solution to this in R, and I would love to hear it, but I was just faced

Hi, I'm analysing some anthropometric data on fifty odd skull bases. We know the gender of each skull, and we are trying to develop a predictor to identify the sex of unknown skulls. Rpart with cross-validation produces two models - one of which predicts gender for Males well, and Females poorly, and the other does the opposite (Females well, and Males poorly). In both cases the error

Dear colleagues, I would be very grateful for your help with the following. I have banged my head off this question several times in the past, and repeatedly over the last week. I have looked in the usual places and found no obvious solution. I fear that this just means I didn't recognize it, but I'd be very grateful for your help. I am scoring 8000 psychometric tests - the SCQ, if