similar to: multiple lines and plot

Please, how to fill the area under the curve? x <- c(1:10) y <- c(rnorm(10)) plot(x,y) lines(x,y) Thanks, Mauricio Cardeal

Hi Does anyone know where is the package: npmc (Nonparametric Multiple Comparisons). I found the reference on R Site Search, but not the package itself on CRAN as suggested. Thanks Mauricio

Um texto embutido e sem conjunto de caracteres especificado associado... Nome: n?o dispon?vel Url: https://stat.ethz.ch/pipermail/r-help/attachments/20060711/1da40382/attachment.pl

Hi all, I'm trying to use lapply on a list with the following command: out<-lapply(mylist,myfun,par1=p,par2=d) (1) where myfun<-function(x,par1,par1) {.....} (2) now this function is in fact a wrapper for some Fortran code I have written so I think this might be the problem. When I call lapply() as in (1) I get the following message: Error in get(x,

Hi everyone ! Could you please help me with this problem ? I??ve trying to write a code that assign to a variable the content from another, but all I??ve got is a message error. For example: if (age <=10) {group == 1} else if (age > 10 & age <= 20) {group == 2} else {group == 3} Syntax error Or if (age <=10) {group == 1} else (age > 10 & age <= 20) {group == 2}

Hello. I asked before and it was great, cause as a beginner I learned a lot. But, if I have this in R (1 and 2 are codes for sex): > sex<-c(1,2,2,1,1,2,2,2) > sex [1] 1 2 2 1 1 2 2 2 I´d like to obtain the proportion according to sex.So I type: > prop.table(sex) [1] 0.07692308 0.15384615 0.15384615 0.07692308 0.07692308 0.15384615 0.15384615 [8] 0.15384615 The result is OK, but I

Hi all # I have this data set and how can I assign NA?s in just one command ? And why the summary(dat) function preserves the value 9 as real. ? x <- c(1,2,3,9,4) y <- c(3,6,9,2,3) z <- c(9,9,2,2,8) w <- c(6,5,3,0,9) dat <- cbind(x,y,z,w) summary(dat) x[x==9] <- NA y[y==9] <- NA z[z==9] <- NA w[w==9] <- NA summary(dat) summary(x) summary(y) summary(z) summary(w)

Please, I??d like to store this sample matrix as a new object. How can I do this ? pulse <- c(67, 67, 68, 68, 68, 69, 69, 69, 69, 69, 70, 70, 70, 70, 71, 71, 72, 72, 73, 74) m <- NULL x <- 0 for (i in 1:5) { x <- sample(pulse,3) m <- mean(x) cat(x,m,"\n") } Thanks, Mauricio

Is it possible to make disappear the id numbers from scatter.dudi (mc analysis) ? a <- as.factor(c(1, 2, 3, 2, 1)) b <- as.factor(c(3, 2, 3, 1, 1)) x <- as.factor(c(1, 2, 2, 1, 3)) y <- as.factor(c(2, 2, 3, 1, 1)) dat <- data.frame(a=a, b=b,x=x,y=y) summary(dat) dat require(ade4) dat.acm <- dudi.acm(dat, scann = FALSE, nf = 2) scatter.dudi(dat.acm) Thank you very much !

Hi everyone! Model: http://pastie.org/511290 Controller: http://pastie.org/511285 View new: http://pastie.org/511286 View _form: http://pastie.org/511284 When I hit submit button there''s an error for usuarios(users in portuguese) instance. It does not seems to save usuarios even though it seems to be posted by the form. Error: http://pastie.org/511292 Need some help here guys plz --

Hi everyone, I?m a student and relatively new to R so apologies in advance if this question seems stupid or obvious to you. I have collected a dataset with about 60 species of diatoms (count data from 19 different sample sites) and environmental variables for each site (salinity, pH, etc.). It?s all in the same dataset but distinct in R through the functions below diat <- diatom [, 1:60]

Hi, This is not a DTrace technical question (so, please don''t flame me for it) :) This is regarding the DTrace (un)conference TShirt. I am curious to know what engine is it. Could not figure it out. Thanks and regards, Sanjeev. -- Solaris Revenue Products Engineering, India Engineering Center, Sun Microsystems India Pvt Ltd. Tel: x27521 +91 80 669 27521

Hi, I wanted to calculate the age of people in my dataset by subtracting the individual's date of birth from their intake into a program. After several hours, searches of help archives, and the downloading of lubiridate, I have had no luck with this. Below is the code I used. > intakeDS$DOB <- as.character(intakeDS$DOB) > intakeDS$DOB <- as.Date(intakeDS$DOB,

Hello everyone, I would like to parse very large xml files from MS/MS experiments and create R objects from their content. (By very large, I mean going up to 5-10Gb, although I am using a 'small' 40M file to test my code.) My first attempt at parsing the 40M file, using the XML package, took more than 2200 seconds and left me quite disappointed. I managed to cut that down to around 40

Dear R-developers, I came to a somewhat unexpected behaviour of read.csv() which is trivial but worthwhile to note -- my data involves a protein named "1433E" but to save space I drop the quote so it becomes, Gene,SNP,prot,log10p YWHAE,13:62129097_C_T,1433E,7.35 YWHAE,4:72617557_T_TA,1433E,7.73 Both read.cv() and readr::read_csv() consider prot(ein) name as (possibly confused by

Dear all, Here is my code which am using to combine 5th column from different data sets. Here is the function to do my job genesymbol.append.file <-NULL gene.column <- NULL readGeneSymbol <- function(files,genesymbol.column=5){ for(i in fnames){ temp <- read.table(fnames,header=T,sep="\t",stringsAsFactors=F,quote="\"")

Hi, I have the following data file to be parsed and captured as a data frame: __DATA__ #GDS_ID GENE_NAME GENE_DESCRIPTION GENE_FUNCTION 1007_s_at | DDR1 | discoidin domain receptor tyrosine kinase 1 | protein-coding 1053_at | RFC2 | replication factor C (activator 1) 2, 40kDa | protein-coding 117_at | HSPA6 | heat shock 70kDa protein 6 (HSP70B') | protein-coding __END__ In particular it is

I have the following senerio. Call file calls 1st party. When connected give called party option to connect to second party. Issue Dial to second party. Caller answers and the two are bridged together. My issue is that 4 out of 5 calls fail to bridge the audio. Am I missing something or is there some kind of bug? Here is my test dialplan ;Dialer Base Code Files. ;Variables

Hi, I have the following data frame: > print(mydatframe) __DATAFRAME__ V1 V2 V3 1 1007_s_at DDR1 discoidin domain receptor tyrosine kinase 1 2 1053_at RFC2 replication factor C (activator 1) 2, 40kDa 3 117_at HSPA6 heat shock 70kDa protein 6 (HSP70B') __END__ Is there a way to create a hash with V2 as Key and V3 as its value? - Gundala Viswanath Jakarta - Indonesia

for the script, please kindly see the script below. At line 10 and line 13, my problems occurs. The first one is I try to retrieve the gene official name from a column of a table. The pattern of official name is something starting with gene_name. For detail problems, please see the according lines. Any suggestions are appreciated example of matching source (extract the Nnat, sometime it would