search for: gardar

...ints out a = 0, not a = 10 [1] ".GlobalEnv" "input" "package:methods" "package:stats" [5] "package:graphics" "package:utils" "Autoloads" "package:base" [1] 0 > Thanks, Gardar ___________________________________________________________________ Gardar Johannesson Lawrence Livermore National Laboratory 7000 East Avenue, L-229 Livermore, CA 94550 johannesson1@llnl.gov Tel: (925) 422-3901, Fax: (925) 422-4141 [[alternative HTML version deleted]]

...x.R' ... OK running code in 'ts-Ex.R' ... OK running code in 'survival-Ex.R' ...bash: survival-Ex.R: No such file or directory mv: cannot access survival-Ex.Rout *** Error code 2 make: Fatal error: Command failed for target `survival-Ex.Rout' Current working directory /home/gardar/local/src/R-1.6.0/tests/Examples *** Error code 1 make: Fatal error: Command failed for target `test-Recommended-Examples' Current working directory /home/gardar/local/src/R-1.6.0/tests/Examples *** Error code 1 make: Fatal error: Command failed for target `test-Recommended-Examples' Curren...

...9; attribute, however, if I don't create the 'n.ch' attribute, it finds it! That is: > tmp <- list(id=1) > attr(tmp,'n') NULL > attr(tmp,'n') <- 1 > attr(tmp,'n') [1] 1 > -- _________________________________________________________ Gardar Johannesson Department of Statistics The Ohio State University 304E Cockins Hall, 1958 Neil Ave. Columbus, OH 43210 Email: gardar at stat.ohio-state.edu Tel: (614) 292-1567 Web: www.stat.ohio-state.edu/~gardar Fax: (614) 292-2096 _______________________________________________...

...ple: Start a new section of R, with no objects: > ls() character(0) > gc() used (Mb) gc trigger (Mb) Ncells 187633 5.1 407500 10.9 Vcells 36061 0.3 786432 6.0 > And the top commands shows: PID USERNAME THR PRI NICE SIZE RES STATE TIME CPU COMMAND 18068 gardar 1 32 0 15M 12M sleep 0:02 0.01% R.bin Creating few objects: > x <- numeric(1e7) > gc() used (Mb) gc trigger (Mb) Ncells 187642 5.1 407500 10.9 Vcells 10036062 76.6 10394092 79.4 > a <- x > gc() used (Mb) gc trigger (Mb) Ncells 1...

...tsSmooth text html latex tsdiag text html latex example R_LIBS= ../../../bin/R CMD INSTALL ERROR: no packages specified *** Error code 1 make: Fatal error: Command failed for target `_tgz.ts' Current working directory /home/gardar/local/src/R-1.6.0/src/library/Recommended *** Error code 1 (ignored) /home/gardar/local/src/R-1.6.0 spatial> That is, the recommened packages are not installed. I traced this down to the solaris version of make (gmake works fine). I think the solaris version has problem with the &q...

...time R-1.4.0-opt: 14.97 sec user CPU time S-6.0 : 13.70 sec user CPU time 4) The result of system.time(B <- A %*% A) where A is 500x500 matrix. R-1.4.0-def: 18.83 sec user CPU time R-1.4.0-opt: 0.37 sec user CPU time I hope this will be of use to somebody... cheers, Gardar _________________________________________________________ Gardar Johannesson Department of Statistics Ohio State University 304E Cockins Hall, 1958 Neil Av. Columbus, OH 43210 Tel: 614-292-1567 Fax: 614-292-2096 e-mail: gardar at stat.ohio-state.edu WWW: www.stat.ohio-state.edu ________________...

...down.stop - Stop a shutdown in progress test.battery.start - Start a battery test test.battery.start.deep - Start a deep battery test test.battery.start.quick - Start a quick battery test test.battery.stop - Stop the battery test An example: upscmd sinus test.battery.start.quick Username (root): gardar Password: Unexpected response from upsd: ERR ACCESS-DENIED I have defined the user in the /etc/ups/upsd.users 74 [gardar] 75 password = hidden 76 allowfrom = local 77 actions = SET 78 instcmds = ALL 79 [local_mon] 80 password = hidden 81 al...

...+ 5 ## the time this loop should end at [block of R commands] while(proc.time()[3] < t.end) {} } This guaranty, at least, 5 second evaluation time for each loop in the above case. Can I delay the loop more efficiently, that is, not using as much CPU time as in the above case? Thanks, Gardar I'm using version 1.1.0 on: Machine hardware: sun4u OS version: 5.7 Processor type: sparc Hardware: SUNW,Ultra-250 _________________________________________________________ Gardar Johannesson Home Addr. Department of Statistics 605 Dennison, #3 Ohio...

Thank you Brian Ripley, Gardar Johannesson, and Marcel Wolbers for your prompt and friendly help! I will share any further learnings as I move through these suggestions. -Bob Abugov Brian Ripley wrote: I just use optim() on the sum of squares with non-negativity constraints. That did not exist in 1999. Gardar Johannesson...

...at package has an equivalent in C (that does handle events during the pause). We do plan to have a version of this for all platforms in due course. ################## By the way; is it 'required' to write a summary of responses and post to r-help? (e.g., as on the s-news list.) Thanks, Gardar _________________________________________________________ Gardar Johannesson Home Addr. Department of Statistics 605 Dennison, #3 Ohio State University Columbus, OH 43215 305A Cockins Hall, 1958 Neil Av. Tel: 614-531-0102 Columbus, OH 43210...

R 1.6.0 has gone into final feature freeze as of today. As a new feature, we'll make interim beta versions available via ftp://cran.us.r-project.org/pub/R/src/base alias http://cran.us.r-project.org/src/base (filename R-1.6.0beta_*.tar.gz, where * is the creation date). If you want to help ensure that the final 1.6.0 works satisfactorily on *your* platform, it might be a good idea to

...ve and have set the problem up with solve.QP and just the non-negativity constraints along the lines of: y as the data vector X as the design matrix D <- t(X) %*% X d <- t(t(y) %*% X) A <- diag(ncol(X)) b <- rep(0,ncol(X)) fit <- solve.QP(D=D,d=d,A=t(A),b=b,meq=0) (as per Gardar Johannesson '01) When I try to add the extra constraint that sum(weights)=1 I get errors owing to incompatibility of matrices. I add the constraint by putting an extra column of all ones to A and setting meq=1. I can work round it I think, by using an intercept and using the extra column...

...know in this case one does only need to carry out the symbolic factorization _once_ and then follow that up with a numerical factorization for each of the Q_i's (re-using the general symbolic factorization each time). Does anybody know if this is possible using the Matrix package? Thanks, Gardar Johannesson Lawrence Livermore National Laboratory

Dear list, I have problems installing R-2.0.1 on SUSE Linux 9.2. I used the following commands in order install R in /usr/local ./configure --with x11 --with-readline make make install When starting R and trying to display a plot on the X-window system output is written to a postscript file. When I try to run x11 with >x11() Error in X11(): X11 is not available. I do have a running

Hello, I am doing mcmc=10000 simulations from a posterior distribution of the parameters of a mixture of K=6 normal densities. I have mcmc by K matrices simMeans, simVars and simWeights containing the simulation output: one row for each simulation, one column for each normal component of the mixture. One thing I would like to do is a plot of the posterior predictive density. In order to do that

Just wondering what information is required for adding a ups to the compatability list... Is upsc output enough? I've got a Inform Sinus SS 210 ( http://www.informups.com/sinus_1_phase_technics.html ) that seems to be working allright with the blazer_ser driver Here is the upsc output: battery.charge: 100 battery.runtime: 1379 battery.voltage: 41.40 battery.voltage.nominal: 36.0

Hi there, I am running Monte Carlo Simulations in R using ordinary "while (condition)" loops. Since the number of iterations is something like 100.000 and within each iteration a given subsample is extended sequentially it takes hours to run the simulation. Does anyone know if there is either a way to avoid using loops in Monte Carlo Simulations or how to include possible faster

Dear R-users I have four simple linear models, which are all in the form of a*X+b The estimated parameter vectors are a <- c(1,2,3,4) b <- c(4,3,2,1) My goal is to draw a plot where x-axis is X (range from -100 to 100) and y-axis is the sum of all four linear models X <- seq(-100,100, length=10000) plot(X, sum of the four linear functions) I started with a function for summing