similar to: height of plots

I have a data frame like this xx<-data.frame(cbind(Sample=c('Ctrl_6h','1+0_6h','1+200_6h','1+5k_6h','Ctrl_5K_6h','ConA_6h'), IFIT1=c(24,25,24.7,24.5,24.2,24.8))) grep('[[:digit:]]h',xx$Sample) yy<-xx$Sample strsplit(yy,"_") I have to extract the time information separated by '_' in the sample

Hello gurus, I have some variables, and i am creating combinations for analysis in the end i need these variables to be displayed like "LEPTIN+SAA+PTH". currently i am using loop to perform this. I would appreciate any pointers to do it without the loop. > mols=c("LEPTIN","SAA","PTH","sEGFR") > samples=mols[1:3] > samples [1]

I would like to know if clogit function can be used as below clogit(group~., data=dataframe) When I try to use in above format it takes a long time, I would appreciate some pointers to get multiple combinations tested. set.seed(100) d=data.frame(x=rnorm(20)+5, x1=rnorm(20)+5, x2=rnorm(20)+5, x3=rnorm(20)+5, x4=rnorm(20)+5, x5=rnorm(20)+5, x6=rnorm(20)+5, x7=rnorm(20)+5, x8=rnorm(20)+5,

so here is the code I need help to fix the line for 'st' when it reaches maximum, that's where i am getting stuck. thanks Sharad set.seed(100) d=data.frame(x=rnorm(20)+5, x1=rnorm(20)+5, x2=rnorm(20)+5, x3=rnorm(20)+5, x4=rnorm(20)+5, x5=rnorm(20)+5, x6=rnorm(20)+5, x7=rnorm(20)+5, x8=rnorm(20)+5) r2=c() cinit=0 for (col in 1:9) { cinit=cinit+1 st=cinit+1

Hello all, I am using 'plot' to create standard curves for elisa data. when I use 'plot' with type='b' i get the points connected with lines and one straigth line from the lowest datapoint to the highest data point. how can i avoid/remove it from the figure. i am using R2.9.1, below is the example of the data. od<-c(10, 8, 6,4,2,1, 10.5,7.8,6.4,3.8,2.1,0.95)

Hello Gurus I have two correlation matrices 'xa' and 'xb' set.seed(100) d=cbind(x=rnorm(20)+1, x1=rnorm(20)+1, x2=rnorm(20)+1) d1=cbind(x=rnorm(20)+2, x1=rnorm(20)+2, x2=rnorm(20)+2) xa=cor(d,use='complete') xb=cor(d1,use='complete') I want to combine these two to get a third matrix which should have half values from 'xa' and half

Hello all, I have arranged my data as per Dennis's suggestion in this post http://www.mail-archive.com/r-help at r-project.org/msg107156.html. the posted code works fine but when I try to apply it to my data, i get "> u2 <- ddply(xxm, .(plateid, cytokine), as.data.frame.function(f)) Error in t.test.formula(conc ~ Self_T1D, data = df, na.rm = T) : grouping factor must have

Hello gurus, I have a dataframe containing two groups viz., 'control' and 'case', each of these groups contains longitudinal data for 100 subjects. I have to plot all these subjects on a single chart and then put a regression line for each of the group for all the subjects. I have written a function to do the chart grpcharts<-function (dat, group,group2,molecule,cutoff){

hello, i try to plot 3 graphs which have the same x.axis underneath each other. i'd like the plots to be aligned without margings between the boxes and draw a single x axis beneath the lowest plot. i managed to get the alignment by setting par(mar), but the middle box is stretched and i cant't figure out how to get around this. par(pin) was my guess, but this doesn't do the job

ear All, I configured Postfix 2.3.3-2 on CentOS 5.2 using SASL Dovecot authentication put I am facing the problem to receive mails . My postconf -n result is: [root at mails ~]# postconf -n alias_database = hash:/etc/aliases alias_maps = hash:/etc/aliases command_directory = /usr/sbin config_directory = /etc/postfix daemon_directory = /usr/libexec/postfix debug_peer_level = 2

Hello, I have data in a dataframe with 139104 rows which is multiple of 96x1449. i have a phenotype file which contains the phenotype information for the 96 samples. the snp name is repeated 1449X96 samples. I haveto merge the two dataframes based on sid and sen. this is how my two dataframes look like dat<-data.frame(snpname=rep(letters[1:12],12),sid=rep(1:12,each=12),

hello gurus, i have a data frame like this HTN HTN_FDR Dyslipidemia CAD t1d_ptype[1:25] 1 Y Y Y T1D 2 T1D 3 Ctrl_FDR 4 T1D 5 Y Ctrl 6 Ctrl 7

If I use the following command to plot points: plot(c(1,2,2,3,3,3), type="p", pch=20, ylab="Y Label", xlab="X Label", xaxt="n") there is a large amount of space between the label "X Label" and the actual x-axis. If I change the xaxt="n" to xaxt="s", the label "X Label" don't move at all. Is there a way to get

I have a data file where each entry is indexed by the time in seconds since epoch (e.g. 1046315697). Is there an easy way to convert this time value into a more friendly time (such as Month-Year) when plotting it? I searched through the manual, mailing lists, and functions like as.POSIXct and strptime, but didn't find what I need. Thanks, Sharad.

Hi, I'm using linux 2.6.21.5 and our kernel is freeze. The problem is, if I create a Software bridge using $brctl command. and add two interfaces say, eth0.0 and eth0.1 using $brctl addbr br-lan $brctl addif br-lan eth0.0 $brctl addif br-lan eth0.1 and when i send traffic from a host connected to one port to host connected at other end, soon all the memory is dried up and and kernel

Hi All, I am new to Ruby and ROR I were trying to create small DSL for conditional validations valid_with_cond :bypass_validation do if self.addresses > 3 errors[:base] << "Can not have more than 3 addresses". end end By this I wanted to create array of method and call them all in custom validation method. this above code I wanted to do attr_accessor

I have written a wrapper for heatmap.2 called heatmap.w.row.and.col.clust which auto-generates breaks using breaks<-round((c(seq(from=(-20 * stddev), to=(20 * stddev))))/20, digits = 2) #(stddev in this case = 2.5) This has always worked well in the past but now I am getting an error that non-finite breaks are being generated. Drilling down, it seems that my wrapper is generating finite

Hi, I am using virConnectOpen(const char *uri) to connect to the local hypervisor. When uri is set to "lxc:///", an error gets logged in /var/log/messages about failed connection. But I do not see any error for "xen:///". The machine does not have these drivers installed, and I am using libvirt 0.6.3. My questions are - 1. Why I see error in /var/log/messages for container

Hi, I have a host with 6 VMs set to autostart. When this host reboots, it starts libvirtd, pegasus cimom and libvirt-cim. It also restarts the VMs. Intermittently we noticed a core during start up that points to the last call made by libvirt-cim to "virInitialize". I am thinking that it is a timing issue where virInitialize is called while libvirt is trying to autostart the VMs.

Hello I have the following data.frame and xyplot. I need this plot without the borders. Does anybody know how to get this xyplot without borders o with white borders? Thank's in advance Juan Hernández my.label <- data.frame(quantil=rep(20,8), my.factor=factor(c('FA','FB','FC','FD','FE','FG','FH','FI'))) library(lattice)