[R] miktex 2.6 and R

Hi, I am new to Latex and in the process of installing miktex to work with R on windows XP professional. From earlier posts, I see that fptex was the preferred version, but it is no longer available. Now miktex is preferred. However, there are known issues with finding the style files in version

[R] miktex 2.6 and R

Hi, I am new to Latex and in the process of installing miktex to work with R on windows XP professional. From earlier posts, I see that fptex was the preferred version, but it is no longer available. Now miktex is preferred. However, there are known issues with finding the style files in versi

[R] factor analysis and pattern matrix

Hi, In a discussion of factor analysis in "Using Multivariate Statistics" by Tabachnick and Fidell, two matrices are singled out as important for interpreting an exploratory factor analysis (EFA) with an oblique promax rotation. One is the "structure matrix". The structure matrix contains the cor

[R] factor analysis and pattern matrix

Hi, In a discussion of factor analysis in "Using Multivariate Statistics" by Tabachnick and Fidell, two matrices are singled out as important for interpreting an exploratory factor analysis (EFA) with an oblique promax rotation. One is the "structure matrix". The structure matrix contains the

[R] strategy for doing an ANOVA on unbalanced data

Hi, I would like some help deciding if and how to average my data before running an ANOVA. Let me first describe the data and what makes it unique. Hopefully, that will generate some ideas because I am not sure what kind of model I need to use. I don't know how to describe this succinctly, so ple

[R] drawing axes on all 4 sides of a trellis plot

Hi, I am creating a 4 x 8 trellis plot. I want the y axis scaled by the data in the rows, and the x axis scaled by data in each column. In March 2006 there was a post describing a way to do this. The link is: tolstoy.newcastle.edu.au/R/help/06/03/22937.html I am using the data from the example

Re: [R] greco-latin square

--- Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1 -----Original Message- From: Richard M. Heiberger [mailto:[EMAIL PROTECTED] Sent: Friday, May 19, 2006 12:16 AM To: Steven Lacey; r-help@stat.math.ethz.ch Subject: Re: [R] grec

[R] greco-latin square

Hi, I am analyzing a repeated-measures Greco-Latin Square with the aov command. I am using aov to calculate the MSs and then picking by hand the appropriate neumerator and denominator terms for the F tests. The data are the following: res

[R] register s3 object in package

Hi, I have a package with s4 methods that match against registered s3 methods (using setOldClass). When I call R CMD INSTALL --build I get warnings that methods are being created for objects with no definition. How do I embed setOldClass(c("pam","partition")) into the package so that when the s

Re: [R] aov contrasts residual error calculation

with respect to aov would be greatly appreciated. Steve -Original Message- From: Jacques Veslot [mailto:[EMAIL PROTECTED] Sent: Friday, April 21, 2006 11:58 AM To: Steven Lacey Cc: r-help@stat.math.ethz.ch Subject: Re: [R] aov contrasts residual error calculation why not using lme() ?

[R] aov contrasts residual error calculation

Hi, I am using aov with an Error component to model some repeated measures data. By repeated measures I mean the data look something like this... subjABC 1 411 15 2 312 17 3 5914 4 610 18 Fo

[R] message posting

Hi, I sent two emails to the R help list and got no reply. While this may be my question, it is unlike the users of this list not to reply (thankfully!). I checked for my message using the archive for April 2006 and found the following where the text of message should have been: An embedded an

[R] more on aov contrasts residual error calculation

Hi, I haven't recieved any replies to my last email, so let me be a bit more specific. I have a dataframe and it has the following structure: Condition Mapping SubjectABC 11 510 15 12

[R] aov contrasts residual error calculation

Hi, I am using aov with an Error component to model some repeated measures data. By repeated measures I mean the data look something like this... subjABC 1 411 15 2 312 17 3 5914 4 610 18 Fo

Re: [R] S4 class slot name 'names' is not allowed

Martin, Thanks! I'll use name for now. Steve -Original Message- From: Martin Maechler [mailto:[EMAIL PROTECTED] Sent: Thursday, April 13, 2006 3:11 AM To: Steven Lacey Cc: [EMAIL PROTECTED]; r-help@stat.math.ethz.ch Subject: S4 class slot name 'names' is not allowed

[R] S4 class slot names

Hi, Why doesn't this work? setClass("test",representation(names="character")) tmp <- new("test") [EMAIL PROTECTED] <- "a" Error in "slot<-"(object, name, TRUE, value) : 'names' attribute [1] must be the same length as the vector [0] Why does this work (replace names with name)? s

Re: [R] documenting s4 methods in package

Why do some functions build chm help and others do not? Is there something in the .Rd file that specifies this? Thanks, Steve -Original Message- From: Duncan Murdoch [mailto:[EMAIL PROTECTED] Sent: Tuesday, April 04, 2006 2:22 PM To: Steven Lacey Cc: 'Thomas Lumley

Re: [R] skipping rows in trellis key

gs$draw <- FALSE key <- do.call("draw.key",drawKeyArgs) vp1<-viewport(x=1,y=1,height=unit(1,"grobheight",key),width=unit(1,"grobwidt h",key),just=c("right","top")) pushViewport(vp1) grid.draw(key) popViewport() What is a good work aro

Re: [R] key position in trellis plotting area

) sets the current viewport to the one for the panel in which I want to position the key. Steve -Original Message----- From: Steven Lacey [mailto:[EMAIL PROTECTED] Sent: Thursday, April 06, 2006 8:55 AM To: 'r-help@stat.math.ethz.ch' Subject: key position in trellis plotting area

[R] key position in trellis plotting area

Hi, I want to do the following: 1) create a trellis plot with 1 x 1 layout 2) add a key in the upper right hand corner of the plotting region (i.e., the panel), but after the initial call to trellis 3) resize the resulting device graphics device without changing the relative position of the ke

Re: [R] skipping rows in trellis key

op")) pushViewport(vp1) grid.draw(key) popViewport() Steve -Original Message- From: Deepayan Sarkar [mailto:[EMAIL PROTECTED] Sent: Wednesday, April 05, 2006 11:53 PM To: Steven Lacey Cc: r-help@stat.math.ethz.ch Subject: Re: [R] skipping rows in trellis key On 4/5/06, Steven Lacey

Re: [R] skipping rows in trellis key

ey and the border? Thanks, Steve -Original Message- From: Deepayan Sarkar [mailto:[EMAIL PROTECTED] Sent: Wednesday, April 05, 2006 10:35 PM To: Steven Lacey Cc: r-help@stat.math.ethz.ch Subject: Re: [R] skipping rows in trellis key On 4/5/06, Steven Lacey <[EMAIL PROTECTED]> wr

[R] skipping rows in trellis key

Hi, I would like to add a key to my trellis plot using draw.key. Here is what I want: 3 x 6 key where the first row is a header. row 1: empty, "S-R Mapping", "R^2" row 2: pch=17, "Color", 0.951 row 2: pch=17, "Shape", 0.934 ect... The problem is that I would like the cell in the upper left

Re: [R] documenting s4 methods in package

es one get rid of it? Thanks, Steve -Original Message- From: Steven Lacey [mailto:[EMAIL PROTECTED] Sent: Tuesday, April 04, 2006 1:54 PM To: 'Thomas Lumley' Cc: 'r-help@stat.math.ethz.ch' Subject: RE: [R] documenting s4 methods in package Thomas, If I delete the \seeal

Re: [R] documenting s4 methods in package

Thomas, If I delete the \seealso section, I still get the same strange behavior. Shouldn't removing that line fix the problem? Thanks, Steve -Original Message- From: Thomas Lumley [mailto:[EMAIL PROTECTED] Sent: Tuesday, April 04, 2006 10:13 AM To: Steven Lacey Cc: r

[R] documenting s4 methods in package

Hi, I have written a package that contains many s4 generic functions and associated methods. I am having a lot of trouble getting R to build the help pages for these generic functions without reporting, "missing link(s): ~~fun~~, which means that it cannot find the appropriate function when code

Re: [R] object size vs. file size

MAIL PROTECTED] Sent: Tuesday, March 28, 2006 3:02 PM To: Steven Lacey Cc: 'Gabor Grothendieck'; r-help@stat.math.ethz.ch Subject: Re: [R] object size vs. file size On 3/28/2006 2:54 PM, Steven Lacey wrote: > Duncan, > > I wrote an R package to process my data. The package was w

Re: [R] object size vs. file size

erstanding of environments and functions. Would you recommend a good reference? Thanks, Steve -Original Message- From: Duncan Murdoch [mailto:[EMAIL PROTECTED] Sent: Saturday, March 25, 2006 9:12 PM To: Steven Lacey Cc: r-help@stat.math.ethz.ch Subject: Re: [R] object size vs. file size

Re: [R] object size vs. file size

[EMAIL PROTECTED] Sent: Saturday, March 25, 2006 5:31 PM To: Steven Lacey Cc: r-help@stat.math.ethz.ch Subject: Re: [R] object size vs. file size On 3/25/2006 7:32 AM, Steven Lacey wrote: > Hi, > > There is rather large discrepancy in the size of the object as it > lives in R and

[R] object size vs. file size

Hi, There is rather large discrepancy in the size of the object as it lives in R and the size of the object when it is written to the disk. The object in question is an S4 of a homemade class "sa". I first call a function that writes a list of these objects to a file called "data.RData". The siz

[R] setting argument defaults in setMethod

Hi, I want to set a default value in a method of a generic function. This seems as though it should be possible. From R help on setMethod... Method definitions can have default expressions for arguments. If those arguments are then missing in the call to the generic function, t

[R] setMethod bug

Hi, I define a generic function that has many parameters (~20). I then want to define methods for the generic function with setMethod. The default behavior for setGeneric is to allow dispatching on any argument in the def function, expect ... For example, # define a generic function setGen

[R] creating objects with a slot of class formula, using new

Hi, This works fine. >setClass("a",representation(b="list",c="list")) >new("a",b=list(7)) >An object of class "a" Slot "b": [[1]] [1] 7 Slot "c": list() But, now suppose you want a slot to accept an object of class formula... >setClass("a",representation(b="list",c="formula")) >new("a",b=li

[R] help with strip.default

Hi, I am creating a multi-conditioned trellis plot. My data look something like this: Factor AFactor BIVDV X 1 X 2 X 3 X 4 Y 1 Y 2 Y 3 Y 5 Z 1 Z 2 Z 3 Z

[R] organizing plot drawing routines; creating complex expressions

Hi, My general goal is to find a coding strategy to efficiently store and retrieve drawing routines for different plots. This is my approach so far. In a single text file I store multiple drawing routines where each routine draws a different plot. A drawing routine may look like this:

[R] slow computation of mixed ANOVA using aov

Dear R-help list, I am trying to do a mixed ANOVA on a 8960 x 5 dataframe. I have 3 factors for which I want to test all main effects and interactions : f1 (40 levels), f2 (7 levels), and f3 (4 levels). I also have a subject factor, subject, and a dependent variable, dv. Some more information

[R] location of key in panels of trellis plot

Hi, In my previous posting I forgot my system information. Sorry. It is listed below. I want to insert a key into each panel of a trellis plot, which I can do with a custom panel function that calles draw.key. The problem arises because I want the top right hand corner of the key to start

[R] location of key in panels of trellis plot

Hi, I want to insert a key into each panel of a trellis plot, which I can do with a custom panel function that calles draw.key. The problem arises because I want the top right hand corner of the key to start in the top right hand corner of the panel. If you run my code below, you can see that t

[R] understanding nlm

Hi, I am using the nlm() function to fit the following exponential function to some data by minimizing squared differences between predicted and observed values: csexponential<- function(x, t1, ti, p){ ti + abs(t1 - ti)*(exp(-(p*(x-1 } As background, the data is performance measured

[R] calling bwplot within a for loop

Hi, I am working with R 1.81. When I call bwplot() it prints the output to the windows device as it should. For example, d<-data.frame(y=c(2,3,4,5,12,14,16,11),x=c(rep("group1",4),rep("group2",4))) bwplot(y~x,data=d) This code results in a parallel boxplot. That is a single plot with

[R] nesting vectors in a dataframe

Hi, I want to create a dataframe where one of the columns does not hold individual elements, but a vector instead. For example, imagine the following dataframe without any vectors as elements . . . ex X1 X2 1 5 15 2 7 12 3