similar to: Segmentation fault with xyplot

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 !)

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 have faced a problem with the merge() function when trying to merge two data frames that have a common index but the second one does not have cases for all indexes in the first one. With usual variables R fills in the missing cases with NA if all=T is requested. But if the variable is a matrix R seems to insert NA only to the first column of the matrix and fill in the rest of the columns by

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") >

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, 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

Hello >From the following script I received a grafic output when I called: - xyplot.test( 'green3' ) - call.xyplot.test( 'blue3' ) I did NOT receive a grafic output when I called: - loop.xyplot.test( 'red3' ) What's the Problem? NB: I am using R 1.4.1 on Linux. --------- START OF SCRIPT ---------------- n <- 1000 x <- seq( 1, n ) y <- rnorm( n )

I know I should not mix base plotting functions with grid/lattice functions, but I have used a "quick-and-dirty" trick of par(new=T) in the past for annotating a trellis-drawn graph in various versions of S-PLUS. The sequence goes something like this: > windows(width = 5, height = 5, pointsize = 10) # open up the device > xyplot(y ~ x) > par(new=T) > xyplot(y2 ~ x) >

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

Hi, Suppose I have a matrix with, say 12 columns. I would like to extract out column 2 ~ 8 and 11 ~ 12 then combine them together. I tried with success to extract out the columns by doing: foo <- test[,2:8] goo <- test[,11:12] However then I am having trouble combining foo and goo. Helps are appreciated! Cheers, Kevin ----------------------------------- Ko-Kang Kevin Wang

Hi R experts, I'm using the xyplot function in lattice to draw a multipanel plot consiting of 5x6 scatterplots. Now I need to link single points in each of those scatterplots (=panel),but the points, that need linking are different for each panel. I tried to use the panel.segments function for that, but I can't address each panel separately. Links right for panel 1, show up in all other

In the hypothetical example below, how do I add two segments() into the two panels, respectively? Say segments(x0=5, y0=10, x1=5, y1=20) on the left and segments(x0=15, y0=-10, x1=15, y1=-2) on the right? Many thanks in advance, Yuelin Li. ps. part of the code came from a solution given by Deepayan Sarkar. ------------------- library(lattice) set.seed(12345) x <- 0:20 y.male.obs <- - 1.2

Full_Name: Gang Liang Version: 1.4.1 OS: mandrake-8.2 Submission from: (NULL) (128.32.81.135) tseries version: 0.9-1 quadprog version: 1.4-4 mva version: 1.4.1 version: platform i686-pc-linux-gnu arch i686 os linux-gnu system i686, linux-gnu status major 1 minor 4.1 year 2002 month 01

I have a dataframe that looks like this: > str(chr) 'data.frame': 84 obs. of 7 variables: $ county: Factor w/ 3 levels "Broome","Nassau",..: 3 3 3 3 3 3 3 3 3 3 ... $ item : Factor w/ 28 levels "Access to healthy foods",..: 21 19 20 18 16 3 2 6 17 8 ... $ value : num 8644 15 3.5 3.9 7.7 ... $ low : num 7897 9 2.5 2.6 7 ... $ high : num 9390