[R] how to resolve the install problem in redhat linux

I try to install a R resource package (in R.2.15 or R 3.0) But when I run the command configure,the Makefile does not generated succefully as error information shown as below --with-readline=yes (default) and headers/libs are not available my linux sysem is Red Hat Enterprise Linux Client

Re: [R] (no subject)

On 07/17/2013 01:41 PM, labib obaid wrote: hi I need to use boxplot for two subsets at the same page Hi labib obaid, Try this: test.df-data.frame(a=rnorm(80)+4,b=rnorm(80)+4, c=rep(LETTERS[1:4],each=20), d=rep(rep(letters[1:4],each=4),5)) boxplot(test.df$a[test.df$c %in% c(A,B)],

Re: [R] how to resolve the install problem in redhat linux

On 17/07/2013 07:02, Jie Tang wrote: I try to install a R resource package (in R.2.15 or R 3.0) There are no such versions of R: see the posting guide. But when I run the command configure,the Makefile does not generated succefully as error information shown as below --with-readline=yes

Re: [R] ERROR: missing value where TRUE/FALSE needed

Hi What is mctp? Why do you think your data are well defined? Did you try str(Datos)? Where you get this error? When you do call mctp or summary? The main problem is that you do not follow posting guide, do not provide reproducible code and/or at least some info about your data. Without that

Re: [R] Setting Derived Class Slots

Dohso simple when you know how! That worked straightaway! Thank you Martin! -- View this message in context: http://r.789695.n4.nabble.com/Setting-Derived-Class-Slots-tp4671683p4671737.html Sent from the R help mailing list archive at Nabble.com.

Re: [R] Setting Derived Class Slots

Hi Simon.many thanks to you for your explanation. - I'll check out setReplaceMethod. Point taken on nomenclature too... -- View this message in context: http://r.789695.n4.nabble.com/Setting-Derived-Class-Slots-tp4671683p4671738.html Sent from the R help mailing list archive at

[R] error message in gev

  Hi r-users,   I would like to use gev and my data (annual rainfall ) is as follows:   head(dat,20) ABCDEFGHI J 1 45.1 41.5 58.5 50.1 46.0 49.1 37.7 49.1 59.8 54.0 2 50.3 39.8 49.4 56.4 49.4 48.8 42.1 49.8 49.4 58.3 3 41.7 39.3 44.6 39.1 35.7 41.5 40.8

[R] extract beta.sim from dist type data (package betapart)

Dear List, This is Elaine. I am using beta.part to calculate the beta diversity index. The function of beta.part can generate three kinds of beta diversity indice at one time. However, it is Simpson index that I want. The calculation result is composed of three rows and N/A columns (dist class).

Re: [R] error message in gev

Hello, You should say which package is gev() coming from. I believe it's from evir, and what follows assumes that. The problem is that 'dat' is a data.frame, a special type of list. It is a list of vectors and gev is expecting just one vector. You could try instead, to fit generalized extreme

Re: [R] Box plot with 5th and 95th percentiles instead of 1.5 * IQR: problems implementing an existing solution...

Hello, you may just first change the boxplot.stats directly: boxplot.stats - function (x, coef = NULL, do.conf = TRUE, do.out = TRUE) { nna - !is.na(x) n - sum(nna) stats - quantile(x, c(.05,.25,.5,.75,.95), na.rm = TRUE) iqr - diff(stats[c(2, 4)]) out

[R] Help: Error when installing R 2.8.1 in IBM AIX system from source code

Hi R-Help, When we tried to install R 2.8.1 in our AIX system, we faced some error, please help. OS level: AIX 7.1 64 bit R 2.8.1 We had a guide about How to install R 2.8.1 in AIX from R source code. Here is what I did: 1. Download the R source file, for example R-2.8.1.tar.gz , from

[R] Some problems with back to back bar plots in ggplot2

Hi, I'm making a simple population pyramid using two back-to-back bar plots for 18 different age groups with totals for males and females. Yesterday I achieved a fairly serviceable plot using the following. library(ggplot2) library(reshape2) library(plyr) #make sample data

Re: [R] error message in gev

Hello, It's better if you Cc the list, the odds of getting more and better answers are greater. Have you tried with different block values? Rui Barradas Em 17-07-2013 13:33, Roslina Zakaria escreveu: Hi Rui, Yes, I am using gev from 'evir' install.packages(evir) library(evir) I have tried

Re: [R] error in predict.gam used with bam

Solved it!! ;) The problem was that the test-data contained factor-levels the training-data didn't. So when trying to fit this new factor-levels to a model which didn't have this levels, the error occurred. When excluding the factor-levels not used when fitting the model or when taking care all

[R] Fw: error message in gev

- Forwarded Message - From: Rui Barradas ruipbarra...@sapo.pt To: Roslina Zakaria zrosl...@yahoo.com Cc: 'r-help' r-help@r-project.org Sent: Wednesday, July 17, 2013 8:49 PM Subject: Re: [R] error message in gev Hello, It's better if you Cc the list, the odds of getting more and

Re: [R] Question on plotting with googleVis

Hello Arnaud, Thank you for your pointer. However I need to more clarification. I want to control the max. and min. values for the x-axis, as well as number of vertical gridlines to be displayed. I tried the following: MyData - data.frame(Names1 = paste(XXX, 1:150), Values1 = 1:150 + 10,

Re: [R] error message in gev

Hi Rui,   I have tried different block values and below are the error message I got:   out - gev(dat, A) Error in as.POSIXlt.default(attributes(data)$times) : do not know how to convert 'attributes(data)$times' to class “POSIXlt”   out.list - lapply(dat, gev, block = 10) Error in

Re: [R] [R-sig-eco] extract beta.sim from dist type data (package betapart)

Replied to on the R-sig-ecology list, with a request not to cross-post. Briefly, the output is a list with named components, so can be extracted as usual. Sarah On Wed, Jul 17, 2013 at 5:50 AM, Elaine Kuo elaine.kuo...@gmail.com wrote: Dear List, This is Elaine. I am using beta.part to

Re: [R] Question on plotting with googleVis

Hi You can do a rotation and use gvisColumnChart instead gvisBarChart plot(gvisColumnChart(MyData, xvar=Names1, yvar=c(Values1, Values2),options=list(width=2500,height=1000))) Michel Le 17/07/2013 15:57, Christofer Bogaso a écrit : Hello Arnaud, Thank you for your pointer. However I need to

Re: [R] Splitting dataframes and cleaning extraneous characters

Hi, YOu could try. ?split() split(ats,ats$Project_NBR) You also mentioned about two columns. split(ats,list(ats$col1, ats$col2)) You should have provided an example dataset using ?dput() ( dput(head(data,10)) ) for testing. Also, gsub(^-[^-]*-,,-005-190) #[1] 190 A.K. Problem: I have a

Re: [R] error message in gev

Hello, And what if you try it without a block value? out.list - lapply(dat, gev) Rui Barradas Em 17-07-2013 14:53, Roslina Zakaria escreveu: Hi Rui, I have tried different block values and below are the error message I got: out - gev(dat, A) Error in

Re: [R] Some problems with back to back bar plots in ggplot2

Thanks, I've now got the plot I was after. No idea what the stacking warning is about, but it does not seem to compromise my output. Gavin. -- View this message in context: http://r.789695.n4.nabble.com/Some-problems-with-back-to-back-bar-plots-in-ggplot2-tp4671744p4671765.html Sent from the

Re: [R] Some problems with back to back bar plots in ggplot2

Hi, I'm making a simple population pyramid using two back-to-back bar plots for 18 different age groups with totals for males and females. I can run both variants without error, though both throw the same warning regarding stacking beein 'not well defined ..' However, the axis scaling in

Re: [R] Help: Error when installing R 2.8.1 in IBM AIX system from source code

Help with outdated versions of R, and particularly with nonstandard patch files, are not really on topic here. Posting in HTML format is not advised (see posting guide) since it actually makes it harder to read R code, so text coloring doesn't come through. It appears to me that your patch

Re: [R] Question on plotting with googleVis

Hi Arnaud, Thanks for your answer. However I prefer to have Horizontal bar chart, because length of the strings for x-axis are quite large, therefore it would be better to put them vertically one-below-another instead horizontally. Therefore I would really appreciate if someone points me how to

[R] R Package License

HI Helpers, How could we use R and R packages licensed under GPL into commercial products? Is it allowed to load a library and get the results from it and using the results for commercial use? Thank you so much! Regards, Yan [[alternative HTML version deleted]]

Re: [R] Masking oceans using polypath

Many thanks, this works beautifully. Louise On Wed, Jul 17, 2013 at 7:42 AM, Paul Murrell p...@stat.auckland.ac.nzwrote: Hi There are a couple of problems: 1. Your 'outline' is much bigger than it needs to be. For example, the following produces just Australia ... outline -

[R] R-squared and GLM

Dear users, I want to compute r-squared values from a glm regression using a gamma distribution and an identity link-function, but find no such thing when using the summary() or names() function. My next guess was to calculate it by hand, i.e. r2 = (sum((estimate - xbar)^2) /sum((x-xbar)^2))

Re: [R] R Package License

On Jul 17, 2013, at 10:56 AM, Li, Yan yan...@ibi.com wrote: HI Helpers, How could we use R and R packages licensed under GPL into commercial products? Is it allowed to load a library and get the results from it and using the results for commercial use? Thank you so much! Regards, Yan

[R] simplify a dataframe

Hi Arun I have two questions always about the question of symplify a dataframe I would like 1) to transform the vector x1 into the vector y1 x1 - c(1,1,1,-1000, 1,-1000, 1,1,1,1,1,1,-1000) y1 - c(1,1,1,1,2,2, 3,3,3,3,3,3,3) 2) to transform the vectors

Re: [R] writing multiple lines to a file

Hi, No problem. You could try: printer = file(out.txt,w)  writeLines(This is line.,con=printer,sep= )  writeLines(The same line.,con=printer)  close(printer) #or cat(sprintf(This is line %d. ,1),file=out.txt,append=TRUE) cat(The same line.,file=out.txt,append=TRUE) A.K. Thank you very much, I

[R] truncation and approximation

What is the function to do the truncation to a certain decimal digit of a number. And the function approximation? Thanks. Francesco M [[alternative HTML version deleted]]

Re: [R] Splitting dataframes and cleaning extraneous characters

HI, One problem with using ?subst() would be it depends upon the number of digits, characters etc.  For eg. substring(-005-190,6) #[1] 190  substring(-0057-190,6) #[1] -190 #whereas gsub(^-[^-]*-,,-0057-190) #[1] 190 Probably, your dataset doesn't have that sort of problem. dat1-

Re: [R] simplify a dataframe

Hello, As for question (1), try the following. y2 - cumsum(c(TRUE, diff(x1) 0)) identical(as.integer(y1), y2) # y1 is of class numeric As for question (2) I'm not understanding it. Hope this helps, Rui Barradas Em 17-07-2013 18:21, Arnaud Michel escreveu: Hi Arun I have two questions

[R] C Interface in R

Dear R Users--- I spent some time implementing the AS75 WLS regression algorithm in C in order to learn the R interface to C. I did not find an easy, brief C+R example of an algorithm that used persistent storage across calls. they probably exist, but I could not find a nice brief template for

[R] How to open .WTG file Extension in R?

Hi R user, I am wondering how we can open the data file which have .WTG extension? Thanks for your help cheers, KG [[alternative HTML version deleted]] __ R-help@r-project.org mailing list

Re: [R] How to open .WTG file Extension in R?

I am wondering what .WTG files have to do with R. Care to elaborate? (Suggestion: Read and follow the posting guide) Cheers, Bert On Wed, Jul 17, 2013 at 11:48 AM, Kristi Glover kristi.glo...@hotmail.com wrote: Hi R user, I am wondering how we can open the data file which have .WTG

Re: [R] How to open .WTG file Extension in R?

On Jul 17, 2013, at 1:48 PM, Kristi Glover kristi.glo...@hotmail.com wrote: Hi R user, I am wondering how we can open the data file which have .WTG extension? Thanks for your help cheers, KG A quick search led me to http://wasptechnical.dk/Forum/viewtopic.php?id=433, which suggests that

Re: [R] How to open .WTG file Extension in R?

Hi Bert, This is the data related to climate. I can open it in Excel but the problem is I do have 190 files . thanks Date: Wed, 17 Jul 2013 12:00:58 -0700 Subject: Re: [R] How to open .WTG file Extension in R? From: gunter.ber...@gene.com To: kristi.glo...@hotmail.com CC:

Re: [R] How to open .WTG file Extension in R?

That reduces the number of possibilities for the file structure to ones that R can probably read, but we need to know more. Try opening a sample file in a text editor to see if it is plain text. If so, you can probably use read.table() or read.csv() to get the file into R. If not, you will need

Re: [R] truncation and approximation

Hello, As for truncation, you can write a one line function: dec_trunc - function(x, digits = 0) trunc(x * 10^digits)/10^digits As for approximation, maybe you're looking for ?signif. Hope this helps, Rui Barradas Em 17-07-2013 18:41, Francesco Miranda escreveu: What is the function to do

Re: [R] How to open .WTG file Extension in R?

Hello, So it's a XML file, like Marc said. See package XML. install.packages('XML', dependencies = TRUE) Hope this helps, Rui Barradas Em 17-07-2013 20:04, Kristi Glover escreveu: Hi Bert, This is the data related to climate. I can open it in Excel but the problem is I do have 190 files .

Re: [R] simplify a dataframe

Thank you for the question (1) Sorry for the imprecision for the question (2) : Suppose the date frame df df1 - data.frame( Debut =c ( 24/01/1995, 01/05/1997 ,31/12/1997, 02/02/1995 ,28/02/1995 ,01/03/1995, 13/03/1995, 01/01/1996, 31/01/1996) , Fin = c ( 30/04/1997, 30/12/1997 ,31/12/1997,

[R] anti off diagonal min max mean - how to rotate matrix by 90 degree

How can id do this efficiently in R ? 1 0 0 0 2 0 0 0 3 rotate right 0 0 1 0 2 0 3 0 0 rotate left 0 0 3 0 2 0 1 0 0 What I want to do is described here: http://stackoverflow.com/questions/13049575/r-min-max-and-mean-of-off-diagonal-elements-in-a-matrix but I want do to it for all

Re: [R] simplify a dataframe

Hi, You could try: df1[,1:2]-lapply(df1[,1:2],as.character)  df2New- data.frame(Deb=unique(with(df1,ave(Debut,INDX,FUN=function(x) head(x,1,Fin=unique(with(df1,ave(Fin,INDX,FUN=function(x) tail(x,1) identical(df2New,df2) #[1] TRUE A.K. - Original Message - From: Arnaud Michel

[R] EWMA --http://www.forecastingfinancialrisk.com/3.html

Hi, I followed the code given in the link http://www.forecastingfinancialrisk.com/3.html EWMA-matrix(nrow=T,ncol=5) S-cov(returns) dim(S) [1] 5 5 EWMA[1,] = c(S)[c(1,4,2)] # extract the variances and covariance *here we will get a vector but we are equating with a matrix. hence the

Re: [R] simplify a dataframe

#or library(plyr) res-ddply(df1,.(INDX),summarize,Debut=head(Debut,1),Fin=tail(Fin,1)) res$INDX-factor(res$INDX,levels=unique(df1$INDX)) res[order(res$INDX),-1] #   Debut    Fin #3 24/01/1995 31/12/1997 #4 02/02/1995 12/03/1995 #1 13/03/1995 30/06/1995 #2 01/01/1996 31/01/1996 A.K.

[R] EWMA error

hi, Could anyone help me in solving the following error: I have 5 stocks returns data (returns) EWMA = matrix(nrow=T,ncol=5) # create a matrix to hold the covariance matrix for each t lambda = 0.94 S-cov(returns) # initial (t=1)

[R] combine select data from 2 dataframes sharing same variables

# The following dataframes are the result of two analyses performed on the same set of numeric data. # The first analysis involved calculations that did not include zero values: StatsUTAH = data.frame(MWtotaleesDue = c(8.428571,2.496256,7,6.604472,1,17,3.593998,4.834573,12.02257),

Re: [R] anti off diagonal min max mean - how to rotate matrix by 90 degree

On Jul 17, 2013, at 3:04 PM, Witold E Wolski wewol...@gmail.com wrote: How can id do this efficiently in R ? 1 0 0 0 2 0 0 0 3 rotate right 0 0 1 0 2 0 3 0 0 rotate left 0 0 3 0 2 0 1 0 0 What I want to do is described here:

Re: [R] simplify a dataframe

Thanks Arun and Rui for your helps Michel Le 17/07/2013 22:20, arun a écrit : #or library(plyr) res-ddply(df1,.(INDX),summarize,Debut=head(Debut,1),Fin=tail(Fin,1)) res$INDX-factor(res$INDX,levels=unique(df1$INDX)) res[order(res$INDX),-1] # DebutFin #3 24/01/1995 31/12/1997 #4

[R] How to capture the printout on the screen?

This is most likely a silly question. First I run the following: require(car) mod.ok - lm(cbind(pre.1, pre.2, pre.3, pre.4, pre.5, post.1, post.2, post.3, post.4, post.5, fup.1, fup.2, fup.3, fup.4, fup.5) ~ treatment*gender, data=OBrienKaiser) phase -

[R] cut into groups of equal nr of elements...

I would like to cut a vector into groups of equal nr of elements. looking for a function on the lines of cut but where I can specify the size of the groups instead of the nr of groups. -- Witold Eryk Wolski __ R-help@r-project.org mailing list

Re: [R] ERROR: missing value where TRUE/FALSE needed

In addition to what Petr said, That message most likely means that you have some missing data where you're not allowed to have missing data in your Datos. Look and see if the mctp() function has an argument that controls what it does with missing data. -Don -- Don MacQueen Lawrence Livermore

[R] robust smoothing - median smoothing.

Looking for robust smoothing methods available in R... a robust version of lowess or smooth.spline. Was searching for median and smoothing but the search engine isn't calibrated . -- Witold Eryk Wolski __ R-help@r-project.org mailing list

Re: [R] How to capture the printout on the screen?

If you are happy with the character strings printed out then you can use the capture.output function. If you want the numbers (without needing to convert) then look at the print method for the class of object that you are working with and see what the code is. It is probably calling another

Re: [R] cut into groups of equal nr of elements...

You could use the quantile function to choose the cut points and pass that to cut. Or you could sort the vector and just take the first n elements as the 1st group, etc. On Wed, Jul 17, 2013 at 3:43 PM, Witold E Wolski wewol...@gmail.com wrote: I would like to cut a vector into groups of

Re: [R] cut into groups of equal nr of elements...

HI, Not sure whether this is what you wanted.  vec1- 1:7  fun1- function(x,nr) {((x-1)%/%nr)+1}  fun1(vec1,2) #[1] 1 1 2 2 3 3 4  fun1(vec1,3) #[1] 1 1 1 2 2 2 3 split(vec1,fun1(vec1,2)) A.K. - Original Message - From: Witold E Wolski wewol...@gmail.com To: r-help@r-project.org Cc:

Re: [R] robust smoothing - median smoothing.

Hello, Maybe ?stats::smooth Hope this helps, Rui Barradas Em 17-07-2013 22:47, Witold E Wolski escreveu: Looking for robust smoothing methods available in R... a robust version of lowess or smooth.spline. Was searching for median and smoothing but the search engine isn't calibrated . --

[R] glmnet on Autopilot

Dear List, I'm running simulations using the glmnet package. I need to use an 'automated' method for model selection at each iteration of the simulation. The cv.glmnet function in the same package is handy for that purpose. However, in my simulation I have p N, and in some cases the selected

Re: [R] cut into groups of equal nr of elements...

Sorry, there was a mistake: fun1 should be: fun1- function(x,nr) {((seq_along(x)-1)%/%nr)+1} vec3- c(4,5,7,9,8,5)  fun1(vec3,2) #[1] 1 1 2 2 3 3 split(vec3,fun1(vec3,2)) A.K. - Original Message - From: arun smartpink...@yahoo.com To: Witold E Wolski wewol...@gmail.com Cc: R help

Re: [R] cut into groups of equal nr of elements...

?qcut On Jul 17, 2013 5:45 PM, Witold E Wolski wewol...@gmail.com wrote: I would like to cut a vector into groups of equal nr of elements. looking for a function on the lines of cut but where I can specify the size of the groups instead of the nr of groups. -- Witold Eryk Wolski

Re: [R] cut into groups of equal nr of elements...

On Jul 17, 2013, at 4:43 PM, Witold E Wolski wewol...@gmail.com wrote: I would like to cut a vector into groups of equal nr of elements. looking for a function on the lines of cut but where I can specify the size of the groups instead of the nr of groups. In addition to the other options,

Re: [R] combine select data from 2 dataframes sharing same variables

Tena koe Without reading your request in detail, I will suggest you look at ?merge. It is often the answer when 'combine' is in the question. Peter Alspach -Original Message- From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org] On Behalf Of bcrombie Sent:

Re: [R] EWMA error

Hello, 1) In the provided example, you have 2 stock returns and EWMA = matrix(nrow=T, ncol=3) See the number of columns (2+1) Please modify the number of columns according to the number of stock returns. 2) As you have 5 stock returns, the following cannot work EWMA[1,] = c(S)[c(1,4,2)]

Re: [R] robust smoothing - median smoothing.

Hello, require(sos) findFn(median and smoothing) Regards, Pascal 2013/7/18 Rui Barradas ruipbarra...@sapo.pt Hello, Maybe ?stats::smooth Hope this helps, Rui Barradas Em 17-07-2013 22:47, Witold E Wolski escreveu: Looking for robust smoothing methods available in R... a robust

Re: [R] How to capture the printout on the screen?

Thanks a lot for the pointer! After I downloaded the source code and saw the innards of the print function, I know what to do now. Thanks again, Gang On Wed, Jul 17, 2013 at 6:00 PM, Greg Snow 538...@gmail.com wrote: If you are happy with the character strings printed out then you can use the

Re: [R] robust smoothing - median smoothing.

Another option is library(locfit) ?locfit There is also a web page from the help which gives a link to Clive Loaders book Regards Duncan Duncan Mackay Department of Agronomy and Soil Science University of New England Armidale NSW 2351 Email: home: mac...@northnet.com.au At 08:23

Re: [R] glmnet on Autopilot

On Jul 17, 2013, at 5:26 PM, Axel Urbiz wrote: Dear List, I'm running simulations using the glmnet package. I need to use an 'automated' method for model selection at each iteration of the simulation. The cv.glmnet function in the same package is handy for that purpose. However, in my

Re: [R] R-squared and GLM

On 18/07/13 02:36, Chris89 wrote: Dear users, I want to compute r-squared values from a glm regression using a gamma distribution and an identity link-function, but find no such thing when using the summary() or names() function. My next guess was to calculate it by hand, i.e. r2 =

Re: [R] error message in gev

Hi Rui,   Yes, it works beatifully this time.   Thank you so much for your help. From: Rui Barradas ruipbarra...@sapo.pt Cc: R help forum R-help@r-project.org Sent: Wednesday, July 17, 2013 11:28 PM Subject: Re: [R] error message in gev Hello, And

Re: [R] combine select data from 2 dataframes sharing same variables

Hi, Not sure if this is what you wanted: #If columns are arranged in the same order in both data.frames. lst1-lapply(seq_len(ncol(StatsUTAH)),function(i) {x1-cbind(StatsUTAH[,i],sStatsUTAH[,i]);row.names(x1)-row.names(StatsUTAH);colnames(x1)-c(zeroNO,zeroYES);x1})  names(lst1)-