[R] kruskal wallis post hoc?

2012-01-11 Thread Iasonas Lamprianou
Dear all, I run a kruskal wallis test and found significant results. Then, I conducted all pairwise comparisons and found no significant results. Could anyone please give me a hint as to why this happens or redirect me towards a specific web page where I can find more info? (I used alpha=5% an

Re: [R] R error in make check

2012-01-11 Thread Prof Brian Ripley
You should ask this on R-devel, and tell them exactly what compilers you used. See the posting guide: this falls under 'Questions likely to prompt discussion unintelligible to non-programmers should go to to R-devel.' But yes, it is something to worry about and it looks as if your compiler

Re: [R] summing NAs in aggregate

2012-01-11 Thread R. Michael Weylandt
Perhaps this would work: spitzSum <- function(x) if(all(is.na(x))) NA else sum(x, na.rm = TRUE) Michael On Wed, Jan 11, 2012 at 11:18 PM, Spitzer, Matthew wrote: > Hello, > I would like to please ask for assistance with aggregate sum.  I have a data > set with consisting of two grouping variab

Re: [R] cannot get tweets from twitteR

2012-01-11 Thread Sachinthaka Abeywardana
Btw I tried this on the terminal in R and it worked, restarted Rstudio and it still doesn't work there. Rstudio bug or something else? Sachin [[alternative HTML version deleted]] __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman

[R] summing NAs in aggregate

2012-01-11 Thread Spitzer, Matthew
Hello, I would like to please ask for assistance with aggregate sum. I have a data set with consisting of two grouping variables (id, visit) and several other variables. I would like to sum the variables for each id and visit, but am having problems with na.rm. na.rm=TRUE seems to replace all

[R] Points inside a polygon

2012-01-11 Thread Hasan Diwan
I have a list of bounds for a series of polygons. I do understand the formula to determine whether point i is within polygon X (X[x1] < i[x] & X[x2] > i[x] & X[y1] < i[y] & X[y2] > i[y]), and I can apply this throughout the dataset. However, this naive algorithm doesn't scale very well. The data se

Re: [R] RGL- Drawing Circle

2012-01-11 Thread Duncan Murdoch
On 12-01-11 9:08 AM, gli wrote: thanks Duncan. I think want my circle to be in user coordinates. I tried the first code u gave me and it gives me a circle. but how can i: 1) change the radius ? 2) place the circle at a given x,y,z coordinate? 3) turn it 90 degree up like these circle plate b

Re: [R] summarizing a complex dataframe

2012-01-11 Thread jim holtman
will this do it for you: > x <- read.table(text = " m1_1 m2_1 m3_1 m1_2 m2_2 m3_2 +111222 + 211222 + 221222", header = TRUE) > # split out the main names of the column > x.names <- do.call(rbind, strsplit(names(x), "_")) > x

[R] cannot get tweets from twitteR

2012-01-11 Thread Sachinthaka Abeywardana
Hi all, I was following the tutorial in: http://www.rdatamining.com/examples/text-mining and using the package twitteR i attempted to execute userTimeline and got the error: rdmTweets<-userTimeline("rdatamining",n=100)Error in .self$twFromJSON(out) : Error: Rate limit exceeded. Clients may not

[R] multcomp two-way anova with interactions within and between

2012-01-11 Thread Schreiber, Stefan
Hi all, I'd like to compare all levels of my interaction with each other. I read the pdf 'Additional multcomp Examples' but even though there is an example with an interaction it doesn't work for me when I want to compare within and between groups. Here is an example: d.fr<-data.frame(id=r

Re: [R] 2D filter in R?

2012-01-11 Thread David Winsemius
On Jan 11, 2012, at 6:21 PM, Michael wrote: Okay thanks guys! I was actually looking for image processing libraries in R. Specifically, I am looking for something like "filter2d" for images(pixels)... (speed is not a concern for me, since my images are really small, 200x200, etc.) And in

Re: [R] plotOHLC(alpha3): Error in plotOHLC(alpha3) : x is not a open/high/low/close time series

2012-01-11 Thread Ted Byers
Thanks Joshua, That did it. Cheers, Ted -- View this message in context: http://r.789695.n4.nabble.com/plotOHLC-alpha3-Error-in-plotOHLC-alpha3-x-is-not-a-open-high-low-close-time-series-tp4283217p4286963.html Sent from the R help mailing list archive at Nabble.com. __

[R] Error in charToDate(x)

2012-01-11 Thread Monchichi
Dear all, I have a problem while working with hourly data of fx rates. I've read from a csv file, the following way: csv-file like: Date,Open,High,Low,Close,Volume 2011-08-11 03:00:00,1.41758,1.42205,1.41625,1.42174,8974 ... 2011-08-12 04:00:00,1.42175,1.42413,1.42067,1.42172,7229 ... 2011-12-30 05

[R] R error in make check

2012-01-11 Thread Anirban Jana
Hi, After building R 2.14.1 on a Linux system (SUSE Linux Enterprise Server 11 (x86_64) ), I ran "make check" but it halted with the foll error: Testing examples for package âbaseâ Error: testing 'base' failed Execution halted I looked in /test/Examples/base-Ex.Rout.fail and found the foll: **

Re: [R] Vegan(ordistep) error: Error in if (aod[1, 5] <= Pin) { : missing value where TRUE/FALSE needed

2012-01-11 Thread Nevil Amos
Whilst the constrained proportion was 0 in the example below I am getting the error in cases where this is not so. See new output pasted below [1] "ORDISTEP RDA1 BOTH WAYS" [1] "" [1] "SFW ALL Adult

Re: [R] 2D filter in R?

2012-01-11 Thread Michael
Thanks a lot Bert! Are z, y and x are of the same length? And after the loess I just "reshape" the "stacked" outputs back into a 2d matrix? Thanks again! On Wed, Jan 11, 2012 at 5:33 PM, Bert Gunter wrote: > You have to set the data up as 3 columns: > > z = Values > y = row > x = column >

Re: [R] 2D filter in R?

2012-01-11 Thread Michael
I looked the "loess" up... But I am confused... what should be my "formula" if I just want to smooth the numbers of my matrix x? And those lots of parameters are daunting... what are the essential ones that I should definitely set up? I remember for 2d Gaussian filter, I just need to set one param

Re: [R] 2D filter in R?

2012-01-11 Thread Bert Gunter
The Bioconductor list/repository will have lots of image processing stuff. -- Bert On Wed, Jan 11, 2012 at 3:21 PM, Michael wrote: > Okay thanks guys! > > I was actually looking for image processing libraries in R. > > Specifically, I am looking for something like "filter2d" for > images(pixels)

Re: [R] 2D filter in R?

2012-01-11 Thread Michael
Okay thanks guys! I was actually looking for image processing libraries in R. Specifically, I am looking for something like "filter2d" for images(pixels)... (speed is not a concern for me, since my images are really small, 200x200, etc.) And in fact, these are not images, these are 2d matrix nu

Re: [R] 2D filter in R?

2012-01-11 Thread Bert Gunter
?loess ## in base R does 2d filtering , though it's not a Gaussian one -- Bert On Wed, Jan 11, 2012 at 3:00 PM, Rolf Turner wrote: > On 12/01/12 09:17, Michael wrote: >> >> Hi all, >> >> I am looking for a command for doing 2D filtering (rectangular or >> Gaussian) in R... >> >> I have looked

Re: [R] 2D filter in R?

2012-01-11 Thread Rolf Turner
On 12/01/12 09:17, Michael wrote: Hi all, I am looking for a command for doing 2D filtering (rectangular or Gaussian) in R... I have looked at ksmooth, filter and convolve but they seem to be 1D... Any thoughts? Thanks a lot! I'm not sure what you really need/want to do, but you might have

Re: [R] Storing/Restoring R objects

2012-01-11 Thread David Winsemius
On Jan 11, 2012, at 5:16 PM, Rich Shepard wrote: On Wed, 11 Jan 2012, David Winsemius wrote: Then R will not find it. (Actually there will be no .RData file until R executes save.image() at the end of the session. R does not create a fresh .RData at the beginning of a session.) David,

Re: [R] Restricting R session

2012-01-11 Thread Thomas Lumley
2012/1/12 Uwe Ligges : > > > On 10.01.2012 20:30, Antonio Rodriges wrote: >> >> Hello, >> >> Is it possible to use R on public server where each user has its own >> restricted R session? > > > This entirely depends on the definition of "restricted", otherwise the > answer is "yes". More precisely,

Re: [R] rjava on FreeBSD

2012-01-11 Thread Rainer Hurling
On 11.01.2012 22:51 (UTC+1), Hasan Diwan wrote: Trying to install Rjava on FreeBSD 9 and am getting the following error: install.packages('rJava') trying URL 'http://cran.cnr.Berkeley.edu/src/contrib/rJava_0.9-3.tar.gz' Content type 'application/x-gzip' length 537153 bytes (524 Kb) opened URL =

Re: [R] summarizing a complex dataframe

2012-01-11 Thread David Winsemius
On Jan 11, 2012, at 3:55 PM, Christopher G Oakley wrote: I need some help summarizing complex data frames (small example below): m1_1 m2_1 m3_1 m1_2 m2_2 m3_2 i1111222 i1211222 i2221222 For an arbitrary number of co

Re: [R] Storing/Restoring R objects

2012-01-11 Thread Rich Shepard
On Wed, 11 Jan 2012, David Winsemius wrote: Then R will not find it. (Actually there will be no .RData file until R executes save.image() at the end of the session. R does not create a fresh .RData at the beginning of a session.) David, That's good to know. Could I then specify which .RDa

[R] Checking dates for entry errors

2012-01-11 Thread Paul Miller
Hello Everyone,   I have a question about how best to check dates for entry errors. I recently discovered that R will read the incorrectly entered date "11/23/21931" without producing a warning or an error message at least under some circumstances.   > as.Date("11/23/21931", format = "%m/%d/%Y")

Re: [R] summarizing a complex dataframe

2012-01-11 Thread Bert Gunter
Well, if I understand what you want to do, it's straightforward, ?"[" (pay attention to the use of column names) and ?grep would pick out the columns you want and you could then use mapply or maybe rowMeans or whatever to get your summaries. HOWEVER ... I think what you should really should do is

Re: [R] summarizing a complex dataframe

2012-01-11 Thread Steve Lianoglou
Hi, On Wed, Jan 11, 2012 at 3:55 PM, Christopher G Oakley wrote: > I need some help summarizing complex data frames (small example below): > >    m1_1 m2_1 m3_1 m1_2 m2_2 m3_2 > i1    1    1    1    2    2    2 > i1    2    1    1    2    2    2 > i2    2    2    1    2    2    2 > > > For an arb

Re: [R] Storing/Restoring R objects

2012-01-11 Thread David Winsemius
On Jan 11, 2012, at 2:56 PM, Rich Shepard wrote: One of my projects has generated quite a few objects (data frames) related to one portion of this project. They can be listed with the ls() function. What I would like to do is move them to another directory so that data frames for other por

[R] rjava on FreeBSD

2012-01-11 Thread Hasan Diwan
Trying to install Rjava on FreeBSD 9 and am getting the following error: > install.packages('rJava') trying URL 'http://cran.cnr.Berkeley.edu/src/contrib/rJava_0.9-3.tar.gz' Content type 'application/x-gzip' length 537153 bytes (524 Kb) opened URL ==

Re: [R] Storing/Restoring R objects

2012-01-11 Thread Jeff Newmiller
?save ?load You may need to create single use environment objects to hold the whole file while you separate multiple objects. --- Jeff NewmillerThe . . Go Live... DCN:Basics:

Re: [R] rjags installation trouble

2012-01-11 Thread Ben Bolker
On 12-01-11 11:45 AM, Martyn Plummer wrote: > It looks like you are cross-linking to an earlier version of the JAGS > library at run time. Check > > /sbin/ldconfig -p | grep jags - > > When compiling rjags, you can hard-code the location of the jags library > using the [GNU-specific] configure op

Re: [R] turning a list of vectors into a data.frame (as rows of the DF)?

2012-01-11 Thread Marc Schwartz
On Jan 11, 2012, at 1:40 PM, Chris Conner wrote: > As a newer R practicioner, it seems I stump myself weekly (at least) with > issues that have spinning my wheels. Here is yet another... I'm trying to > turn a list of numeric vectors (of uneual length) inot a dataframe. Each > vector held in

Re: [R] turning a list of vectors into a data.frame (as rows of the DF)?

2012-01-11 Thread William Dunlap
Perhaps the following does what you want. It extends each element of your list to a common length, converts that to a matrix, then to a data.frame: f <- function(data) { nCol <- max(vapply(data, length, 0)) data <- lapply(data, function(row) c(row, rep(NA, nCol-length(row data <

Re: [R] Confidence Interval from Moments?

2012-01-11 Thread R. Michael Weylandt
Assuming a distribution defined solely by those moments it is possible (e.g., z- or t-test confidence intervals) but this isn't really the place to discuss such things since there's no R content to your question: try stats.stackexchange.com Michael On Wed, Jan 11, 2012 at 4:56 AM, lambdatau wrot

[R] summarizing a complex dataframe

2012-01-11 Thread Christopher G Oakley
I need some help summarizing complex data frames (small example below): m1_1 m2_1 m3_1 m1_2 m2_2 m3_2 i1111222 i1211222 i2221222 For an arbitrary number of columns (say m1 …. m199) where the column names have

Re: [R] New to R, Curious about Project Idea

2012-01-11 Thread Daniel Nordlund
> -Original Message- > From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org] > On Behalf Of arbeaupg > Sent: Wednesday, January 11, 2012 6:36 AM > To: r-help@r-project.org > Subject: [R] New to R, Curious about Project Idea > > Good morning, > > I am a student whom is c

Re: [R] RGL- Drawing Circle

2012-01-11 Thread Uwe Ligges
On 11.01.2012 15:08, gli wrote: thanks Duncan. I think want my circle to be in user coordinates. I tried the first code u gave me and it gives me a circle. but how can i: 1) change the radius ? 2) place the circle at a given x,y,z coordinate? 3) turn it 90 degree up like these circle plate

Re: [R] Storing/Restoring R objects

2012-01-11 Thread Uwe Ligges
On 11.01.2012 20:56, Rich Shepard wrote: One of my projects has generated quite a few objects (data frames) related to one portion of this project. They can be listed with the ls() function. What I would like to do is move them to another directory so that data frames for other portions of the

Re: [R] Storing/Restoring R objects

2012-01-11 Thread R. Michael Weylandt
R is the natural tool to operate on a .RData file: you can name it however you wish either with your OS or with the save() command in R. You can load any .RData files with the load() command but the startup routine only looks for .RData (to my knowledge) unless you put specific instructions in you

Re: [R] turning a list of vectors into a data.frame (as rows of the DF)?

2012-01-11 Thread R. Michael Weylandt
Most methods take the rows of data.frame()s to be very significant (indicating multiple values from a single observation) so what you're doing seems like it may be "against the spirit of R", but if you want a simple NA padding at the end, this should do it: listToDF <- function(inputList, fill = N

Re: [R] Restricting R session

2012-01-11 Thread Antonio Rodriges
Thank you, Uwe, below are my comments >> In particular, how to prohibit some set of functions, for example, >> from "base" package? > > > You can't: R is free software. This does not imply it must be inflexible and unsuitable for cloud services > Well, of course you could build your own > versi

Re: [R] New to R, Curious about Project Idea

2012-01-11 Thread R. Michael Weylandt
It sounds quite possible but you'll probably get more specialized help if you ask on the r-sig-geo mailing list. Michael On Wed, Jan 11, 2012 at 9:36 AM, arbeaupg wrote: > Good morning, > > I am a student whom is currently working on a term project for my GIS > Program. I am looking for a softwa

Re: [R] stacked barplot colour coding

2012-01-11 Thread Richard M. Heiberger
This graph would be easier under lattice graphics. biomass <- data.frame(bg=c(0.41, 0.37, 0.31, 0.32), ag=c(2.81, 2.91, 2.06, 2.39)) b2 <- stack(biomass) names(b2) <- c("mass", "type") b2$type <- factor(b2$type, levels=c("bg","ag")) b2$AB <- rep(c("A","A","B","B"), 2) b2$loca

Re: [R] Problem with segmented

2012-01-11 Thread Sarah Goslee
We really need the small reproducible example requested in the posting guide, including sample data, the actual R commands you used, the libraries required, and your OS and version of R. Sarah On Wed, Jan 11, 2012 at 9:08 AM, Filoche wrote: > Hi there. > > Here's the error message. > > Error in

[R] 2D filter in R?

2012-01-11 Thread Michael
Hi all, I am looking for a command for doing 2D filtering (rectangular or Gaussian) in R... I have looked at ksmooth, filter and convolve but they seem to be 1D... Any thoughts? Thanks a lot! [[alternative HTML version deleted]] __ R-help@r-

Re: [R] 2 sample wilcox.test != kruskal.test

2012-01-11 Thread Daniel Nordlund
> -Original Message- > From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org] > On Behalf Of syrvn > Sent: Wednesday, January 11, 2012 2:36 AM > To: r-help@r-project.org > Subject: Re: [R] 2 sample wilcox.test != kruskal.test > > Hi, > > thanks for your answer. Unfortuna

Re: [R] plotOHLC(alpha3): Error in plotOHLC(alpha3) : x is not a open/high/low/close time series

2012-01-11 Thread Joshua Ulrich
On Wed, Jan 11, 2012 at 11:10 AM, Ted Byers wrote: > Hi Joshua, > > Thanks. > > I had used irts because I thought I had to.  The tick data I have has some > minutes in which there is no data, and others when there are hundreds, or > even thousands.  If xts supports irregular data, the that is one

[R] Storing/Restoring R objects

2012-01-11 Thread Rich Shepard
One of my projects has generated quite a few objects (data frames) related to one portion of this project. They can be listed with the ls() function. What I would like to do is move them to another directory so that data frames for other portions of the project can be more easily seen and used

Re: [R] Restricting R session

2012-01-11 Thread Uwe Ligges
On 10.01.2012 20:30, Antonio Rodriges wrote: Hello, Is it possible to use R on public server where each user has its own restricted R session? This entirely depends on the definition of "restricted", otherwise the answer is "yes". In particular, how to prohibit some set of functions, fo

[R] Getting a "root edge" error when trying to read the phylocom megatree into R

2012-01-11 Thread Megan Bartlett
Hi, I'm trying to use the picante package in R to build phylogenetic trees, based on a list of taxa I have data for and the phylocom APG3 megatree (version R20091110; http://www.phylodiversity.net/phylomatic/). However, trying to read the tree in R using: read.tree("R20091110.new.txt") Gives the

[R] Confidence Interval from Moments?

2012-01-11 Thread lambdatau
Hi all, I'm wondering whether it is possible to construct a confidence interval using only the mean, variance, skewness and kurtosis, i.e. without any of the population? If anyone could help with this it'd be much appreciated (even if just a confirmation of it being impossible!). Thanks. -- Vie

Re: [R] general question on Spotfire

2012-01-11 Thread Spencer Graves
Roughly 5 years ago, a Spotfire rep at the Joint Statistical Meetings told me they routinely interfaced with both R and S-Plus. I'm not 100% certain, but I believe they have many customers who use that facility today. Spencer On 1/11/2012 10:37 AM, Bert Gunter wrote: Peter et. al: 1. I ag

[R] turning a list of vectors into a data.frame (as rows of the DF)?

2012-01-11 Thread Chris Conner
As a newer R practicioner, it seems I stump myself weekly (at least) with issues that have spinning my wheels.  Here is yet another... I'm trying to turn a list of numeric vectors (of uneual length) inot a dataframe.  Each vector held in the list represents a row, and there are some rows of uneq

Re: [R] rjags installation trouble

2012-01-11 Thread Martyn Plummer
It looks like you are cross-linking to an earlier version of the JAGS library at run time. Check /sbin/ldconfig -p | grep jags - When compiling rjags, you can hard-code the location of the jags library using the [GNU-specific] configure option --enable-rpath. Martyn On Tue, 2012-01-10 at 10:40

Re: [R] fix and edit don't work: unable to open X Input Method->segfault

2012-01-11 Thread Patrick Connolly
On Sun, 08-Jan-2012 at 03:32PM -0600, Paul Johnson wrote: |> I can't run fix() or edit() anymore. Did I break my system? |> |> I'm running Debian Linux with R-2.14.1. As far as I can tell, the R |> packages came from Debian's testing "wheezy" repository. I would |> like to know if users on othe

[R] meta-analysis normal quantile plot metafor

2012-01-11 Thread Ricc
Hello, I once used the metawin software to perform a meta-analysis (see metawinsoft, Rosenberg et al.) and produced normal qqplot to test for a potential bias in the dataset. I now want to re-use the same dataset with the package metafor by W. Viechtbauer (great package btw). I run the qqnorm.rma

[R] New to R, Curious about Project Idea

2012-01-11 Thread arbeaupg
Good morning, I am a student whom is currently working on a term project for my GIS Program. I am looking for a software package which can aid me in my project, and I was curious if R would be able to address my goals. My project includes power outage data from a hydro company (point data, with

Re: [R] Problem with segmented

2012-01-11 Thread Filoche
Hi there. Here's the error message. Error in seg.lm.fit(y, XREG, Z, PSI, weights, offs, opz) : (Some) estimated psi out of its range I have tried many ways to specify the arguments, but apparently the error message is related to the estimated break point being invalid. However, my estimation

Re: [R] Constructing a data.frame from csv files

2012-01-11 Thread Vincy Pyne
Dear Sir, Thanks a lot for your guidance. I have understood my mistake. It was naming the columns viz.   colnames(price_returns) = entity_returns which was creating the problems. Code is running excellently once I got rid of this particular line. I will use melt from reshape etc to get the requi

Re: [R] odfWeave UTF-8 error and latin characters

2012-01-11 Thread staffan7s
Hello, I am using R and Libreoffice on Ubuntu 11.10 (64-bit) and have been experiencing similar problems with character encoding (Swedish utf8) in odfWeave. Here is an example of what it looks like: Should be: "Hör Ärland dåligt?" Appears as: "Hör Ärland dÃ¥ligt?" I found a (pretty clumsy) sol

Re: [R] Accomplishing a loop on multiple columns

2012-01-11 Thread iliketurtles
Lists are the answer. LIST<-list() for(i in 1:ncol(results6)) { LIST[[i]]<-lm(results6[,i]~data$observed) } You'll now have a 91 entry list of lm(). You can then do something like this: LIST2<-list() for(i in 1:length(LIST)) { LIST2[[i]]<-LIST[[i]]$r.squared } This should now be a li

Re: [R] RGL- Drawing Circle

2012-01-11 Thread gli
thanks Duncan. I think want my circle to be in user coordinates. I tried the first code u gave me and it gives me a circle. but how can i: 1) change the radius ? 2) place the circle at a given x,y,z coordinate? 3) turn it 90 degree up like these circle plate bus stop? http://www.geocities.co.

Re: [R] Getting Sphericity Tests for Within Subject Repeated Measure Anova (using "car" package)

2012-01-11 Thread maximilian.mueller
Dear John, thanks for your help and sorry for answering this late. My question is a follow up question of an older thread posted several mongths ago, but your statement helped a lot. Thanks, Max - M a x i m i l i a n M ue l l e r PhD-Student Department of Business Studies Leuphana Univer

[R] Problem concerning withRestarts and R2WinBUGS

2012-01-11 Thread Jongerling, J. (Joran)
Dear R-users, I have a question regarding the withRestarts function in R. I´m running a simulation code in which I analyse data using both lme and R2WinBUGS. Now, I want to run this code for 1000 replications, however the model I´m using is a little ´sensitive´, so sometimes the WinBUGS analy

Re: [R] Max value of an integer

2012-01-11 Thread Mikko Korpela
On 01/11/2012 12:01 PM, Rui Esteves wrote: > Is there any constant that represents the maximum value of an integer? Yes, there is (assuming you refer to the 'integer' type). See ?.Machine. > .Machine$integer.max [1] 2147483647 > as.integer(2147483647) [1] 2147483647 > as.integer(2147483648) [1]

Re: [R] 64bit R under 32bit winxp

2012-01-11 Thread 孟欣
Thanks At 2012-01-11 16:55:32,"Jeff Newmiller" wrote: >You cannot install 64-bit R on 32-bit OS, but you can install a 32-bit R on a >64-bit OS, and you can later install 64-bit R as well. That is, installing >32-bit R does not interfere with your option to later install a 64-bit R. >--

[R] Accomplishing a loop on multiple columns

2012-01-11 Thread Nerak
Hello, I have a question concerning ‘for loops’ on multiple columns. I made 91 columns with results (all made together with a for loop) and I want to us lm to fit the model. I want to compare the results of all these calculated columns (91) with one column with observed values. I use the function

[R] Assist me on how I can arrange trend data of rainfall and temperature for analysis

2012-01-11 Thread MukhalaC
I am a student doing my MSc Research Methods, i am working on my thesis research on analysing and modelling of crop failure risks due to drought in selected districts in Malawi. The analysis and modelling will focus on two crop stages of development: just after planting and flowering stages. I have

Re: [R] Plot maps with R

2012-01-11 Thread hwright
Dear Alex Two other packages that create maps are: maps mapproj alaios wrote > > Dear all I would like to use R and make some maps. > I want to have strict control, over the details of the produced map, like > remove borders, city names, add markers, add labels. > > Is there any package apart

Re: [R] plotOHLC(alpha3): Error in plotOHLC(alpha3) : x is not a open/high/low/close time series

2012-01-11 Thread Ted Byers
Hi Joshua, Thanks. I had used irts because I thought I had to. The tick data I have has some minutes in which there is no data, and others when there are hundreds, or even thousands. If xts supports irregular data, the that is one less step for me to worry about. Alas, your suggestion didn't h

Re: [R] general question on Spotfire

2012-01-11 Thread Louis Bajuk-Yorgan
Hello, I am a Product Manager at Spotfire, focused on integrating statistical capabilities from R & S+ into Spotfire, so I will make a few comments: 1. We have a quite a few customers who use Spotfire and R side-by-side for doing ad hoc data analysis. Sometimes by the same user, sometimes by di

Re: [R] stacked barplot colour coding

2012-01-11 Thread Jean V Adams
acacia21 wrote on 01/09/2012 07:01:28 PM: > Hi all, > i'm fairly new to R and its graphing, but having unsuccessfully 'googled' > and checked this forum to find answer to my problem, i'm posting my question > here. > > I'm trying to plot stacked barplot. I have simple data that looks like this

Re: [R] general question on Spotfire

2012-01-11 Thread Bert Gunter
Peter et. al: 1. I agree with Duncan: wrong list. 2. AFAIK, Spotfire **already** can interface with R. -- Bert On Wed, Jan 11, 2012 at 8:17 AM, peter dalgaard wrote: > > On Jan 11, 2012, at 16:28 , Duncan Murdoch wrote: > >> On 12-01-11 10:13 AM, John Smith wrote: >>> Dear R users, >>> >>> I h

Re: [R] 2 sample wilcox.test != kruskal.test

2012-01-11 Thread Nordlund, Dan (DSHS/RDA)
> -Original Message- > From: r-help-boun...@r-project.org [mailto:r-help-bounces@r- > project.org] On Behalf Of syrvn > Sent: Tuesday, January 10, 2012 10:28 AM > To: r-help@r-project.org > Subject: [R] 2 sample wilcox.test != kruskal.test > > Hello, > > > I think I am right in saying th

Re: [R] Generating unque patient IDs

2012-01-11 Thread Chris Campbell
Dear Ayyappa Unique identifiers can be created from numbers using factor. These are coded as integers in R which you could use to relabel your dataset. > x <- rep(16:18, each = 2) > x [1] 16 16 17 17 18 18 > y <- factor(x) > levels(

Re: [R] Odp: Generating unque patient IDs

2012-01-11 Thread Ayyappa Chaturvedula
Dear all, I am sorry if I misstated the problem. The roundig issue is with NONMEM software not with R. But the suggestions are helpful. Regards,Ayyappa Chaturvedula On Jan 11, 2012, at 12:14 PM, David Winsemius wrote: Unfortunately the "rounding effect" (which I assumed was related to

Re: [R] Odp: Generating unque patient IDs

2012-01-11 Thread jim holtman
One of the reasons that I specified the 'seq' command as it was was to make sure it used numerics: > x <- seq(123456789012.0, length = 10, by = 1.0) > x [1] 123456789012 123456789013 123456789014 123456789015 123456789016 123456789017 123456789018 [8] 123456789019 123456789020 123456789021 > str

Re: [R] CairoPDF and greek letter spacing

2012-01-11 Thread Chris Campbell
As a workaround you could use escape characters, then adjust the font style as necessary. cairo_pdf(file = "zend.pdf") print(xyplot(y ~ x, data = data.frame(x = 1:10, y = 1:10), main = "Length (\u03BCm)")) dev.off() Regards Chris Campbell MANGO SOLUTIONS Data Analysis that Delivers +44 124

[R] rgl/ x11 problem

2012-01-11 Thread jack hietpas
Hello- I am having problems with plot3d... I keep receiving the follow messages: when I attempt to load the package: library(rgl, pos=4), I get this error message: [6] WARNING: Warning in rgl.init(initValue) : Warning in rgl.init(initValue) : Warning in rgl.init(initValue) : Warning in rgl.init(i

Re: [R] Odp: Generating unque patient IDs

2012-01-11 Thread David Winsemius
Unfortunately the "rounding effect" (which I assumed was related to the automatic conversion from integer to numeric) is only going to show up above 2147483647L, so I question whether you really demonstrated a solution to what I understood was the fundamental problem. -- David. On Jan 11

Re: [R] Help with speed (replacing the loop?)

2012-01-11 Thread Dimitri Liakhovitski
Thanks a lot, Steve! match sounds very promising - that means I only need a loop across predictors. As far as "get more memory" advice is concerned: I already have "more" memory :) On Wed, Jan 11, 2012 at 11:14 AM, Steve Lianoglou wrote: > Hi, > > On Wed, Jan 11, 2012 at 10:50 AM, Dimitri Liakho

Re: [R] Generating unque patient IDs

2012-01-11 Thread David Winsemius
On Jan 11, 2012, at 11:12 AM, Ayyappa Chaturvedula wrote: Dear group, I am trying to prepare a NONMEM friendly dataset for population PK analysis. My patient IDs are 10 digit long and NONMEM is losing precison and rouding the last couple of digits. Are you sure? I need to generate unique

Re: [R] Odp: Generating unque patient IDs

2012-01-11 Thread jim holtman
Does this do it for you: > sprintf("%010.0f", seq(10.0, length = 250, by = 1.0)) [1] "10" "11" "12" "13" "14" "15" "16" [8] "17" "18" "19" "100010" "100011" "100012" "100013" [15] "1000

[R] test for condition during whole r session or after each command

2012-01-11 Thread Jannis
Dear List, is there any way to test for certain conditions during the whole r session or after the execution of each command? I am debugging my code and sometimes a certain logical error causes a program error much later in the script/function so especially with loops etc it is hard to backt

[R] Odp: Generating unque patient IDs

2012-01-11 Thread Petr PIKAL
Hi > > Dear group, > > I am trying to prepare a NONMEM friendly dataset for population PK > analysis. My patient IDs are 10 digit long and NONMEM is losing precison > and rouding the last couple of digits. I need to generate unique Patient > IDs fromt he current 10-digit IDs. Ihave total 250 s

Re: [R] Help with speed (replacing the loop?)

2012-01-11 Thread Jeff Newmiller
It is common that performance problems are addressed by using more memory. If your algorithm needs to join those tables and do calculations, then you can either pay the piper in memory (usually the most appropriate answer) or you can reinvent those optimized algorithms in a compiled language and

Re: [R] general question on Spotfire

2012-01-11 Thread peter dalgaard
On Jan 11, 2012, at 16:28 , Duncan Murdoch wrote: > On 12-01-11 10:13 AM, John Smith wrote: >> Dear R users, >> >> I have been using R for 10 years, and I love it very much. But in my daily >> job for drug discovery, some people use Spotfire. I tried Spotfire on >> couple of data sets. It sounds

Re: [R] Help with speed (replacing the loop?)

2012-01-11 Thread Steve Lianoglou
Hi, On Wed, Jan 11, 2012 at 10:50 AM, Dimitri Liakhovitski wrote: > Thanks a lot, Steve. I have one question (below): > > >> >> library(data.table) >> >> ## your data >> xx <- data.table(group=c(rep("group1",5),rep("group2",5)), >>                 a=1:10, b=seq(10,100,by=10), key="group") >>

[R] Generating unque patient IDs

2012-01-11 Thread Ayyappa Chaturvedula
Dear group, I am trying to prepare a NONMEM friendly dataset for population PK analysis. My patient IDs are 10 digit long and NONMEM is losing precison and rouding the last couple of digits. I need to generate unique Patient IDs fromt he current 10-digit IDs. Ihave total 250 subjects so I apprec

Re: [R] general question on Spotfire

2012-01-11 Thread John Smith
I am struggling whether I should learn Spotfire or not. I just want some statisticians inputs. Thanks On Wed, Jan 11, 2012 at 10:28 AM, Duncan Murdoch wrote: > On 12-01-11 10:13 AM, John Smith wrote: > >> Dear R users, >> >> I have been using R for 10 years, and I love it very much. But in my

Re: [R] CairoPDF and greek letter spacing

2012-01-11 Thread Walker, John Stephen
Thanks for the input, but I can confirm the either form of the paste command i.e. expression(paste("Length (", mu*m, ")")) or my original which had expression(paste("Length (", mu, "m)")), and the expression command alone have the same effect. It may be limited to Linux machines (which I did st

Re: [R] get the percentage rank of a value based on an empirical data vector

2012-01-11 Thread Martin Batholdy
> > findInterval(6.2, sort(x)) > [1] 704 > > xecdf <- ecdf(x) > > xecdf(6.2) > [1] 0.704 thanks, that helped a lot! On 11.01.2012, at 14:58, David Winsemius wrote: > > On Jan 11, 2012, at 8:12 AM, Martin Batholdy wrote: > >> Hi, >> >> I have a vector with values: >> >> x <- rnorm(1000, 5

Re: [R] Help with speed (replacing the loop?)

2012-01-11 Thread Dimitri Liakhovitski
Thanks a lot, Steve. I have one question (below): > > library(data.table) > > ## your data > xx <- data.table(group=c(rep("group1",5),rep("group2",5)), >                 a=1:10, b=seq(10,100,by=10), key="group") > yy <- data.table(group=c("group1","group2"), a=c(10,20), b=c(2,3), >          

Re: [R] plotOHLC(alpha3): Error in plotOHLC(alpha3) : x is not a open/high/low/close time series

2012-01-11 Thread Joshua Ulrich
Hi Ted, On Tue, Jan 10, 2012 at 1:59 PM, Ted Byers wrote: > R version 2.12.0, 64 bit on Windows. > > Here is a short script that illustrates the problem: > library(tseries) > library(xts) > setwd('C:\\cygwin\\home\\Ted\\New.Task\\NKs-01-08-12\\NKs\\tests') > x = read.table("quotes_h.2.dat", heade

Re: [R] general question on Spotfire

2012-01-11 Thread Duncan Murdoch
On 12-01-11 10:13 AM, John Smith wrote: Dear R users, I have been using R for 10 years, and I love it very much. But in my daily job for drug discovery, some people use Spotfire. I tried Spotfire on couple of data sets. It sounds I still need do some data manipulation before plot figures. For ex

Re: [R] Help with speed (replacing the loop?)

2012-01-11 Thread Steve Lianoglou
Hi, On Wed, Jan 11, 2012 at 9:57 AM, Dimitri Liakhovitski wrote: > Dear R-ers, > > I have a loop below that loops through my numeric variables in data > frame x and through levels of the factor "group" and multiplies (group > by group) the values of numeric variables in x by the corresponding > g

[R] general question on Spotfire

2012-01-11 Thread John Smith
Dear R users, I have been using R for 10 years, and I love it very much. But in my daily job for drug discovery, some people use Spotfire. I tried Spotfire on couple of data sets. It sounds I still need do some data manipulation before plot figures. For example, I can not plot figure with data arr

Re: [R] get the percentage rank of a value based on an empirical data vector

2012-01-11 Thread R. Michael Weylandt
If performance is an issue, I think mean(x < y) will be as quick as it can be done in R alone (you could do it in C in a single pass if needed which might be a good first exercise in using compiled code) Michael On Jan 11, 2012, at 8:58 AM, David Winsemius wrote: > > On Jan 11, 2012, at 8:12

  1   2   >