Re: [R] Displaying Compositional Data With 46 Parts

2015-07-17 Thread Rich Shepard
On Fri, 17 Jul 2015, Aaron Mackey wrote: One immediate question is how independent you believe the 46 components to be, and whether certain components could be reduced or otherwise coordinately-modeled; a heatmap of your 46x46 pairwise correlations should be informative. Also consider

Re: [R] Displaying Compositional Data With 46 Parts

2015-07-17 Thread Rich Shepard
On Fri, 17 Jul 2015, John Kane wrote: I think this is more a technical question for the subject matter experts than for R-help if I am understanding the question correctly. John, I agree completely. Unfortunately, there is no R SIG devoted to CoDA, nor any other mail list or Web forum that

Re: [R] Displaying Compositional Data With 46 Parts

2015-07-17 Thread John Kane
Hi Rich, I think this is more a technical question for the subject matter experts than for R-help if I am understanding the question correctly. That said, there seems to be an R package called compositional and a corresponding book http://www.springer.com/us/book/9783642368080 that may help.

[R] Displaying Compositional Data With 46 Parts

2015-07-17 Thread Rich Shepard
The compositional data have been divided into two data frames: 46 response variables (the compositional components) and 5 explanatory variables. There are 209 observations of each. With no experience analyzing large compositions with so many parts your advice on how to plot and report results

Re: [R] TukeyHSD troubles

2015-07-17 Thread tdiver
Bart, I want to thank you for your code. I was having similar problems as Amy, even after setting my numeric variable as a factor using as.factor(). I used is.factor() to confirm and received the answer as TRUE from R; however after running the TukeyHSD() my set factor in my aov() was not read

[R] Hausman Test trouble - plm

2015-07-17 Thread TDix
Hi there. I am a student / very fresh R user who is currently having some issues running the procedure for a Hausman test in R. The head for my data sheet named data looks like this: BirdSeason Gully Grouping Food Habitat.Type 1 83 111 0.152

Re: [R] User defined function within a formula

2015-07-17 Thread Kunshan Yin
Thanks Bill for your quick reply. I tried your solution and it did work for the simple user defined function xploly. But when I try with other function, it gave me error again: OPoly-function(x,degree=1,weight=1){ weight=round(weight,0)# weight need to be integer

Re: [R] installation of package ‘devtools’ had non-zero exit status in a POWERPC

2015-07-17 Thread aidaph
I'm stuck on this and i don't know how to install Curl without any problem. Anyone? -- View this message in context: http://r.789695.n4.nabble.com/installation-of-package-devtools-had-non-zero-exit-status-in-a-POWERPC-tp4709687p4709994.html Sent from the R help mailing list archive at

Re: [R] Hausman Test trouble - plm

2015-07-17 Thread TDix
Might have just solved my own problem team! I assumed that the issue here was the replicated samples, and so added a column and gave a number to each replicate. R seemed to like this and was happy to run the test! A significant result tells me that the fixed effects model is the most preferable

[R] gridPackage not available in r

2015-07-17 Thread Issoufou Ouedraogo
Dear Responsible, Hello! I am a PhD student at Université Catholique de Louvain, I contact you because I have some difficulties to install neuralnet package in R for my research. However, to install this package, we need two packages MASS and grid. The grid package is not available in my R

Re: [R] User defined function within a formula

2015-07-17 Thread Kunshan Yin
Thank you very much. It worked. I think I need to digest this further later. Thanks again for the help. On Thu, Jul 16, 2015 at 4:51 PM, William Dunlap wdun...@tibco.com wrote: This might do what you want: OPoly - function(x, degree=1, weight=1, coefs=NULL, rangeX=NULL){ weight -

Re: [R] installation of package ‘devtools’ had non-zero exit status in a POWERPC

2015-07-17 Thread Jeff Newmiller
That is an operating-system configuration question, not a question about R. There are many OSs out there... please find a forum with users of your OS in which to pursue this question. --- Jeff Newmiller

Re: [R] Displaying Compositional Data With 46 Parts

2015-07-17 Thread John Kane
See in-line John Kane Kingston ON Canada -Original Message- From: rshep...@appl-ecosys.com Sent: Fri, 17 Jul 2015 07:33:43 -0700 (PDT) To: r-help@r-project.org Subject: Re: [R] Displaying Compositional Data With 46 Parts On Fri, 17 Jul 2015, John Kane wrote: I think this is

Re: [R] powerTransform warning message?

2015-07-17 Thread Brittany Demmitt
Thank you so much for the explanation. That was very helpful! :-) Thanks! Brittany On Jul 16, 2015, at 6:16 PM, John Fox j...@mcmaster.ca wrote: Dear Brittany, On Thu, 16 Jul 2015 17:35:38 -0600 Brittany Demmitt demmi...@gmail.com wrote: Hello, I have a series of 40 variables

Re: [R] removing the columns with 0 or NA or 1or NA or 2 or NA

2015-07-17 Thread Jim Lemon
Hi Lida, I think that your matrix is actually a data frame, so try this: mat[,sapply(mat,function(x) var(x,na.rm=TRUE)0)] Jim On Fri, Jul 17, 2015 at 12:58 AM, Lida Zeighami lid.z...@gmail.com wrote: I have ma matrix which its elements are NA,0,1,2 ! I got my answer bout removing the columns

[R] Release of R 3.2.2 scheduled for August 14

2015-07-17 Thread Peter Dalgaard
We intend to have a patch release on August 14, nickname will be Fire Safety. The detailed schedule will be made available via developer.r-project.org as usual. -- Peter Dalgaard, Professor, Center for Statistics, Copenhagen Business School Solbjerg Plads 3, 2000 Frederiksberg, Denmark Phone:

Re: [R] Displaying Compositional Data With 46 Parts

2015-07-17 Thread Rich Shepard
On Fri, 17 Jul 2015, Bert Gunter wrote: I believe John Aitchison's book and papers are the authoritative basic resources. Have you read them? Bert, Yes, I have. The problem is that the support of the distributions are (hyper)simplexes, not Euclidean space, due to the requirement that the

Re: [R] Displaying Compositional Data With 46 Parts

2015-07-17 Thread Rich Shepard
On Fri, 17 Jul 2015, John Kane wrote: Then it sounds like you are one of the experts. Do whatever you think appropriate and either set the standard for future research or get enough feedback to do even better next time. :) John, Far from an expert, but becoming more capable with each

Re: [R] gridPackage not available in r

2015-07-17 Thread David Winsemius
On Jul 17, 2015, at 7:21 AM, Issoufou Ouedraogo wrote: Dear Responsible, Hello! I am a PhD student at Université Catholique de Louvain, I contact you because I have some difficulties to install neuralnet package in R for my research. However, to install this package, we need two

[R] OPTIMX: non-finite finite-difference value [26] and scaling problem

2015-07-17 Thread Olu Ola via R-help
Hello,I am running a nonlinear GMM using the optimx wrapper. I am trying to estimate 37 variables however and my code for the optimx is: nlgmm = optimx(par=b0, fn=obj,method = BFGS, itnmax=1, control=list(follow.on = TRUE,kkt=FALSE,starttests=TRUE,save.failures=TRUE, trace=0)) My staring

[R] Nagelkerke Pseudo R-squared

2015-07-17 Thread varin sacha
Dear R-Experts, I have fitted an ordinal logistic regression with just 1 explanatory variable for the reproducible example here below. Everything is working, now I try to calculate the Nagelkerke Pseudo R-squared. I have found a package BaylorEdPsych providing many Pseudo R-squared, but the

[R] Warning message with maxLik()

2015-07-17 Thread Maram SAlem
Dear All, I'm trying to get the MLe for a certain distribution using maxLik () function. I wrote the log-likelihood function as follows: theta -vector(mode = numeric, length = 3) r- 17 n -30 T-c(7.048,0.743,2.404,1.374,2.233,1.52,23.531,5.182,4.502,1.362,1.15,1.86,1.692,11.659,1.631,2.212,5.451)

Re: [R] matching strings in a list

2015-07-17 Thread tryingtolearn
Thank you all very much! -- View this message in context: http://r.789695.n4.nabble.com/matching-strings-in-a-list-tp4709967p4710015.html Sent from the R help mailing list archive at Nabble.com. __ R-help@r-project.org mailing list -- To

[R] modifying a package installed via GitHub

2015-07-17 Thread Steve E.
Hi Folks, I am working with a package installed via GitHub that I would like to modify. However, I am not sure how I would go about loading a 'local' version of the package after I have modified it, and whether that process would including uninstalling the original unmodified package (and,

Re: [R] Displaying Compositional Data With 46 Parts

2015-07-17 Thread Jim Lemon
Hi Rich, Being in a position of relative ignorance on this topic, I'll offer some suggestions that may well be useless. You mention ternary diagrams, which use position to represent compositional proportions. These will not scale up to 46 values in any way that I can imagine. If you want to

Re: [R] Nagelkerke Pseudo R-squared

2015-07-17 Thread David Winsemius
On Jul 17, 2015, at 4:33 PM, varin sacha wrote: Dear R-Experts, I have fitted an ordinal logistic regression with just 1 explanatory variable for the reproducible example here below. Everything is working, now I try to calculate the Nagelkerke Pseudo R-squared. I have found a