Re: [R] Simple Error Bar

On 12/06/2013 04:16 PM, mohan.radhakrish...@polarisft.com wrote: Hi, Basic question with basic code. I am simulating a set of 'y' values for a standard 'x' value measurement. So here the error bars are very long because the number of samples are very small. Is that correct ? I am

[R] Simple Error Bar

Hi, Basic question with basic code. I am simulating a set of 'y' values for a standard 'x' value measurement. So here the error bars are very long because the number of samples are very small. Is that correct ? I am plotting the mean of 'y' on the 'y' axis. Thanks, Mohan x

Re: [R] Test ADF differences in R and Eviews

On Dec 5, 2013, at 3:18 PM, nooldor wrote: > Hi, > > > In attachment you can find source data on which I run adf.test() and > print-screen with results in R and Eviews. > > Results are very different. Did I missed something? Yes. You missed the list of acceptable file types for r-help. -- D

Re: [R] How to generate a smoothed surface for a three dimensional dataset?

On 13-12-05 7:02 PM, Jun Shen wrote: Thanks again, Duncan. Please allow me to ask one more question. Is it possible to generate a contour plot overlaying with the plot3d() plot? Yes, the contourLines function from grDevices can calculate the lines, then you can use rgl::lines3d to draw them.

[R] help with the nested anova formulas

I am modeling grade as a function of membership in various cohorts. There are four "cohorts". (NONE, ISE07,ISE08,ISE09) and two times of cohorts coded as ISE = TRUE (ISE0#) or FALSE (NONE). There is clear co-linearity but that is to be expected. running the following code CutOff <-0 fit.base <

Re: [R] convert factor to numeric in bulk

Hi, colnames(x1) <- gsub("NULL\\.","",colnames(x1)) ##not needed. A.K. On Thursday, December 5, 2013 5:04 PM, arun wrote: Hi, Try: library(XML) x1 <- readHTMLTable('http://www.fns.usda.gov/pd/15SNAPpartPP.htm',skip.rows=c(1:2,4,58:60),colClasses=c("character",rep("FormattedNumber",5)),strings

Re: [R] How to graph categorical data by percent

Hi, Try: dat1 <- structure(list(Site = c("Downstream", "Downstream", "Downstream", "Downstream", "Downstream", "Downstream", "Downstream", "Downstream", "Downstream", "Midstream", "Midstream", "Midstream", "Midstream", "Midstream", "Upstream", "Upstream", "Upstream", "Upstream", "Upstream", "Up

[R] Test ADF differences in R and Eviews

Hi, In attachment you can find source data on which I run adf.test() and print-screen with results in R and Eviews. Results are very different. Did I missed something? Best, T.S. __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listin

Re: [R] convert factor to numeric in bulk

Hi, Try: library(XML) x1 <- readHTMLTable('http://www.fns.usda.gov/pd/15SNAPpartPP.htm',skip.rows=c(1:2,4,58:60),colClasses=c("character",rep("FormattedNumber",5)),stringsAsFactors=FALSE)[[1]] colnames(x1) <- gsub("NULL\\.","",colnames(x1))  str(x1) #'data.frame':    53 obs. of  6 variables: # $ S

Re: [R] How to generate a smoothed surface for a three dimensional dataset?

On Dec 5, 2013, at 4:02 PM, Jun Shen wrote: > Thanks again, Duncan. Please allow me to ask one more question. Is it > possible to generate a contour plot overlaying with the plot3d() plot? > This is not going to give you the kewl rotation capabilities that `rgl` offers but you might be intere

Re: [R] help with the nested anova formulas

I am modeling grade as a function of membership in various cohorts. There are four "cohorts". (NONE, ISE07,ISE08,ISE09) and two times of cohorts coded as ISE = TRUE (ISE0#) or FALSE (NONE). There is clear co-linearity but that is to be expected. running the following code CutOff <-0 fit.

Re: [R] How to generate a smoothed surface for a three dimensional dataset?

Thanks again, Duncan. Please allow me to ask one more question. Is it possible to generate a contour plot overlaying with the plot3d() plot? Jun On Thu, Dec 5, 2013 at 11:08 AM, Duncan Murdoch wrote: > On 05/12/2013 10:33 AM, Jun Shen wrote: > >> Hi Federico/Duncan/David/Bert, >> >> Thank you f

[R] How do I print predicted effect sizes in forest plot?

Hi, I am struggling a bit with creating a forest plot containing the predicted effect size. As seen in other studies these effect sizes are shown per study usually as a light grey diamond - which is what I want to achieve. The calls I use are: iat_result = rma(yi=Mean, vi=Variance_rounded, ni=N

Re: [R] Can I find the number of subscribers to r-help and r SIG mailing lists?

On Dec 5, 2013, at 10:16 AM, Andrew Hoerner wrote: > Dear Folks-- > Is there a reasonably straightforward way for me to find the number of > subscribers That's not easily accessed by any publicly avaiable method I know but I can give the current count: 7783 Non-digested Members of R-help: 7563

[R] r SIG birth, life and death

I would like to correspond or talk on the phone with anyone who has founded a SIG or participated in numerous SIGs and has thoughts or lessons learned about what makes them thrive, stagnate or die. Also, if anyone has written anything useful on this topic (SIG health, SIG growth, SIG moderation (in

[R] ASA John M. Chambers Statistical Software Award - 2014

John M. Chambers Statistical Software Award - 2014 Statistical Computing Section American Statistical Association The Statistical Computing Section of the American Statistical Association announces the competition for

[R] ML factor analysis

Does anyone know how to get standard errors for factor loadings from factanal? Thanks, Greg __ 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 an

[R] ASA Stat. Computing & Stat. Graphics Student Paper Competition 2014

Dear all, Please see the Student Paper Competition/Travel Award announcement to Attend the Joint Statistical Meetings to be held in Boston in 2014. Call info below: Please send your inquiries, entries, and application material to me at (aarti.mun...@ucdenver.edu

[R] Can I find the number of subscribers to r-help and r SIG mailing lists?

Dear Folks-- Is there a reasonably straightforward way for me to find the number of subscribers or distinct posters to the r-help list and the various r SIG lists? Ideally I would like a time series from founding, but unadorned current membership numbers would still be useful, as would founding dat

Re: [R] Double Infinite Integration

On Dec 4, 2013, at 12:07 PM, Aya Anas wrote: > Hello all, > > I need to perform the following integration where the integrand is the > product of three functions: > f(x)g(y)z(x,y) Since f(x) does not depend on y, could you not do: Int( f(x) * Int( Int( g(y)*z(x,y).dy)).dx) # not R code Then

Re: [R] how to link different versions of R

On Dec 5, 2013, at 5:15 AM, Charles Thuo wrote: > Hi, > > I had downloaded the package zoo under R version 3.0.0. Afterwards i > downloaded R version 3.0.2 on the same workstation. Is it possible to use > the zoo package from 3.0.0. in 3.0.2. > Should be not only possible but trivial. > Is i

Re: [R] How do I extract Random Forest Terms and Probabilities?

Gents: This discussion is now off-topic here, I believe. Please take it private or move to somewhere else more appropriate (SO, maybe). Cheers, Bert On Thu, Dec 5, 2013 at 8:32 AM, Lopez, Dan wrote: > Hi Andy, > > I have used predict before and in fact when I do that to the train set I get > a

Re: [R] GAM Assumption Tests

Hi Mike, I recently had this issue and didn't find any package that implemented these tests directly for the gam object. I found it simplest just to pull the residuals from it and run tests like shapiro.test directly. Best, Collin. On Thu, 5 Dec 2013, Mike.lang wrote: > Dear all

Re: [R] How do I extract Random Forest Terms and Probabilities?

Hi Andy, I have used predict before and in fact when I do that to the train set I get a perfect model (i.e. ROC right angle curve in upper left quadrant) which just looks like it overfit the data. This is not the case with the test set were I get auc of .77. I wanted to attempt a couple of cal

Re: [R] How to generate a smoothed surface for a three dimensional dataset?

On 05/12/2013 10:33 AM, Jun Shen wrote: Hi Federico/Duncan/David/Bert, Thank you for your thoughtful comments and it's a great learning experience. I can see the critical point here is to find a right function to make the prediction. So I was thinking to start with "loess". However the predict.l

Re: [R] Merging two data frames, but keeping NAs

Hi, Try ?join() library(plyr) y$ref <- y$id > join(x,y,by="ref")    ref id val 1   NA NA  NA 2   NA NA  NA 3   NA NA  NA 4   10 NA  NA 5    9  9 106 6    8  8 107 7    7  7 105 8    6  6 104 9    5 NA  NA 10  NA NA  NA 11   1 NA  NA 12   2  2 101 13   3  3 102 14   4  4 103 15   5 NA  NA A.K.

Re: [R] how to store multiple data frame name and print in R?

Hi, May be this helps: i<-1 li<-list() while(i<3) { li[[i]]<-get(paste("L", i, sep = "")) print(li[[i]]) i<-i+1 } li[[1]] #   A  B #1 V1 v2 #2 V2 V3 #you can also store it by: l2 <- list(L1=L1,L2=L2,L3=L3) A.K. I have 3 data frame L1<-data.frame(A=c("V1","V2"),B=c("v2","V3")) L2<-d

Re: [R] How to generate a smoothed surface for a three dimensional dataset?

Jun: You need to: 1. Read the docs 2. Consult a local statistician, as you seem to be out of your depth statistically here. Re: 1: ?predict.loess explicitly says: "If newdata was the result of a call to expand.grid, the predictions (and s.e.'s if requested) will be an **array** of the appropri

Re: [R] How to generate a smoothed surface for a three dimensional dataset?

Hi Federico/Duncan/David/Bert, Thank you for your thoughtful comments and it's a great learning experience. I can see the critical point here is to find a right function to make the prediction. So I was thinking to start with "loess". However the predict.loess gave me an error as follows Error in

Re: [R] Double Infinite Integration

> -Original Message- > I need to perform the following integration where the integrand is the > product of three functions: > f(x)g(y)z(x,y) > > the limits of x are(0,inf) and the limits of y are(-inf,inf). > > Could this be done using R? I tried using the function integrate 2 times, but

Re: [R] Merging two data frames, but keeping NAs

Adding the argument all.x=TRUE to merge() will retain the NA values, but the only reliable way I've found to preserve order with NA values in a merge is to add an index column to x, merge the data, sort on the index column, then delete it. Sarah On Thu, Dec 5, 2013 at 9:56 AM, Rainer M Krug wrot

Re: [R] how to link different versions of R

> -Original Message- > I had downloaded the package zoo under R version 3.0.0. Afterwards i > downloaded R version 3.0.2 on the same workstation. Is it possible to use > the zoo package from 3.0.0. in 3.0.2. Usually a slightly older package version will work with a more recent version of

[R] Merging two data frames, but keeping NAs

-BEGIN PGP SIGNED MESSAGE- Hash: SHA1 Hi My brain is giving up on this... I have the following two data.frames: x <- data.frame(ref=c(NA, NA, NA, 10:5, NA, 1:5)) y <- data.frame(id = c(2, 3, 4, 6, 7, 9, 8), val = 101:107) Which look as follow: > x ref 1 NA 2 NA 3 NA 4

[R] how to link different versions of R

Hi, I had downloaded the package zoo under R version 3.0.0. Afterwards i downloaded R version 3.0.2 on the same workstation. Is it possible to use the zoo package from 3.0.0. in 3.0.2. Is is possible to upgrade from version 3.0.1 to 3.0.2 for example without downloading afresh and retain the

[R] p value for mu: anova()

Hello I have run an anova analysis for the fallowing model: H_obs=mu+REGION+MANAGEMENT + e When I run it in ASRelm I get the p-value for mu, and, of course also for the two dependent variables (REGION and MANAGEMENT) When I run it in R, I do not get the pvalue for mu. Can some one help me to

[R] GAM Assumption Tests

Dear all, currently I set up a GAM for my dataset (~32k records). I assume a normal distribution, constant variance and no correlation effects. With gam.check() it is possible to check those assumptions graphically. But is there also any option to do quantitative tests like the Wald-Test, shapi

[R] plot of a bagging tree

by considering this general example ##start code library(ipred) data("Ionosphere", package = "mlbench") Ionosphere$V2 <- NULL # constant within groups iono<-bagging(Class ~ ., data=Ionosphere, coob=TRUE) print(iono) ##end code does anybody knows any possibility to plot the (average) plot of th

Re: [R] O-ring statistic

This is the R-help mailing list, not Prof. Baddeley's personal email. If you want to direct a question to Prof. Baddeley (as is indeed appropriate in this instance since your question is about spatstat) then sent an email to Prof. Baddeley, not to this list. I *think* that the following will d

[R] O-ring statistic

Dear Prof. Baddeley, According to your suggestion (https://stat.ethz.ch/pipermail/r-help/2005-December/084274.html), I can run O-ring in library spatstat in R. Could you please told me how to draw the confidence envelop lines of O-ring function in R? Thanks! —————————