similar to: R-beta: NaN & Inf

Thanks to Ted and Guido for your quick help, I think I know what happened. ernesto ______________________________________________________ Get Your Private, Free Email at http://www.hotmail.com -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html Send "info", "help", or

Hi I'm writing a function that uses four parameters (scalars) and I need to run it in an iterative process (the parameters vary to find the minimum RSS). I don't want to use loops and so tried the do.call function. However it didn't work. My understanding is that the do.call simple runs the function replacing the arguments (scalars by vectors), instead of runing the function for

Hi I'm working on a remote computer with R (SuSE 8.1) and I want to do some ploting. Does anyone knows if it is possible to export the graphics window in to my display and how I can do it ? Thanks EJ -- Ernesto Jardim <ernesto at ipimar.pt> Marine Biologist Research Institute for Agriculture and Fisheries Lisboa, Portugal Tel: +351 213 027 000 Fax: +351 213 015 948

Hi I'm computing a bca interval using bca.ci from the boot package. When I try to use this I get an error > library(boot) > boot(logglm.data,boot.fishpower,2500,coef.vec=coeflm.vec)->blm8901 > bca.ci(blm8901,index=29) Error: cannot allocate vector of size 456729 Kb However my machine has 2GB of memory and without R running I only have 112M of memory used. Is there something

Hi I have a list with several data.frames, all with the same number of colunms but different number of rows, and I'd like to transform this list into a single dataframe. I need to mimic an rbind of all dataframes ... Transform doesn't seem to work :-( Thanks EJ -- Ernesto Jardim <ernesto at ipimar.pt> Bi?logo Marinho/Marine Biologist IPIMAR - Instituto Nacional de Investiga??o

Hi I'm doing a bootstrap of a linear model using: boot.fishpower <- function(data, i){ data <- data[i,] fplm <- lm(log(U)~Q+S+P+B+D, data=data) fp <- coef(fplm) exp(fp) } > boot(logglm.data,boot.fishpower,100) Error in "[<-"(*tmp*, r, , value = statistic(data, i[r, ], ...)) : number of items to replace is not a multiple of replacement length

Try factor(vec2) in your bwplot() call. -----Original Message----- From: r-help-bounces at stat.math.ethz.ch [mailto:r-help-bounces at stat.math.ethz.ch]On Behalf Of Ernesto Jardim Sent: Friday, July 09, 2004 9:41 AM To: Mailing List R Subject: [R] Problem with bwplot Hi, I'm ploting some box-and-whisker plots with bwplot but I'm not getting any box-and-whiskers ... just dots.

Full_Name: Ernesto Jardim Version: 1.2.2 OS: SuSE Linux 6.4 Submission from: (NULL) (193.137.98.92) I wanted to import an array, dim=c(28,12,17), from splus into R.On splus I did dump(c("array"),file="dump.out) I checked and Splus could source it. Than, in R I tryied to source it with: source("dump.out") and got and error message: > source("dump.out")

Hi everyone I'm trying to use sigma2.mat<-apply(ri.mat,2,'-n/2*log(sum(ri.mat^2)/n)') but it's not working. Is it wrong to use that kind of FUN in apply ? or am I doing something wrong ?! Thank's ernesto -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html Send

Hi I've just submitted to the NEdit development team a R highligth patterns for NEdit (www.nedit.org). NEdit is a text editor for LINUX which I use to write small R functions. In my opinion is a very good tool for small scripts, maybe for huge projects is not the best. I'm attaching the R-5.1.pats file so you can try it. As usual contributions and comments are welcome. Regards EJ

Hi I'm exporting some graphs from R to M$Word. I used png, jpeg and bmp and the quality is poor when comparing with the postscript. What is the best way to export a chart to be included in a M$Word file ? Thanks EJ

Hi, I want to use the first digit of the elements of a vector. I've tried grep but didn't work. Any help is welcome. Thanks EJ > grep("^[0-9]",as.character(runif(100,0,2))) [1] 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 [19] 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 [37] 37 38 39 40 41 42 43 44 45

Hi I'm doing a xyplot and I wand to reduce the number of tick marks in the x axis. My x axis are month and I want to reduce the 12 tick marks to 4. I used the scales argument but it doesn't seem to work, althougth it works on y axis if I use scales=list(tick.number=4).

Hi How can I do the same has interp in s+ with R ? I've looked to "approxfun" but didn't figured out how to do it. Thanks EJ -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- 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",

Hi I have a matrix (4 x 950) and I want to remove 3 rows, where the values from the first column are 713, 714 and 715. I can select the rows, one by one, with mat[mat$first==713,] mat[mat$first==714,] ... but I'm unable to (i) select the 3 rows at once, (ii) select the matrix excluding those rows. How can I do it ? Thanks EJ

Hi I'm getting the following error that do not make sense to me, what am Idoing wrong ? > acf(Recsim[1,], lag.max=1) Error in acf(Recsim[1, ], lag.max = 1) : 'lag.max' must be at least 1 Regards EJ

Hi I'm trying to make an histogram using the plot() function. Mainly I do so because I need to make the histogram of the values in the vector and not of the "counts" or "probabilities" of it. The main problem is that when using the flags type="h", lwd=10 I get a histogram with rounded bars ... and I didn't found out how I can change this. All help will be

Hi, I'm ploting a portuguese map with "plot" but I'm getting a strange vertical line. To avoid posting binary files, I've rolled the data file, an R script and the resulted ps file in this tgz file http://ernesto.freezope.org/ploterror.tgz. Can someone give me an help on this ? Thanks EJ

Hi! I 've a problem to have a lognorm distribution with mean=1 and var (or sigma)=1. rlnorm(1000,0,0) rlnorm(1000,1,1) rlnorm(1000,0,1) .... ? Can you help me?

Hi, I'm using the function krige.conv from package geoR but I'm getting a warning that I do not understand: Warning messages: 1: NaNs produced in: sqrt(variance) 2: NaNs produced in: sqrt(variance) I'm performimg an ordinary kriging with a log transform (lambda=0). Does anyone knows where this comes from ? Thanks and regards EJ