Re: [R] advise for modeling a linear mixed model

2003-08-23 Thread Jason Turner
Frank Mattes wrote: Dear R help-list reader, I'm trying to investigate my data with linear mixed model and are seeking advise how to write the model in R. I was trying to get hold of the recommended book from Bates et al, but neither the major bookshop nor our university library had the book.

[R] Thank you

2003-08-23 Thread governor
Thank you for your email. I appreciate you taking the time to keep me informed of your opinions and the issues that are important to you. Your comments help keep me informed as we strive to make California a better place in which to work and live. Your email has been directed to members of my

Re: [R] advise for modeling a linear mixed model

2003-08-23 Thread Douglas Bates
Jason Turner [EMAIL PROTECTED] writes: As for the book, Bates Pinheiro, I'm working my way slowly through Actually, it's Pinheiro Bates. __ [EMAIL PROTECTED] mailing list https://www.stat.math.ethz.ch/mailman/listinfo/r-help

[R] help--kernel distribution dynamics

2003-08-23 Thread Eugene Salinas
Deall all, I'm just learning R, but unfortunately I need to urgently do a rather more complex task so I need some help. I have just learnt the very basics a few days ago and am not ready yet to deal with panels and kernel densities, so a soft guidance would be most appreciated. I have a (very)

[R] Majordomo results: Re: Thank you!

2003-08-23 Thread Majordomo-owner
See the attached file for details Illegal command! No valid commands processed. Ignoring part of type application/octet-stream. __ [EMAIL PROTECTED] mailing list https://www.stat.math.ethz.ch/mailman/listinfo/r-help

Re: [R] help--kernel distribution dynamics

2003-08-23 Thread Deepayan Sarkar
On Saturday 23 August 2003 10:36, Eugene Salinas wrote: Deall all, I'm just learning R, but unfortunately I need to urgently do a rather more complex task so I need some help. I have just learnt the very basics a few days ago and am not ready yet to deal with panels and kernel densities, so

Re: [R] Substitute in legend

2003-08-23 Thread Uwe Ligges
Peter Dalgaard BSA wrote: [EMAIL PROTECTED] writes: I tried to use substitute in legend as follows: pval - 0.04 plot(0) legend(1,0.5,substitute(hat(theta) == p, list(p = pval))) For some reason the legend is repeated 3 times. Any suggestions or is this a bug? It's a bug. The code is

Re: [R] help--kernel distribution dynamics

2003-08-23 Thread Eugene Salinas
--- Deepayan Sarkar [EMAIL PROTECTED] wrote: [...] As a first step, you could create a matrix (with 50 rows, one for each time point) where each row holds the kernel density estimate for that time point. e.g. (with a grid of size 100 for each estimated density), foo - matrix(0, 50,

Re: [R] help--kernel distribution dynamics

2003-08-23 Thread Deepayan Sarkar
On Saturday 23 August 2003 13:02, Eugene Salinas wrote: [...] Hi, Thanks a lot. This seems like what I want to do. I don't know all the syntax yet so just a clarification... Is the []n = 100)$y there in order to condition on y which is the year and derive the conditional kernel

Re: [R] help--kernel distribution dynamics

2003-08-23 Thread Prof Brian Ripley
As a possible enhancement, I would think about using the same bandwidth at all the time points --- indeed I would probably start by looking at a few time points, playing with the bandwidth and then using e.g. persp on density estimates at all 50 time points with that bandwidth. As a first

RE: [R] filling a matrix who's entries are a function of the ind

2003-08-23 Thread Ted Harding
On 23-Aug-03 Douglas G. Scofield wrote: What's the best way in R to fill a matrix who's entries depend on some function of the indices? I'm currently doing: Q - matrix(0, k, k) for (A in 1:k) { for (B in 1:k) { Q[A,B] - my.function(A,B) } } but I wonder

[R] explanation of lm's coefficients

2003-08-23 Thread Justin Fay
I don't understand the coefficients returned from the lm function. I expected these to be the mean values for each factor in the model. Given this data and model: data-c(rnorm(10,mean=0,sd=1),rnorm(10,mean=1,sd=1),rnorm(10,mean=-.5,sd=1)) ftr-as.factor(rep(1:3,each=10)) fit-lm(data ~ ftr) the