similar to: Status of "data.table" package

Hi All. (This is probably an incredibly stupid question :)) I have just noticed that every time I start an R2.3 session and request a random number (say with rnorm) that I receive the same number. AFAIK I have not set a seed at any point. Returning to an older version of R(v1.8) I get the behaviour that I have come to expect, where a new seed is automatically generated for each new session

Dear R-Help List, I am using optim, with method=L-BFGS-B, to maximize a likelihood inside a large simulation exercise. This runs fine for most simulated data sets, but for some reason, about 1 out of 100 times, optim will just hang. Using a dumb approach to the problem (i.e. printing the parameter values each time the function being maximized is evaluated), I tracked down when this happens,

R-Devel, I note from the CHANGES log accompanying the Windows version of R-2.2.0 that the behaviour with respect to the Windows registry has changed. It says: "If the user chooses to register R during installation, a registry entry HKEY_LOCAL_MACHINE\Software\R-core\R\{version}\InstallPath will be added. Users require administrative privileges to create this key. For others, the same key

I was having some problems persuading optim() to give me the answers I wanted, & simplified down to: sqr<-function(x){(x+1)^2} optim(1,sqr) I accept this is a hammer to crack a nut, but was still expecting the answer -1. I got: $par [1] -0.8 $value [1] 0.04 $counts function gradient 12 NA $convergence [1] 0 $message NULL so I've

Hi, I've got the following command that works fine; it prompt for a password then runs okay: "rsync -az -e ssh myserver.ac.uk:/var/lib/mysql /www/rsynctagsassets/databases" I wanted to run it as a script so put it into a file called rsync.sh and chmod +x 'd it. It runs as a script okay but still prompts for a password. I read "man ssh" and the bit about

Apologies if this is a naive beginners question. I am trying to create a dotplot with the lattice dotplot function in which my dots are colored differently depending on if positive or negative and sized by sp.nc.bdrs.data$mwZ I have tried... dotplot(sporder ~ cvarorder | direct, data=sp.nc.bdrs.data, cex=abs(sp.nc.bdrs.data$mwZ * 0.05), xlab="climate variables",

Hello list, I was hoping I could get some help on something which is really giving me a headache. I am using R version 2.14.1 (2011-12-22) (Platform: x86_64-pc-mingw32/x64 (64-bit)) An object which is supposed to have times has a few elements listed as not times but NA's > ls() [1] "todelete" > class(todelete) [1] "POSIXlt" "POSIXt" > todelete [1]

G'Day, I'm a little confused why the c function has the followng effect on classes - is this a feature ? My workround [ class(cc) <- c("POSIXt", "POSIXct") ] seems to do the job. Many thanks Bernie McConnell "R version 1.3.1, 2001-08-31" on NT > aa <- as.POSIXct("2001-09-23") > bb <- as.POSIXct("2001-09-24") > cc

Is there a way to add residuals to plots produced by plot.gam in the mgcv package? I'm looking for something like what you get using resid=T in Splus plot.gam. Thanks in advance Toby -----Original Message----- From: Simon Wood [mailto:snw at mcs.st-and.ac.uk] Sent: 23 January, 2002 8:14 PM To: Toby.Patterson at csiro.au Cc: r-help at stat.math.ethz.ch Subject: Re: [R] multiple surfaces in

Hi All. I have R code whose functionality is being replicated within a C+ program. The outputs are to be compared to validate the conversion somewhat - however (as is always the case) I have stuffed my code with random number calls. Random uniform numbers in C+ are being produced using the (Boost) mersenne-twister generators (mt11213b & mt19937) - which is the default type of generator

Hello All, I am looking for some help with was is probably very simple: I want the y-axis title '$B&D(B13C$B"s(B VPDB' and similarly for nitrogen on a secondary axis. I am able to get the delta notation and per mil notation separately but one messes up the other with the inclusion of quotations: ylab=expression("{delta}^13*C \211 VPDB") gives the per mil just fine but

Hi All, Thought I would let those that are interested know that I had success in running 24G on an Asus P6T with 24G kit of Kingston DDR3. While I was putting this together I saw lots of forum posts asking if anyone had tried it. Well we did here at our work and all looks great including running "memtest86" overnight. I have a fluid dynamics simulation running on it with 90% memory

Is Samba (1.19) incompatible with FoxPro? We have tried to network a database package (FoxPro is the underlying DBMS) by putting it in a Samba share. The database starts up and you can browse the tables, but it looks as though FoxPro's record locking system causes the entire share to be locked, so that further access to the shared drive (both in FoxPro and in the Windows98 Explorer) is

Hi, Is it possible to compare two GAM objects created with the gam() function from the mgcv package. I use a slightly modified version of anova.glm() named anova.gam(), modified from John Fox (2002). It often gives me some aberant responses, especially with "F" test. I use a quasibinomial model and scale (dispersion) is calculated and used in the calculation of the F value. Does someone

Dear all, I'm trying to use the Pspline add-on package to fit a quintic spline (norder =3), but I keep running into a Singularity error. > traj.spl <- smooth.Pspline(time, x, norder=3 ) Error in smooth.Pspline(time, x, norder = 3) : Singularity error in solving equations > Playing around with the other parameters produces an "unused arguments" error: > traj.spl

Hello, Once again, I posted a message without a subject line. Sorry.... here is the question again. Is there a simple way to modify the circ.mean function in the CircStats package to include a vector of weights to obtain a weighted average angle? Thanks! Martin -- Martin Biuw Sea Mammal Research Unit Gatty Marine Laboratory, University of St Andrews St Andrews, Fife KY16 8PA Scotland Ph:

Hello all, Does anyone know of any already written functions for carrying out stochastic dynamic programming in R? Thanks! Martin -- Martin Biuw Sea Mammal Research Unit Gatty Marine Laboratory, University of St Andrews St Andrews, Fife KY16 8PA Scotland Ph: +44-(0)1334-462637 Fax: +44-(0)1334-462632 Web: http://smub.st.and.ac.uk

Hi Guys, We are having the weirdest problem with all SAMBA version 2.05 to 2.07pre 3 SAMBA is running on DEC UNIX 4.0F For no reason that we can see some users are no longer being allowed to connect to their shares. We see the problem with WIN95, 98 2000 (all plain passwd) Let me explain more clearly If I map netowrk drive //joequimby/naismith my own user and passwd are not recognised

Hello, This is probably very simple but I am trying to select a region from the OIv2 SST dataset that crosses 0 longitude. The longitude variable is arranged from 0 to 360 and I get an error message when trying to go from 330 to 30 for example. The data are from a NOAA netCDF file. *temp = open.ncdf("c:/documents and settings/theoni/desktop/sstoi_v2.nc") print(temp) y1 = get.var.ncdf(

Hello, I'm using the nls function and would like to increase the number of iterations. According to the documentation as well as other postings on R-help, I've tried to do this using the "control" argument: nls(y ~ SSfpl(x, A, B, xmid, scal), data=my.data, control=nls.control(maxiter=200)) but no matter how much I increase "maxiter", I get the following error