[R] Box plot without original data

2014-04-03 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-help@r-pr

Re: [R] Doing box plots for part of the data??

2014-04-03 Thread Frede Aakmann Tøgersen
Or using the original data and the subset argument of boxplot: boxplot(Barium~Well.ID,data=mydata, main="Barium", ylab="mg/L", subset = Well.ID %in% c("MW-1", "MW-2")) Yours sincerely / Med venlig hilsen Frede Aakmann Tøgersen Specialist, M.Sc., Ph.D. Plant Performance & Modeling Technology &

Re: [R] NHPP with Cyclic Rate

2014-04-03 Thread Andrew Wachira
Yes. Infact i'm looking at the documentation right now - though there is not explicit mention of polynomial-trigonometric rate. I also found SAPP which comes closer to what i was looking for but only does polynomial rate. Thanks.. anymore recommendations will be highly appreciated. On 4 April 2

[R] MANOVA post hoc testing

2014-04-03 Thread Nicholas.Surawski
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 not end up being statistically significant, but the single factor's effect on the multiva

Re: [R] figure margins too large

2014-04-03 Thread 张以春
Dear Jim, Yes, what you said is right. I have resolved this question by your method. Thank you very much. Best, Yichun > -原始邮件- > 发件人: "Jim Lemon" > 发送时间: 2014年4月3日 星期四 > 收件人: "张以春" > 抄送: r-help@r-project.org > 主题: Re: [R] figu

Re: [R] grouping explanatory variables into "sets" for GLMM

2014-04-03 Thread Don McKenzie
Reading the Intro, as Bert suggests, would likely solve some of your problems. If you think about how many combinations it would take, using only one variable from each group in any one model, you would see that the number of individual models (12) is not so onerous that you couldn’t specify the

Re: [R] Doing box plots for part of the data??

2014-04-03 Thread Jim Lemon
On 04/04/2014 02:20 PM, David Doyle wrote: HI folks I'm wanting to do box plots of just 2 of my wells. I can do it for all of them by using boxplot(Barium~Well.ID,data=mydata, main="Barium", ylab="mg/L") to do for all 4 wells Anyone have a sugestion on how to do it for only 2 wells like MW-1

Re: [R] grouping explanatory variables into "sets" for GLMM

2014-04-03 Thread Bert Gunter
Unless there is reason to keep the conversation private, always reply to the list. How will anyone else know that my answer wasn't satisfactory? 1. I don't intend to go through your references. A minimal reproducible example of what you wish to do and what you tried would help. 2. Have you read A

[R] Doing box plots for part of the data??

2014-04-03 Thread David Doyle
HI folks I'm wanting to do box plots of just 2 of my wells. I can do it for all of them by using boxplot(Barium~Well.ID,data=mydata, main="Barium", ylab="mg/L") to do for all 4 wells Anyone have a sugestion on how to do it for only 2 wells like MW-1 and MW-2? The data can be found? http://www.

Re: [R] unzip error with odfWeave and OpenOffice

2014-04-03 Thread Pascal Oettli
Hello, Probably a path problem: https://stat.ethz.ch/pipermail/r-help/2010-March/232135.html Regards, Pascal On Fri, Apr 4, 2014 at 11:30 AM, Bob Kelly wrote: > Hello r-help mailing list readers, > > I am trying to create a report with OpenOffice (ver 4.0.1) using the > odfWeave (ver 0.8.4) pac

[R] unzip error with odfWeave and OpenOffice

2014-04-03 Thread Bob Kelly
Hello r-help mailing list readers, I am trying to create a report with OpenOffice (ver 4.0.1) using the odfWeave (ver 0.8.4) package. This is completely new to me so I am using part of an example from the book R in Action. I have used 7-zip on the file and this opens it correctly, am I doing s

Re: [R] Icelandic Characters in Mac

2014-04-03 Thread David Winsemius
On Apr 3, 2014, at 12:24 PM, peter dalgaard wrote: > > On 03 Apr 2014, at 00:59 , David Winsemius wrote: > >> >> On Apr 2, 2014, at 2:18 PM, Stefán Hrafn Jónsson wrote: >> >>> Dear R community >>> >>> I have few students that use Mac. When creating graphs they inform me that >>> when they u

[R] Nested Logit Model

2014-04-03 Thread Tim Marcella
Hi All, I am using the package 'mlogit' I am having trouble constructing a nested logit model. I am headed down this path as I am violating the IIA assumption when running a multinomial logistic regression. I am analyzing the choice a seabird floating on the surface of the water makes when appro

[R] Computing standardised beta in linear regression

2014-04-03 Thread Ruijie
Hi all, I​ know that there is a package QuantPsyc that has the function lm.beta to calculate standardised coefficients. However, the function is not able to handle categorical variables with more than 2 levels. I coded an alternative function and hope to know if I have done anything wrongly. ​

Re: [R] Time Series Decomposition On Zoo Objects: Errors

2014-04-03 Thread Rich Shepard
On Thu, 3 Apr 2014, Roy Mendelssohn wrote: The state-space approach has the advantage in the appropriate situations that you can model the trends and seasonals and cycles in a way that doesn't assume stationarity and provides a lot of flexibility. To me a lot of it depends on if the nature of th

Re: [R] Time Series Decomposition On Zoo Objects: Errors

2014-04-03 Thread Roy Mendelssohn
The state-space approach has the advantage in the appropriate situations that you can model the trends and seasonals and cycles in a way that doesn't assume stationarity and provides a lot of flexibility. To me a lot of it depends on if the nature of the irregularity is an inherent property of

Re: [R] Time Series Decomposition On Zoo Objects: Errors

2014-04-03 Thread Rich Shepard
On Thu, 3 Apr 2014, Roy Mendelssohn wrote: How irregular is irregular. kalman filter based methods, such as those in KFAS and DLM, can handle missing data, and often "irregular" data can be thought of as regular data with missing values, A lot depends on how irregular and how big the gaps, to th

Re: [R] Time Series Decomposition On Zoo Objects: Errors

2014-04-03 Thread Roy Mendelssohn
HI Rich: How irregular is irregular. kalman filter based methods, such as those in KFAS and DLM, can handle missing data, and often "irregular" data can be thought of as regular data with missing values, A lot depends on how irregular and how big the gaps, to the point where the analysis can

Re: [R] ff package: problem with reading in numerical values

2014-04-03 Thread Qiang Kou
Can you paste the line of data which caused the error? It will be much easier for other people to help you. Best, KK On Thu, Apr 3, 2014 at 6:35 AM, Srdjan Santic wrote: > I'm trying to read in a fairly large (3.5 Gb) csv file into R, using the > read.csv.ffdf() function from the ff package.

Re: [R] A question on IRanges package

2014-04-03 Thread Qiang Kou
As I as I know, the type should be one of "any", "start", "end", "within", "equal". You are using a undefined type. Best, KK On Thu, Apr 3, 2014 at 2:48 PM, Yuan Luo wrote: > Hi All, > Sorry for possible spam, but I am trying to customize IRanges package > locally. For what I am doing, I int

Re: [R] Time Series Decomposition On Zoo Objects: Errors

2014-04-03 Thread Rich Shepard
On Thu, 3 Apr 2014, arun wrote: Not sure if this helps you. http://stackoverflow.com/questions/12623027/how-to-analyse-irregular-time-series-in-r A.K., Yes, it does. I've read all the zoo docs I can find and have been searching for more information on irregular time series data. Environment

Re: [R] Equation of a curve

2014-04-03 Thread Greg Snow
The mean value theorem of integration (I have a cross-stitch of this theorem hanging on my wall (between cross-stitches of the central limit theorem and Bayes theorem)) tells us that the area under a curve is equal to the width of the area of interest times the average height of the curve. Often w

Re: [R] NHPP with Cyclic Rate

2014-04-03 Thread Rolf Turner
On 03/04/14 23:17, Andrew Wachira wrote: Is there an R-package that computes the parameters of a nonhomogenous Poisson process with cyclic rate (the rate function has polynomial and trigonometric components)? Something more like NPPMLE ( http://www.ise.ncsu.edu/jwilson/files/johnson94orl.pdf ) bu

[R] Time Series Decomposition On Zoo Objects: Errors

2014-04-03 Thread Rich Shepard
I have irregular time series as zoo objects; one example: structure(c(6, 5, 14, 9, 8, 9, 8, 5, 5, 5, 3, 3, 4, 3, 9, 6.94, 7.44, 3.09, 0.84, 5.35, 4.76, 4.21, 1.58, 2.6, 3.41, 9.59, 7.1, 5, 5, 5, 3, 1.5, 2.4, 3.9, 5.8, 2.3, 3.6, 4.1, 15.4, 7.8, 4.2, 5.8, 3, 4.5, 8.1, 9.6, 9.3, 7.9, 3.8, 3.2,

Re: [R] nearest neighbour

2014-04-03 Thread Rolf Turner
On 04/04/14 01:16, Monaly Mistry wrote: Hi, If I have a data frame of the location of individuals (x and y coordinate), how do I find the 5 nearest neighbours for each individual. The nnwhich() function from the spatstat package will do this for you. E.g.: require(spatstat) # You need to hav

[R] Stats courses in Australia

2014-04-03 Thread Highland Statistics Ltd
We are planning to run a series of four courses in Australia (and/or New Zealand) in August 2014. Potential courses: 1. Data Exploration, Regression, GLM & GAM with introduction to R 2. Introduction to Linear Mixed Effects Models and GLMM with R 3. Zero Inflated Models & GLMM with R 4. Beginner’

Re: [R] rpart and randomforest results

2014-04-03 Thread Schillo, Sonja
Hi, the random forest should do that, you're totally right. As far as I know it does so by randomly selecting the variables considered for a split (but here we set the option for how many variables to consider at each split to the number of variables available so that I thought that the random

Re: [R] Icelandic Characters in Mac

2014-04-03 Thread peter dalgaard
On 03 Apr 2014, at 00:59 , David Winsemius wrote: > > On Apr 2, 2014, at 2:18 PM, Stefán Hrafn Jónsson wrote: > >> Dear R community >> >> I have few students that use Mac. When creating graphs they inform me that >> when they use Icelandic characters in title or xlab they get some wrong >> re

[R] A question on IRanges package

2014-04-03 Thread Yuan Luo
Hi All, Sorry for possible spam, but I am trying to customize IRanges package locally. For what I am doing, I introduced another option to type parameters to the findOverlaps method. In the file findOverlaps-methods.R, I modified every instance of type = c("any", "start", "end", "within", "equal"),

Re: [R] Equation of a curve

2014-04-03 Thread Boris Steipe
You may have a different use in mind, but I think integration does not make sense for growth curves. And there is no simple, general equation that I'm aware of: When you determine the area under the curve (integration), you are essentially multiplying bacterial mass by time. Imagine that you w

Re: [R] Equation of a curve

2014-04-03 Thread Frede Aakmann Tøgersen
Hi Frances The short answer is YES, you can do anything in R if you bother. Do you have any figures of the growth curves (something close to exponential growth until they suddenly commit suicide or the environment kills)?. We NEED more information, thank you please. Yours sincerely / Med venli

Re: [R] How to store interim print results

2014-04-03 Thread Katherine Gobin
Dear Sir, Thanks a lot for your guidance and efforts. Appreciate it. Thanks again. Katherine On Thursday, 3 April 2014 6:55 PM, jim holtman wrote: This will get you close: > settlement = as.Date("2013-11-25") > maturity   = as.Date("2015-10-01") > coupon     = 0.066 > yield      = 0.1040 > b

[R] Equation of a curve

2014-04-03 Thread Frances Cheesman
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 [[alternative HTML version deleted]]

Re: [R] How to store interim print results

2014-04-03 Thread Katherine Gobin
Dear Jean Thanks a lot for this solution. Its very useful. I did only one small change and defined  cashflow.tenure <- numeric(0) instead of character(0). This helps me in further numerical calculations using these dates like finding the difference between two dates etc. Thanks again, Regar

Re: [R] text minning script

2014-04-03 Thread Frede Aakmann Tøgersen
And Anan, PLEASE make sure that you respond to the whole list and not only to me (answer to all in your preferred email client). Also PLEASE make sure that you post to this list in plain text and not HTML format (see below for your email formatting that I received). Happy text mining ;-) Yours

Re: [R] text minning script

2014-04-03 Thread Frede Aakmann Tøgersen
Hi Anan Don't expect people on this list to do YOUR work. If you have some specific problems with R then perhaps some people are devoted to help you. Now is the time when you help yourself. Google for "r text mining" or "r text mining example" and you will get some help. Certainly not any help

Re: [R] mzR and Rcpp version bug

2014-04-03 Thread Qiang Kou
Hi, I think error message is clear enough "mzR has been built against a different Rcpp version". mzR on bioconductor is built using Rcpp 0.10, and you are using Rcpp 0.11. The solution I suggest is to install the mzR from source not using the binary package. Best, KK On Wed, Apr 2, 2014 at 6:

[R] text minning script

2014-04-03 Thread Anan Halabi
I need simple script to work with text minning [[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 http://www.R-project.org/posting-guide.html a

Re: [R] Mixed Effects MANOVA

2014-04-03 Thread Bert Gunter
This is R-help, not a list that provides statistical help (primarily; they do intersect at times). Post to the r-sig-mixed-models list instead. You're likely to do better there anyway for this sort of thing. Cheers, Bert Bert Gunter Genentech Nonclinical Biostatistics (650) 467-7374 "Data is not

Re: [R] grouping explanatory variables into "sets" for GLMM

2014-04-03 Thread Bert Gunter
Have you read "An Introduction to R" (or other online tutorial)? If not, please do so before posting further here. It sounds like you are missing very basic knowledge -- on factors -- which you need to learn about before proceeding. ?factor gives you the answer you seek, I believe. Cheers, Bert

Re: [R] Plotting Satellite

2014-04-03 Thread Dao_De
Ok i figured it out. Here is the script that gives you the map. The second map should express with a two-color scale or code [indicating significant/non significant] whether your trends are truly significant at a 95% confidence level. All i got is this. And this time i am really stuck.. I

[R] Mixed Effects MANOVA

2014-04-03 Thread Nick Negovetich
Greetings, I have a question regarding data analysis of habitat use of animals. These animals were radio collared and tracked periodically throughout the year. When they were sighted/detected, the habitat type was marked. Our dataset recorded the sex of the animal, and we know the data whe

[R] NHPP with Cyclic Rate

2014-04-03 Thread Andrew Wachira
Is there an R-package that computes the parameters of a nonhomogenous Poisson process with cyclic rate (the rate function has polynomial and trigonometric components)? Something more like NPPMLE ( http://www.ise.ncsu.edu/jwilson/files/johnson94orl.pdf ) but in R? [[alternative HTML version

[R] nearest neighbour

2014-04-03 Thread Monaly Mistry
Hi, If I have a data frame of the location of individuals (x and y coordinate), how do I find the 5 nearest neighbours for each individual. [[alternative HTML version deleted]] __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/l

[R] Converting corpus into dataframe in r tm package

2014-04-03 Thread vikrant
HI , I would like to do stemming operation on vector of words using tm package. I am using R 2.8.1 and tm package for same. For doing any operations i tm package the data first needs to be converted to corpus and then use various commands in tm package. No problem is I am not able to convert these

[R] grouping explanatory variables into "sets" for GLMM

2014-04-03 Thread Maria Kernecker
Dear all, I am trying to run a GLMM following the procedure described by Rhodes et al. (Ch. 21) in the Zuur book Mixed effects models and extensions in R . Like in his example, I have four "sets" of explanatory variables: 1. Land use - 1 variable, factor (forest or agriculture) 2. Location - 1

Re: [R] 'rms' package error

2014-04-03 Thread Göran Broström
On 04/03/2014 04:29 AM, Lucy Leigh wrote: Hi everyone, I am attempting to use the R package 'rms' http://biostat.mc.vanderbilt.edu/wiki/Main/Rrms to implement a PH weibull model, using the pphsm() function. However, I get the following error, f.ph <- pphsm(f) Warning message: In pphsm(f) :

Re: [R] How to store interim print results

2014-04-03 Thread Adams, Jean
Katherine, One easy way to do this for small data is by using the append() function (see code below). But, if you have a lot of data, it may be too slow for you. In that case, you can gain some efficiency if you determine in advance how long the vectors will be, then use indexing to fill in the

Re: [R] trouble using readOGR() function

2014-04-03 Thread Roger Bivand
remiss free.fr> writes: > > Hi all, > I got some trouble trying to open a .kml file into R. Usually, the readOGR > package works great for it but here I get a message error that I can't > understand. > > When I'm typing this : > myfile <-readOGR(dsn="/windows/landuse.kml",layer="agricultura

Re: [R] How to store interim print results

2014-04-03 Thread jim holtman
This will get you close: > settlement = as.Date("2013-11-25") > maturity = as.Date("2015-10-01") > coupon = 0.066 > yield = 0.1040 > basis = 1 > frequency = 2 > redemption = 100 > > # __ > > add.months = function(date,

Re: [R] trouble using readOGR() function

2014-04-03 Thread Michael Sumner
On 1 Apr 2014 22:25, "remiss" wrote: > > Hi all, > I got some trouble trying to open a .kml file into R. Usually, the readOGR > package works great for it but here I get a message error that I can't > understand. > > When I'm typing this : > myfile <-readOGR(dsn="/windows/landuse.kml",layer="a

Re: [R] trouble using readOGR() function

2014-04-03 Thread Roger Bivand
remiss free.fr> writes: > > Hi all, > I got some trouble trying to open a .kml file into R. Usually, the readOGR > package works great for it but here I get a message error that I can't > understand. > > When I'm typing this : > myfile <-readOGR(dsn="/windows/landuse.kml",layer="agricultura

Re: [R] 'rms' package error

2014-04-03 Thread Kevin E. Thorpe
On 04/02/2014 10:29 PM, Lucy Leigh wrote: Hi everyone, I am attempting to use the R package 'rms' http://biostat.mc.vanderbilt.edu/wiki/Main/Rrms to implement a PH weibull model, using the pphsm() function. However, I get the following error, f.ph <- pphsm(f) Warning message: In pphsm(f) : at

[R] ff package: problem with reading in numerical values

2014-04-03 Thread Srdjan Santic
I'm trying to read in a fairly large (3.5 Gb) csv file into R, using the read.csv.ffdf() function from the ff package. Whenever the function comes up to a numerical value (such as 40.0, or 2 - doesn't mater if there is a decimal point or not), it throws the following error: Error in scan(file, wh

[R] How to store interim print results

2014-04-03 Thread Katherine Gobin
Dear R forum, Following is an customized extract of a code I am working on. settlement = as.Date("2013-11-25") maturity   = as.Date("2015-10-01") coupon     = 0.066 yield      = 0.1040 basis      = 1   frequency = 2 redemption = 100 # _

Re: [R] data.frame(1)*1:4 = 1?

2014-04-03 Thread Spencer Graves
Hi, Josh & Fabian: Thanks for the replies. On 4/3/2014 12:07 AM, Joshua Wiley wrote: > Hi Spencer, > > One piece is that a data frame of the same dimensions as went in comes > out. The second piece is that the vector is recycled. > > So in your first example: > > data.frame(1) * 1:4 >

Re: [R] figure margins too large

2014-04-03 Thread Prof Brian Ripley
On 03/04/2014 06:32, 张以春 wrote: Dear R experts, I tried to plot some figures in R using postscript(), but it always shows that the fugures margin is too large. I don't know how to change it. The following is my example: postscript("All.eps",width=3.27,height=1.416,pointsize=12,family="Aria

Re: [R] data.frame(1)*1:4 = 1?

2014-04-03 Thread Joshua Wiley
Hi Spencer, One piece is that a data frame of the same dimensions as went in comes out. The second piece is that the vector is recycled. So in your first example: data.frame(1) * 1:4 you only end up with the first element: data.frame(1) * 1 If you try: data.frame(1) * 4:1 you get a data fr

Re: [R] figure margins too large

2014-04-03 Thread Jim Lemon
On 04/03/2014 04:32 PM, 张以春 wrote: Dear R experts, I tried to plot some figures in R using postscript(), but it always shows that the fugures margin is too large. I don't know how to change it. The following is my example: postscript("All.eps",width=3.27,height=1.416,pointsize=12,family="A