similar to: Stochastic dynamic programming

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, I''m trying to fit a linear mixed effects model of the form: lme(y ~ x * Sex * Year, random=x|subject) where Sex and Year are factors with two and three levels respectively. I want to compare the fixed effects for each level to the overall mean, but the default in R is to compare to the first level. This can be changed by adding the term -1 to the righthand side of the model

Hi, Sorry if this is an obvious one, but is there a simple way to modify the lm function to test whether a slope coefficient is significantly different from 1 instead of different from 0? Thanks, Martin Martin Biuw SEA MAMMAL RESEARCH UNIT Gatty Marine Laboratory School of Environmental and Evolutionary Biology University of St Andrews Fife, Scotland KY16 8LB UK phone +44 (0)1334 462630 fax

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

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

G'Day again, I am attempting to read a large MSAccess file into R, but get memory problems. With the first 100 rows of the table ("Macca99") things, as shown below, are fine and the resulting object is 33,780 bytes. But when I read the entire table ("MaccaDiv99") which is 218,000 rows R falls over with the message: Rgui.exe - Application error The instruction at

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

it seems that word can not read encapsupalted postscripts generated by R I used this command postscript("output.eps",horizontal=F,onefile=TRUE) since onefile=TRUE produces an encapsualted postscript actually what I'm trying to do is to insert the postsript file into a word document since other formats like jpeg and bmp do not reproduce the same quality like postscript formats any

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(

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

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

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

A non-text attachment was scrubbed... Name: not available Type: text/enriched Size: 905 bytes Desc: not available Url : https://stat.ethz.ch/pipermail/r-help/attachments/20010420/61ada248/attachment.bin

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

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

Hi All, When a couple of EXT4 filesystems are mounted in a server I get the message Oct 1 18:49:42 sraid3 kernel: EXT4-fs (sdb): mounted filesystem without journal Oct 1 18:49:42 sraid3 kernel: EXT4-fs (sdc): mounted filesystem without journal in the system logs. My confusion is why are they mounted without a journal? They were both created with mkfs -t ext4 /dev/sdb mkfs -t ext4