similar to: anova subhypotheses

I simply put, > NEVER=subset(infants$bwt,ISNO1) > UNTILPREGNANT=subset(infants$bwt, ISNO2) > ONCENOTNOW=subset(infants$bwt, ISNO3) and I wanna combine those three. I do it like ISNO=NEVER&UNTILPREGNANT&ONCENOTNOW and R tells me 1: In NEVER & UNTILPREGNANT : longer object length is not a multiple of shorter object length 2: In NEVER & UNTILPREGNANT & ONCENOTNOW

I'm attempting the heruclean task of teaching myself Introductory Statistics and R at the same time. I'm working through Peter Dalgaard's Introductory Statistics with R, but don't understand why the answer to one of the exercises works. I'm hoping someone will have the patience to explain the answer to me, both in the statistics and R areas. Exercise 6.1 says: The zelazo data

Hi all, for the data below I would like to 1. generate a dummy variable for each group "gr" of the same composition by people, then save each portion in a separate file, 2. compute the frequency of "1"'s in "x" for each person by group "gr". So, "mike" will have freq=2/3, as he has two "1" and one "0" in 3 groups.

When I create a lattice/Trellis type graph, I typically write a function that returns the graph, as in do.graph <- function(x, y, ...) { require(lattice) return(xyplot(y~x, ...)) } My question today is this: If I want two graphs on one page, one way of achieving it is to print the objects into defined areas, as in gr1 <- xyplot(rnorm(111) ~ runif(111)) gr2 <-

Dear R-users, I would like to maximize the function g above which depends on 4 parameters (2 vectors, 1 real number, and 1 matrix) using optim() and BFGS method. Here is my code: # fonction to maximize g=function(x) { x1 = x[1:ncol(X)] x2 = x[(ncol(X)+1)] x3 = matrix(x[(ncol(X)+2):(ncol(X)+1+ncol(X)*ncol(Y))],nrow=ncol(X),ncol=ncol(Y)) x4 = x[(ncol(X)+1+ncol(X)*ncol(Y)+1):length(x)]

Folks, I'm running into something which looks like a bug in linearscan allocator. Of course I can't be 100% sure it's not some unobvious mistake on my part, so I'd like to hear your opinion. First, I attach two files -- LLVM asm and the asm for my target. The problem with assembler is: on line 171 it uses register gr2, which is copied from gr6 above, on line 161. The only

Hi, looking at the fresh CVS state I see: class Register<string n> : RegisterBase<n> { list<RegisterBase> Aliases = []; } while previously the Register class did not require any parameters. The change log is just: * Target.td: Revamp the Register class, and allow the use of the RegisterGroup class to specify aliases directly in register definitions. and I

On Tue, 22 Jun 2004, Vladimir Prus wrote: > First, I attach two files -- LLVM asm and the asm for my target. The problem > with assembler is: on line 171 it uses register gr2, which is copied from gr6 > above, on line 161. The only predecessor of this basic block is jump on line > 90. The problem is that gr6 is not initialized in the interval from the > function entry till the

Hello, In a study exploring transgenerational transmission of anxiety disorder we investigate whether infants react to experimentally induced mood changes of their mothers. We measured the time that an infant needed to cross a cliff (=crossing time) depending on whether his mother had previously undergone a mood induction (treatment) or not (control). The treatment is thus a

On Tue, 16 Nov 2004, Vladimir Prus wrote: > and I could not find any discussions in the archives. > > Why the change was necessary? Writing: > > def gr0 : Register<"gr0">; > def gr1 : Register<"gr1">; > def gr2 : Register<"gr2">; > def gr3 : Register<"gr3">; > def gr4 :

Dear expeRts, Why does na.blank=TRUE not replace the NA's in the following LaTeX table? x <- matrix(1:72, ncol=4, nrow=8) colnames(x) <- c("gr1.sgr1", "gr1.sgr2", "gr2.sgr1", "gr2.sgr2") rn <- apply(expand.grid(beta=c(0.25, 0.75), n=c(100, 500), d=c(10, 100))[, 3:1], 2, rmNames) x <- cbind(rn, x) # append columns containing the row labels

Dear RHelp-list, ?? Using dplyr and the group_by approach on a dataframe, I want to be able to apply a specific action according to the group name. The code bellow works, but I am not able to write it in a more esthetic way using dplyr. Can somebody help me to find a better solution ? Thank you Best regards Laurent df_test <- data.frame( x1=1:9, x2=1:9,

I was starting to write a note to a prospective R-user and came to the point of explaining how to get useful introductory information on R. After mentioning the Into and the FAQs I went on to try to explain how to use a lot of the contributed information. However I realised that there seems to be no direct way to get to Other Publications or Contributed Documenation. The best I have seen is

Hi everyone. I want to solve the following problem. I have a data.frame and I create a dotplot using lattice. Then I want to use the grid-package to create a combined graphic which contains the dotplot as well as a textplot() (using package gplots) of the data.frame next to the dotplot. Example code: library(lattice) library(grid) library(gplots) xx <- data.frame(f=factor(rep(1:5, each=5)),

Laurent: As I don't use dplyr, this won't help you, but I hope you and others may find it entertaining anyway. If I understand you correctly (and ignore this if I have not), there are a ton of ways to do this in base R, including using switch() along the lines you noted in your post. However, when the functions get sufficiently complicated or numerous, it may be useful to store them in a

Although there may well be many ways to do what is being asked for with the tidyverse, sometimes things are simple enough to do the old-fashioned way. The request seems to have been to do something to all rows in ONE specific group but was phrased in the sense of wanting to know which group your functionality is being called in. What grouping gains you is more worthwhile if you are interested in

Hi I have three steps like this: Given ''a user called "$name" of type $user_type and status $user_state in organisation "$organisation_name"'' do |name, user_type, user_state, organisation_name| $runs ||= 1 puts "MOO #{$runs}" $runs += 1 # ... end Given ''a user called "$name" with username

Hi, I wish to run the application "Roger Wilco" on my box. Is there anybody out there allready successfully running it ? Or are there negatives to report ? Please tell me if or if not it is possible Peter -- Ours is a world of nuclear giants and ethical infants. -- General Omar N. Bradley

TASK: problem with "glm" with binomial errors STATUS: Open FROM: p.dalgaard@kubism.ku.dk in glm(,binomial) it's possible that loss of significant digits make expected values 0 or 1 even though there's no divergence of the fit. (Happened to me with menarche data, infants and grown-ups included) [ Need the example data. Glm needs a

Hi, perhaps this is a stupid question, but i need some help about Helmert contrasts in the Cox model. I have a survival data frame with an unordered factor `group' with levels 0 ... 5. Calculating the Cox model with Helmert contrasts, i expected that the first coefficient would be the same as if i had used treatment contrasts, but this is not true. I this a error in reasoning, or is it