Re: [R] barplot with more than 1 variable

2012-02-19 Thread Jim Lemon
On 02/19/2012 11:31 PM, Francesco Sarracino wrote: Dear R listers, I am trying to produce a simple (for a stata user) barplot with 4 countries on the x axis, each country observed in 2 subsequent years and 3 variables. Basically, I should have three bars for each year for each country. I am attac

Re: [R] random numbers

2012-02-19 Thread Petr Savicky
On Sun, Feb 19, 2012 at 08:39:13PM -0800, Ryan Murphy wrote: > Hello, > > Is there a way to create random numbers that fit a certain specified > requirement other than distributional characteristics. > In particular, I would like to create simulated income distributions with > certain gini coeffic

Re: [R] R

2012-02-19 Thread Daniel Nordlund
> -Original Message- > From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org] > On Behalf Of Steve Lianoglou > Sent: Sunday, February 19, 2012 8:08 PM > To: elliot.we...@virgin.net > Cc: r help > Subject: Re: [R] R > > Hi, > > When replying to mails from r-help, please h

[R] random numbers

2012-02-19 Thread Ryan Murphy
Hello, Is there a way to create random numbers that fit a certain specified requirement other than distributional characteristics. In particular, I would like to create simulated income distributions with certain gini coefficient Ryan -- Ryan Murphy 2012 B.A. Economics and Mathematics

Re: [R] R

2012-02-19 Thread Steve Lianoglou
Hi, When replying to mails from r-help, please hit "reply all" so that people can follow along so that everyone can benefit from the help, as well as getting more eyeballs on your questions so you get better help. Now: On Sun, Feb 19, 2012 at 8:54 PM, wrote: > Brilliant, thanks for the help. D

[R] repeating or looping within an apply statement to handle multiple variables

2012-02-19 Thread Matt Spitzer
Dear R experts, I would like to please ask for your help with repeating steps in an apply statement. I have a dataframe that lists multiple variables for a given id and visit, as well as drug treatment. > head(exp) id visit variable1 variable2 variable3 variable4 drug 1 3 113

Re: [R] help updating package rJava (on ubuntu)

2012-02-19 Thread Simon Urbanek
On Feb 19, 2012, at 4:31 PM, Karl Brand wrote: > Hi Hasan, > > Success. For myself and FWIW to other useR's here's how i spent the sunny > half of my sunday to achieve it :/ > > Many thanks for your and Simon's input, > FWIW you should not need to set any custom settings if you system is prop

[R] Reporting Kaplan-Meier / Cox-Proportional Hazard Standard Error, km.coxph.plot, survfit.object

2012-02-19 Thread Paul Johnston
What is the best way to report the standard error when publishing Kaplan-Meier plots? In my field (Vascular Surgery), practitioners loosely refer to the "10% error" cutoff as the point at which to stop drawing the KM curve. I am interpreting this as the *standard error of the cumulative hazard*,

Re: [R] svm with GRASS GIS

2012-02-19 Thread etiennebr
Look at ?predict.svm, you'll see that you need to provide a Matrix, not a data.frame. Etienne -- View this message in context: http://r.789695.n4.nabble.com/svm-with-GRASS-GIS-tp4388006p4403005.html Sent from the R help mailing list archive at Nabble.com. ___

Re: [R] R

2012-02-19 Thread Steve Lianoglou
Hi, On Sun, Feb 19, 2012 at 7:47 PM, wrote: > Hello, just a couple of short questions that would be much appreciated. > > Is there a way of putting a vectors inputs in size order from low to high? See: ?order, ?sort > And if I had a random list or TRUE and FALSE, is there a way of finding the

Re: [R] installing the package Rcplex

2012-02-19 Thread zheng wei
I did not know this before. I installed it as you suggested. what to do next?   Best, Wei From: David Winsemius Cc: "r-help@r-project.org" Sent: Sunday, February 19, 2012 6:44 PM Subject: Re: [R] installing the package Rcplex Just a guess (based on eviden

Re: [R] time-series trend classification

2012-02-19 Thread John Kohr
Thanks for your responses and sorry if it wasn't clear what I am trying to do. Yes the points are 5-10-15 no more than these and let's say we are interested to capture only the increase or decrease (are more classes than these but let's start simple). Applying linear regression and comparing th

[R] R

2012-02-19 Thread elliot.welch
Hello, just a couple of short questions that would be much appreciated. Is there a way of putting a vectors inputs in size order from low to high? And if I had a random list or TRUE and FALSE, is there a way of finding the 100th TRUE? Thank you very much, Elliot Welch elliot.we...@virgin.net

Re: [R] installing the package Rcplex

2012-02-19 Thread David Winsemius
On Feb 19, 2012, at 7:45 PM, zheng wei wrote: > I did not know this before. I installed it as you suggested. what to > do next? Read the Installation Manual? -- David. > > Best, > Wei > > From: David Winsemius > To: zheng wei > Cc: "r-help@r-project.org" > Sent: Sunday, February 19, 2012

Re: [R] installing the package Rcplex

2012-02-19 Thread David Winsemius
Just a guess (based on evidence of failing to read other instructions that were made available), but have you perhaps forgotten to install Rtools for Windows? http://cran.r-project.org/bin/windows/Rtools/ -- David. On Feb 19, 2012, at 5:13 PM, zheng wei wrote: Dear All, I have difficulty

[R] installing the package Rcplex

2012-02-19 Thread zheng wei
Dear All,   I have difficulty in installing the package of Rcplex. Its description web is at http://cran.r-project.org/web/packages/Rcplex/index.html   First, I tried -- url <- "http://cran.r-project.

[R] mixed anova models with aov

2012-02-19 Thread Nicolas Schuck
Dear R Users, I am trying to calculate a mixed effect Anova with 2 nested between subject factors and two nested within subject factors. Specifically, I have data from two Agegroups (between subject variable old vs young) that are characterized by their genotype (between subjects, NESTED withi

Re: [R] time-series trend classification

2012-02-19 Thread Bert Gunter
If I understand you correctly, 5-10 points is not near enough to fit anything beyond a simple linear trend unless you know a priori what the form of the "trend" is. Period. Admitting that the data are insufficient is better than performing misleading analyses. * Cheers, Bert * My personal view.

Re: [R] time-series trend classification

2012-02-19 Thread Roy Mendelssohn
On Feb 19, 2012, at 9:46 AM, John Kohr wrote: > > Hello everyone, > > I was looking for a way to classify time-series based on the curve-fit. I try > to campute several trends so i was thinking to link each trend with a > function. increase with exponential for example, increase and decrease

Re: [R] help updating package rJava (on ubuntu)

2012-02-19 Thread Karl Brand
Hi Hasan, Success. For myself and FWIW to other useR's here's how i spent the sunny half of my sunday to achieve it :/ Many thanks for your and Simon's input, Karl Since: $ javac -version returned nothing i believe you (and Simon) were right, i.e, it (and JDK) were missing on my system. F

Re: [R] parse .ps files with R?

2012-02-19 Thread Nick Matzke
Thanks very much!! Nick On 2/18/12 7:39 PM, baptiste auguie wrote: Hi, the grImport package provides some tools for this. HTH, b. On 19 February 2012 16:06, Nick Matzke wrote: Hi, Is there a way to parse a postscript (*.ps) file with R (or perhaps with some other command-line utility)? E

Re: [R] R help

2012-02-19 Thread Rolf Turner
On 20/02/12 09:55, Bert Gunter wrote: Folks: Perhaps I am naive, ignorant, or foolish, but this whole discussion seems rather ridiculous. What possible relation to reality could a multivariate normal of the size requested have? Even for simulation, it seems like nonsense. Right on, bro'!!!

Re: [R] Non-parametric test for repeated measures and post-hoc single comparisons in R?

2012-02-19 Thread peter dalgaard
Repeated measures means that you have multiple measurements on the same individual. Usually, the same person measured at different time points. So if you have N individuals and T times, then you can place your observations in an N*T layout. In this layout, you can have 1 observation per cell o

Re: [R] Basic Model Setup Question from a Beginner

2012-02-19 Thread Bert Gunter
You need statistical consultation. The R list is not meant for this. Any advice you receive here -- or indeed, anywhere online-- is likely to be flawed because it is difficult to communicate essential aspects of the issues remotely. I would strongly recommend that you seek local statistical help.

Re: [R] R help

2012-02-19 Thread Bert Gunter
Folks: Perhaps I am naive, ignorant, or foolish, but this whole discussion seems rather ridiculous. What possible relation to reality could a multivariate normal of the size requested have? Even for simulation, it seems like nonsense. Cheers, Bert On Sun, Feb 19, 2012 at 11:35 AM, Petr Savicky

[R] GTK+

2012-02-19 Thread Idielle Walters
Hi I am struggling to install GTK+ for Windows 7. RGtk2 needs this package to load. Does anybody know of a installation file that works? __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide

[R] coxme: model simplification using LR-test?

2012-02-19 Thread Simon Tragust
Hi I'm encountering some problems with coxme My data: I'm looking at the survival of animals in an experiment with 3 treatments, which came from 4 different populations, two of which were infected with a parasite and two of which were not. I'm interested if infected animals differe from uninfect

[R] time-series trend classification

2012-02-19 Thread John Kohr
Hello everyone, I was looking for a way to classify time-series based on the curve-fit. I try to campute several trends so i was thinking to link each trend with a function. increase with exponential for example, increase and decrease with a gaussian etc. The possiblities are endless though an

Re: [R] Non-parametric test for repeated measures and post-hoc single comparisons in R?

2012-02-19 Thread Tal Galili
The following post would not answer your question at full, but might give some good code/information: http://www.r-statistics.com/2010/02/post-hoc-analysis-for-friedmans-test-r-code/ Contact Details:--- Contact me: tal.gal...@g

[R] Basic Model Setup Question from a Beginner

2012-02-19 Thread All Day
Hello all! I would like to start off by saying that I am still really new to the vast world of R so please excuse my very limited vocabulary in the program. I have collected data from monkey videos and would like to setup some model(s) in R that can help with my hypotheses. I am having trouble fi

Re: [R] R help

2012-02-19 Thread Petr Savicky
On Sat, Feb 18, 2012 at 06:00:53PM -0500, li li wrote: > Dear all, > I need to generate numbers from multivariate normal with large dimensions > (5,000,000). Hi. I am replying to your first email, since the other did not arrive to my folder, possibly filtered out by a spam filter. I see them at

Re: [R] barplot with more than 1 variable

2012-02-19 Thread ilai
You could have found the solution in http://had.co.nz/ggplot2/geom_bar.html yourself since all help pages for ggplot refer you to the web site. But to speed things up for you, try this for starters: TUSE2 <- data.frame(country = rep(c("United States","Italy","Germany","Netherlands"),each=2), year3

Re: [R] Non-parametric test for repeated measures and post-hoc single comparisons in R?

2012-02-19 Thread saschaview
On 19.02.12 19:31, Tal Galili wrote: The following post would not answer your question at full, but might give some good code/information: http://www.r-statistics.com/2010/02/post-hoc-analysis-for-friedmans-test-r-code/ Thank you, Tal! As you already mentioned, your interesting post supports

[R] Non-parametric test for repeated measures and post-hoc single comparisons in R?

2012-02-19 Thread saschaview
Some attribute x from 17 individuals was recorded repeatedly on 6 time points using a Likert scale with 7 distractors. Which statistical test(s) can I apply to check whether the changes along the 6 time points were significant? set.seed( 123 ) x <- matrix( sample( 1:7, 17*6, repl=T ), nrow =

Re: [R] R help

2012-02-19 Thread David Winsemius
On Feb 19, 2012, at 11:46 AM, li li wrote: Actually I still get an error for the case of equal correlation. No. You need to read the error message for meaning. Below is the code m <- 10 n <- 500 m1 <- 0.5*m mu <- c(rep(2, m1), rep(0, m-m1)) rho <- 0.5 x.mod1 <- matrix(rnorm(n*m, sd=sq

Re: [R] R help

2012-02-19 Thread peter dalgaard
On Feb 19, 2012, at 16:53 , li li wrote: > Petr, > Thanks for the help. That certainly makes sense and solves my current > problem. But, in general, for arbitrary covariance matrix (instead of the > special equi-correlation case), it there a way to generate numbers from > multivariate normal wh

Re: [R] World Heat/Rainbow Map

2012-02-19 Thread ilai
Figure 13.10 in "Lattice: Multivariate Data Visualization with R" might give you some ideas http://lmdvr.r-forge.r-project.org/figures/figures.html On Sun, Feb 19, 2012 at 9:05 AM, Lorenzo Isella wrote: > Dear All, > I would like to do the following: make a plot of the world and color a few > se

Re: [R] R help

2012-02-19 Thread li li
Actually I still get an error for the case of equal correlation. Below is the code > m <- 10 > n <- 500 > m1 <- 0.5*m > mu <- c(rep(2, m1), rep(0, m-m1)) > rho <- 0.5 > x.mod1 <- matrix(rnorm(n*m, sd=sqrt(1-rho)), nrow=n, ncol=m)+rnorm(n, sd=sqrt(rho))+t(replicate(n,mu)) Error: cannot allocate

[R] World Heat/Rainbow Map

2012-02-19 Thread Lorenzo Isella
Dear All, I would like to do the following: make a plot of the world and color a few selected states. Some states have an associated scalar (i.e a simple number) and I would like to paint them using a heat map and show legend for the color map in the plot. One or two states do not have any numbe

Re: [R] Multiple planes in a scatterplot3d

2012-02-19 Thread Uwe Ligges
On 18.02.2012 06:42, kpiratla wrote: I am relatively new to R and scatterplot3d function. I need to draw a 3d scatterplot and then add three different planes (each parallel to xy, yz and zx with varying intercepts) to a 3d scatterplot. I got the plot running and i was also able to produce one p

Re: [R] R help

2012-02-19 Thread li li
Petr, Thanks for the help. That certainly makes sense and solves my current problem. But, in general, for arbitrary covariance matrix (instead of the special equi-correlation case), it there a way to generate numbers from multivariate normal when the dimension is very high? Thanks.

Re: [R] how to delete all classes & variables in current & previously saved history?

2012-02-19 Thread Bert Gunter
On Sun, Feb 19, 2012 at 5:33 AM, Jim Holtman wrote: > the real question is why are you reloading the former session?  I always > start with a clean workspace and build what I need. Because that's the default? See ?Startup and especially the part about "--no-restore" (sans quotes) for how to

Re: [R] Plot OctTree

2012-02-19 Thread Uwe Ligges
On 18.02.2012 23:34, Jaimin Dave wrote: Hi Everyone, I have csv file which is in following format xmin,xmax,ymin,ymax,zmin,zmax I want to plot 3d graph with all the all octants being displayed as well. Any idea? I have used scatterplot3d package but it does not seem to have anything by whic

Re: [R] JGR Segmentation fault

2012-02-19 Thread Simon Urbanek
On Feb 19, 2012, at 2:48 AM, dberringer wrote: > I have been trying to get JGR up and running but so far no luck. I installed > R from source using the homebrew package manager and everything works fine > and I was able to install JGR and all its dependencies. Every time I try to > start JGR I

Re: [R] changing names with different character sets

2012-02-19 Thread Prof Brian Ripley
On 19/02/2012 12:43, peter dalgaard wrote: On Feb 19, 2012, at 08:49 , Prof Brian Ripley wrote: On 19/02/2012 07:30, Erin Hodgess wrote: Dear R People: I'm trying to replicate something that I saw on an R blog. The first step is to load in the .rda file, which is fine. However, some of the

Re: [R] how to delete all classes & variables in current & previously saved history?

2012-02-19 Thread Jim Holtman
the real question is why are you reloading the former session? I always start with a clean workspace and build what I need. Sent from my iPad On Feb 19, 2012, at 6:51, Alfredo Alessandrini wrote: > rm(list=ls()) > > 2012/2/19 sagarnikam123 > >> i have variables stored in previous history,n

Re: [R] how to delete all classes & variables in current & previously saved history?

2012-02-19 Thread Alfredo Alessandrini
rm(list=ls()) 2012/2/19 sagarnikam123 > i have variables stored in previous history,next time i log in,& load > previous history,i want to use same variable name for using other values, > & i don't know exactly how many are they ? & their names also.. > > How to delete all variable created in cu

Re: [R] changing names with different character sets

2012-02-19 Thread peter dalgaard
On Feb 19, 2012, at 08:49 , Prof Brian Ripley wrote: > On 19/02/2012 07:30, Erin Hodgess wrote: >> Dear R People: >> >> I'm trying to replicate something that I saw on an R blog. >> >> The first step is to load in the .rda file, which is fine. >> >> However, some of the names of the columns in

Re: [R] read.csv error: invalid multibyte string

2012-02-19 Thread Milan Bouchet-Valat
Le samedi 31 décembre 2011 à 07:05 -0800, Dennis Fisher a écrit : > R version: 2.13.1 > OS X > > Colleagues, > > I am working with a CSV file; for testing purposes, I created an XLS version > of the file. > When I read these files using read.xls (gdata) or read.csv, I encounter an > error: >

[R] how to delete all classes & variables in current & previously saved history?

2012-02-19 Thread sagarnikam123
i have variables stored in previous history,next time i log in,& load previous history,i want to use same variable name for using other values, & i don't know exactly how many are they ? & their names also.. How to delete all variable created in current & previous session i.e. (known & unknown) wh

[R] how to identify which value is missing from array/variable

2012-02-19 Thread sagarnikam123
i have variable C as > > c [1] 1 2 3 6 7 8 12 13 14 15 16 17 18 want to check variable sequentially for missing values/integers & reporting their positions, latter on,i want to fill with appropriate /missed values increasingly. -- View this message in context: http://r.7

Re: [R] R help

2012-02-19 Thread Petr Savicky
On Sat, Feb 18, 2012 at 06:00:53PM -0500, li li wrote: > Dear all, > I need to generate numbers from multivariate normal with large dimensions > (5,000,000). > Below is my code and the error I got from R. Sigma in the code is the > covariance > matrix. Can anyone give some idea on how to take car

Re: [R] > 14 points for lines(y~x)

2012-02-19 Thread Jim Lemon
dardisco wrote: > Problem: lines plots 1st 14x points then stops. > > boxplot(y1~lx) > lines(y1~lx) > > Note tilda appearing as -(minus) on this thread. > Variables: y1, lx; both=continuous > Goal: display boxplots for y1 for each value of lx, then join with ine > (Loess are nice, any will do). >

Re: [R] > 14 points for lines(y~x)

2012-02-19 Thread Jeff Newmiller
Time for you to read the posting guide mentioned at the bottom of every R-help message. Your assertions make no sense, so we will need a reproducible example to understand you. --- Jeff NewmillerThe

Re: [R] changing names with different character sets

2012-02-19 Thread Jeff Newmiller
Refer to the columns by their position (numerically). --- Jeff NewmillerThe . . Go Live... DCN:Basics: ##.#. ##.#. Live Go... Live:

[R] > 14 points for lines(y~x)

2012-02-19 Thread dardisco
Problem: lines plots 1st 14x points then stops. boxplot(y1~lx) lines(y1~lx) Note tilda appearing as -(minus) on this thread. Variables: y1, lx; both=continuous Goal: display boxplots for y1 for each value of lx, then join with line (Loess are nice, any will do). Have tried with other versions o

[R] JGR Segmentation fault

2012-02-19 Thread dberringer
I have been trying to get JGR up and running but so far no luck. I installed R from source using the homebrew package manager and everything works fine and I was able to install JGR and all its dependencies. Every time I try to start JGR I get the error below. Any idea why I'm getting the Segmen

Re: [R] Counting value changes

2012-02-19 Thread peter dalgaard
On Feb 19, 2012, at 04:25 , jim holtman wrote: > For completeness, if you want to count all possible four transitions: > >> x <- c(0,1,0,1,0,0,0,1,1,1,0,0,0,1) >> # lets keep count of the 4 different transitions that can happen >> indx <- cbind(head(x, -1), tail(x, -1)) %*% c(2, 1) >> table(indx