[R] Box plot without original data

2014-04-04 Thread Mª Teresa Martinez Soriano
Hi to everyone! I need to plot a box plot but I don't have the original data, is it possible to make it just with the median, mean, sd and range of values? Any idea is welcome. Thanks in advance. __

[R] Simulating data

2014-04-04 Thread Mª Teresa Martinez Soriano
Hi to everyone To simulate data I only know this command: rnorm(n , mean, sd)   It exists another way to specify Median and range as well? The point is that I have  this information about the variable: Median: 4.3 Mean: 4.2 SD: 1.8 Range: 0-8 and I need a boxplot, but I don't have the

Re: [R] Box plot without original data

2014-04-04 Thread Jim Lemon
On 04/04/2014 05:36 PM, Mª Teresa Martinez Soriano wrote: Hi to everyone! I need to plot a box plot but I don't have the original data, is it possible to make it just with the median, mean, sd and range of values? Any idea is welcome. Hi Teresa, You can get the return value from the boxplot

Re: [R] Equation of a curve

2014-04-04 Thread Keith Jewell
On 03/04/2014 16:26, Frances Cheesman wrote: Hi all, I have a number of bacterial growth curves I would like to find the equations for these and then integrate them to find the area under the curves for me to do stats on later. Is there any way I can do this in R? Thanks, Frances

[R] [R-pkgs] rDVR package

2014-04-04 Thread John Harrison
I would like to announce the release of version 0.1.1 of the rDVR package on CRAN. rDVR is an R package that allows cross Platform (Lin/Win/OSx) video recording from R. It is light weight and allows up to 10 minutes of recording by default. rDVR has a project page at

Re: [R] Equation of a curve

2014-04-04 Thread Frances Cheesman
Thanks everyone for al your help, I don't think it's necessarily as easy as I first thought. I'm going to have a think about it and try some things out. And I'll be back if I get stuck! Thanks very much, Frances On 4 April 2014 09:11, Keith Jewell keith.jew...@campdenbri.co.uk wrote: On

Re: [R] Simulating data

2014-04-04 Thread Frede Aakmann Tøgersen
Hi Teresa Try this: Median - 4.3 Mean - 4.2 SD - 1.8 RangeLower - 0 RangeUpper - 8 par(mfcol = c(1,2)) (bpstats - boxplot(rnorm(100, mean = Mean, sd = SD), outl = FALSE)) ## assuming normality of original data since we don't know the distribution ## get the 25% and 75% quantile q25 -

Re: [R] MANOVA post hoc testing

2014-04-04 Thread John Fox
Dear Nicholas, On Fri, 4 Apr 2014 04:59:09 + nicholas.suraw...@csiro.au wrote: Greetings, I'm interested in performing some post hoc tests after conducting a multivariate analysis of covariance (MANCOVA) which I performed using the Anova function in the car package. The covariate did

[R] simulation data

2014-04-04 Thread thanoon younis
dear sir i want to simulate multivariate ordinal data matrix with categories (1,4) and n=1000 and p=10. thanks alot thanoon [[alternative HTML version deleted]] __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help

[R] How to plot data using ggplot

2014-04-04 Thread drunkenphd
Hi, I have a list of cities and their coordinates, and also for each city I have a variable varA which I want to represent in a map using ggplot. For example : CityA lat 22.93977 lon 46.70663varA 545 CityB lat 23.93977 lon 46.70663varA 122 VarA values begin from 0 to 3000. I

[R] Lognormal AR(0,1) model

2014-04-04 Thread Chris89
Hi everyone! I am trying to make two log-normal AR(0,1) model using R with a given correlation between them, \rho, on the form: X_t = \alpha X_{t-1} + a_t Y_t = \beta Y_{t-1} + b_t At the moment I have been making n values of correlated log-normal data, called a_t and b_t, and generated a

[R] MICE, POST-PROCESSING imputations with two conditions

2014-04-04 Thread Leah Schinasi
Dear R listserve, I am writing with the hope or receiving assistance using the MICE package.My troubles relate to post-processing imputations. I am attempting to impute variables for which I want to specify two conditions using post-processing imputation: 1) if the value of a variable

[R] Mistakes in date conversion for future date/time (POSIXct)

2014-04-04 Thread Winkler, Matthias
Dear R-users, I'm working on datasets which contain data from the years 1960 to 2100 with a timestep of one hour. Every year has 365 days, leap years are ignored. After reading the dataset with R I convert the column which contains date/time to POSIXct: as.POSIXct(strptime(MyData [,1],

Re: [R] simulation data

2014-04-04 Thread Charles Determan Jr
Hi Thanoon, How about this? # replicate p=10 times random sampling n=1000 from a vector containing your ordinal categories (1,2,3,4) R - replicate(10, sample(as.vector(seq(4)), 1000, replace = T)) Cheers, Charles On Fri, Apr 4, 2014 at 7:10 AM, thanoon younis thanoon.youni...@gmail.comwrote:

Re: [R] Mistakes in date conversion for future date/time (POSIXct)

2014-04-04 Thread Rmh
high probability you are in a daylight savings time problem. see the archives for repair strategies. probably it will be enforcing standard time on all measurements. Rich Sent from my iPhone On Apr 4, 2014, at 10:55, Winkler, Matthias matthias.wink...@ibp.fraunhofer.de wrote: Dear

Re: [R] Mistakes in date conversion for future date/time (POSIXct)

2014-04-04 Thread Duncan Murdoch
On 04/04/2014 10:55 AM, Winkler, Matthias wrote: Dear R-users, I'm working on datasets which contain data from the years 1960 to 2100 with a timestep of one hour. Every year has 365 days, leap years are ignored. After reading the dataset with R I convert the column which contains date/time to

[R] Base within reverses column order

2014-04-04 Thread Dan Murphy
I just noticed this annoyance, but I'm not the first one, apparently -- see http://lists.r-forge.r-project.org/pipermail/datatable-help/2012-May/001176.html The thread never answered the OP's question Is this a bug? so I assume the answer, unfortunately, is No. If not a bug, do users of within

Re: [R] Base within reverses column order

2014-04-04 Thread Duncan Murdoch
On 04/04/2014 1:32 PM, Dan Murphy wrote: I just noticed this annoyance, but I'm not the first one, apparently -- see http://lists.r-forge.r-project.org/pipermail/datatable-help/2012-May/001176.html The thread never answered the OP's question Is this a bug? so I assume the answer,

[R] par(mfrow)

2014-04-04 Thread Emanuele Belli
Hi, I have some problems using the par function: I want to split the screen into 2 rows and 4 col and I type par(mfrow=c(2, 4)) but when I do that, instead of setting a graphical parameter, it creates a white Quarz. I'm currently using the R base version for Mac Os, 3.0.3 . Could you give me

[R] average of rows of each column

2014-04-04 Thread eliza botto
Dear useRs, I have a matrix of 120 row and 1000 columns.What I want is to get an average of a set of 12 rows starting from 1 till 120 for each column. Precisely, for column 1 the average of 1:10 rows, 11:20 rows 111:120. similarly for column 2, 3, 4 1000. So in the end i should have a

Re: [R] How to plot data using ggplot

2014-04-04 Thread Dennis Murphy
1. Look into the ggmap package if you want to overlay your data onto a map. 2. Re your color scale representation, define 'appropriately'. Do you mean a continuous range expressible in a colorbar or a discrete range, and if the latter, what intervals did you have in mind? Dennis On Fri, Apr 4,

[R] Revolutions blog: March roundup

2014-04-04 Thread David Smith
Revolution Analytics staff write about R every weekday at the Revolutions blog: http://blog.revolutionanalytics.com and every month I post a summary of articles from the previous month of particular interest to readers of r-help. In case you missed them, here are some articles related to R from

Re: [R] par(mfrow)

2014-04-04 Thread Franklin Bretschneider
Hi, re: I have some problems using the par function: I want to split the screen into 2 rows and 4 col and I type par(mfrow=c(2, 4)) but when I do that, instead of setting a graphical parameter, it creates a white Quarz. I'm currently using the R base version for Mac Os, 3.0.3 . Could

Re: [R] average of rows of each column

2014-04-04 Thread David Carlson
Something like (only 20 columns here): x - matrix(rnorm(120*20), 120, 20) xagg - aggregate(x, list(rep(1:12, each=10)), mean) - David L Carlson Department of Anthropology Texas AM University College Station, TX 77840-4352 Original Message- From:

[R] solicit help to read in 384 plate color image

2014-04-04 Thread Waverley @ Palo Alto
Hi, I am doing an experiment which results different colors of different intensities in the 384 micro titer plate. I took a picture of the plate by scanning in the image as a jpeg file and now I want to 1. read in the image file 2. grid the content 3. need to extract the intensity and color of

Re: [R] average of rows of each column

2014-04-04 Thread Rui Barradas
Hello, Try the following. m - 120 n - 10 # in your case this is 1000 mat - matrix(rnorm(n*m), nrow = m) fun - function(x, na.rm = TRUE){ tapply(x, rep(1:12, each = 10), mean, na.rm = na.rm) } apply(mat, 2, fun) apply(mat, 2, fun, na.rm = FALSE) # alternative Hope this helps, Rui

[R] how to change annotations in contour function from rsm package

2014-04-04 Thread John Smith
Dear All, I am using R3.0.3 and rsm 2.04. Since I want do log transformation on my data, so I want change annotation of Slice at ... For example, I want change (Slice at W = 1.25, L=2) in the first graph into (Slice at W = log(1.25), L=log(2)). Can anyone tell me how can I change it? Thanks

Re: [R] Mistakes in date conversion for future date/time (POSIXct)

2014-04-04 Thread David Winsemius
On Apr 4, 2014, at 9:54 AM, Duncan Murdoch wrote: On 04/04/2014 10:55 AM, Winkler, Matthias wrote: Dear R-users, I'm working on datasets which contain data from the years 1960 to 2100 with a timestep of one hour. Every year has 365 days, leap years are ignored. After reading the dataset

Re: [R] average of rows of each column

2014-04-04 Thread eliza botto
Dear David and Rui, Thankyou very much. :D Eliza From: dcarl...@tamu.edu To: eliza_bo...@hotmail.com; r-help@r-project.org Subject: RE: [R] average of rows of each column Date: Fri, 4 Apr 2014 15:14:32 -0500 Something like (only 20 columns here): x - matrix(rnorm(120*20), 120, 20) xagg

[R] Sweave files into LaTex

2014-04-04 Thread Axel Urbiz
Hi, I'm writing a thesis in Latex (say master.tex). I'd like to include R code/results from an .Rwd file. I've naively tried: 1) Add ONLY the code below in Rcode.Rnw file: \section{Exploratory data analysis} eval=TRUE, echo=FALSE= library(ggplot2) data(diamonds) head(diamonds) @ 2) Then, in

Re: [R] Sweave files into LaTex

2014-04-04 Thread Duncan Murdoch
On 04/04/2014, 7:10 PM, Axel Urbiz wrote: Hi, I'm writing a thesis in Latex (say master.tex). I'd like to include R code/results from an .Rwd file. I've naively tried: 1) Add ONLY the code below in Rcode.Rnw file: \section{Exploratory data analysis} eval=TRUE, echo=FALSE= library(ggplot2)

[R] rehh package for iHS and Rsb on NGS data

2014-04-04 Thread James Abugri
Hello expert community. Is there anybody whose used rehh package for iHS and rsb. I have gone through the tutorial a dozen times, I need expert help to enable me apply it to my p.falciparum snp data from the illumina platform. Thanks James Abugri -- * The information contained in this email

Re: [R] How to plot data using ggplot

2014-04-04 Thread drunkenphd
djmuseR wrote Do you mean a continuous range expressible in a colorbar or a discrete range, and if the latter, what intervals did you have in mind? I mean diskrete values like : 54508 25 101 420 95 928 24 1656 108 18 213 70. The problem is that range is from 0-4000

[R] Subsetting data by ID with different constraints

2014-04-04 Thread Lib Gray
Hello, I have a data set with many individuals all with multiple timed observations, and I would like to subset the data to exclude later timed observations. However, I would like to exclude different amounts of data for each individual. These individuals have two types of data: DV and dose. What

Re: [R] rehh package for iHS and Rsb on NGS data

2014-04-04 Thread David Winsemius
On Apr 4, 2014, at 11:46 PM, James Abugri wrote: Hello expert community. Is there anybody whose used rehh package for iHS and rsb. I have gone through the tutorial a dozen times, I need expert help to enable me apply it to my p.falciparum snp data from the illumina platform. This sounds

Re: [R] Removing White spaces with NA

2014-04-04 Thread arun
Hi, Check if this works: n - 1e7 dat - data.frame(Col1=c(A, , B,C, ,,D,rep('',n)), stringsAsFactors=FALSE) dat2 - dat dat[dat==''] - NA which(dat$Col1=='') #integer(0) #or dat2$Col1[dat2$Col1==''] - NA  which(dat2$Col1=='') #integer(0) A.K. On Wednesday, April 2, 2014 2:20 PM, arun

[R] Aggregate time series from daily to monthly by date and site

2014-04-04 Thread Zilefac Elvis
Hi, I have daily data arranged by date and site. Keeping the number of columns as there are, I will like to aggregate (FUN=mean) from daily to monthly the following data (only part is shown here) which starts in 1971 and ends in 1980.     Year Month Day Site Sim001 Sim002 Sim003 Sim004 1   1971

Re: [R] Aggregate time series from daily to monthly by date and site

2014-04-04 Thread Jeff Newmiller
You have been around long enough that we should not have to tell you how to provide data in a reproducible manner... read ?dput. --- Jeff NewmillerThe . . Go Live...

Re: [R] Aggregate time series from daily to monthly by date and site

2014-04-04 Thread Zilefac Elvis
Hi, I have daily data arranged by date and site. Keeping the number of columns as there are, I will like to aggregate (FUN=mean) from daily to monthly the following data (only part is shown here) which starts in 1971 and ends in 1980. structure(list(Year = c(1971, 1971, 1971, 1971, 1971,