similar to: ' R ' - General Question (newbie)

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

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,

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

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

Could someone help me on the following: SAS has DUD (Does not Use Derivatives) for nonlinear regression. Does "R" has a similar capability? I am not good at derivatives and may get my derivative wrong before feeding it to a nonlinear regression procedure. Any help would be much appreciated. Liu Hancock Forest Management NZ Tokoroa, New Zealand DDI: 07-8850387 Mobile: 021-1576178

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

Hi All, I would like to output the results of a function into a text file, legible as a such. The function produces a summary quite like: summary(lm(x ~ y + w * z)) [for instance] and I am not clear how to save this summary into a text file 'automagically', because I need to be able to do it in a for() loop. Cheers, Federico -- Federico C. F. Calboli Department of Epidemiology

Dear R-users, I would like to write Greek letters followed by subscripts in a graph (on the X-axis and in a legend). I would appreciate any help. thank you, delphine courvoisier ************************************ Delphine Courvoisier Clinical Epidemiology Division University of Geneva Hospital +4122 37 29029

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 All I am pretty new to R but saw stata and sas's macro facilities and am looking for how such things work in R. I am trying to piece together a series of statements: n = 5 #want to have it dynamic with respect to n for (j in 1:n) { eval(paste("x", j, "=x[", j, "]", sep="")) } I want the created statements 'x1=x[1]' immediately executed

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

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

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

My institute has been heavily dependent on SAS for the past while, and SAS is starting to charge us a very deep amount for license renewal. Since we are a non-profit organization that is definitely not sustainable. The team is brainstorming possibility of switching to R, at least gradually. I am talking about the entire institute with considerable number of analysts using SAS their entire

Hello all, My statistical analysis training up until this point has been entirely done in SAS. The code I frequently used was: *Yield Champagin; data yield; set stress; if field='YV' then delete; if field='HB' then delete; if barcode='16187DD4015' then delete; if barcode='16187DD6002' then delete; if barcode='16187DD2007' then delete; if

Regarding point 3, as a moderator I have been helping Andrew get this post out to the list over the past week. His previous attempts were encoded in some way that the listserv rejected. He sent me the post via his gmail account and viewing the source I saw it had at least both plain test and HTML an I said it was worth a try to post it. Certainly on my mail client his post displays acceptably

something like ... > x<-c(1,2,3,0,5,6,0) > is.na(x[x==0])<-T > x [1] 1 2 3 NA 5 6 NA ~~~~~~~~~~~~~~~~ Robert Kinley ~~~~~~~~~~~~~~~~ [[alternative HTML version deleted]]

Hi, I came across a case where there's a discrepancy between minimum and maximum values reported by 'summary' and the 'min' and 'max' functions: ---<---------------cut here---------------start-------------->--- R> str(tt) num [1:1397] 1952 1970 1976 1967 1946 ... R> summary(tt) Min. 1st Qu. Median Mean 3rd Qu. Max. 1920 1960 1970