[R] Unable to load R

2009-05-28 Thread anupam sinha
Dear all, I have recently installed R on a Red Hat Enterprise Linux (RHEL5) system. But I am unable to load R and it is giving the following error : */usr/local/lib/R/bin/exec/R: error while loading shared libraries: libreadline.so.5: cannot open shared object file: No such file

Re: [R] Unable to load R

2009-05-28 Thread Zeljko Vrba
On Thu, May 28, 2009 at 11:42:41AM +0530, anupam sinha wrote: I have checked for the presence of the above mention library and found that the library is present. I have run out of ideas. Can anyone help me out??? I will be greatly indebted. First: how did you check that the library is

Re: [R] Interaction plots as lines or bars?

2009-05-28 Thread Dieter Menne
Michael Kubovy wrote: An editor has suggested that I use bar plots to capture an interaction of two 2-level factors and an interaction of a 2 by 3 factorial experiment. (It would seem that there's a fear that someone might try to interpolate between, e.g., 'male' and 'female'.) In

Re: [R] question about using a remote system

2009-05-28 Thread Ted Harding
On 28-May-09 00:58:17, Erin Hodgess wrote: Dear R People: I would like to set up a plug-in for Rcmdr to do the following: I would start on a Linux laptop. Then I would log into another outside system and run a some commands. Now, when I tried to do system(ssh e...@xxx.edu) password

Re: [R] alternative to built-in data editor

2009-05-28 Thread Dieter Menne
urlwolf wrote: I often have to peek at large data. While head and tail are convenient, at times I'd like some more comprehensive. I guess I debug better in a more visual way? I was wondering if there's a way to override the default data editor. I have never seen the data editor. The

Re: [R] question about using a remote system

2009-05-28 Thread Mark Wardle
Hi. Do you need an interactive session at the remote machine, or are you simply wanting to run a pre-written script? If the latter, then you can ask ssh to execute a remote command, which conceivably could be R CMD x If you explain exactly why and what you are trying to do, then perhaps

Re: [R] How do I get removed from this mailing list?

2009-05-28 Thread Peter Dalgaard
Andrey Lyalko wrote: How do I get removed from this mailing list? [[alternative HTML version deleted]] __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide

[R] 1L and 0L

2009-05-28 Thread bogaso.christofer
Hi, Recently I come through those R-expressions and understood that 1L means 1 and 0L means 0. Why they are so? I mean, what the excess meanings they carry, instead writing simple 1 or 0? Is there any more this kind of expressions in R? Regards [[alternative HTML version deleted]]

Re: [R] Object-oriented programming in R

2009-05-28 Thread Mark Wardle
Hi. I remember considering these options myself but concluded that for most analyses a strictly procedural approach was satisfactory. Although I may re-run multiple analyses, the data manipulation (and subsequent analysis - the former always more complex than the latter IMHO) is fairly project-

Re: [R] How do I get removed from this mailing list?

2009-05-28 Thread Peter Dalgaard
Duncan Murdoch wrote: On 27/05/2009 8:25 PM, Andrey Lyalko wrote: How do I get removed from this mailing list? Like most lists nowadays, it gives the instructions in each message header: List-Unsubscribe: https://stat.ethz.ch/mailman/options/r-help,

Re: [R] How do I get removed from this mailing list?

2009-05-28 Thread Wacek Kusnierczyk
Peter Dalgaard wrote: Duncan Murdoch wrote: On 27/05/2009 8:25 PM, Andrey Lyalko wrote: How do I get removed from this mailing list? Like most lists nowadays, it gives the instructions in each message header: List-Unsubscribe: https://stat.ethz.ch/mailman/options/r-help,

Re: [R] ggplot2 adding vertical line at a certain date

2009-05-28 Thread Ivan Alves
check out geom_vline + geom_vline(xintercept=as.numeric(as.Date(2002-11-01))) [you may not need to convert the date to numeric in the most recent ggplot2 version] On 27 May 2009, at 20:31, stephen sefick wrote: library(ggplot2) melt.updn - (structure(list(date = structure(c(11808,

[R] optima in unimode

2009-05-28 Thread ms.com
Dear all i could not estimate the optima value or range value in unimodal plot in glm please help me out thanking you regard madan _ Hotmail® has a new way to see what's up with your friends. orial_WhatsNew1_052009

Re: [R] 1L and 0L

2009-05-28 Thread Gavin Simpson
On Thu, 2009-05-28 at 12:50 +0530, bogaso.christofer wrote: Hi, Recently I come through those R-expressions and understood that 1L means 1 and 0L means 0. Why they are so? I mean, what the excess meanings they carry, instead writing simple 1 or 0? Is there any more this kind of expressions

Re: [R] How do I get removed from this mailing list?

2009-05-28 Thread Gavin Simpson
On Thu, 2009-05-28 at 09:33 +0200, Wacek Kusnierczyk wrote: Peter Dalgaard wrote: Duncan Murdoch wrote: On 27/05/2009 8:25 PM, Andrey Lyalko wrote: How do I get removed from this mailing list? Like most lists nowadays, it gives the instructions in each message header:

Re: [R] optima in unimode

2009-05-28 Thread Gavin Simpson
On Thu, 2009-05-28 at 07:50 +, ms.com wrote: Dear all i could not estimate the optima value or range value in unimodal plot in glm please help me out thanking you regard madan You're going to have to give us more to go on than that, but... Are you an ecologist and you want to

Re: [R] Unable to load R

2009-05-28 Thread anupam sinha
Actually I did a search using locate but could not find the file. But when I do yum install readline (package containing the dependency libreadline.so.5) I get the following error: yum install readline Loading security plugin Loading rhnplugin plugin rhel-x86_64-client-5 100%

Re: [R] Unable to load R

2009-05-28 Thread Zeljko Vrba
On Thu, May 28, 2009 at 02:02:52PM +0530, anupam sinha wrote: Actually I did a search using locate but could not find the file. But when Locate reports useful results only if its database is up-to-date. Try running one of ldd /usr/lib/R/bin/exec/R ldd /usr/lib64/R/bin/exec/R and see what it

[R] Delta Kronecker

2009-05-28 Thread FMH
Hi, Could some give some ideas on how to compute the spatiotemporal covariance matrix by a sum of Kronecker products in R. Is there any special function that can be used? Cheers. Firdaus [[alternative HTML version deleted]] __

Re: [R] R help

2009-05-28 Thread Richard . Cotton
I am interested in modeling hydrological extreme events. I found MSClaio2008 very interesting function. In this function four criterions for choosing distributions. Can we call these criterions as model selection techniques or goodness of fit techniques or both? Because goodness of fit

[R] r-plot problem resolved

2009-05-28 Thread durden10
Thanks a lot for all the helpful comments! It finally works :handshake: (I settled with the code of Gavin) Best, Durden -- View this message in context: http://www.nabble.com/r-plot-tp23739356p23757458.html Sent from the R help mailing list archive at Nabble.com.

[R] R connectivity with Oracle DB

2009-05-28 Thread Madan Mohan
Hi Friends, I am working on R-2.9.0 and i want to connect oracle through R, but i could not find ROracle (zip file for windows) on CRAN Packages. Can anyone suggest me how to connect Oracle from R in windows platform? I would be grateful if u can provide me code aswell. Thanks Regards, Madan

Re: [R] How do I get removed from this mailing list?

2009-05-28 Thread Wacek Kusnierczyk
Gavin Simpson wrote: The View message source would be a more direct way of viewing the full email, not just the bits TBird shows you in the preview pane. I forget how it is named exactly and under which menu it is found as it has been quite a while since I used TBird. indeed; it's C-U,

Re: [R] no internal function int.unzip in R 2.9.0 for Windows

2009-05-28 Thread Uwe Ligges
The most recent version on CRAN is already fixed. We had a non-maintainer update end of March that fixed the R2HTML issues - the maintainer was unresponsive to our requests, unfortunately. Hence please update from CRAN. Best, Uwe Ligges Romain Francois wrote: Hi, I'll try to fix this

Re: [R] How do I get removed from this mailing list?

2009-05-28 Thread Patrick Connolly
On Thu, 28-May-2009 at 09:20AM +0200, Peter Dalgaard wrote: Duncan Murdoch wrote: On 27/05/2009 8:25 PM, Andrey Lyalko wrote: How do I get removed from this mailing list? Like most lists nowadays, it gives the instructions in each message header: List-Unsubscribe:

Re: [R] no internal function int.unzip in R 2.9.0 for Windows

2009-05-28 Thread Romain Francois
Hello Uwe, I have recently taken over maintenance of the package, which is now version controlled at r-forge (http://r-forge.r-project.org/projects/r2html/), and I am planning these changes: - rework the HTML function so that the html code is generated from a brew template, which in time

[R] sample unique pairs from a matrix

2009-05-28 Thread jos matejus
Dear R users, I have a matrix of both negative and positive values that I would like to randomly sample with the following 2 conditions: 1. only sample positive values 2. once a cell in the matrix has been sampled the row and column of that cell cannot be sampled from again. #some dummy data

Re: [R] How do I get removed from this mailing list?

2009-05-28 Thread Duncan Murdoch
On 28/05/2009 5:03 AM, Gavin Simpson wrote: On Thu, 2009-05-28 at 09:33 +0200, Wacek Kusnierczyk wrote: Peter Dalgaard wrote: Duncan Murdoch wrote: On 27/05/2009 8:25 PM, Andrey Lyalko wrote: How do I get removed from this mailing list? Like most lists nowadays, it gives the instructions

Re: [R] Object-oriented programming in R

2009-05-28 Thread Gabor Grothendieck
p.plot is an S3 object and you are attempting to define an S4 class so you if you really want to do that then you would need this to let S4 know about the ggplot class: # S4 class with S3 slot setClass(ggplot) # make ggplot visible to S4 setClass(AClass, representation(mFirst = numeric, mSecond

Re: [R] Labeling barplot bars by multiple factors

2009-05-28 Thread Jim Lemon
Thomas Levine wrote: I want to plot quantitative data as a function of three two-level factors. How do I group the bars on a barplot by level through labeling and spacing? Here http://www.thomaslevine.org/sample_multiple-factor_barplot.png's what I'm thinking of. Also, I'm pretty sure that I

Re: [R] R connectivity with Oracle DB

2009-05-28 Thread Prof Brian Ripley
See http://cran.r-project.org/bin/windows/contrib/2.9/@ReadMe (ReadMe files are there to be read!) and consider the use of RODBC, which is used by lots of Oracle users as an R client on Windows. On Thu, 28 May 2009, Madan Mohan wrote: Hi Friends, I am working on R-2.9.0 and i want to

Re: [R] 1L and 0L

2009-05-28 Thread Duncan Murdoch
On 28/05/2009 3:20 AM, bogaso.christofer wrote: Hi, Recently I come through those R-expressions and understood that 1L means 1 and 0L means 0. Why they are so? I mean, what the excess meanings they carry, instead writing simple 1 or 0? Is there any more this kind of expressions in R? Gavin

Re: [R] R connectivity with Oracle DB

2009-05-28 Thread Duncan Murdoch
On 28/05/2009 6:03 AM, Madan Mohan wrote: Hi Friends, I am working on R-2.9.0 and i want to connect oracle through R, but i could not find ROracle (zip file for windows) on CRAN Packages. Can anyone suggest me how to connect Oracle from R in windows platform? I would be grateful if u can

Re: [R] sample unique pairs from a matrix

2009-05-28 Thread Richard . Cotton
I have a matrix of both negative and positive values that I would like to randomly sample with the following 2 conditions: 1. only sample positive values 2. once a cell in the matrix has been sampled the row and column of that cell cannot be sampled from again. #some dummy data

Re: [R] consultation

2009-05-28 Thread David Winsemius
http://cran.r-project.org/web/views/Survival.html On May 11, 2009, at 1:46 PM, Nancy Tejerina wrote: *Dear R Users;* * * *I´m writing to ask you how can I do Survivals Curves using Time- dependent covariates? Which packages I need to Install?* *Thanks for your help, I look foward you

Re: [R] sample unique pairs from a matrix

2009-05-28 Thread David Winsemius
On May 28, 2009, at 6:33 AM, jos matejus wrote: Dear R users, I have a matrix of both negative and positive values that I would like to randomly sample with the following 2 conditions: 1. only sample positive values 2. once a cell in the matrix has been sampled the row and column of that

Re: [R] no internal function int.unzip in R 2.9.0 for Windows

2009-05-28 Thread Uwe Ligges
Will respond privately, since the thread becomes uninteresting for R-help, I guess. Uwe Romain Francois wrote: Hello Uwe, I have recently taken over maintenance of the package, which is now version controlled at r-forge (http://r-forge.r-project.org/projects/r2html/), and I am planning

[R] Read name multiple excel sheets using RODBC

2009-05-28 Thread simeon duckworth
I'd like to be able to read multiple sheets from an excel workbook and use the sheet name to name the resulting dataframe using RODBC. at the moment i've figured out how to do it the long way (see below) but feel sure that there is a speedier possibly automatic way to do it in R. i've tried to

Re: [R] Full likelihood from survreg

2009-05-28 Thread Terry Therneau
How can I find out what survreg generates: the full likelihood or a likelihood with unnecessary constants dropped? Survreg returns a log-likelihood, not a deviance (LL with parameters dropped); and the result is labeled as such. It turns out that the deviance is not as easily defined for

Re: [R] How do I get removed from this mailing list?

2009-05-28 Thread Peter Dalgaard
Wacek Kusnierczyk wrote: Gavin Simpson wrote: The View message source would be a more direct way of viewing the full email, not just the bits TBird shows you in the preview pane. I forget how it is named exactly and under which menu it is found as it has been quite a while since I used TBird.

Re: [R] Labeling barplot bars by multiple factors

2009-05-28 Thread Thomas Levine
Both of those worked, but hierobarp looked a bit easier, so I used that. The one annoying thing is that it sorts alphabetically. Tom On Thu, May 28, 2009 at 6:46 AM, Jim Lemon j...@bitwrit.com.au wrote: Thomas Levine wrote: I want to plot quantitative data as a function of three two-level

Re: [R] sample unique pairs from a matrix

2009-05-28 Thread jos matejus
Dear Ritchie and David, Thanks very much for your advice. I had thought of this potential solution, however it doesn't really fullfill my second criteria which is that once a particular cell has been sampled, the row and column of that cell can't be sampled from subsequently. In other words, the

Re: [R] sample unique pairs from a matrix

2009-05-28 Thread jim holtman
This is one way of sampling all the values in the matrix that meet your criteria: set.seed(1) (x - matrix(rnorm(100), 10)) [,1][,2][,3][,4] [,5] [,6][,7] [,8] [,9] [,10] [1,] -0.6264538 1.51178117 0.91897737 1.35867955

Re: [R] split strings

2009-05-28 Thread Wacek Kusnierczyk
(diverted to r-devel, a source code patch attached) Wacek Kusnierczyk wrote: Allan Engelhardt wrote: Immaterial, yes, but it is always good to test :) and your solution *is* faster and it is even faster if you can assume byte strings: :) indeed; though if the speed is immaterial

Re: [R] Re ad name multiple excel sheets using RODBC

2009-05-28 Thread Dieter Menne
simeon duckworth wrote: I'd like to be able to read multiple sheets from an excel workbook and use the sheet name to name the resulting dataframe using RODBC. In Microsoft theory, something like the below should be ok (note the $), but never managed to get this to work. The same method

Re: [R] Still can't find missing data

2009-05-28 Thread Dieter Menne
Farley, Robert wrote: I can't get the syntax that will allow me to show NA values (rows) in the xtabs. lengthy non-reproducible example removed If you want a reproducible answer, prepare a reproducible result. And check that the syntax is na.action=na.pass Dieter -- View this

Re: [R] question about using a remote system

2009-05-28 Thread Erin Hodgess
My goal is for a user to sit down at a Linux laptop, get to an Rcmdr type screen, submit jobs on a remote system and then get the results back in R. We will assume that the user is naive, and the only thing he/she can do is get to the Rcmdr screen. The Rcmdr plugin will have a submit jobs menu.

[R] How could I terminate a script (without leaving R)?

2009-05-28 Thread Mario Valle
Good morning! Which is the preferred method to leave a sourced script and returning back to the '' prompt? For example I search for certain files to be processed, but nothing should be done if they are not present. Normally I do in my script: f-dir(pattern=qq) if(length(f) 0) { ...process

Re: [R] R-help Digest, Vol 75, Issue 28

2009-05-28 Thread Ista Zahn
From: Jeff Newmiller jdnew...@dcn.davis.ca.us To: R Heberto Ghezzo, Dr heberto.ghe...@mcgill.ca Date: Wed, 27 May 2009 09:00:44 -0700 Subject: Re: [R] R in Ubunto R Heberto Ghezzo, Dr wrote:    Hello , I do not know anything abount Ubunto, but I found a Portable Ubunto for Windows

[R] can you help me please :)

2009-05-28 Thread NOUF AL NUMAIR
hi there :) i want to use barplot with if else but i dont know how to do it ? i tried this but it is not working with me SNP - read.table(my.txt) SNP[,2] [1] 1175 483 240 170 99 79 76 45 38 35 21 16 14 19 16 [16]333 10216868

Re: [R] How could I terminate a script (without leaving R)?

2009-05-28 Thread Duncan Murdoch
On 5/28/2009 8:58 AM, Mario Valle wrote: Good morning! Which is the preferred method to leave a sourced script and returning back to the '' prompt? For example I search for certain files to be processed, but nothing should be done if they are not present. Normally I do in my script:

Re: [R] Re ad name multiple excel sheets using RODBC

2009-05-28 Thread Hans-Peter Suter
2009/5/28 Dieter Menne dieter.me...@menne-biomed.de: If you only have the sheet names, you should use package xlsReadWrite which is rather fast, but has some limitations in the non-commercial version. what limitations, i.e. features do you miss? Cheers, Hans-Peter

Re: [R] [Rd] split strings

2009-05-28 Thread William Dunlap
Bill Dunlap TIBCO Software Inc - Spotfire Division wdunlap tibco.com -Original Message- From: r-devel-boun...@r-project.org [mailto:r-devel-boun...@r-project.org] On Behalf Of Wacek Kusnierczyk Sent: Thursday, May 28, 2009 5:30 AM Cc: R help project; r-de...@r-project.org; Allan

Re: [R] Re ad name multiple excel sheets using RODBC

2009-05-28 Thread Dieter Menne
Hans-Peter Suter wrote: If you only have the sheet names, you should use package xlsReadWrite which is rather fast, but has some limitations in the non-commercial version. what limitations, i.e. features do you miss? Reading of named ranges. Dieter -- View this message in

Re: [R] Re ad name multiple excel sheets using RODBC

2009-05-28 Thread Erich Neuwirth
Using the package rcom (by Thomas Baier) you can do everything you want since you have the full power of COM at your disposal. Dieter Menne wrote: Hans-Peter Suter wrote: If you only have the sheet names, you should use package xlsReadWrite which is rather fast, but has some limitations

[R] re gularly spaced points on multi-segement line

2009-05-28 Thread PLAFF
Dear all ... I need your help for the following purpose : I would like to create regularly spaced points on a multi-segment line (a 'psp class' object). A function of the spatstat library( = pointsOnLines() ) is dedicated to that objective but the problem is that the probability of falling on

[R] Changing data point size in quilt.plot

2009-05-28 Thread Wilson, R.C.
Hi, I am new to R and am using quilt.plot (from fields package) to plot gas concentrations in Europe. I am using the following code to plot my data library(maps) library(fields) test = read.csv(%change 1996_2005.txt, sep=\t) colnames(test) = c(Station, Measurement_BaseO3, Model_BaseO3,

[R] Urgent Statistician Position near Philadelphia

2009-05-28 Thread José Ignacio Bustos Melo
Dear R-list, I know is going to be weird, but I need to do something about and I know the R-list is big group of good people. I’m moving to Philadelphia in August 2009 (my fiancé is there) and I need to find a job. I know to do so many things, but my skill is in science and research. I would like

[R] Plot error

2009-05-28 Thread NatsS
Hello, I am an R amateur. I want to plot data such that the 3 time points(a,b,c) lie on the X-axis and the values of these times points are on Y-axis for n samples (e.g.100). So, I have an object x, dim 100 4, it is a dataframe (when checked the class) x = name a b c 1

Re: [R] Plot error

2009-05-28 Thread Richard . Cotton
I want to plot data such that the 3 time points(a,b,c) lie on the X-axis and the values of these times points are on Y-axis for n samples (e.g.100). So, I have an object x, dim 100 4, it is a dataframe (when checked the class) x = name a b c 10.11 1.11 0.86 2

Re: [R] question about using a remote system

2009-05-28 Thread Greg Snow
Will R be able to send one batch of commands to the remote machine, then wait for the output? Or will there be back and forth with new commands based on the output from the first commands? You may want to look at the nws package for one way to have jobs run on a remote machine. -- Gregory

[R] maximum over one dimension of a 3-dimensional array

2009-05-28 Thread eric lee
Hi, I'm running R 2.7.2 on windows XP. I'd like to find the maximum of a 3-d array over it's third index to create a 2-d array. For example: x - array(c(1,2,3,10,11,12,3:8),c(2,3,2)) x , , 1 [,1] [,2] [,3] [1,]13 11 [2,]2 10 12 , , 2 [,1] [,2] [,3] [1,]3

Re: [R] sample unique pairs from a matrix

2009-05-28 Thread David Winsemius
Your last step will either be a single number (not really a sampling operation) or a non-positive number. So at best you really only have an number that depends entirely on the prior sequence of draws. -- David. On May 28, 2009, at 8:21 AM, jos matejus wrote: Dear Ritchie and David,

Re: [R] maximum over one dimension of a 3-dimensional array

2009-05-28 Thread David Winsemius
On May 28, 2009, at 11:25 AM, eric lee wrote: Hi, I'm running R 2.7.2 on windows XP. I'd like to find the maximum of a 3-d array over it's third index to create a 2-d array. For example: x - array(c(1,2,3,10,11,12,3:8),c(2,3,2)) x , , 1 [,1] [,2] [,3] [1,]13 11 [2,]2

[R] String replacement in an expression

2009-05-28 Thread Caroline Bazzoli
Dear R-experts, I need to replace in an expression the character Cl by Cl+beta But in the following case: form-expression((Cl-(V *ka) ) +(V *Vm *exp(-(Clm/Vm) *t))) gsub(Cl,(Cl+beta),as.character(form)) We obtain: [1] ((Cl+beta) - (V * ka)) + (V * Vm * exp(-((Cl+beta)m/Vm) * t))

Re: [R] String replacement in an expression

2009-05-28 Thread Gabor Grothendieck
Try matching on word boundaries as well: gsub(\\bCl\\b,(Cl+beta),as.character(form)) [1] ((Cl+beta) - (V * ka)) + (V * Vm * exp(-(Clm/Vm) * t)) See ?regexp On Thu, May 28, 2009 at 11:41 AM, Caroline Bazzoli caroline.bazz...@inserm.fr wrote: Dear R-experts, I need to replace in an expression

[R] getTree visualization (randomForest)

2009-05-28 Thread Károly Kovács
Dear All, I would like to visualize a tree extracted from a random forest using getTree {randomForest}. I'm wondering if there is any way to do it directly or to convert my tree to any other class of tree repesentation, which then can be plotted? Thank you in advance, Karoly

Re: [R] can you help me please :)

2009-05-28 Thread Stefan Grosse
On Thu, 28 May 2009 15:47:18 +0300 NOUF AL NUMAIR noufalnum...@hotmail.com wrote: NAN barplot(zz,width = 4,names.arg= xx,axes = TRUE, axisnames = TRUE, NAN main=title,xlab=xlab,ylab=ylab,ylim=c(0,1175),xlim=c(0,226),col = NAN c(for (i in zz){if i70 col= lightblue else col= mistyrose)) NAN NAN

Re: [R] Neural Network resource

2009-05-28 Thread jude.ryan
The package AMORE appears to be more flexible, but I got very poor results using it when I tried to improve the predictive accuracy of a regression model. I don't understand all the options well enough to be able to fine tune it to get better predictions. However, using the nnet() function in

Re: [R] String replacement in an expression

2009-05-28 Thread William Dunlap
From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org] On Behalf Of Caroline Bazzoli Sent: Thursday, May 28, 2009 8:41 AM To: r-help@r-project.org Subject: [R] String replacement in an expression Dear R-experts, I need to replace in an expression the character Cl by

Re: [R] boxplot

2009-05-28 Thread Dieter Menne
amor Gandhi wrote: Hi gues, This should read: Hi, guess what I want   amor Gandhi wrote: Is there any function in R for boxplot with different time points? t1 - c(rep(1,20),rep(2,20)) t2 - c(rep(1,10),rep(2,10),rep(1,10),rep(2,10)) x - rnorm(40,5,1) dat - data.frame(t1,t2,x)

Re: [R] Hmisc package: deff() command's formula for the design effect

2009-05-28 Thread jjh21
One additional question: It seems that in practice the design effect is often calculated on the response variable or on individual predictor variables. I have some OLS models using observational data that are clustered. Does it make sense to calculate the design effect on the residuals of one of

[R] ggplot2 legend

2009-05-28 Thread Felipe Carrillo
Hi: I need some help with the legend. I got 14 samples(Muestreo) and I am trying to plot a smooth line for each sample. I am able to accomplish that but the problem is that the legend only displays every other sample. How can I force the legend to show all of my Muestreos? Thanks in

Re: [R] Changing point color/character in qqmath

2009-05-28 Thread Kevin W
Apologies once more, there was a slight error in our example. Here is a correct version of how to color the points in qqmath according to another variable. qqmath(~ yield | variety, data = barley, groups=year, auto.key=TRUE, prepanel = function(x, ...) { list(xlim =

[R] plotting time series with data gap using type line- but do not want to connect gap with line

2009-05-28 Thread Josef . Kardos
I have a time series of about 1500 measurements. There are sporadic data gaps. I would like to plot the time series with type line. I only want a line to connect the periods with data; I don't want a line to connect the points across a data gap. Is there a function or recommended method to

Re: [R] plotting time series with data gap using type line- but do not want to connect gap with line

2009-05-28 Thread Gabor Grothendieck
Just place a point with an NA value between the two segments. Here is one way to do it in zoo: set.seed(123) library(zoo) # create sample data tt - c(1:1000, 1200:2000) z - zoo(rnorm(length(tt)), tt) # this will fill in omitted values with NAs z - as.zoo(as.ts(z)) plot(z) On Thu, May 28,

[R] installing gRain package on ubuntu 8.10

2009-05-28 Thread Alex Biedermann
Dear all I am trying to install the gRain package on ubuntu 8.10, running R Vers. 2.9. Unfortunately, I get a bad Exit-Status. What could cause that problem? Any help is appreciated. Many thanks in advance. A. Biedermann install.packages(gRain) Warnung in install.packages(gRain) :

[R] how to avoid add 'X' before numeric colnames when read.table

2009-05-28 Thread Zheng, Xin (NIH) [C]
Hi all, Is there any way to keep numeric colnames as is? Any hint will be appreicated. Xin Zheng [[alternative HTML version deleted]] __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the

Re: [R] String replacement in an expression

2009-05-28 Thread Wacek Kusnierczyk
Caroline Bazzoli wrote: Dear R-experts, I need to replace in an expression the character Cl by Cl+beta But in the following case: form-expression((Cl-(V *ka) ) +(V *Vm *exp(-(Clm/Vm) *t))) gsub(Cl,(Cl+beta),as.character(form)) We obtain: [1] ((Cl+beta) - (V * ka)) + (V * Vm *

Re: [R] how to avoid add 'X' before numeric colnames when read.table

2009-05-28 Thread jim holtman
?read.table check.names=FALSE On Thu, May 28, 2009 at 3:25 PM, Zheng, Xin (NIH) [C] zheng...@mail.nih.gov wrote: Hi all, Is there any way to keep numeric colnames as is? Any hint will be appreicated. Xin Zheng [[alternative HTML version deleted]]

Re: [R] how to avoid add 'X' before numeric colnames when read.table

2009-05-28 Thread Greg Snow
This really depends on where the 'X' is coming from, which you did not tell us (see the posting guide). For example, if the 'data.frame' function is the one adding the 'X', then you can use the 'check.names' argument to prevent the addition (or you can use data.frame with the argument to

Re: [R] how to avoid add 'X' before numeric colnames when read.table

2009-05-28 Thread Zheng, Xin (NIH) [C]
Thank you guys. But I also want 'check.names' to help keeping colnames unique. It seems there's no proper option for that in 'read.table'. I have to alter colnames with strsplit or other similar functions. From: jim holtman [mailto:jholt...@gmail.com] Sent: Thursday, May 28, 2009 3:40 PM To:

[R] How this addition works?

2009-05-28 Thread bogaso.christofer
I have following addition : 1:2 + 1:10 [1] 2 4 4 6 6 8 8 10 10 12 I could not understand how R adding those two unequal vector? Any help? [[alternative HTML version deleted]] __ R-help@r-project.org mailing list

Re: [R] how to avoid add 'X' before numeric colnames when read.table

2009-05-28 Thread Zheng, Xin (NIH) [C]
substring(header, 2,) works for the purpose perfectly. Xin Zheng [[alternative HTML version deleted]] __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide

[R] filter function speed

2009-05-28 Thread Bishara, Anthony J
I am trying to use a linear filter to reduce loops and thereby increase the speed of an existing program. However, while the filter function (stats package) should have reduced the looping by about 30-fold, the time to complete the program remained about the same. This surprised me, because I

Re: [R] String replacement in an expression

2009-05-28 Thread William Dunlap
Bill Dunlap TIBCO Software Inc - Spotfire Division wdunlap tibco.com -Original Message- From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org] On Behalf Of Wacek Kusnierczyk Sent: Thursday, May 28, 2009 12:31 PM To: Caroline Bazzoli Cc: r-help@r-project.org

Re: [R] Labeling barplot bars by multiple factors

2009-05-28 Thread Thomas Levine
Ah, that makes sense. But now another two issues have arisen. Firstly, the error bars look like confidence intervals, and I'm pretty sure that they are but does some document verify this? I suppose I could check the code too. Secondly, I just read about how dynamite plots should be avoided. It's

Re: [R] How this addition works?

2009-05-28 Thread Rolf Turner
On 29/05/2009, at 8:00 AM, bogaso.christofer wrote: I have following addition : 1:2 + 1:10 [1] 2 4 4 6 6 8 8 10 10 12 I could not understand how R adding those two unequal vector? Any help? Look at the help for ``+'' (?+) and look at ``Value''. There you will see: These

Re: [R] How this addition works?

2009-05-28 Thread baptiste auguie
recycling rule: repeat the shorter element as many times as necessary, all.equal(1:2 + 1:10 , rep(1:2, length=10) + 1:10) # TRUE HTH, baptiste On 28 May 2009, at 22:00, bogaso.christofer wrote: I have following addition : 1:2 + 1:10 [1] 2 4 4 6 6 8 8 10 10 12 I could not

Re: [R] How this addition works?

2009-05-28 Thread Sarah Goslee
R recycles the shorter one to match the longer one: 1 2 3 4 5 6 7 8 9 10 + 1 2 1 2 1 2 1 2 1 2 = 2 4 4 6 6 8 8 10 10 12 R does this recycling in many cases, and it can sometimes trap the unwary. Sarah On Thu, May 28, 2009 at 4:00 PM, bogaso.christofer

Re: [R] ggplot2 legend

2009-05-28 Thread Mike Lawrence
First, your example didn't work because fish_ByMuestreo didn't build properly (deleting the first structure(list(data = bit solves this). To solve your plotting problem, note that as constructed, the Muestreo column is numeric, whereas you seem to want to treat it as a factor. Solution: convert

Re: [R] How this addition works?

2009-05-28 Thread Luc Villandre
bogaso.christofer wrote: I have following addition : 1:2 + 1:10 [1] 2 4 4 6 6 8 8 10 10 12 I could not understand how R adding those two unequal vector? Any help? [[alternative HTML version deleted]] __

[R] max.col weirdness

2009-05-28 Thread Daryl Morris
Hi, I think there's some rounding issue with returning the max column. (running 2.9.0 on an Apple, but my buddy found it on his PC) x - matrix(c(1234.568,1234.569,1234.567),1) max.col(x) [1] 2 x - matrix(c(12345.568,12345.569,12345.567),1) max.col(x) [1] 3 x -

Re: [R] ggplot2 legend

2009-05-28 Thread Felipe Carrillo
Thanks for your help Mike, it works like a charm now!! --- On Thu, 5/28/09, Mike Lawrence mike.lawre...@dal.ca wrote: From: Mike Lawrence mike.lawre...@dal.ca Subject: Re: [R] ggplot2 legend To: Felipe Carrillo mazatlanmex...@yahoo.com Cc: r-h...@stat.math.ethz.ch Date: Thursday, May 28,

Re: [R] sample unique pairs from a matrix

2009-05-28 Thread jos matejus
Dear All, Many thanks for all your useful suggestions. Much appreciated. Jos 2009/5/28 David Winsemius dwinsem...@comcast.net: Your last step will either be a single number (not really a sampling operation) or a non-positive number. So at best you really only have an number that depends

Re: [R] max.col weirdness

2009-05-28 Thread Daryl Morris
Hi, Sorry, I was too hasty. I see that the tolerance is in the documentation (and that random is the default tie-breaker). Would it be possible to allow tolerance to be a parameter? thanks, Daryl Daryl Morris wrote: Hi, I think there's some rounding issue with returning the max column.

Re: [R] max.col weirdness

2009-05-28 Thread Bert Gunter
Try reading the man page, which says: Details When ties.method = random, as per default, ties are broken at random. In this case, the determination of a tie assumes that the entries are probabilities: there is a relative tolerance of 1e-5, relative to the largest (in magnitude, omitting

[R] custom sort?

2009-05-28 Thread Steve Jaffe
Sounds simple but haven't been able to find it in docs: is it possible to sort a vector using a user-defined comparison function? Seems it must be, but sort doesn't seem to provide that option, nor does order sfaics -- View this message in context:

[R] why warning appears when set both 'scale' and 'breaks' meanwhile in heatmap.2 of gplots

2009-05-28 Thread Zheng, Xin (NIH) [C]
The warning I met with is attached here: Warning message: In heatmap.2(as.matrix(val.sf.ns[-1:-3]), Rowv = F, Colv = T, dendrogram = column, : Using scale=row or scale=column when breaks arespecified can produce unpredictable results.Please consider using only one or the other. What does

[R] avoid a loop

2009-05-28 Thread KARAVASILIS GEORGE
Hello, R users. I have the following code: a=1:10 b=-3:15 n=5 x - rep(0,n) for (i in 1:n) x[i] - sum( outer(a,b, function(s,t) abs(a-b-i)==0) ) Can someone tell me if I could avoid the for command? Thank you in advance. __ R-help@r-project.org

  1   2   >