similar to: R glm

I have a monthly price index series x, the related return series y = diff(log(x)) and a POSIXlt date-time variable dp. I would like to apply annual blocks to compute for example annual block maxima and mean of y. When studying the POSIX classes, in the first stage of the learning curve, I computed the maximum drawdown of x: > mdd <- maxdrawdown(x) > max.dd <- mdd$maxdrawdown > from

I have R-0.62.3 and RedHat 5.0 with fort77. At the office R is running but when I install a package, eg. 'akima' I get a message make: f77: Command not found 'R INSTALL library' does not accept fort77, but R installation checks f77/g77/fort77. Is there an easy way to edit some file and replace f77/g77 with fort77? At home I have the same R and Linux versions, but when

Hi, I have a graph with seven parallel horizontal lines. Is it possible to shade the area between two adjacent lines? Thank you in advance, Hannu

Hi, I am trying to maximize a utility function using optim. I have a problem with optim, since some of the values in the caw, mom, me and btm matrixes in the code bellow are missing. Is there a handy way just to skip the missing values in the loop? g <- 5 retp <- NULL object <- function (x) { b1 <- x[1] b2 <- x[2] b3 <- x[3] for(i in 1:nrow(ret)){ for(j

Hi, when using a laptop and trying to install packages for R for Windows 2.2.1 I get the following error > utils:::menuInstallPkgs() --- Please select a CRAN mirror for use in this session --- Warning: unable to access index for repository http://cran.dk.r-project.org/bin/windows/contrib/2.2 Warning: unable to access index for repository

Hello All, I would like to optimize a (log-)likelihood function subject to a number of linear constraints between parameters. These constraints are equality constraints of the form A%*%theta=c, ie (1,1) %*% 0.8,0.2)^t = 1 meaning that these parameters should sum to one. Moreover, there are bounds on the individual parameters, in most cases that I am considering parameters are bound between zero

Hi, I have adviced my students to install R and an editor on a USB drive for working in the computer class. With R everything works fine following these instructions: http://personal.bgsu.edu/~mrizzo/Rmisc/usbR.htm. But several editors (e.g., Tinn-R and WinEdt) require administrator rights. I have found out that Emacs (Vincent Goulet's Emacs for Windows Modified

My apologies in advance for my thickness but I can't seem to solve the following, seemingly simple, data manipulation problem: I have a data frame that contains multiple factors and multiple continuous response variables, but duplicates of some factor combinations. The duplicates contain bad data, so I would like to eliminate the duplicates. I would like to retain the entire rows

R v. 1.7.1, Windows 2000. A particular journal wants me to provide scatter plots with no box, but with axes that meet in the lower left corner. It seems as though there must be an easy way of doing this, but my reading the help on plot.default, axis, and box have not provided any clues. I would be most appreciative of any feedback. Thank you, Hank Stevens Dr. Martin Henry H. Stevens,

Hi All, I am running R 2.2.0 on Mac OS 10.4.2, dual G5 processors with 8 Gig RAM. I am running a simulation with lsoda that requires ~378 s to complete one set of time intervals. I need to optimize the parameters, and so need to considerably speed up the simulation. I have tried to figure out how to change the appropriate memory allocation and have search R help and Introductory

Hi - I am getting different results when I run the numerical integrator function lsoda (odesolve package) on a Mac and a PC. I am trying to simulating a system of 10 ODE's with two exogenous pulsed inputs to the system, and have had reasonably good success with many model parameter sets. Under some parameter sets, however, the simulations fail on the Mac (see error message below). The

Hi Folks, How do I get red bold font on my y axis and black standard font on my x axis? plot(runif(10), ylab="Red, Bold?", xlab="Black, standard?") Any pointers or examples would be great. Thanks! Hank Dr. Hank Stevens, Assistant Professor 338 Pearson Hall Botany Department Miami University Oxford, OH 45056 Office: (513) 529-4206 Lab: (513) 529-4262 FAX: (513) 529-4243

How do I cite a package (not R itself - I know how to do that)? Any thoughts or links? Many thanks in advance? Hank Stevens Dr. Martin Henry H. Stevens, Assistant Professor 338 Pearson Hall Botany Department Miami University Oxford, OH 45056 Office: (513) 529-4206 Lab: (513) 529-4262 FAX: (513) 529-4243 http://www.cas.muohio.edu/botany/bot/henry.html http://www.muohio.edu/ecology/

Hi, I just found out that Tinn-R_2.3.6.0 is now available on http://sourceforge.net/projects/tinn-r/. When configuring Tinn-R: R > Configure > Permanent (Rprofile.site) I get error C:\Program Files\R\R-2.12.0\bin\etc\RProfile.site The above file was not found. Please try to repeat the procedure! The correct path for RProfile.site is C:\Program Files\R\R-2.12.0\etc\Rprofile.site I was

All, In the function below I have 24 individuals and 6 calculations per individual. The 6 calculations are collected each time in a 1:24 loop when calculating "delta". I'd like to collect all 144 = 24*6 calculations in one vector ("delta.patient.comb"). The function works as is via indexing, but is there an easier way to collect the measurements via appendinng the 6

hi i'm looking for a gibbs sampling algorithm for R for the case of mixture of K normals, and in particular for the case of bivariate normals. i'd be grateful if anyone could send its own R-routine, at least for the univariate case. thank you in advance matteo

I recently formatted and reinstalled my laptop. Since then html-help function is not working properly. I can open the page, however, searching for any term or opening any link does not work. The Internet explorer just mentions an error on the page if I try to open a links and instead of searching it reloads the page. Reinstalling R did not help and now I am not sure if it is the IE, or R, or

Using R v. 1.7.0 on Windows 2000 I would like to plot the fitted values of a model as a function of a continuous covariate, augmented with data (e.g., augPred) grouping by combinations of fixed effects. I have not been able to use augPred effectively, and am wondering if it does not handle unbalanced data (3 out of 192 missing). I include below the model and an xyplot that almost does the

Respected Sir, I am trying to import excel file into R, but I need to truncate some columns from the original file. How to delete unwanted columns when I import data from excel file. How to use cor.test for the data when I want the output rowwise. How to do grouping and use cor.test on that data I need some help regarding how to calculate the correlation. I don't know whether you

I have data on number of vines per tree for ~550 trees. Over half of the trees did not have any vines and the data is fairly skewed (median = 0, mean = 1.158, 3rd qu. = 1.000). I am attempting to investigate whether plot location (four sites), species (I'm using only the four most common species), or tree dbh has a significant influence on the number of vines per tree. When I