Re: [R] uniform integer RNG 0 to t inclusive

that I need. Thanks again to both Prof Ripley and Duncan Murdoch, Sean O'Riordain affiliation - NULL On 19/09/06, Prof Brian Ripley [EMAIL PROTECTED] wrote: On Tue, 19 Sep 2006, Sean O'Riordain wrote: Hi Duncan, Thanks for that. In the light of what you've suggested, I'm now using

Re: [R] uniform integer RNG 0 to t inclusive

-.Machine$double.eps)) } } cheers and Thanks, Sean On 18/09/06, Duncan Murdoch [EMAIL PROTECTED] wrote: On 9/18/2006 3:37 AM, Sean O'Riordain wrote: Good morning, I'm trying to concisely generate a single integer from 0 to n inclusive, where n might be of the order of hundreds

[R] uniform integer RNG 0 to t inclusive

rev38247 language R version.string Version 2.3.1 (2006-06-01) Many thanks in advance for your help. Sean O'Riordain affiliation - NULL __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read

Re: [R] subsetting a data set

Hi Graham, Try creating a new column with the two levels that you want... something along the lines of (warning untested!!!) GQ1[(GQ1$Status == Expert) | (GQ1$Status == Ecol),]$newColumn - AllEcol GQ1[GQ1$Status == Stake,]$newColumn - Stake and then do the by(GQ1[,Max], list(GQ1$NewColumn),

Re: [R] how to the p-values or t-values from the lm's results

Hi there Zhang, While there might be a better way... an ugly but generic way of accessing this type of information is to use str() and a little experimentation... here is a little history() of what I did to find it... a str(a) str(logr) a[[1]] a[[2]] a[[3]] a[[4]] a[[4]][[1]] a[[4]][1,]

Re: [R] String manipulation and formatting

Does it have to be a stop char ., or could it be a separate parameter, i.e. put in a comma, then the 10 becomes an integer... s/ On 18/07/06, Prof Brian Ripley [EMAIL PROTECTED] wrote: On Tue, 18 Jul 2006, Bashir Saghir (Aztek Global) wrote: Thanks to Richard, Gabor and Marc for some nice

Re: [R] References verifying accuracy of R for basic statisticalcalculations and tests

Please don't shoot! q: would it be a good idea to use these datasets as a basis for some regression tests? Sean On 14/07/06, Rau, Roland [EMAIL PROTECTED] wrote: Hi, [mailto:[EMAIL PROTECTED] On Behalf Of Corey Powell Do you know of any references that verify the accuracy of R for

Re: [R] Generating random normal distribution with mean 0 and standard deviation 1

?rnorm On 14/07/06, Neuro LeSuperHéros [EMAIL PROTECTED] wrote: Hello, This must be really simple, but I can't find it on R Site search. I need to generate a random normally distributed series with mean 0 and sd 1. In Matlab, this code is randn(n). The closest I found is runif(20,-1,1)

Re: [R] PowerPoint

or try win.metafile() if i'm in a hurry, (in windows) i just right click on the graph and select Copy as Metafile and paste directly into powerpoint... Sean On 23/06/06, Gabor Grothendieck [EMAIL PROTECTED] wrote: Note that jpg, bmp and png are in less desirable bit mapped formats whereas

Re: [R] Qurey : How to add trendline( st. line) in Graph

lines(predict(lm(g~f))~f,lty=2) or if its a curvy line you're looking for... lines(predict(loess(g~f))~f,lty=3) On 19/06/06, priti desai [EMAIL PROTECTED] wrote: How to add trendline (i.e. straight line passing through maximum points) in graph. I have worked on the data given below. Please

Re: [R] can I call user-created functions without source() ?

Indeed, some folk say that the documentation should be written before the code... Sean On 19/06/06, Duncan Murdoch [EMAIL PROTECTED] wrote: 2. If you use the package system, you will be encouraged to create man pages for your functions. This is work, but I think it pays off in the end. I

Re: [R] how to successfully remove missing values for a repeated measures analysis

Spencer, A number of times I've been wrecking my head trying to solve a problem, and when I thought I'd tried everything possible (all the docs that I can think of, web searches etc...) I start re-reading the posting guide (I've no idea how many times I've read it now)... and 9 times out of 10,

Re: [R] a question about a simply figure

Hi Zhang Jian, If I say plot(t$no,t$leiji), then the lower bound is neatly set at about 25... (I'm not sure how I can measure the bounds on the current plot - but I'm sure it can be found!) You can set the bounds on the y-axis to be between 0 and 100 by saying

Re: [R] plot xy data with error bars

Hi Petr, I just needed to add an 'attach/detach' to get it going. A-data.frame(x=1:10, y=1:10, erx=rnorm(10)^2, ery=rnorm(10)^2) attach(A) plot(x,y) arrows(x,y,x+erx/10,y, angle=90, length=.1) arrows(x,y,x-erx/10,y, angle=90, length=.1) arrows(x,y,x,y-ery/10, angle=90, length=.1)

[R] win2k memory problem with merge()'ing repeatedly (long email)

requirements for this merge just explode... cheers, (and thanks in advance!) Sean O'Riordain == version _ platform i386-pc-mingw32 arch i386 os mingw32 system i386, mingw32 status Patched major 2 minor

Re: [R] win2k memory problem with merge()'ing repeatedly (long email)

, 32, 64... rows for that date). -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf Of Sean O'Riordain Sent: Monday, May 22, 2006 10:12 AM To: r-help Subject: [R] win2k memory problem with merge()'ing repeatedly (long email) Good afternoon, I

Re: [R] Unreadable labels

Hi Maciej, fyi - on winxp-pro 2.3.0-patched and up to date MASS - I can clearly read the labels, except the start of the word 'black' is clipped, but the font size is fine. cheers, Sean On 21/05/06, Maciej Bliziński [EMAIL PROTECTED] wrote: Playing around with examples from MASS4, I found a

[R] merge problem... extra lines appear in the presence of NAs

Patched (2006-05-11 r38037) win-xp-pro sp2 - binary installs from CRAN it works in a similar way if I say atot - merge(a1,a2,by.x=mdate,by.y=mdate,all=T) or even atot - merge(a1,a2,by=mdate,all=T) also tested on versions 2.2.1, 2.3.0 cheers, Sean O'Riordain (ps. ctrl-v paste wouldn't work

Re: [R] (Nothing to do with) merge problem... extra lines appear in the presence of NAs

on 'mdate'. You example is not reproducible, of course, since you used random values. Perhaps you intended a1[floor(runif(nacount)*count), value] - NA On Sat, 20 May 2006, Sean O'Riordain wrote: Good morning! [Or afternoon in Europe, ] I've searched the docs etc... Am I doing

Re: [R] Precision in estimating log

Anthony, in the same way that we're not allowed to say if(x==0) if x is a real number, we can't say that 0.05=1-0.95... as 1-0.95 is not represented as a base 10 number on the computer, but in some base 2^i depending on your computer...and the representation is not necessarily exact... i.e.

Re: [R] uniform and clumped point plots

Perhaps you're looking for something along the lines of Sobol sequences - refer Section 7.7 of Numerical Recipes in C by Press et al. Sean On 18/05/06, Ben Bolker [EMAIL PROTECTED] wrote: Beutel, Terry S Terry.Beutel at dpi.qld.gov.au writes: I am trying to generate two dimensional random

Re: [R] Please help me to combine two datasets.

merge(data.set.1, data.set.2) but how do I know which rows to drop in ds1? if data.set.2 had dates, you could say merge(data.set.1, data.set.2, by.x=Date, by.y=Date, all=T) or do you just want to drop the first length(data.set.1)-length(data.set.2) from data.set.1? do you have NA values?

[R] data input strategy - lots of csv files

each) like that which are timeseries with missing data... which doesn't go through the above awkwardness (and language manipulation) but still ends up with a nice data.frame with NA values correctly coded etc. Many thanks, Sean O'Riordain __ R-help

Re: [R] Please help me to combine two datasets.

Kumar Saha [EMAIL PROTECTED] wrote: Dear Sean, Thanks for this reply. But my problem is not solved. Actually I want to drop first two rows from dataset-1 and then combine them. Can you give me any idea? Thanks and regards, Arun On 5/11/06, Sean O'Riordain [EMAIL PROTECTED] wrote

Re: [R] Boxplot and range of x-axis

Hi Rainer, boxplot(liste1,ylim=c(1,6)) boxplot(liste2,ylim=c(1,6)) or would you prefer boxplot(liste1,liste2) cheers, S/ On 11/05/06, Rainer Hahnekamp [EMAIL PROTECTED] wrote: Hello there, could somebody help me with this: I have to create a lot of boxplots but with same range concerning

Re: [R] data input strategy - lots of csv files

On 5/11/06, Sean O'Riordain [EMAIL PROTECTED] wrote: Good morning, I have currently 63 .csv files most of which have lines which look like 01/06/05,23445 Though some files have two numbers beside each date. There are missing values, and currently the longest file has 318 rows

Re: [R] Newbie question about read.table

~age,binomial)) Error in eval(expr, envir, enclos) : object low not found The read.table function works fine and when look at the data it shows the variable names across the top. The data come from SPSS which I read out into a .dat file. Thanks - Original Message - From: Sean

Re: [R] R binaries, platform independent and Design of Experiments

On 02/09/05, Nam-Ky Nguyen [EMAIL PROTECTED] wrote: b) You do NOT want to do numerical computations on software available in Java byte code. You do not want to do heavy numerical computations with R either. Most statistical calculation using R requires a fraction of a second and I cannot

Re: [R] Reference manual is not available in the help menu of the rgui

Actually, I've started reading the reference manual... :-) I printed it out 2-to-a-page and I'm working my way through it, in order to learn about the full capabilities of the base system... I know I'm not going to remember everything, but when I bump into a particular problem, I'll know what

Re: [R] The Perils of PowerPoint

I can't lay my hands n it at the moment - its around here somewhere, but in Numerical Methods That Work by Forman Acton, the author points out that the result of computation should be insight, not numbers ps. an excellent book if you haven't seen it.

Re: [R] memory

Ferran, What are you trying to do with such a large matrix? with 7e9 cells and a linear algorithm which is quite unlikely, your problem solution is likely to take a very long time(tm)... just quickly... at one micro-second per operation (very optimistic?) and 7e9 operations, thats 7e9/1e6/60

Re: [R] [newbie] Want to perform some simple regressions.

Hi Thomas, I'm not an expert - so I might use incorrect terminology, but hopefully you'll get the picture! Assuming that you've got your data in a .CSV file, you'd first read in your data, where the first three lines might look like... x,y 0,2.205954909440447 1,8.150580118785099 # load the info

Re: [R] writing to a fixed format (fortran) file

?sprintf more C than fortran, but you get the idea :-) cheers! Sean On 27/08/05, Duncan Golicher [EMAIL PROTECTED] wrote: Could anyone help with what should be a simple task? I have data as a fixed format (fortran) table. I have no trouble getting it into R using read.table. Each column is

Re: [R] plotting GAM

Hi Emil, can you give us a working example of what you're trying to do? cheers! Sean ps. as per the posting-guide... :-) On 24/08/05, Tkadlec Emil [EMAIL PROTECTED] wrote: Dear colleagues, I would like to have GAM regression lines (package gam) thicker than the default setting does. Is

Re: [R] plotting GAM

Emil, I've no clue what gam() does nor how it works, but saying that is it possible to use lines(,lwd=???) in someway and plot on top of the original line? cheers, Sean On 24/08/05, Tkadlec Emil [EMAIL PROTECTED] wrote: Sear Sean, here is the more working example (I use the latest version of

Re: [R] read a table ignoring specific rows ?

Can you read in the entire file as a data.frame; and then construct a new data.frame which excludes some rows? perhaps something along the lines of... fred.file - data.frame() fred.file - edit(fred.file) fred.file colA colB colC 1142 2233 3354 443

Re: [R] A. Mani : Avoiding loops

Hi, I'm not sure what you actually want from your email (following the posting guide is a good way of helping you explain things to the rest of us in a way we understand - it might even answer your question! I'm only a beginner at R so no doubt one of our expert colleagues will help me... fred

Re: [R] Console

Hi Daniela, Which platform are you working on? If you're working within a console on windows-98, then the answer is entirely different to working under linux or RGui on windows. This is why the Posting Guide says to give platform details :-) cheers! Sean On 18/08/05, Daniela Salvini [EMAIL

Re: [R] How to repeat code snippet for several variables in a data frame?

Sander, consider writing a function to do your plotting, then pass in the dataframe name so... something along the lines of... # create a function which takes two arguments # mydf - a dataframe of a particular format... eg. the abc column is number 4 # i the column we want to aggregate and plot

Re: [R] How to get a list work in RData file

you know about ?history Sean On 16/08/05, Liaw, Andy [EMAIL PROTECTED] wrote: If you want to keep track of the function call that produced an object, usually you need to do that inside the function that's being called, e.g., test.xy - function(x,y) { + xy - x+y + attr(xy, Call) -

Re: [R] Problems runing R CMD check

hi! start - control panel - system - advanced - environmental variables probably better to be working in a directory with no spaces in the path... s/ On 13/08/05, Aleš Žiberna [EMAIL PROTECTED] wrote: Hello! I have a problem checking the package. Firstly, I do not know how to specify the

Re: [R] evaluating string variables

use paste() to construct the file name e.g. fn - paste(i,.mat,sep=) Sean On 12/08/05, Leonardo Sepulveda Durán [EMAIL PROTECTED] wrote: Hello!!! I have a folder (C:/R/) with matrix files, named by number, i.e. 0.mat, 1.mat,...,1250.mat. In this case, they are 5x5 simetric matrices. I

Re: [R] connexion problem getHdata (HMisc)

Hi Anne, does update.packages() work, or does it just hang? cheers! Sean On 11/08/05, [EMAIL PROTECTED] [EMAIL PROTECTED] wrote: Sorry about the missing info, I'm just trying to do too many things at the same time! I installed R on my Windows2000 computer (no choice there) . version

[R] solution: package install probs. win-xp. file move

Hi, This is just an FYI documenting a conflict between R and Google-Desktop. The solution was to click on the Google-Desktop icon in the systray and click on Pause Indexing. I also temporarily suspended anti-virus scanning before successfully install.packages(VR) many times without getting an

Re: [R] help regarding functions loops

Hi Rangesh, Perhaps I mis-understand your question, but it could be as simple as... p - 1-exp(-exp(s)) In R, this is vectorized such that a new vector is calculated - one value for each value of s, so p will have the same length as s. In the Introduction to R read up on vectors and how to

Re: [R] nnet package for OS X

Hi Lisa, as far as I know, nnet is part of the VR bundle which is a recommended package... refer http://www.stats.bris.ac.uk/R/src/contrib/Descriptions/VR.html so just try library(nnet) and see what happens! cheers! Sean On 02/08/05, Lisa Schweitzer [EMAIL PROTECTED] wrote: Hello every one--

Re: [R] R: graphics devices

Alan, I'm not sure what you mean... perhaps plot(y~x) # to the screen pdf(myplot.pdf) plot(y~x) # write the plot to the file dev.off() # close the file dev.off() # close the graphics window s/ On 29/07/05, Clark Allan [EMAIL PROTECTED] wrote: a simple question how does one produce plots on

Re: [R] R: graphics devices

like to do the following: say plot(y~x) on one screen and then plot(q~w) on another screen so that one can see both of them together but not on the same graphics device. Sean O'Riordain wrote: Alan, I'm not sure what you mean... perhaps plot(y~x) # to the screen pdf(myplot.pdf

Re: [R] How to save the object of linear regression in file and load it later

Hi Bahoo! I've found the R/Rpad Reference Card quite good at helping me find this sort of information... i.e. towards the bottom of the first column of the first page it says... save(file,...) saves the specified ojects (...) in the XDR platform independent binary format if I then say ?save it

Re: [R] The steps of building library in R 2.1.1

Just a few thoughts... Good documentation helps everybody - the beginners and the experts (less beginner questions if there is thorough and accessible documentation. I fully appreciate that this is a volunteer effort - I'm just trying to pin down some places where we have documentation issues.

Re: [R] convert to chron objects

are those dates in m/d/y or d/m/y ? ?chron and watch out for format = c(dates = d/m/y, times = h:m:s) On 13/07/05, Gabor Grothendieck [EMAIL PROTECTED] wrote: [I had some emails problems so I am sending this again. Sorry if you get it twice.] On 7/13/05, Gabor Grothendieck [EMAIL

Re: [R] R: to the power

(-8+0i)^(1/3) [1] 1+1.732051i ie complex... On 12/07/05, [EMAIL PROTECTED] [EMAIL PROTECTED] wrote: hi all why does R do this: (-8)^(1/3)=NaN the answer should be : -2 a silly question but i kept on getting errors in some of my code due to this problem. i solve the problem as

Re: [R] Problem compiling R 2.1.* on SUSE 9.2

Hi! what point exactly in the command sequence does it fail? during make? or make install? You might try saving the output and the error from the build into a file? Something along the lines of ./configure 21 filename or even ./configure 21 filename | tee configure_output.txt etc. cheers!

Re: [R] Tempfile error

Hi Harold, I know nothing about Sweave (though it certainly sounds like a great idea!). Does Sweave hold these files open simultaneously - many operating systems wouldn't be able to cope with 20,000 simultaneously open files. Could you be running out of disk space? or inodes - if you're on a

[R] finding out more about an object, e.g. lm

] 0.267 so now I can say my.r.sq - summary(lm(y~x))$r.squared but is there a better way? on my small lm model unlist(lm(y~x)) is a pretty long list! :-) Many thanks in advance, Sean O'Riordain __ R-help@stat.math.ethz.ch mailing list https

Re: [R] hellow

Hi, I suggest that you read the 'posting guide' - refer the bottom of all emails on the list. http://www.R-project.org/posting-guide.html (il ne faut que d'aller chez http://www.r-project.org/ vous meme). Sean karima alem wrote: can i ask you for logiciel R

Re: [R] access R remotely

This likely means that either R is not installed on that machine or the R binary is not in your PATH. : try asking your IT guy where R is installed... say Mr.IT.Guy says that R is installed in /usr/local/bin... then /usr/local/bin/R to start R... cheers, Sean. array chip wrote: Hi, not

Re: [R] R commands

Hi Andy, I know it doesn't answer all of your questions but have you read the Posting guide (see the bottom of all posts) ? http://www.R-project.org/posting-guide.html and in particular, try reading An Introduction to R, chapter 9? Note that in R, looping is relatively slow - if you can do

Re: [R] R programming

Hi Henrick, if you checked the posting guide (refer the bottom of all emails) you'd have spotted a link # R Language Definition (aka `R Language Manual') describes the R language, data objects, etc. http://cran.r-project.org/doc/manuals/R-lang.pdf cheers, Sean ps. R is indeed a full and

Re: [R] LCG with modulo 2^30

Hi Raheem, Firstly - fair warning...I'm not an R expert at all! However it is my understanding that the expression for i in 1:m creates a full vector in memory of all the consecutive numbers 1 to m... (i presume these are 4-byte ints here otherwise it would have fallen over before 2^29), but

Re: [R] R killed on Solaris

Hi Ott, could it be a per process limit or ulimits set? i don't have a solaris box here so i can't remember what the appropriate commands are... cheers, Sean Ott Toomet wrote: Hello, I am running a preliminary data-processing job on solaris. Basically, the program reads a large number of

Re: [R] R and RS232-Interface?

Hi, Personally, I would avoid Java for hardware handling - in my experience this is tricky and very hardware specific. Sascha, is this for a Unix type machine or under windows? If under Unix/Linux, it should not be particularly difficult to send characters and receive characters in C

Re: [R] Ignoring Errors in Simulations

?try Ken Kelley wrote: Hello all. I'm doing a simulation study and every so often I get an error that stops the simulation. I would like to ignore the errors *and* identify the particular iterations where they occurred. I have tried: options(error = expression(NULL)) which I thought would

Re: [R] Visualising Vectors

Hi Laura, you should find some useful information in the latest R news Volume 3/2, October 2003 http://cran.r-project.org/doc/Rnews/ refer page 8 where Paul's figure 2 shows some novel symbols showing China Sea Wind Speed, Direction and Temperature... plotted by lat.long. s. Laura Quinn

[R] r-ish ? how can i improve my code?

someone help me make this a bit more R-ish please? Or is there a completely different approach I should be taking? Many thanks in advance, Sean O'Riordain seanpor AT acm.org n - 900; # number of valid items required... x - numeric(n); y - numeric(n); z

Re: [R] r-ish ? how can i improve my code?

} Patrick Burns Burns Statistics [EMAIL PROTECTED] +44 (0)20 8525 0696 http://www.burns-stat.com (home of S Poetry and A Guide for the Unwilling S User) Sean O'Riordain wrote: Hi Folks, I'm trying to learn R. One of my intentions is to do some Monte-Carlo type modelling of road accidents. Below