search for: cullens

...aph of this over the range [-3,3] so I used: plotit <-function(x=seq(-3,3,0.01),...){ y<-sapply(x,mdevstdev) plot(x,y,...) } > plotit() This produces the graph, but some discontinuities appear on it. I've produced the same graph in Mathematica 5 (http://econserv2.bess.tcd.ie/cullens/R/DOverDelta.eps), and it was smooth over this range (it takes ages to run, though). Is this a numerical precision problem? Any suggestions on how to improve the precision? I'm running R 1.9 on WinXP, PIII. I haven't changed any R parameters that I know of. -- SC Simon Cullen Room...

...t;- table(cut2(z.trun, cuts=breaks)) Having done this I've found that the code takes much longer to run - up to 10x as long. Is there a more effecient way of doing this? Anyone have any thoughts? -- SC Simon Cullen Room 3030 Dept. Of Economics Trinity College Dublin Ph. (608)3477 Email cullens at tcd.ie

Hi all, I've been trying to run a model using R2winBUGS, and recurrently I get the message: "Error in FUN(X[[3L]], ...) : invalid to change the storage mode of a factor" My model is the following: sink("GLMM_Poisson.txt") cat(" model{ mu~dnorm(0,0.01) beta1~dnorm(-1,1) for(j in 1:nsite){ alpha[j]~dnorm(mu.alpha,tau.alpha) } mu.alpha~dnorm(0,0.01)

Hi all I want to conditionally operate on certain elements of a matrix, let me explain it with a simple vector example > z<- c(1, 2, 3) > zz <- c(0,0,0) > null <- (z > 2) & ( zz <- z) > zz [1] 1 2 3 why zz is not (0, 0, 3) ????? the null <- assignment is to keep the console silent in the other hand, it curious that null has reasonable values > null [1]

Hi, I installed the "kwatch-mysql-ruby" gem on my Debian server but when I try to run the "db:migrate" task I get the following error: ** Invoke db:migrate (first_time) ** Invoke environment (first_time) ** Execute environment ** Execute db:migrate rake aborted! uninitialized constant MysqlCompat::MysqlRes

Dear all, I need to plot an cumulative distribution plot of a variable and then to fit a distribution to that, probably a weibull or lognormal. I have plotted the ecdf as > plot(ecdf(x)) but I haven't managed to fit the distribution. I have as well attached the data. I would appreciate if you could help me on that. Thank you. Kind regards Maria -------------- next part --------------

Hi, Well, i need some help, practical and theoretical. I am wondering why the fitdistplus (mle function) is returning an error for this code: [code] x1 <- c(100,200,140,98,97,56,42,10,2,2,1,4,3,2,12,3,1,1,1,1,0,0); plotdist(x1); descdist(x1, boot =1000); y<- sum(x1); d= as.vector(length(x1)); for(i in 1:length(x1)){ d[i] = x1[i]/y; } fitdist(d, "beta") [/code] Error:

Hi everyone, Here's a follow-up on the two BoFs we've had about doing royalty-free codecs at the IETF. Well, the good new is that the proposal is now in IETF last call until January 20th (see below). No final decision has been made, so it's important to get as much support as possible for the Working Group proposal. You can see the ongoing discussion on the mailing list archive:

Hi everyone, Here's a follow-up on the two BoFs we've had about doing royalty-free codecs at the IETF. Well, the good new is that the proposal is now in IETF last call until January 20th (see below). No final decision has been made, so it's important to get as much support as possible for the Working Group proposal. You can see the ongoing discussion on the mailing list archive:

Hi, I'm trying to label the horizontal axis of a plot with a symbol that is the equivalent of \mathbb{R} in LaTeX. I've had a look through the help pages for plotmath and for Hershey and haven't found the symbol. Could someone give me a pointer, please? Using R 1.9.1 on Win32. -- SC

Hello List, I have a data frame object (wa2) that I am exporting for use in another statistics package. Using library(foreign) write.foreign(wa2, choose.files(), choose.files(), package='SPSS') I noticed that there were several differences between the data sets as seen within R (View(wa2)) and what was produced in SPSS. Examining the data file produced by write.foreign (before running

Hi, I have a strange problem viewing data in real time. I wonder can anyone help? I run tail -f on a file which is served by a windows 2000 machine. I start a simple program running on windows 2000 which writes numbered lines of text to this file. I notice that tail -f is missing some of the lines. If I edit the file afterwards the data is there it's just that tail seems to miss bits of it

Hi, I have just switched from FreeBSD 4.11 using 3.0.21b to NetBSD 3.0.1 using 3.0.22. With the same configuration, and the same clients, I am now seeing what appears to be some rather large memory leakage [ I am not sure if this is the fault of Samba or NetBSD ]. Eventually the process will run out of memory (machine only has 256MB of RAM) and the Windows clients (2000, SP4 and all the

Hi, I'm trying to get C code working with R. This is my first time writing C on Windows and I'm making a mess of it. Help! I'm following the example in Roger Peng's "An Introduction to the .C interface to R". The C code is: #include <R.h> void hello(int *n){ int i; for(i=0; i < *n; i++) { Rprintf("Hello, world!\n"); } } I seem to be unable

Hello, I'm trying to figure out how to keep a mirror and 7 day incremental backup between to 2 mount points. I want to rsync everything from /mnt/production/ to /mnt/backup/production/ on the same server. Nothing fancy. But , I Want to utilize the --delete option for files that no longer exist on /mnt/production/ that have been on /mnt/backup/production/ for longer then 7 days. I

We asked Henning Schulzrinne (the author of RTP and currently the RTP/AVG maintainer) about getting a static RTP payload type code assigned for Speex, since we meet the criteria in rfc1890. He said no. Here is his answer: On Fri, 22 Nov 2002, Henning Schulzrinne wrote: > Sorry, this is not going to happen, regardless of the codec or its > merit. See the new draft which will replace RFC

Tiger's rsync anomalies I'm trying to sort out some anomalies using rsync locally. Booted to Volumes/LaCie_Disk_A I create a full direct clone with: sudo asr -erase -source /Volumes/Mirrored_HD_Set -target /Volumes/LaCie_Disk_B then boot to it and check it out. Then booted back to Mirrored_HD_Set I do a little other work (i.e. some files are changed). Then for a differential

Hi all, It's been a while since we've given an update on scalable vector support in LLVM. Over the last 12 months a lot of work has been done to make LLVM cope with scalable vectors. This effort is now starting to bear fruit with LLVM gaining more capabilities, including an intrinsics interface for AArch64 SVE/SVE2, LLVM IR Codegen for scalable vectors, and several loop-vectorization

...060807 at stat.auckland.ac.nz> > Content-Type: text/plain; charset=us-ascii; format=flowed > > Hi > > > Deepayan Sarkar wrote: > > On Thursday 26 August 2004 11:41, Martin Maechler wrote: > > > >>>>>>>"Simon" == Simon Cullen <cullens at tcd.ie> > >>>>>>> on Thu, 26 Aug 2004 15:30:06 +0100 writes: > >>>>>> > >> Simon> On Thu, 26 Aug 2004 15:24:33 +0200, Trenkler, > >> Simon> Dietrich <dtrenkler at nts6.oec.uni-osnabrueck.de> > >> >...

Our cluster is has recently been logging messages like this: Apr 18 08:05:39 lv2 kernel: (26235,1):o2hb_do_disk_heartbeat:776 ERROR: Device "blk_shared": another node is heartbeating in our slot! where blk_shared is the name of the shared disk that we have ocfs2 mounted on. Can someone tell me what this means?