search for: gilescran

...ing, then one could write raw= " x y z 1 2 3 4 5 6 " df = read.table(raw,head=TRUE) Of course, one can cat() to write raw into a temporary file, and read.table() from that file. However, direct reading might be a good option? Hope this is useful. Cordially Giles Crane gilescrane at verizon.net

Greetings, The aov() function may mislabel the random effects as in the example below: Has anybody else noticed this? Cordially, Giles Crane, MPH, ASA, NJPHA gilescrane at verizon.net > m2 Call: aov(formula = y ~ ap + pe + Error(ju), data = d) Grand Mean: 77.50667 Stratum 1: ju Terms: ap Sum of Squares 4322.538 Deg. of Freedom 12 13 out of 25 effects not estimable Estimated effects may be unbalanced Stratum 2: Within Term...

...Mandel obtrain the following adjusted means mandel = c(1.728, 1.817, 1,588, 1.785, 1.512, 1. 678, 1.720, 1.658) Can anyone shed light on how to obtain in R the correct treatment means, adjusted for both row and column effects (blocks) in these designs? Cordially, Giles Crane gilescrane at verizon.net

Full_Name: Giles L Crane Version: 2.8.0 2008-10-20 OS: Windows XP Submission from: (NULL) (71.168.215.186) paste(rep(" ",10),sep="") [1] " " " " " " " " " " " " " " " " " " " " As I read the documentation for paste(), this should have produced " " that is,

Colleagues: I am having trouble loading data from within .Rdata file within the environment of a function. That is, the following always loads to the global environment: f1 <- function() { load("mydata.Rdata") # compute with objects from mydata.Rdata. } I wish to load mydata.Rdata only within the function f1. Perhaps I have misunderstood

Greetings. How can I insert gridlines in the following density plot call? Must one compose a panel function? 'data.frame': 46 obs. of 2 variables: $ fallrates: num 5.2 7.1 7.1 9.8 3.7 7.5 5 6.2 1.5 2.9 ... $ prepost : Factor w/ 2 levels "post","pre": 2 2 2 2 2 2 2 2 2 2 ... library(lattice) densityplot(~fallrates, groups = prepost,

Surprising behavior: Most R users are aware that print() must be used inside functions to gain output on the console. Apparently, print() is sometimes required when interactively using the console. For example, the followingmay be entered into the R console with different results. sample(1:8,8) #prints a permutation of 1 to 8 for(i in 1:5) #does not sample(1:8,8)

There are instances when assignment should be indicated by "<-" rather than "=". One example is the function Sys.time(). sytem.time( x=1:10000) # x is interpreted as an argument and is an error here. system.time( x<-1:1000) # times the creation of x with numbers 1 to 10000 However, grouping the statement with curly brackets also gives the correct system time.