similar to: Question on mixed effect models with LME

Hello Running R2.9.2 on Windows XP I am puzzled by the performance of LME in situations where there are missing data. As I understand it, one of the strengths of this sort of model is how well it deals with missing data, yet lme requires nonmissing data. Thus, m1.mod1 <- lme(fixed = math_1 ~ I(year-2007.5)*TFC_, data = long, random =

Good afternoon Running R2.10.0 on Windows I have a data frame that includes (among much else) a factor (In_2006) and a continuous variable (math_3_4). I would like to find the 2 cases for In_2006 = 0 that are closest to each case where In_2006 = 1. My data looks like In_2006 math_3_4 0 55.1 1 51.6 1 18.1 1 26.6 1 14.1

Hi, We could use qqplot to see how two distributions are different from each other. To show better how they are different (departs from the straight line), how is it possible to plot the straight line that goes through them? I am looking for some thing like qqline for qqnorm. I thought of abline but how to determine the slope and intercept? Best wishes, Carol

Good morning I just got a new computer with Windows 7. R works fine, but the editor I am used to using "RWinEdt" does not. I did find one blog post on how to get RWinEdt to work in Windows 7, but I could not get those instructions to work either. Is there a patch for RWinEdt? If not, is there another good R editor that works under Windows 7? I tried RSiteSearch with various

Good morning I am learning about NLME and LME4, using Pinheiro and Bates and other materials from Douglas Bates, but I have not seen anything on how to do variable selection sensibly in this type of model. In OLS regression, I frequently use the lasso, but googling did not reveal a method for lasso with mixed models. Most of the material I've seen on these packages is about models with very

Hi everyone, I am wondering if there exists a stepwise regression function for the Bayesian regression model. I tried googling, but I couldn't find anything. I know "step" function exists for regular stepwise regression, but nothing for Bayes. Thanks -- View this message in context: http://www.nabble.com/Bayesian-regression-stepwise-function--tp26013725p26013725.html Sent from

Dear All, Is there some way of drawing a boxplot, with R, when one does not have the original continuous data, but only the data grouped in classes? The function boxplot() can only deal with original data. Thanks in advance, Paul

Hello, just a simple question How can I do to invert the sign of a number? (for example: -4 to 4, 2 to -2 and so on..) I was looking for a specific function in R but I didn't found it... thank you Francesco -- View this message in context: http://www.nabble.com/Invert-the-sign-of-a-number-tp26110267p26110267.html Sent from the R help mailing list archive at Nabble.com.

Hi everybody, well, I definitely don't understand anything. Why the hist function with freq=FALSE gives such a strange result??? R <- c(-1.10, 0.79, -1.17, -0.53, -0.26, -0.22, 0.29, -0.26, -0.26, 0.39) hist(R, freq=FALSE, breaks=10) Thanks everybody [[alternative HTML version deleted]]

Hi I can install RWinEdt if I start R with administrator rigths, but it does not paste my code to the console. I found advice in the link below how to manage the problem, but it did not work, any other idea? http://yusung.blogspot.com/2009/01/rwinedt-and-windows-vistawindow-7.html Thanks a lot,Johannes >From: Uwe Ligges <ligges_at_statistik.tu-dortmund.de> >Date: Sun, 08 Nov 2009

Dear R Community, I am attempting to transpose a dataset from rows to columns but am stuck. I have tried using reshape() with little luck, possibly due to the categorical nature of the data. For example: id<-c(1,2,2,3,3,3) author<-c("j","k","k","l","l","l")

I have a nested ANOVA, with a fixed factor "tmt" nested within "site" (random). There are missing values in the data set. aeucs, tmt and site have been defined as objects I have tried: model1=lme(aeucs~tmt,random=~1|tmt/site) I get the following error message Error in na.fail.default(list(aeucs = c(0.833333333, 1, 1, 1, 1, 1, 0.666666667, : missing values in object I

Hi All, I'm working with some SAS code to analyze an experiment set up as follows: 66 subjects (colonies) treated with a random treatment (1-8) and measured at three time points. The data structure looks like: input colony tmt y1 y2 y3; y=y1; date=*1*; output; y=y2; date=*2*; output; y=y3; date=*3*; output; datalines; 1

I have a local library which depends on the expm library. The expm library is loaded into my personal space and I have the environment variable R_LIBS_USER set appropriately. The command "library(expm)" works just fine from the command line, and in fact the package works if I do the source() and dyn.load() commands by hand. The following sequence works: tmt% R CMD build

Martin, I am aware of --vanilla; I use it myself for some testing. In this case R_LIBS_USER was set externally (part of my login) and does not involve any of the R scripts. That means it is inherited by any subprocess. For example: tmt1495% R --vanilla --no-environ R version 3.3.1 (2016-06-21) -- "Bug in Your Hair" Copyright (C) 2016 The R Foundation for Statistical Computing

A user asked me about this and I can't figure it out. tmt% R R Under development (unstable) (2018-04-09 r74565) -- "Unsuffered Consequences" Copyright (C) 2018 The R Foundation for Statistical Computing Platform: x86_64-pc-linux-gnu (64-bit) > library(survival) > data(cgd0) Warning message: In data(cgd0) : data set ?cgd0? not found ---- The data set is present and can be

I am getting failure of build and check, for an Rd file that has a long argument list. Guess diagnosis: a quoted string beyond a certain point in the argument list is fatal. Example: Use the function below, create an Rd file for it with prompt(). Move the .Rd file to the man directory (no need to edit it) and try building dart.control <- function(server=c("production",

I have a local library with functions that interrogates an institution-specific web API, so is not of interest to anyone outside of Mayo. For some reason the R CMD INSTALL command fails. See below: Build the library, then install it. tmt-local2127% R CMD build dart * checking for file ?dart/DESCRIPTION? ... OK * preparing ?dart?: * checking DESCRIPTION meta-information ... OK * installing

This is likely yet another instance of round off error, but it caught me by surprise. tmt% R --vanilla (headers skipped, version 3.0.2 on Linux) > load('qtest.rda') > length(temp) [1] 3622 > max(temp) >= quantile(temp, .98) 98% FALSE I can send the file to anyone who would like to understand this more deeply. The top 3% of the vector is a single repeated value. Terry

Hi, I would like to know how (if) I can extract some of the information from the summary of my nlme. at present, I get a summary looking something like this: > summary(fit.nlme) Nonlinear mixed-effects model fit by maximum likelihood Model: MLKYLD ~ W4(DIM, logA, B, C) Data: ADHIS.x0 AIC BIC logLik 265314 265401.6 -132647 Random effects: Formula: list(logA ~ 1 , B ~