Re: [R] Warning message with maxLik()

Dear Maram - Please do not start a new thread for the same issue but reply to previous messages in this thread [1]. - Please read my previous responses [1] more carefully, e.g. to use "theta <- exp( param )" which guarantees that all elements of "theta" are always positive. [1] http://r.789695.

Re: [R] matching strings in a list

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 UNSUBSCRIBE

[R] modifying a package installed via GitHub

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, convers

[R] Warning message with maxLik()

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

Re: [R] Nagelkerke Pseudo R-squared

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

[R] Nagelkerke Pseudo R-squared

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 ex

Re: [R] Displaying Compositional Data With 46 Parts

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 displa

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

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

Re: [R] Displaying Compositional Data With 46 Parts

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

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 projec

[R] Release of R 3.2.2 scheduled for August 14

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

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

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

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] powerTransform warning message?

Thank you so much for the explanation. That was very helpful! :-) Thanks! Brittany > On Jul 16, 2015, at 6:16 PM, John Fox wrote: > > Dear Brittany, > > On Thu, 16 Jul 2015 17:35:38 -0600 > Brittany Demmitt wrote: >> Hello, >> >> I have a series of 40 variables that I am trying to trans

Re: [R] "grid"Package not available in r

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 tw

Re: [R] Hausman Test trouble - plm

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] "grid"Package not available in r

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

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

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 Nabble.co

Re: [R] User defined function within a formula

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 wrote: > This might do what you want: > > OPoly <- function(x, degree=1, weight=1, coefs=NULL, rangeX=NULL){ > weight <- round(weight,0)# weig

Re: [R] TukeyHSD troubles

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 p

Re: [R] Displaying Compositional Data With 46 Parts

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

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 log-scalin

Re: [R] Displaying Compositional Data With 46 Parts

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

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 of

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

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 wrote: > I have ma matrix which its elements are NA,0,1,2 ! I got my answer bout > removing the columns with 0 or NA

[R] Hausman Test trouble - plm

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.15

Re: [R] User defined function within a formula

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 if(length(weight)!=length(x)