similar to: A question about digest mode

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

Good afternoon Using R 2.9.2 on a machine running Windows XP I have a longitudinal data set, with data on schools and their test scores over a four year period. I have centered year, and run the following m1.mod1 <- lme(fixed = math_1 ~ I(year-2007.5)*TFC_, data = long, random = ~I(year-2007.5)|schoolnum, na.action = "na.omit") where

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 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

I wrote As a beginner, I agree .... the for loop is much clearer to me. Wacek Kusnierczyk <Waclaw.Marcin.Kusnierczyk at idi.ntnu.no> replied > >well, that's quite likely. especially given that typical courses in >programming, afaik, include for looping but not necessarily functional >stuff -- are you an r beginner, or a programming beginner? > Both. My PhD is in

Proposal That a new mailing list be established that pertains exclusively to R documentation. The purpose of the list would be to discuss weak sections of the documentation and establish fixes for those weak spots. Pro If it works, there would be better documentation. It would be an excellent opportunity for newish and/or less technical people to contribute to R. In some respects such people

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

I deal with a huge amount of Biology data stored in different databases. The databases belongig to Bioconductor organization can be accessed through Bioconductor packages. Unluckily some useful data is stored in databases like, for instance, miRDB, miRecords, etc ... which offer just an interactive HTML interface. See for instance http://mirdb.org/cgi-bin/search.cgi,

Christophe Genolini <cgenolin at u-paris10.fr> wrote >Thanks for yours answers. So if I understand: > - Trajectories are continuous, the other are discrete. > - The difference between time series and longitudinal is that time >series are made at regular time whereas longitudinal are not ? > - Repeated measures are over a short period of time. > > >So if I measure

Hello all I have a feeling this is very simple......but I am not sure how to do it My boss has two variables, one is an average of 4 numbers, the other is an average of 3 of those numbers i.e var1 = (X1 + X2 + X3 + X4)/4 var2 = (X1 + X2 + X3)/3 all of the X variables are supposed to be measuring similar constructs not surprisingly, these are highly correlated (r = .98), the question is how

A time series text with a title that seems designed to hide its wide scope is: Forecasting with Exponential Smoothing The State Space Approach Hyndman, R.J., Koehler, A.B., Ord, J.K., Snyder, R.D. Springer 2009. This book is actually an excellent overview of time series theory, ARIMA as well as state space. It is of course, in part, a manual for the forecast and other packages in what has been

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

Pursing my earlier question, when I tried loading Lindsey's gnlm, I got a message Loading required package: rmutil Warning message: There is no package called 'rmutil' in: library(package, character.only = TRUE, logical = TRUE, warn.conflicts = warn.conflicts, According to the R documentation http://finzi.psych.upenn.edu/R/doc/html/packages.html rmutil is in the standard

When I try citypro <- sas.get('c:/ndri/cvar/data', member = 'citypro2') I get the following errors and warnings: 'sas' is not recognized as an internal or external command, operable program or batch file. 'less' is not recognized as an internal or external command, operable program or batch file. Error in sas.get("c:/ndri/cvar/data", member =

Does anyone know of research on the effects of rounding on regression? e.g., when you ask people "How often have you _______?" you are more likely to get answers like 100, 200, etc. than 98, 203, etc. I'm interested in investigating this, but don't want to reinvent the wheel. thanks Peter Peter L. Flom, PhD Assistant Director, Statistics and Data Analysis Core Center for

Thanks for all the help on my earlier questions..... How do you plot a simple time series with unequal intervals? I have the following dateofpoll <- as.ts(c("6/1/02", "7/1/02", "10/1/02", "1/4/03", "1/25/03", "6/7/03", "7/16/03", "8/17/03", "9/4/03")) reelect <- c(51, 47, 49, 51, 49, 49, 46, 45, 40)

Hello Using R 2.2.1 on a Windows machine. Has anyone programmed the Welch test for equality of variances? I tried RSiteSearch, but this gave references to t test and oneway.test, which are not quite what I need.....I need the Welch test itself, for use in a meta-analysis (to determine if variances are equal). TIA Peter Peter L. Flom, PhD Assistant Director, Statistics and Data Analysis

Hello I am running R 1.8.1 on a Windows platform I am attempting to fit an ordinal logistic regression model, using the polr function, as described in Venables and Ripley. But when I try model4 <- polr(ypsxcat~committed + as.factor(sex) + as.factor(drugusey) + anycsw + as.factor(sex)*committed + as.factor(sex)*as.factor(drugusey)+as.factor(sex)*anycsw, data = duhray) I get a message

Running R 1.8.1 on a Windows machine In dotchart, I would like to shrink the labels on the tick marks (that is, the numbers) without shrinking anything else. I could not find this in either the Rhelp archives or in ?dotchart, which recmmends cex to avoid 'label overlap', but cex shrinks all the characters in the plot. Is there a way to do this? Thanks Peter Peter L. Flom, PhD