similar to: vector of center of

Hello everybody, I'm testing the randomForest package in order to do some simulations and I get some trouble with the prediction of new values. The random forest computation is fine but each time I try to predict values with the newly created object, I get an error message. I thought I was because NA values in the dataframe, but I cleaned them and still got the same error. What am I

Dear members, I'm not very used with Unix systems, but I got an old Sun/UltraSparc workstation running Solaris 2.5.1 and I tried to install R on it to able able to do some batch R processing while working on my pc. I downloaded and installed following required and recommanded programs before installing R : perl 5.8.0, readline 4.3, gzip 1.3.5, zlib 1.1.4, make 3.80, jpeg 6b, libpng 1.2.4

Hi, I'm using R to plot some data and associated linear models, and I was wondering if there was a parameter to fix the line width in the graphs (for axis, lines, ...) because I explored the associated help pages and could not find it. Thx in advance, -- Ir. Yves Brostaux - Statistics and Computer Science Dpt. Gembloux Agricultural University 8, avenue de la Facult? B-5030 Gembloux

Dear list, Here's a little problem I already solved with my own coding style, but I feel there is a more efficient and cleaner way to write it, but had no success finding the "clever" solution. I want to produce a factor from a subset of the combination of two vectors. I have the vectors a et b in a data-frame : > df <- expand.grid(a=c(0, 5, 10, 25, 50), b=c(0, 25, 50,

Hi, I have this table (BTW, published by FBI and representing age distribution of the WTC tragedy victims): "LABEL" "FREQ" "MIDPOINT" "1" "Infant (under 1)" 0 "2" "1 to 4" 5 "3" "5 to 8" 1 "4" "9 to 12" 3 "5" "13 to 16" 0 "6" "17 to

Hello ! I'm struggling again against lattice graprhics. ;) I'm trying to produce a conditionnal xyplot with two conditionning factors (let's say A and B). I want the levels of those factors (A1, A2, etc) to show in the margins of the lattice plot, not in the strips between the panels. A1 A2 A3 plot11 plot12 plot13 B1 plot21 plot22 plot23 B2 I managed to remove the

Hello I want to compare the slope (let's say 'b') of a linear model obtained with lm() to a theoretical value (let's say 'th'). To do so, I think I should compute a 't value' using something like : (b - 'th')/standard.deviation(b) and then look at the p-value of this computed t. I don't understand how to do this in a simple way, just using lm()

Hi, a questionaire data-set with more than one row for one Person make me problems. (1) i want create a function which count the used row for the relevant id ! example: id c.row 1 1 1 2 1 3 2 1 2 2 4 1 4 2 4 3 4 4 4 5 (2) i have got a problem to filter & analyse , because the filter criteria is unlucky a row above or under the ratings which i would like to

What is the recommended/optimal way to perform aggregates on data frames with chron objects? Here is an example: >raw.data 1 07/09/01 4000 2 07/09/01 2000 3 07/11/01 1000 4 07/13/01 800 5 07/13/01 700 6 07/16/01 600 7 07/17/01 500 I'm trying to construct a function that would first aggregate the data (second column) by day (grouping by the first column) according to a

Dear members, I'm using R to do some test computation on a set of parameters of a function. This function is included in three for() loops, first one for replications, and the remaining two cycling through possible parameters values, like this : for (k in replicates) { data <- sampling from a population for (i in param1) { for (j in param2) { result <- function(i,

Hello again, First I want to thank all the people who answered my question about line width in graphs. I promise I will learn the 'par' help page by heart for the end of the month ! I now want to trace some ellipses to emphasize groups of data. I found how to trace circles with 'symbols()', but no ellipse. I'm planning on writing my own function based on

Hello, I'm trying to plot a set of 3 dependant variables (y) against 4 predictors (x) in a matrix-like plot, sharing x- an y-axis for all the plot on the same column/line : y1/x1 y1/x2 y1/x3 y1/x4 y2/x1 y2/x2 y2/x3 y2/x4 y3/x1 y3/x2 y3/x3 y3/x4 In fact, this plot is a rectangular selection of the result of pairs(), limited to the relations between x's and y's

Dear list, I'm stuck with a little graphical problem. I'm generating several lattice plots which are printed in a single postcript device opened by > trellis.device(postscript, theme=canonical.theme("postscript", color=F), file="an_phase2_graph.ps", paper="a4", pointsize = 10, onefile=TRUE, horizontal=TRUE) Everything works fine,but some of these

Hello, I'm trying to analyse data from an incomplete design with four factor : - fr : number of the batch - op : ID of operator - meth : method used - mat : nature of the material used and one variable - mv : mesure - trmv : transformed mesure str(matvol) `data.frame': 120 obs. of 6 variables: $ fr : Factor w/ 30 levels "1","2","3","4",..:

Using R-2.12.1 and latticeExtra-0.6-14, I would like to understand why a rootogram displaying samples from the Poisson distribution looks like I expected it, whereas a rootogram using the normal distribution does not: library(latticeExtra) rootogram(~rpois(1000, lambda = 50), dfun = function(x) dpois(x, lambda = 50)) rootogram(~rnorm(1000), dfun = function(x) dnorm(x,mean(x),sd(x))) I

Hello again I'm trying to put error bars onto a barplot. I've tried something that Bill Simpson suggested a while ago, ie: x<-c(1,2,3,4,5) y<-c(1.1, 2.3, 3.0, 3.9, 5.1) ucl<-c(1.3, 2.4, 3.5, 4.1, 5.3) lcl<-c(.9, 1.8, 2.7, 3.8, 5.0) plot(x,y, ylim=range(c(lcl,ucl))) arrows(x,ucl,x,lcl,length=.05,angle=90,code=3) #or segments(x,ucl,x,lcl) but I can't get it to work on a

Full_Name: Gerald Guglielmo Version: 2.8.1 (2008-12-22) OS: OSX Leopard Submission from: (NULL) (131.225.103.35) When I attempt to use the hist() function to bin integers the behavior seems very odd as the bin boundary seems inconsistent across the various bins. For some bins the upper boundary includes the next integer value, while in others it does not. If I add 0.1 to every value, then the

Dear R-users, I need to produce a histogram where for every breaks there are the mean of the data. I tried tu use the function >hist(x, break=20 ... ) but this return the numerosity for every breaks, not the mean. Any hint? Thanks in advance, francesco --------------------------------------- Francesco Nutini CNR-IREA Ist. per il Rilevamento Elettromagnetico dell'Ambiente Via

Hi, I got two data point vectors. Now I want to make a ks.test(). I you print both vectors you will see, that they fit pretty fine. Here is a picture: http://www.jochen-bauer.net/downloads/kstest-r-help-list-plot.png As you can see there is one histogram and moreover there is the gumbel density function plotted. Now I took to bin-mids and the bin-height for vector1 and computed the

Is there a simple way to plot a histogram with colors? For example, suppose I generate random points in the N(2,1) distribution: x <- rnorm(100000, mean = 2, sd = 1) Now I would like to plot the histogram: hist(x) but I would like to show the bars with x < 0 in red, and the bars with x >= 0 in lightgreen. Is there any simple way to do it? I think I can do it in two steps: