similar to: Histogram/Bar plot graph

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

I want to use gamm to generate smoothed trend line for three groups identified by dummy variable genea and geneb. My question is how to add in an interaction term between the time and another dummy variable such as gender? fitm<-gamm(change_gfr~ genea+geneb+s(timea_n,bs="ps")+s(timeb_n,bs="ps")+s(timec_n,bs="ps"),data=mm,random=list(time_n=~1|PID)) -- View this

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

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

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() #

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

-----BEGIN PGP SIGNED MESSAGE----- Hash: SHA1 - -----BEGIN PGP SIGNED MESSAGE----- Hash: SHA1 The error codes for' doveadm user' should be improved to avoid error handling workarounds when scripting commands: $ doveadm user doveadm user [-a <userdb socket path>] [-x <auth info>] <user mask> [?] $ echo $? 1 $ doveadm user foo userdb lookup: user foo doesn't exist $

> -----Original Message----- > From: William Dunlap > Sent: Tuesday, November 25, 2008 9:16 AM > To: 'johannes_graumann at web.de' > Subject: Re: [R] Efficient passing through big data.frame and > modifying select fields > > > Johannes Graumann johannes_graumann at web.de > > Tue Nov 25 15:16:01 CET 2008 > > > > Hi all, > > > >

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''m seeing these in my logs on my puppetmaster: puppetmasterd[4111]: Could not store configs: SQLite3::SQLException: SQL logic error or missing database: DELETE FROM fact_values WHERE "id" = 8 currently this is 0.24.0 relect from some previously collected resource that wasn''t cleaned up? upgrade issue? nothing to worry about? that''s the only

Greetings, I have a question that I'd like to get input on. I have a classic toxicology study where I artificially fertilized and exposed embryos to a chemical and counted defects. In addition, I kept track of male-female pairs that I used to artificially fertilize and generate embryos with. I need to use logistic regression to model the response, but also check that the genetics of the

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., ...)

I forgot the answer for this: in http://wiki.dovecot.org/AuthDatabase/PasswdFile, we have the passwd file defined as user:password:uid:gid:(gecos):home:(shell):extra_fields What should the uid and gid be associated with? In my setup I told postfix to use Dovecot LDA's deliver, which then is saving files as uid:gid of 1500:1500. How would that affect the uid:gid used in the passwd file?

Hello, I use the function 'aggregate' a lot. One small annoyance is that it is necessary to name the factors in the 'by' list to get the names in the resulting data.frame (else, they appear as Group.1, Group.2...etc). For example, I am forced to write: aggregate(y,list(f1=f1,f2=f2),mean) instead of aggregate(y,list(f1,f2),mean) (for two factors with short names, it is not such

I''m running into a very odd occurrence. I started with a new puppet server and was able to add a few classes, however as I try and add new ones I get odd parsing error messages, the most common being Could not match ''}'' at whatever file. Here is the exact error I get err: Could not retrieve catalog: Could not parse for environment production: Could not match

3 years ago (!) I aksed about this and finally I implemented it today: >> # NEW >> passdb { >> driver = passwd >> # firstname.lastname at charite.de:password:uid:gid:(gecos):home:(shell):user=realloginname >> args = /usr/local/etc/dovecot.login_via_email >> } >If you're only using it as a passdb, you don't need the uid/gid/home.

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

All, I'm having problems w/ a simple attempt to subset an xyplot. The first plot below is a plot of y versus x for certain values of a third categorical variable z. Now I'd like to further restrict this to certain values of variable y. Neither of the two attempts below work. Any suggestions much appreciated. (note: I don't want to merely use ylim since I have a loess plot and I

All, How does one replace plot symbols with say subject IDs when using xyplot? Or superimpose them next to plot symbols? I searched the archives under various key words but haven't had much. Any suggestions or links much appreciated. Sample code below. David junk.frm = data.frame(ID = rep(1:16, each = 2), x, y, z = rep(c("D", "P"), 16)) y = c( 0.4, 0.6, -0.1,

Hi, I'm having a problem with dovecot and postfix. I'm using a dovecot version 2.0.13. Postfix (Version: 2.8.5) config is changed according to the http://wiki2.dovecot.org/LDA/Postfix: I have thought to build a virtual user and virtual mailbox, and use the functionality of the Sieve. However, when the test set, the mail not received the mailbox. Record that there was a "temporary