similar to: Linear models (lme4) - basic question

I do not know if I have a problem or not. The R script at the end of this email seems to run properly and a I get a boxplot that looks proper but I get the long string of messages during execution of the script looking like: ... The following object(s) are masked from dat ( position 8 ) : sugar The following object(s) are masked from dat ( position 9 ) :

Hi All, My apologies if this is a totally newbie question. I want to calculate the probability that a particular set of genes is picked repeatedly for 3 samplings. For example, if from a total of 'N' genes, I pick 'A' number of genes in the first pick, 'B' number of genes in the second pick, and 'C' number of genes in the third pick, I would want to calculate

Dear group, I have these 2 following data frame: allcon <- structure(list(DESCRIPTION = structure(1:17, .Label = c("COFFEE C Jul/10", "COPPER May/10", "CORN Jul/10", "CORN May/10", "COTTON NO.2 Jul/10", "CRUDE OIL miNY May/10", "GOLD Jun/10", "HENRY HUB NATURAL GAS May/10", "ROBUSTA COFFEE (10)

Dear group, Here is my data frame: futures <- structure(list(DESCRIPTION = c("CORN Jul/10", "CORN Jul/10", "CORN Jul/10", "CORN Jul/10", "CORN Jul/10", "LIVE CATTLE Aug/10", "LIVE CATTLE Aug/10", "SUGAR NO.11 Jul/10", "SUGAR NO.11 Jul/10", "SUGAR NO.11 Jul/10", "SUGAR NO.11

Dear group, I have 3 data frames I would like to merge. Here they are: pose16 <- structure(list(DESCRIPTION = structure(c(1L, 2L, 3L, 4L, 5L, 6L, 7L, 8L, 9L, 12L), .Label = c("COPPER May/10", "COTTON NO.2 Jul/10", "CRUDE OIL miNY May/10", "GOLD Jun/10", "ROBUSTA COFFEE (10) Jul/10", "SOYBEANS Jul/10", "SUGAR NO.11

Dear R users, i am using the following code to produce barcharts with lattice: Compound<-c("Glutamine", "Arginine", "Glutamate", "Glycine", "Serine", "Glucose", "Fructose", "Raffinose", "Glycerol", "Galacglycerol", "Threitol", "Galactinol", "Galactitol")

Dear group, Here is the kind of data.frame I obtain every day with my function : futures <- structure(list(DESCRIPTION = c("CORN Jul/10", "CORN Jul/10", "CORN Jul/10", "CORN Jul/10", "CORN Jul/10", "LIVE CATTLE Aug/10", "LIVE CATTLE Aug/10", "SUGAR NO.11 Jul/10", "SUGAR NO.11 Jul/10", "SUGAR

Dear group, Here is my df (obtained with a read.csv2()): df <- structure(list(DESCRIPTION = c("COTTON NO.2 Jul/10", "COTTON NO.2 Jul/10", "PALLADIUM Jun/10", "PALLADIUM Jun/10", "SUGAR NO.11 Jul/10", "SUGAR NO.11 Jul/10"), CREATED.DATE = c("13/05/2010", "13/05/2010", "14/05/2010",

Dear all, I am a total beginner in R, so sorry if this is the wrong place. I am using R 2.10.1 on a Mac (Mac OS 10.6.2). I have this small dataset : growth sugar 75 C 72 C 73 C 61 F 67 F 64 F 62 S 63 S I have no problem reading the table, or getting the summary, but if I try boxplot(growth~sugar, ylab="growth", xlab="sugar", data=Dataset), I have the following error :

Dear list, I have the following function, sugar = function(fun, id = "id"){ ff <- formals(fun) if( id %in% names(ff)) stop("id is part of args(fun)") formals(fun) <- c(unlist(ff), alist(id=)) fun } which one may use on a function foo, foo = function(x){ x } sugar(foo) # results in the extended closure, function (x, id) { x } Its limitation (other

Hi, Is anyone able to help Mr Stoner? -Ivo ---------- Forwarded message ---------- From: Ed Stoner <ed@whsd.net> Date: Oct 9, 2007 5:32 PM Subject: [sugar] video encoding To: sugar@lists.laptop.org I would like to provide the functionality of a "video export" from within an activity I am working on. I'm starting with a png for each frame of the video and would like to do

Dear list, I'm facing the following problem : A cox model with my sex variable interacting with several continuous variables : cph(S~sex*(x1+x2+x3)) And I'd like to make a nomogram. I know it's a bit tricky and one mights argue that nomogram is not a good a choice... I could use the parameter interact=list(sex=("male","female"),x1=c(a,b,c))... but with rcs or pol

Dear list, I am still struggling with something that should be easy: I compare two data frames with a lot of common rows and want to keep only rows that are NOT in both data frames, unique. Here are an example of these data frame. reported <- structure(list(Product = structure(c(1L, 1L, 1L, 1L, 2L, 2L, 3L, 4L, 5L, 5L), .Label = c("Cocoa", "Coffee C", "GC",

Dear group, I know this issue has been already covered, and before you reply I must say I have read the R-FAQ and search the mailing list archive. I still can't manage to change my factor to numeric as I couldn't find any clear answer. Here is my df : Pose1 <- structure(list(DESCRIPTION = structure(c(1L, 2L, 3L, 4L, 5L, 8L), .Label = c(" SUGAR NO.11 May/10 ", "COTTON

Hi, Just built the gem from trunk for the first time and ran into an issue with the generated spec.opts file (rev. 1493). Moving the file out of the way fixed the problem. I''m not on windows, so this is not a win32console gem issue as per the README. /usr/lib/ruby/gems/1.8/gems/rspec-0.8.0/lib/spec/expectations/sugar.rb:13:in `call'': undefined local variable or method

A bit new to wxRuby and have 2 questions regarding wxSugar: 1. I only find one sample application for wxSugar, but as I understand from the docs, wsSugar is the preferred approach. Are there any more wxSugar samples available somewhere? 2. When I run sugar-sample.rb, I get a couple of warnings: "listen is deprecated, use evt_xxx directly". Just wanted to check: Is this sample

So, it's advisable to grow cane as near to the factory as possible. He is the fourth person to resign at NFA. They sold it to Russia at virtually give away prices. uuuu WATCH LITL ON WEDNESDAY DEC 13! Company: L INTL COMPUTERS INC (Other OTC:LITL.PK) Symbol: LITL Mon Close: $0.99 5-day Target: $5 L Reveals 8-Way Processing Professional Mobile Computing Design CALL YOUR BROKER TO PLACE

So, it seems there are enough people on the plus side, I just wanted to make sure we evaluate all sides before taking a decision to add syntactic sugar to LLVM assembler. My main concern is still the same as earlier this year: the integrated assembler for ARM is still not complete, and the more extensions we add to the back-end, the harder it'll be to get it into production quality. That

Dear group, Here is my df : df3 <- structure(list(DESCRIPTION = c("COPPER May/10", "COTTON NO.2 Jul/10", "CRUDE OIL miNY May/10", "GOLD Jun/10", "ROBUSTA COFFEE (10) Jul/10", "SOYBEANS Jul/10", "SUGAR NO.11 Jul/10", "SUGAR NO.11 May/10", "WHEAT Jul/10", "SPCL HIGH GRADE ZINC USD",

Hi, I wanted to check the difference in results (using lme4) , if I treated a particular variable (beadchip) as a random effect vs if I treated it as a fixed effect. For the first case, my formula is: lmer.result <- lmer(expression ~ cancerClass + (1|beadchip)) For the second case, I want to do: lmer.result2 <- lmer(expression ~ cancerClass + beadchip) However, I get an error in