similar to: boxplot help

After checking the original data in Excel for blanks and running Summary(cm3) to identify any null values in my data, I'm unable to identify an instances. Yet when I attempted to use the data in Random Forest, I get the following error. Is there something that Random Forest is reading as null which is not actually null? Is there a better way to check for this? > library(randomForest) >

Dear R users, I have the following two problems, related to the function sub, grep, regexpr and similia. The header of the file(s) I have to import is like this. c("y (m)", "BD (g/cm3)", "PR (Mpa)", "Ks (m/s)", "SP g./g.", "P (m3/m3)", "theta1 (g/g)", "theta2 (g/g)", "AWC (g/g)") To get rid of spaces and

Please I need some help using R to analyze my data. What I would like to do is to repeat the same basic process (e.g. linear regression between wood density and distance from pith) for at least 240 data subsets within the main data-frame. Within the main data-frame, these data subsets will be defined by three variables namely,  species, individual and core (i.e. 20 species, at least 6

Hi, dear all, I want to show numerical values with decimal points on the boxplots. Here is what I did: For example; x1<-rnorm(100,2,2); x2<-rexp(100); label<-rep(1:2, rep(100,2)) median<-round(c(median(x1), median(x2)),3) boxplot(c(x1, x2)~label, medpch=paste(median), medcex=1.2) It only shows the integers at the median position in the boxplots. How to make it show more decimal

Dear all, I am a total beginner in R, so sorry if this is the wrong place. I am using R 2.10.1 on a Mac (Mac OS 10.6.2). I have this small dataset : growth sugar 75 C 72 C 73 C 61 F 67 F 64 F 62 S 63 S I have no problem reading the table, or getting the summary, but if I try boxplot(growth~sugar, ylab="growth", xlab="sugar", data=Dataset), I have the following error :

Dear Wiza[R]ds, I have a data.frame header that looks like this: v2FfaPre15 v2FfaPre10 v2FfaPre5 v2Ffa2 v2Ffa3 v2Ffa4 I need it to look like this, 15 10 5 2 3 4 i.e., with v2FfaPre and v2Ffa stripped off Any suggestions, Thanks in advance! -- Oscar Oscar A. Linares, MD Translational Medicine Unit LaPlaisance Bay, Bolles Harbor Monroe, Michigan 48161

Some of you may remember back in May the thread on using Asterisk as a voicemail server for a Cisco Callmanager system. My own Callmanager system is integrated into an Asterisk server for voicemail (and other things). Back in May I was using H323 for integration, but since I've upgraded to CCM 4.1 I have switched over to SIP. The integration with H323 required using Call forwarding to send

Dear Rexperts, Given, aar <-function(command) { switch(command, {scrn = cat("scrn :Screening","\n")} {dx = cat("dx :Diagnosis","\n")} {df = cat("df :Don't Forget","\n")} ) } I want to be able to do: aar("dx") # function does cat("dx :Diagnosis","\n")

Dear Wiza[R]ds, I have the following data in a data.frame. I need to do an anova with multiple comparison but I don't know how to factor the groups for analysis. There are 3 groups, 1,2 and 3 labelled in column 1. Help appreciated with thanks in advance. Group SI Sif SG Io I2 lol2 1 9.08e-05 9.08e+00 0.060842287 1.798556446 32.574500

Dear Wiz[R]ds, I am deeply grateful for the help from Duncan Murdoch, Gray Calhoun, and others. We are almost there. For whatever reason, I can't change the symbol from a circle to a triangle in the upright posture plots. Any ideas? I have included the problem in full. # tritiated (3H)-Norepinephrine(NE) disappearance from plasma # concentrations supine and upright # supine datasu <-

I have a data set for different time intervals. The data has three comment lines before data for each time interval. For each time interval there are 500 data points. I want to change the dataset such that I have the following format: t1 t2 t3 ................ 0.00208 0.00417 0.00625 ................. a1 a2 a3 ...................

No objections from Martin or elsewhere, so I have now committed this. Thanks, Uwe. On Tue, 24 Feb 2009, Uwe Ligges wrote: > [CCing Martin and Brian who had both done most svn commits of boxplot.R so > far] > > > A very minor wishlist item that I should have already reported years ago: > > All the time when I need presentation/publication quality boxplots, I add >

Dear Wiza[R]ds, # I have the following experimentally observed data: csdata <- data.frame( time=c(0,1,3,9,20), conc=c(638.697,395.69,199.00,141.58,112.16) ) # weighting resp means response wt.MM<- function(resp, time,A1,a1,A2,a2) { pred <- A1*exp(-a1*time)+A2*exp(-a2*time) (resp - pred) / sqrt(pred) } # Fit using nls cs.wt <- nls( ~ wt.MM(conc, time,A1,a1,A2,a2),

>[Rd] boxplot() defaults {was "boxplot in extreme cases"} >Martin Maechler maechler at stat.math.ethz.ch >Mon Nov 8 10:36:42 CET 2004 > > AndyL> Try: > > AndyL> x <- list(x1=rep(c(0,1,2),c(10,20,40)), > x2=rep(c(0,1,2),c(10,40,20))) > AndyL> boxplot(x, pars=list(medpch=20, medcex=3)) > > AndyL> (Cf ?bxp, pointed to from

Dear Customer! This is an automatic response from our support system, because you have sent a message to our support system e-mail address, which does not refer to an existing support case tracking number. We regret to inform you that we can no longer process such messages and your e-mail has, therefore, NOT been delivered. We have recently implemented a new support case management system in

Dear Wiza[R]ds, I am very grateful to Duncan Murdoch for his assistance with this problem. His help was invaluable. However, the problem has become a little more complicated for me. Now, in each plot, I need to plot the observed and fitted values of a supine and upright posture experiment. Here is what I have and how far I got. # tritiated (3H)-Norepinephrine(NE) disappearance from plasma #

Hello, I have some trouble using step() and stepAIC() functions. I'm predicting recruitment against several factors for different plant species using a negative binomial glm. Sometimes, summary(step(model)) or summary(stepAIC(model) does not select the best model (lowest AIC) but just stops before. For some species, step() works and stepAIC don't and in others, it's the opposite.

Hi, I am trying to find m breakpoints of a linear regression model. I used the segmented package. It works fine for small number of predicators and breakpoints.(3 r.v. 3 points). However, my model has 14 variables it even would not work even for just one breakpoints!. The error message is always estimated breakpoints are out of range. Since my problem is time related problem. So I

Full_Name: Karl W Broman Version: 1.7.1 OS: Linux Submission from: (NULL) (162.129.44.18) There appears to be a slight problem with round(); for rounding off a 5, it should go to the even digit, but there appear to be some cases where it goes to an odd digit; see below. Thanks! karl ~[1007]$ R --vanilla R : Copyright 2003, The R Development Core Team Version 1.7.1 (2003-06-16) R is free

Dear all, I'm trying to use the fold function as described here: http://cran.r-project.org/doc/contrib/Fox-Companion/appendix-cox-regression.pdf Page9 It does say that you can use this when you have more than one time varying covariate: in the description of the argument cov it says: "cov: A vector giving the column numbers of the time-dependent covariate in data, or a list of