similar to: how to add in interaction terms in gamm

Hi, I have the following data: > Myvalues Gene ES MEF Embryo ESHyp 1 GeneA -0.38509507 0.00 1.6250 1.7039921 2 GeneB 0.06262914 0.00 1.6250 -0.272033 and so on... I want to plot the expression values of GeneA and GeneB in the different cell/embryo/conditions (columns 2:5 above). Now, if I do: >library(ggplot2) > qplot(x=Gene, Embryo, geom =

Hi, I have a dataframe "gexp_1" with 115 rows of samples and 27000 columns of gene expression measurements with each column corresponding to one gene. I now have a smaller vector of genes "tp" for which I need to pull out the data form the data frame. # first gene from the list with 16 genes x <- tp[1,] x [1] geneA 16 Levels: geneA, geneB... #unsuccessful gexp_1$x NULL

Dear All, I have a following for-loop code which is basically intended to read in many excel files (each file has many columns and rows) in a directory and extract the some rows and columns out of each file and then combine them together into a dataframe. I use for loop which can do the work but quite slow. How to make it faster using lapply function ? Thanks in advance! temp.df<-c() #

Hi Stephen, I am not sure that the "for loop" is the source of slowness. You seem to be doing a lot of unnecessary work each time through the loop. e.g. no need to check if it's the last file, just move that section outside of the loop. It will be executed when the loop finishes. As it is you are calling list.files() each time through the loop which could be slow. In any case

Eric's approach seems reasonable to me, and I agree that it's probably not the use of a "for" loop that makes the original version slow. As Eric mentioned, there are lots of unnecessary things happening in the loop. For example, list.files() was called twice inside the loop, which is unnecessary, and will definitely slow the loop down (though probably not by much). Call it

It seems that confint.default returns an empty data.frame for objects of class mlm. For example: ``` nobs <- 20 set.seed(1234) # some fake data datf <- data.frame(x1=rnorm(nobs),x2=runif(nobs),y1=rnorm(nobs),y2=rnorm(nobs)) fitm <- lm(cbind(y1,y2) ~ x1 + x2,data=datf) confint(fitm) # returns: 2.5 % 97.5 % ``` I have seen proposed workarounds on stackoverflow and elsewhere, but

>>>>> steven pav >>>>> on Thu, 19 Jul 2018 21:51:07 -0700 writes: > It seems that confint.default returns an empty data.frame > for objects of class mlm. For example: > It seems that confint.default returns an empty data.frame for objects of > class mlm. Not quite: Note that 'mlm' objects are also 'lm' objects, and so it is

Dear Users, I am new to both of things, so do not blame me too much... I am busy with semiparametric regression and use WinBUGS to sample posteriors. The code to call Winbugs is as follows: data <- list("y","X","n","m") #My variables inits.beta <- rep(0,K) inits.beta0 <- 0 inits <-

Following on from my previous mail! plot(1:10,1:10) mtext("title", side=3, adj=0, font=3, cex=1.5) This works as expected and puts the font in italics. tag <- "A)" suffix <- "genea::" plot(1:10,1:10) mtext(bquote(.(tag) ~ Delta * .(suffix)), side=3, adj=0, font=3, cex=1.5) Here, the font isn't in italics, it's normal. I presume this is some

I have 2 files containing data analysed by 2 different methods. I would like to find out which genes appear in both analyses. Can someone show me how to do this? _________________________________________________________________ [[trailing spam removed]] [[alternative HTML version deleted]]

I try to display the grouping variable in coplot. It work, but it's special solution and rather ugly. Any better idea? # Simulate my data frame data(state) x77<-data.frame(state.x77) x77$region<-state.region coplot(Life.Exp ~ Income | region, data=x77, show.given=F, subscripts=T, panel = function(x, y,subscripts, ...) { panel.smooth(x, y, span = 1., ...)

Hello list, I was going to buy this weekend a Wisip from http://www.pulverinnovations.com/, but jeff got out of stock and he wont have Wisip for the next 3 to 4 weeks. So I start searching for other wifi phones because I was really upset about it and I found IPC5000 from http://www.fahdtel.com/sip_phones.htm, I liked so much the pic that I email the guy and he send me the PDF with all the details