similar to: error constraints in lme

I am learning about sparse matrices and wonder if R can create them from a full matrix. Can anyone tell me how I might be able to accomplish this. ------ Harold C. Doran Director of Research and Evaluation New American Schools 675 N. Washington Street, Suite 220 Alexandria, Virginia 22314 703.647.1628 <http://www.edperform.net/> [[alternate HTML version deleted]]

Dear R: Am I having trouble downloading the LME4 library. I am using Windows and am using ver 1.7 I have tried the following: 1) Install package from CRAN, but LME4 is not listed 2) Downloaded LME4 from http://cran.us.r-project.org/, however, I cannot open the file when I try install from local drive. I get the following error: Error in file(file, "r") : unable to open connection

Dear List: I have a data set with over 7000 students with about 4 observations over time per student. I want to examine the within-group fits of a random sample of this group as it takes forever to compute and draw all 7000 regressions. Here is the code I have used so far. >group<-groupedData(math~year|childid, data=scores) >group.list<-lmList(group)

Dear list In HLM, one can specify a covariate at one of the "levels". For example, if the data structure are repeated observations nested within students nested within schools, school size might be a covariate that is used at level 3, but not at the other levels. In HLM this is rather easy to do. However, how can one specify a covariate in R for only one of the levels? I have a

Dear List I have math test scores for male and female students where gender is a dummy code (female =1). I also have a variety of other demographic variables. However to begin, I want to create a very simple stripchart where female math scores are a blue circle and male scores are a red triangle. I am having difficulty using conditional statements to accomplish this. Thank you. ------

Dear List: I am trying to find the area between two ECDFs. I am examining the gap in performance between two groups, males and females on a student achievement test in math, which is a continuous metric. I start by creating a subset of the dataframe male<-subset(datafile, female="Male") female<-subset(datafile, female="Female") I then plot the two CDFs via

Dear R: I am trying to fit a doubly multivariate LME (DM) where I have two response variables measured on two occasions per person. Specifically, reading and math scores measured at the beginning and ending of a school year. The response variables have a correlation of r = .85. The response variables in the data matrix are stacked in a vector with a dummy code flagging each outcome and with

I am trying to import a file using R Commander. It was working a few days ago, but now I get the following message when I try to import from SPSS. Any thoughts? Error in parse(file, n, text, prompt) : parse error ------ Harold C. Doran Director of Research and Evaluation New American Schools 675 N. Washington Street, Suite 220 Alexandria, Virginia 22314 703.647.1628

I am new to R and would like to get some practice analyzing multilevel data. I wonder if anyone can point me to a sample data set and command lines that I might replicate for a sample session. I would then compare my output with HLM output. Any help is appreciated. ------ Harold C. Doran Director of Research and Evaluation New American Schools 675 N. Washington Street, Suite 220 Alexandria,

I have loaded the foreign package and am still having problems with an import. I get a message that reads, unable to open file. Whe I try different files I get the same message. Here is the code I used. Am I missing something? I am using 1.7 and have also tried this in 1.6 with the same problem. hsb<-read.spss("C:\HLM504_Student\Examples\AppendxA\HSB1.SAV", use.value.labels=TRUE,

Dear List: I am having difficulties with the fitted values at different levels of a multilevel model. My data set is a series of student test scores over time with a total of 7,280 observations, 1,720 students nested witin 60 schools. The data set is not balanced. The model was fit using eg.model.1<-lme(math~year, random=~year|schoolid/childid, data=single). When I call the random

Have a look at ?optim. I don't think it has the BHHH algorithm as an option, though. =========================================== David Barron Jesus College University of Oxford -----Original Message----- From: r-help-bounces at stat.math.ethz.ch [mailto:r-help-bounces at stat.math.ethz.ch]On Behalf Of Harold Doran Sent: 10 July 2003 15:43 To: Fohr, Marc [AM]; R-help at stat.math.ethz.ch

Dear listers, I can't find the variance or se of the coefficients in a multilevel model using lme. I want to calculate a Chi square test statistics for the variability of the coefficients across levels. I have a simple 2-level problem, where I want to check weather a certain covariate varies across level 2 units. Pinheiro Bates suggest just looking at the intervals or doing a rather

If these things bother you, could you please supply fixes? For some reason you omitted the part of the help page dealing with one of these (which is as documented), and I also think you are confusing `plot' and `device' sizes. The latter will be the displayed device surface, I believe. Nothing on the help page suggested that the fine details of the R read-only parameters will be as you

Dear List: Earlier this week I posted a question and received no response, and I continue to struggle with my model. My original question is pasted below. I am using lme and want to fix the variance of the within group residual at 1 (e~n(0,1). I think the varFixed function should be used to accomplish this, but I am struggling to figure out how to do this. Can anyone offer suggestions on how

Part 1) is basically that the initial rescaling is done if the window is too large to fit on the screen. That's true and obvious (no scrollbars) and intentional and on the help page. As ever, you need to report bugs on what is documented, not what you would like. This will be changed for 1.2.1, but was several hours' work. I have also cleaned up the scrollbar behaviour. Part 2) I just

Hi, I have a consultants nightmare -- I was given a project that another consultant did and I was told to do the same calculations, but there's no documentation on what he did. Basically, I have yes/no answers to survey questions about the effectiveness of product attributes by brands. There are 44 attributes and 13 brands. The other guy scaled the proportion of respondents who said

Wayne, Hello, I saw a post of yours out on samba.org and was wondering if you might be able to shed a little light on a problem I'm having with rsync. My rsync source files are all off of nfs mounts and the destination directories are all nfs mounts. The owner and groups all end up as nobody/nobody on the destination, despite being all sorts of owners and groups on the source. Have you

Dear useRs, I have a data frame and I want to plot all rows. Each row is represented as a line that links the values in each column. The plot looks like this: dfr <- data.frame(A=sample(1:50,10),B=sample(1:50,10), C=sample(1:50,10),D=sample(1:50,10)) xa <- 10*1:4 plot(c(10,40),c(0,50)) for (i in 1:nrow(dfr)) { lines(xa,dfr[i,],pch=20,type="o") } Things get more complicated

Good Morning to All, I have a Red Hat 7.3 box with Samba 2.2.5 successfully integrated with a Windows 2000 domain. I am using a Win2K Domain Controller to hold all of the user accounts. I use WebMin 1.0 for the administration of my Linux box, and am able to pull up domain user accounts when putting permissions on files on the RH machine. I used a combination of resources to do this.