similar to: piecing together statements (macro?)

Léandre BASSOLE PhD Student CNRS-CERDI 65 Bd Francois Mitterrand Boite Postale 320 63009 Clermont-Ferrand CEDEX 1 FRANCE Tel : +33 4 73 17 74 45 Fax : +33 4 73 17 74 28 ----- Forwarded Message ---- From: Leandre Bassole <leandrebassole@yahoo.co.uk> To: r-help@r-project.org Sent: Saturday, 27 October, 2007 8:41:05 PM Subject: Hi all!! I am a new user of R. I am very familar to Stata, but

Can lme or lmer fit a plain regular fixed effects anova? Ie a model without a random effect, or have there be at least one random effect in order for these functions to work? Trying to run such, (1) without specifying a random effect produces an error, (2) specifying that there is no random effect does not produce the same output as an anova run in lm(); (2b) specifying that there is no

In working through material on p.272 of MASS (4th ed.), I came across the following model formula: pet1.lm <- lm(Y ~ No/EP - 1, Petrol) I was at a loss to understand the use of "/" until I looked in "An Introduction [!] to R," where I found the explanation. My request is that more complete material on model formulae be lifted from "Introduction to R" (or

Dear R users, I am using "subset" function to filter out specific conditions and would like to use the value of subsetted argument as a name of an object. Specifically, from the following statement: a <- subset(dat, dat$x == "A" & dat$xx == 1 & dat$xxx == "AB" & dat$y == "B" & dat$yy == 2) I would to assign the value of the 3rd

R 2.4.1 on Windows XP. Question: In traditional graphics, is it possible to find out the height of a line of text in units that can be used in arithmetic and then in calls to text()? Context: I have written a function that draws a plot and then, depending on whether some arguments are TRUE or FALSE, draws various lines of text in the plot. The text lines may be turned on or off individually

Hi All I dont know if anyone else has noticed the same thing, but with 2 subsequent runs of the same syntax, I am getting exactly the same results. I was expecting that results differ slighlty, say in the 4th or 5th decimal place. Is this a specialty with R2WinBUGS? Does it have something to do with the seed value? Isnt the seed value reset everytime I restart winbugs? Thanks Toby

Dear List, Is there a possibility to change the font-face for a part of the title of a plot? For example I have the following... plot(nirs, type="l", xlab="Wellenl?nge [nm]", col="darkslategray", main = "Spektrum Deschampsia caespitosa") ...and I would like to change the part of the title-string "Deschampsia caespitosa" to italics? Is

I hope to use the plotmath facility to print titles that mix math and values of R variables. The help for "plotmath" has an example, which after repeated reading, I find baffling. Likewise, I have read the help file for "substitute" (wqhich seems to be needed) without ever understanding what it does, other than being used in some magic incantations. I would like to do

I need help with replacing NaN with zero (the value '0') in my dataset. The reason is that I can't get it to graph because of the NaN in the dataset. I have tried: data[is.nan(data)] <- 0 that others have suggested in the help archives but this does nothing so I am not sure what I am doing wrong. Thanks in advance for your help. [[alternative HTML version deleted]]

I was reading the help but just did not get how to specify starting values for varIdent() of the lme() function, although I managed to do it for corSymm(). Do I specify the values just as they are printed out in an output, like c(1, 1.3473, 1.0195). Or do I need to take the residual and multiply it with these like c(0.2235, 0.2235*1.3473, 0.2235*1.0195) or any other form that I dont know of?

hello, today while trying to extract data from a list for subsequent analysis, i stumbled upon this funny behavior on my system: > x<-c(0.1,0.9) > 1-x[2] [1] 0.1 > x[1] [1] 0.1 > x[1]==1-x[2] [1] FALSE > x[1]>1-x[2] [1] TRUE > x<-c(0.3,0.7) > x[1] [1] 0.3 > x[2] [1] 0.7 > 1-x[2] [1] 0.3 > x[1]==1-x[2] [1] FALSE but: >

Hi, I am thinking about IWLS vs ML estimation. When I use glm() for a 2-parameter distribution (e.g., Weibull), I can otain the MLE of scale parameter given shape parameter through IWLS. Because this scale parameter usually converges to the MLE. In this point, I am wondering: i) can you say that the direct MLE, which is obtained by maximizing a likelihood function, is equalvant to the indirect

Hello R users -- If I have a dataframe such as the following, named "frame" with the columns intended to be named col1 through col6, > frame col1 col2 cmlo3 col4 col5 col6 [1,] 3 10 2 6 5 7 [2,] 6 8 4 10 7 1 [3,] 7 5 1 3 1 8 [4,] 10 6 5 4 9 2 and I want to correct or otherwise change the

Hi All I was wondering if I can avoid a time-consuming for loop on my 600000 obs dataset. school_id y 8 9.87 8 8.89 8 7.89 8 8.88 20 6.78 20 9.99 20 8.79 31 10.1 31 11 There are, say, 143 different schools in this 600000 obs dataset. I need to thave sequential identifiers, 1,2,3,4,5,...,143. I was using an

library(StreamMetabolism) day <- (structure(c(10.08, 10.08, 10.05, 10.03, 10, 9.98, 9.96, 9.95, 9.95, 9.96, 9.96, 9.98, 10.01, 10.05, 10.06, 10.09, 10.11, 10.11, 10.13, 10.13, 10.15, 10.15, 10.13, 10.14, 10.11, 10.13, 10.14, 10.13, 10.12, 10.13, 10.14, 10.16, 10.18, 10.19, 10.23, 10.27, 10.33, 10.37, 10.45, 10.57, 10.6, 10.66, 10.73, 10.77, 10.84, 10.86, 10.87, 10.94, 10.98, 11.01, 11.05,

Do I miss here something? dtaa = read.table("http://www.ats.ucla.edu/stat/mplus/examples/ma_snijders/mlbook1.dat", sep=",") head(dtaa) # shows the data as it should be save(dtaa,"dtaa",file="c:/dtaa") d = load("c:/dtaa") head(d) # all data is lost, it only shows [1] "dtaa" "dtaa" Thanks for your hint on this.

Hi All Maybe I dont understand it, but I would have expected that the design matrix has as many rows as there were observations available to fit the model. Below a small artificial dataset created, then one model fitted and the design matrix outputted, having 27 rows. Then I delete 6 obs, and fit the model on these 21 obs, but the design matrix that comes out has 26 rows? Thanks for your

How on Earth can I know what are the arguments of any of the functions of the tcl/tk package? I tried hard to find, using all search engines available, looking deep into keywords of R, python's tkinter and tcl/tk, but nowhere I found anything remotely similar to a help. For example, what are the possible arguments to tkgetOpenFile? I know that this works: library(tcltk) filename <-

Gurus, Can I rely on the rownames() function, when applied to a matrix, always returning either NULL or an object of type character? It seems that row names can be entered as integers, but as of now (R 2.4.1 on Windows) the function returns character vectors, not numbers (which is my desired result). I am using elements of the returned vector to index the matrix. e.g., nams <-

Hi All I was wondering how to get the same results with gls and lme. In my lme, the design matrix for the random effects is (should be) a identity matrix and therefore G should add up with R to produce the R matrix that gls would report (V=ZGZ'+R). Added complexity is that I have 3 levels, so I have R, G and say H (V=WHW'+ZGZ'+R). The lme is giving me the correct results, I am