similar to: (no subject)

Thanks to John Fox, Brian Ripley, and Peter Dalgaard for responding. The short answer (as in Peter Dalgaard's reply, already posted to the list) is that the models I'm concerned with can in fact be compared using ancova. The key fact is that while the parameters may not be nested, the subspaces I'm examining are. An additional note from Prof. Ripley on AIC and BIC (which I quote in

Hello, all. I've got what is probably a simple question about comparison of models using anova, specifically about the situations in which it's valid. I understand, I think, what's going on when the models are strictly nested (as most are in the demo(lm) examples). My question involves what happens when the models aren't strictly nested. In my particular case, I'm doing

Dear Dr. Sheth, This spring, Dr. Cottler and the EPRG submitted a renewal for the T32 NIMH Training Program in Psychiatric Epidemiology and Biostatistics. We are submitting a letter to Council for the competing continuation for NIMH which better showcases the accomplishments of our past trainees. It is vital that you provide us with an update that includes grants received since the completion

The IP addresses are assigned via DHCP - a Linksys router. David David I. Sommers, Ph.D. Scientific Review Administrator Interventions Review Committee (ITV) National Institute of Mental Health 6001 Executive Boulevard Room 6144, MSC 9608 Bethesda, MD 20892 Rockville, MD 20852 (for express or courier service) 301-443-7861 - voice 301-443-4720 - fax dsommers@mail.nih.gov - email -----Original

Yes you can name the linux boxes anything you want, no you don't need a domain... probably want to add the names and their respective IP addresses to your hosts files. > -----Original Message----- > From: Sommers, David (NIH/NIMH) [mailto:dsommers@mail.nih.gov] > Sent: Thursday, December 26, 2002 3:34 PM > To: 'samba@lists.samba.org' > Subject: [Samba] (no subject) >

hello, I'm trying to use R to take an image matrix and paint it into a tcltk canvas. I'm using R-1.9.0 within ess-5.2.0 on a linux machine running the 2.4.30-31.9smp kernel. I'm using the ActiveTcl8.4.6.1-linux-ix86 tcltk libraries. When I run the following comands, however, I get an R segmentation fault: ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~

How do I add to a trellis plot the best fit line from a robust fit? I can use panel.lm to add a least squares fit, but there is no panel.rlm function. -- Alan S Barnett <asb at mail.nih.gov> NIMH/CBDB

I'm using xyplot to generate a trellis plot with each panel containing a scatterplot and a best fit line. Is it possible to write the slope of the best fit line in each panel? -- Alan S Barnett <asb at mail.nih.gov> NIMH/CBDB

1) I have a data that includes some "bad" data. I want to make a trellis plot where each panel contains a) A scatter plot of the "good" data b) A scatter plot of the "bad" data in a different color c) A best fit line of all the data, and d) A best fit line of the "good" data. I tried using xyplot and setting the "group" argument, but I'm

Matthew Wiener <mcw@ln.nimh.nih.gov> writes: > t1 <- data.frame(matrix(rnorm(16), nc=4)) > t1$X1 <- 1 > t1$X2 <- 2 > print(t1) > Error: dim<- length of dims do not match the length of object Well, it is prototype-compatible. Splus 5.3 does likewise. A way out is t1<-data.frame(unclass(t1)) However, we do seem to have a bug in the area: > t1 <-

Matthew Wiener <mcw@ln.nimh.nih.gov> writes: > t1 <- data.frame(matrix(rnorm(16), nc=4)) > t1$X1 <- 1 > t1$X2 <- 2 > print(t1) > Error: dim<- length of dims do not match the length of object Well, it is prototype-compatible. Splus 5.3 does likewise. A way out is t1<-data.frame(unclass(t1)) However, we do seem to have a bug in the area: > t1 <-

Matthew Wiener <mcw@ln.nimh.nih.gov> writes: > > parameters (xxdefun in gram.y calls lang4). The arg count got changed > > from 0.64.2 to 0.65.0 (from 3 to 4) to accommodate the stored source > > attribute. I suppose that a function stored in a saved workspace might > > cause some kind of mess when restored. > > > > If I understand this correctly, it

================================================= This is an automated summary of the status of the R-bugs repository. Note that this may be neither complete nor perfectly correct at any given instance: Not all bugs are reported, and some reported bugs may have been fixed, but the repository not yet updated. Some bug fixes are difficult to verify because they pertain to specific hardware or

I would like to add an extension to the current name of a variable to create a new variable that is its sqrt transform. Each piece of the equation below works independently, but the left side definition fails on run. I also tried creating the variable name first, but ended up with an object that toString() did not fix. Better ideas? Example 1

Hi everyone, I'm a new R user (if this is a really basic question, please do excuse me...) and I'm having some questions regarding a deciles problem. I have a variable which I need to categorize according to its deciles (X). However, this categorization should be made into another variable (call it NewVar). Ex. for the quartiles case (just for the sake of exposition, since I need

I need to look through a dataset with two factor variables, and depending on certain criteria, create a new variable containing the data from one of those other variables. The problem is, R keeps making my new variable an integer and saving the data as a 1 or 2 (I believe the levels of the factor). I've tried using as.factor in the IF output statement, but that doesn't seem to work. Any

Hello list, how can I go from CALLINGout to just CALLING ? I've tried : exten => s,n,Set(newVAR=${CUT(CALLINGout,,3)}) or exten => s,n,Set(newVAR=$[CUT(CALLINGout,,3)]) But no result : [Jan 4 11:10:12] -- Executing [s at from-S:34] NoOp("SIP/s2-0000003b", "newVAR=") in new stack Asterisk 1.6.10 here. Kind regards, Jonas. -------------- next part

I have several newbie questions - any help much appreciated: 1. I run a combination Windows 98 and XP home network - have been looking at Linux on two boxes at home - one desktop - one laptop - RH 7.3 2 I can get samba to run okay and in each of the two Win98SE machines can see them as localhost 3. I cannot access any resources on the Linux machines - Win98 asks for a password - \\localhost\IPC$

Hi, there I'm currently trying to figure out how to keep my "factor" levels for a variable when moving it from one data frame or matrix to another. Example below: vec1<-(rep("10",5)) vec2<-(rep("30",5)) vec3<-(rep("80",5)) vecs<-c(vec1, vec2, vec3) resp<-rnorm(2,15) dat<-as.data.frame(cbind(resp, vecs))

Hi, all. I'm trying to use some linear regression models in which both the dependent and independent variables are measured with some error. To make things worse, while the errors in the dependent variable are uniform, the errors in the independent (or explanatory, or "x") variables can be heteroskedastic. I've been looking at the book _Measurement Error Models_ by Fuller