similar to: boxplot

Hi everyone,   I have got the following problem:   x1 <- rnorm(10,5,1) x2 <- runif(10,0,1) nam1 <- paste("A",1:4,sep=".") nam2 <- paste("A",6:9,sep=".") nam <- c(nam1,"A.4",nam2,"A.9") mydata <- data.frame(x1,x2) rownames(mydata) <- nam Error in `row.names<-.data.frame`(`*tmp*`, value = c("A.1",

Hi, I have gote the following data x1 <- c(rep(1,6),rep(4,7),rep(6,10)) x2 <- rnorm(length(x1),6,1) data <- data.frame(x1,x2) and I would like to compute the mean of the x2 for each individual of x1, i. e. x1=1,4 and 6? Thank you very much in advance, Amori [[alternative HTML version deleted]]

Hi everyone,   Is there any counter function in R for the following purpose: x <- matrix(c(1,1,0,2,1,0,0,2,0,1,2,1,2,1,0,1),nrow=4) As I would like to know how many zeros, ones, and twos in each row of x?   Many thank in advance, Amor     [[alternative HTML version deleted]]

Hello everybody,   I wonder if there is any swap function in R that does the following: x <- seq(1,10,12) x1 <- swap(x) x1 1 8 3 4 5 6 7 2 10 Thank you very much!   Amor __________________________________________________ Schutz gegen Massenmails. [[alternative HTML version deleted]]

Hi,   I have the following data and I would like to delete douple names, it is almost similar to SAS PROC SORT nodupkey! Is there any function in R does this?   x1 <- rnorm(11,5,1) x2 <- runif(11,0,1) nam <-paste("A", c(1:4,4,5:9,9), sep=".") mydata <- data.frame(x1,x2) crownames(mydata) <- nam   Many thanks in advance, Amor [[alternative HTML version

Hello,   I am using mod1 <- lrm(y~x1+x2,na.action=na.pass,method="lrm.fit") summary(mod1) and I've got the following error: Error in summary.Design(mod1) : adjustment values not defined here or with datadist for x1 x2   Many thank, Amor [[alternative HTML version deleted]]

Hi, R members. I want to save plots. So I use the R as fallows : -------------------------------- a <- c(1,10) b <- c(1,10) c <- c(2,20) d <- c(2,20) require(lattice) trellis.device("bmp", file = "test.bmp", bg = "white") print(plot(a,b)) print(plot(c,d)) dev.off() -------------------------------- However, I got a image of plot(c,d). I want to get a

Dear all, I am trying to use SAMr-library(samr), it gives me p-value = 1.001, any idea why? Many thanks, Amor __________________________________________________ z gegen Massenmails. [[alternative HTML version deleted]]

Hi, Can you please tell me what does the the function logged2 in R do im list or..? As I have > ?logged2 No documentation for 'logged2' in specified packages and libraries: you could try '??logged2' > ??logged2 No help files found with alias or concept or title matching ‘logged2’ using fuzzy matching. Thank you in advance, Amor [[alternative HTML version

Hi R-helpers. #I start with the reproducible example: firm<-c(rep(1,10),rep(2,10),rep(3,10),rep(4,10),rep(5,10)) year<-c(rep(1998:2007,5)) industry<-c(rep(1,20),rep(5,10),rep(7,10),rep(9,10)) X1<-rnorm(50) X2<-rnorm(50,mean=0.5,sd=0.1) Y<-rnorm(50,mean=0,sd=0.5) data<-data.frame(firm, industry,year,X1,X2,Y) data #I need to calculate for all the industries the following

I am making some changes to the dial plan at the request of the company president and have run into some problems. I have a couple of layers of menu's and I am not sure how to handle them. Here is how it should work (sorry for the crappy diagram) main menu --------Dial 1 for support | Dial 2 for special | Dial 3 sales

Hi, A'm about to set up a asterisk for 5000 users, and the customer had a 64bit server - can asterisk compile on that ? I will use a digium X100P for timing use will that do on a 64bit ? (I'm using SUSE91 kernel 2.6) What else ? Is it posible to have only one server for 5000 users ? I gues that it will be 5-700 sim. users only talking sip, and IAX2 to my PSTN-Gateway. The system is

Hi Group, I have a data frame below. Within this data frame there are samples (columns) that are measured more than once. Samples are indicated by "idx". So "id1" is present in columns 1, 3, and 5. Not every id is repeated. I would like to create a new data frame so that the repeated ids are averaged. For example, in the new data frame, columns 1, 3, and 5 of the original

Hi all, I downloaded latest CentOS 5.5 DVD i386 image from one of FTP's in a list. I've burned that image to DVD and created new DVD to use for fresh installations. Now when I install fresh CentOS 5.5 (in VM) I am getting info that there are 50 packages updates. This is ok when I have good internet speed (@work) but when I am home, this update takes a lot of time. I guess I can skip

Dear all, I am trying to shove a number of cmdscale() results into a single plot (k=1 so I'm trying to get multiple columns in the plot). From ?par I learned that I can/should set new=TRUE in either par() or the plot function itself. However with the following reduced code, I get only a plot with a column of data points with x==2. plot(1,10, xlim=range(0,3), ylim=range(0,10),

What SoftPhone working very well with *? S.O. is Debian Linux Thanks for your comments. JRR _________________________________________________________________ MSN Amor: busca tu ? naranja http://latam.msn.com/amor/

I have some code as shown below. Basically, I would like three boxplots to be set next to each other with no ylabels on the two "inner" plots, and I want the same x axis range on all three. However, it seems like boxplot does not respect the xlim setting. I've tried the various ways I thought would work (par, boxplot(...xlim=)) but none of them seem to work. I then tried plot.window,

hello, I try to run simulation of levene's test to find the p-value but the error of replacement has length zero occur, could anyone help me to fix this problem? asim <- 1000 pv<-rep(NA,asim) for(i in 1:asim) {print(i) set.seed(i) g1 <- rnorm(20,0,2) g2 <- rnorm(20,0,2) g3 <- rnorm(20,0,2) x <- c(g1,g2,g3) group<-as.factor(c(rep(1,20),rep(2,20),rep(3,20))) library(Rcmdr)

Dear List, I am trying to find each time a value changes in a dataset. The numbers are variables for day vs. night values, so what I am really getting is the daily sunrise and sunset. A simplified example is the following: x<-seq(1:100) y1<-rep(1,10) y2<-rep(2,10) y<-c(y1,y2,y1,y1,y1,y2,y1,y2,y1,y2) xy<-cbind(x,y) I would like to know each time the numbers change. Correct

Hi, Another one. I got notransfer=yes i iax.conf for both 2109 and dialout, but I still get this in my logfile Attempting native bridge of IAX2[2109@2109]/5 and IAX2[dialout]/6 Asterisk Version is CVS-04/19/04-22:17:41 What's wrong ? I gues it has somethnig to do withe my bilsec-problem as well. /HHA