similar to: merging data frames with matrix objects when missing cases

Hi R Users, I have been trying to work out how to rename column names using grep, basically I have generated these column names using tapply: [1] "NAME" "X1.1" "X2.1" "X3.1" "X4.1" "X5.1" "X6.1" "X7.1" "X8.1" [10] "X1.2" "X2.2" "X3.2" "X4.2"

First I have a matrix called stocks1: > class(stocks1) [1] "matrix" Here are the first 5 rows of the last 4 columns: > stocks1[1:5,2:5] [,1] [,2] [,3] [,4] [1,] 80.73 31.95 25.4 25.69 [2,] 83.66 31.95 27.12 25.2 [3,] 83.27 32.93 28.74 26.29 [4,] 83.9 34.07 29.77 26.6 [5,] 82.74 35.18 30.24 27.41 Now, why can't I convert this into a dataframe?

Is there an existing function to compute Cook-Weisberg confidence curves for non-linear model parameters in R? Kari Ruohonen -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html Send "info", "help", or "[un]subscribe" (in the "body", not the subject !)

How do I calculate the confidence interval for the value x given by the intersection of two quadratics (i.e. parabolas)? I fit two quadratics of the form: y = C1 + B1*x + A1*x^2 y = C2 + B2*x + A2*x^2 to two sets of points N1 and N2. I test for whether they intersect, if they do then I calculate the roots of: 0 = (C1 - C2) + (B1 - B2)*x + (A1 - A2)*x^2 to determine where they intersect

Hi, I wonder if predict.gam is supposed to work with family=negbin() definition? It seems to me that the values returned by type="response" are far off the observed values. Here is an example output from the negbin examples: > set.seed(3) > n<-400 > dat<-gamSim(1,n=n) > g<-exp(dat$f/5) > dat$y<-rnbinom(g,size=3,mu=g) >

I am trying to explore nlmer by running some nlme examples from Pinheiro & Bates (2000). I do not seem to find information how to specify fixed effects covariates to nlmer models. Specifically, I tried to run the "Carbon Dioxide Uptake" example from p. 368 onwards in the PB200 book. The model without fixed effects covariates runs well but how to tell nlmer to include Type and

Hi, I wonder if this is something on my machine locally or R in general. When I do the following: > plot(c(0,1),c(0,1),main=expression(paste(symbol("D"),"D",sep=""))) I get a plot with a title having uppercase delta followed by "D". But in the following > pdf(file="deltaTest.pdf") >

Hi - Are there any known bugs or other issues that may cause R to crash when trying to use xyplot()? For example, > x<-1:100 > y<-rnorm(100) > library(lattice) Loading required package: grid > xyplot(y~x) causes this: Process R segmentation fault at Wed Apr 3 16:56:42 2002 I am running linux debian unstable on i386. R says it is R 1.5.0 in the header text when starting but

Hi, I wonder if someone could help. I needed to transfer (copy) a workspace file that had been generated in linux (R 2.11) to windows running the same version of R 2.11 (but of course windows binary). Usually, there is no problem in doing this and all objects work as expected. I am often doing this to be able to produce wmf or emf graphic files that I need. This time I had some spectra that I

Hi, I would appreciate of someone could explain how the residual standard error is computed for rlm models (MASS package). Usually, one would expect to get the residual standard error by > sqrt(sum((y-fitted(fm))^2)/(n-2)) where y is the response, fm a linear model with an intercept and slope for x and n the number of observations. This does not seem to work for rlm models and I am wondering

Le jeudi 09 f?vrier 2012 ? 15:31 +0200, Kari Ruohonen a ?crit : > Hi, > I am trying to replicate a training of AttributeSelectedClassifier with > CFsSubsetEval, BestFirst and NaiveBayes that I have initially done with > Weka. Now, I am trying to use RWeka in R. > > I have a problem of passing arguments to the CfsSubsetEval, BestFirst > and NaiveBayes. I have first created an

Hi, I get the following error message when trying to use the display function on the ARM package: > display(model) Error in .Internal(round(x, digits)) : no internal function "round" Looks like some kind of mismatch between the ARM package and some others? Can I somehow get around it? I have learned to like the display function to print model summaries. Here is my sessionInfo():

This is maybe not directly an R problem but I have used R to try to solve it so I think somebody may be able to help. I have a mixture model with three components and a quadratic Scheffe polynomial p1x1+p2x2+p3x3+p12x1x2+p13x1x3+p23x2x3 fitted to the response. Now I'd like to compute the mixture corresponding the maximum response. Model for Y1 has the parameters p1=124.02 p2=60.973 p3=41.479

I have a fresh installation of R from debian unstable packages. The html index found in /usr/lib/R/doc/html/index.html works in Mozilla and under the link of 'packages' on this page I have a list and corresponding links. However, it appears that not all packages I have installed from available deb files via apt-get, will get their link updated to this package index. Specifically, from

i have a data frame(dat) which has many variables.and i use the following script to get the crosstable. >danx2<-c("x1.1","x1.2","x1.3","x1.4","x1.5","x2","x4","x5","x6","x7","x8.1","x8.2","x8.3","x8.4","x11",

I am new to generalized linear models and studying McCullagh & Nelder (1989). Especially, I have a problem resembling the \"cheese taste\" example (5.3.1. p. 109) of the book. I tried to analyse the cheese example with R but failed to do so because R allowed me to use logit link function only with binary family that supposes 0 <= y <= 1. Do I need to scale the y\'s or

Hi all, I have a question regarding the boxplot function. The data I am working on has 1 grouping variable (G) and it has many numerical variables (V1, V2, V3, V4, Vx, etc). What I would like to do is create a boxplot where the Y-axis represents the numerical values of variable V1...Vx (all the variables have the same range). The X-axis needs to represent the G-V combination. So suppose the

Full_Name: Giovanni Bonafe' Version: 2.2.0 OS: Linux Submission from: (NULL) (195.62.164.225) If the file "example.dat" is like this: aaa bbb ccc 3.4 1.2 5.6 4.6 10 32 667 343 1.7 With the older 1.9.1, as expected: > data<-read.fwf(file = "example.dat",widths=c(3,4,4),header=TRUE) > data aaa bbb ccc 1 3.4 1.2 5.6 2 4.6 10.0 32.0 3 667.0 343.0

Hi, You might need to check library(boot).? I have never used that before.? So, I can't comment much.? It is better to post on R-help list.? I had seen your postings on Nabble in the past.? Unfortunately those postings were not accepted in R-help.? You have to directly post at ? r-help at r-project.org after registering at: https://stat.ethz.ch/mailman/listinfo/r-help ?

I've sometimes thought it would be useful to have another mail list called "no I haven't read the documentation but someone can save me a lot of time," or perhaps just "newbie" for short. This can be very useful for people just getting started in R and with school starting soon I expect the help list may become overwhelmed. Of course, I think that people should read all