similar to: S4 slot containing either aov or NULL

Hi, In an effort to learn S4 objects, I am trying to port some c based tree code to S4 object. My immediate goal is to use .Call() interface for calling my c code from R. My long term goal is to understand how to write c structs that follows S4 classes and not the other-way-around. The c struct for the node is : typedef struct node{ int c; int n; inode **nodes; //length = n } inode; I

Dear Madame, Dear Sir, I am able to obtain the coefficients from a 'summary' of 'lm', but NOT from a 'summary' of 'aov'. The following example shows my steps. ## Initialize rm(list = ls()) # remove (almost) everything in the working environment utils::data(npk, package="MASS") # get data model <- yield ~ block + N*P*K ## Using lm npk.lm <-

Hi, I am trying to plot an interaction.plot with different color for each level of a factor. It has an erratic behavior. For example, it works for the first interaction.plot below, with the example from the ALDA book, but not with the other plots, from the NPK dataset: # from http://www.ats.ucla.edu/stat/r/examples/alda/ch2.htm tolerance <-

On 03/19/2016 01:22 AM, Michael Lawrence wrote: > > > On Sat, Mar 19, 2016 at 12:10 AM, Herv? Pag?s <hpages at fredhutch.org > <mailto:hpages at fredhutch.org>> wrote: > > On 03/18/2016 03:28 PM, Michael Lawrence wrote: > > > On Fri, Mar 18, 2016 at 2:53 PM, Herv? Pag?s > <hpages at fredhutch.org <mailto:hpages at

Hi all, Perhaps this in an inappropriate post, but I've found a bug in proj I'd like to track down a bit further before making a formal bug report. The example below shows the problem, if you change the rownames proj fails. The problem seems to be that there is a mismatch in the rownames in the qr objects constructed by aov and the rownames that proj is expecting them to have.

On 03/18/2016 03:28 PM, Michael Lawrence wrote: > > On Fri, Mar 18, 2016 at 2:53 PM, Herv? Pag?s <hpages at fredhutch.org > <mailto:hpages at fredhutch.org>> wrote: > > Hi, > > Short story > ----------- > > setClassUnion("ArrayLike", "array") > > showClass("ArrayLike") # no slot > >

Hi, This might be a minor bug, or it could be that I'm just mis-using the tools (again). If you modify the example for model.tables to introduce an NA, N <- c(0,1,0,1,1,1,0,0,0,1,1,0,1,1,0,0,1,0,1,0,1,1,0,0) P <- c(1,1,0,0,0,1,0,1,1,1,0,0,0,1,0,1,1,0,0,1,0,1,1,0) K <- c(1,0,0,1,0,1,1,0,0,1,0,1,0,1,1,0,0,0,1,1,1,0,1,0) yield <-

Hello Everyone, I am trying to excess the inbuit .Fortran and .C codes of R. Can any one help me in that. For example in kmeans clustering the algorithms are written in .Fortran I want to access them and see the .Fortran syntax of the codes. Can any one help me how can I do that? Thanx, Nitin Kumar On Thu, Nov 27, 2008 at 12:00 PM, <r-help-request@r-project.org> wrote: > Send R-help

I'm baffled. When I run make check after installing from source, I get a Error 2. From my understanding of how these things work, it would appear to be coming from this (as at the end of base-Ex.Rout.fail: > has.VR <- require(MASS, quietly = TRUE) Attaching package 'MASS': The following object(s) are masked from package:base : confint confint.lm nclass.FD nclass.scott

Dear statistitians / R users, I was told to analyse the effects of the four factors using lm or update for aov. The following is the result from lm. As you can see that in the last few lines enclosed coefficients. 1)Does this mean I can write a formula, like yield =0.26 + 0.03*H + 1.48*T + 0.04*L + 0.004*C 2) in the two levels design, is there any difference between I use (-1, 1) to

dear all, I try to update my Rversion fro 1.7.1 to 1.8.1 on linux. I have a problem when I try to install pacakges from Bioconductors with : source("http://www.bioconductor.org/getBioC.R") and getBioC(relLevel="release") somme errors append and when I make a library(affy) for example I obtain : > library(affy) Error in setIs("character",

The following R commands were written: >help.search("factorial") >data(npk) >npk >coef(npk.aov) In the output of coef command, please explain me the interpretation of coefficients of block1 to block 6 in this 2^3 confounded factorial experiment. Thanks. [[alternative HTML version deleted]]

Dear all, I have a question regarding setIs and method dispatch in S4 classes: Consider the following example: ##################################################### ## A02 "contains" A01 by setIs: setClass("A01", representation(a="numeric",b="numeric")) setClass("A02", representation(a="numeric",d="numeric"))

Hello, I noticed that model.tables fails when applied to an aov() fit if called inside a function. The problem seems to occur when as.formula is used inside a function on a string containing "<formula> + Error( x / y )" The reason I tried to use as.formula is to generate dynamic calls to aov(). Here is a minimal example illustrating the problem: ## Example test <-

Hi, I've been looking for Image Processing packages. Thresholding, Edge Filters, Dct, Segmentation, Restoration. I'm aware, that Octave, Matlab etc. would be a good address but then I'm missing the "statistical power" of R. Does anybody know of packages, projects etc. Comments on wether the use of R for such matters is useful are welcome. Greetings Thomas Kaliwe

Hi Martin, Thankyou for your response. I suspect that we're not going to agree on the main point. Making it trivially simple (as say Java) to set slots to NULL. So, I'll move on to the other points here. ***Note that cited text uses excerpts only.*** > setClassUnion("character_OR_NULL", c("character", "NULL")) > A = setClass("A", slots =

Dear colleagues, there is a behaviour with S4 (virtual) classes that I find very hard to understand: Depending on the position of setAs(), the tree of inheritance changes. This is my baseline example that defines the classes "grandma", "mother", "daughter" and a virtual class "mr_x". For a new instance if "daughter", "mr_x" is betweeen

Hi, Short story ----------- setClassUnion("ArrayLike", "array") showClass("ArrayLike") # no slot setClass("MyArrayLikeConcreteSubclass", contains="ArrayLike", representation(stuff="ANY") ) showClass("MyArrayLikeConcreteSubclass") # 2 slots!! That doesn't seem right. Long story ----------

Dear members, Is there a way to turn a strwidth of a string into a number of lines that ist needed to display the string when using par(mar = c(?,4,4,2)) #### x = 1:5 names(x) = c("ZZZZZ","ZZZZZzzzzz","ZZZlllll","TTTTT","Zzhtsddfg" ) par(mfrow = c(1,2)) par(mar = c(8,4,4,2)) barplot(x, las = 3, main = "8 is to much") par(mar =

Hi all, I have been trying to reproduce an analysis from Douglas Montgomery?s book on design and analysis of experiments. Table 6.10 of example 6.2 on page 246, gives a table as follows: > NPK <- expand.grid(A=mp,B=mp,C=mp,D=mp) > Rate <- c(45,71,48,65,68,60,80,65,43,100,45,104,75,86,70,96) > filtration <- cbind(NPK,Rate) > filtration A B C D Rate 1 - - - - 45 2