Re: [R] Fitting a curve to weibull distribution in R using nls

2015-10-13 Thread Aditya Bhatia
Yes. I do.I'm trying to repeat the methodology of a paper. They have fitted their data to a weibull curve and so I want to do the same too, but I'm unable to figure out how.. On Wed, Oct 14, 2015, 9:44 AM David Winsemius wrote: > > On Oct 13, 2015, at 2:42 PM, Aditya Bhatia wrote: > > > I am try

Re: [R] unsuscribe

2015-10-13 Thread Jeff Newmiller
We cannot do that. Please read the footer of any email on the list, such as quoted here. >R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see >https://stat.ethz.ch/mailman/listinfo/r-help >PLEASE do read the posting guide >http://www.R-project.org/posting-guide.html >and provide com

Re: [R] Fitting a curve to weibull distribution in R using nls

2015-10-13 Thread David Winsemius
On Oct 13, 2015, at 2:42 PM, Aditya Bhatia wrote: > I am trying to fit this data to a weibull distribution: > > My y variable is:1 1 1 4 7 20 7 14 19 15 18 3 4 1 3 1 1 1 > 1 1 1 1 1 1 > > and x variable is:1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 > 19 20 21 22 23 24

[R] unsuscribe

2015-10-13 Thread Dorota Buczek
Please unsuscribe me from the mailing list. Thank you Dorota __ R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html

Re: [R] heritability parameter

2015-10-13 Thread Peter Alspach
Tena koe Ramendra There are various packages (not libraries) to fit mixed models, and hence the individual animal (or plant) model - perhaps the best known being lme4. There's also MCMCglmm which takes a Baseyian approach, and asreml which is licence managed. This is nothing like an exclusive

[R] Fitting a curve to weibull distribution in R using nls

2015-10-13 Thread Aditya Bhatia
I am trying to fit this data to a weibull distribution: My y variable is:1 1 1 4 7 20 7 14 19 15 18 3 4 1 3 1 1 1 1 1 1 1 1 1 and x variable is:1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 The plot looks like this:http://i.stack.imgur.com/FrIKo.png and

[R] contrasts among simple effects - 2

2015-10-13 Thread James Henson
Greetings R Community Apologize for previously sending a csv file. My goal is to make orthogonal contrasts among simple effects in analysis of repeated measures data. The SAS publication, on page 1224, shows how to make this type of contrasts in SAS. But, my search of books about repeated measu

Re: [R] k-nearest neighbour classification

2015-10-13 Thread Adams, Jean
Ah, I see what you're after. I don't know of a built in function to search for the best number of nearest neighbors. You may have to run the code for each k separately, then compare the resulting errors. Jean On Tue, Oct 13, 2015 at 9:13 AM, Neverstop wrote: > I know that knn.cv(train=predict

Re: [R] how to do away for loop using functionals?

2015-10-13 Thread Michael Hannon
The answer to "another note" is: mapply(rep, w, 5:1) I'll try to look at the first part in more detail later today. -- Mike On Mon, Oct 12, 2015 at 5:55 PM, Annie Hawk via R-help wrote: > HI R-experts, > > > I am trying to speed up my calculation of the A results below and replace the > for

Re: [R] quantile regression: warning message

2015-10-13 Thread T.Riedle
Thank you very much. So, the results are correct and the differences between the solutions are pretty small. Thus, I do not need to worry about the warning message? Yes? -Original Message- From: Roger Koenker [mailto:rkoen...@illinois.edu] Sent: 13 October 2015 21:59 To: T.Riedle Cc: R

Re: [R] R help

2015-10-13 Thread Jeff Newmiller
In general you should read the Posting Guide mentioned at the bottom of every email, which points out that you should post in plain text. Figuring out how to send plain text email in your email client will help you make sure that we see what you see. You should also try to insert reproducible

Re: [R] quantile regression: warning message

2015-10-13 Thread Roger Koenker
see the output from the quantreg FAQ: FAQ() especially point 2. url:www.econ.uiuc.edu/~rogerRoger Koenker emailrkoen...@uiuc.eduDepartment of Economics vox: 217-333-4558University of Illinois fax: 217-244-6678Urba

[R] quantile regression: warning message

2015-10-13 Thread T.Riedle
Greetings R Community, I am trying to run a quantile regression using the quantreg package. My code looks as follows: RegressionUtilitiesUK<-rq(ReturnUtilities~yield.spread.change+ReturnFTSE, tau=0.01,data=State_variables_UK_calm) Unfortunately, the summary() function returns the results but als

Re: [R] R help

2015-10-13 Thread Sarah Goslee
Hi, First of all, you almost certainly want to get your data into R with stringsAsFactors=FALSE rather than creating a factor with 345196 levels. mydata <- read.table("whateverfile", stringsAsFactors=FALSE) I'm just guessing, but you probably also don't really want filter(), but subset() instead

Re: [R] contrasts among simple effects

2015-10-13 Thread David Winsemius
On Oct 13, 2015, at 11:16 AM, James Henson wrote: > Greetings R Community > > My goal is to make orthogonal contrasts among simple effects in analysis of > repeated measures data. The SAS publication, on page 1224, shows how to > make this type of contrasts in SAS. But, my search of books abou

Re: [R] Parsing XML File

2015-10-13 Thread Lorenzo Isella
Dear Jim, Thanks for your reply. What you did is 100% what I need -- I now have a data frame with the relevant data and I can take up from there. Regards Lorenzo On Sun, Oct 11, 2015 at 03:54:10PM -0400, jim holtman wrote: Not sure exactly what you want since you did not show an expected output

[R] class-specific Gini metrics for Random Forest?

2015-10-13 Thread Matt Fagan
library(randomForest) data(iris) fit <- randomForest(Species ~ ., data=iris, importance=TRUE); fit.imp<-importance(fit) fit.imp columns 1-3 of fit.imp show the class-specific variable importance for the Mean Decrease Acuracy measure (MDA). Is there a way to calculate class-specific Gini metrics ra

[R] heritability parameter

2015-10-13 Thread Ramendra Sarma
Is there in any library in R to calculate gcv, pcv and heritability parameters in R? -- *Dr R N SarmaProfessorDepartment of Plant Breeding and GeneticsAss

Re: [R] Problem with binomial gam{mgcv}

2015-10-13 Thread Simon Wood
Following up on Duncan's answer... You can see that the model doesn't fit with a little bit of residual checking. For example gam.check(Model) shows that the binomial assumption is just wrong here: there is massive overdispersion, for example. Even the crudest response to this, of using qua

Re: [R] k-nearest neighbour classification

2015-10-13 Thread Neverstop
I know that knn.cv(train=predictors.training, cl=classes.training, k=3, prob=TRUE) works but by doing so I fix the tuning paramer k to be 3. Isn't cross validation a technique to choose the optimal tuning parameter trying a range of different values for the tuning parameter? -- View this messag

[R] R help

2015-10-13 Thread Clayton Samples
Hello R community, I need a bit of direction. I am new to R, so please forgive any mistakes or confusion. Here is an example of the columns and data contained in my data frame. Number: Factor w 345196 levels "1.001", "1.002", "1.003", "2.001", "2.002", "2.003" Name: factor w 2 levels ("X", "Y")

[R] NaNs produced running icc()

2015-10-13 Thread Chad Danyluck
Hello, I want to determine the inter-rater reliability of ratings made from a random selection of 11 raters and 71 subjects using the irr package. Each rater was randomly assigned to rater only a handful of subjects. So there will be some cells with NAs for raters who did not rate a particular per

[R] contrasts among simple effects

2015-10-13 Thread James Henson
Greetings R Community My goal is to make orthogonal contrasts among simple effects in analysis of repeated measures data. The SAS publication, on page 1224, shows how to make this type of contrasts in SAS. But, my search of books about repeated measures analysis using R, and on-line has not yiel

Re: [R] Reformatting dataframe for use with icc()

2015-10-13 Thread Chad Danyluck
Indeed, I've worked with the solutions provided by others and the reshape package has been very fluid to use. Thanks for the suggestion. Chad On Sat, Oct 10, 2015 at 4:56 PM, Frans Marcelissen < fransiepansiekever...@gmail.com> wrote: > I think this is what reshape is made for... > Frans > -

Re: [R] 3D matrix columns messed up _ looking for your help

2015-10-13 Thread Kristi Glover
Hi Jim, Thank you very much for the message. Sorry for the email that was not clear. Yes, you are right, in A$date, should be t2> ENM, t3> ENL. but using the code it gave t2>ENL, t3>ENM. I am struggling to fix it. If the data set was small, I could do it manually. would you mind to try this e

[R] var.gamma in the prop.odds funtion (timereg package)

2015-10-13 Thread lspirk
Hi R users! I am trying to calculate the variance of the parameter under the null (Ho) using the "prop.odds" function in the timereg package. I know var.gamma will give the variance of the parameter, but how do I find this under Ho? For the Cox PH model, I used iter = 0 and the "vcov" function:

Re: [R] k-nearest neighbour classification

2015-10-13 Thread Adams, Jean
The argument k should be a scalar, not a vector. So, for example, this works: knn.cv(train=predictors.training, cl=classes.training, k=3, prob=TRUE) Jean On Tue, Oct 13, 2015 at 3:59 AM, Neverstop wrote: > Hi, I'm trying to perform a cross validation to choose the optimal k in the > k-neares

Re: [R] Windows release R 3.2.2 : internet.dll module

2015-10-13 Thread Duncan Murdoch
On 13/10/2015 4:06 AM, Philippe DESSEN wrote: > Hi > > With the Windows release R-3.2.2, I have a problem of internet.dll module > which prevents the use of automatic package loading in applications and > other access to internet (ex: BrB Array Tools ) > There was no problem with the previous r

[R] Windows release R 3.2.2 : internet.dll module

2015-10-13 Thread Philippe DESSEN
Hi With the Windows release R-3.2.2, I have a problem of internet.dll module which prevents the use of automatic package loading in applications and other access to internet (ex: BrB Array Tools ) There was no problem with the previous release 3.2.0 on the same DELL Windows 7 machine. Is some

[R] k-nearest neighbour classification

2015-10-13 Thread Neverstop
Hi, I'm trying to perform a cross validation to choose the optimal k in the k-nearest-neighbors algorithm for classification. I'm using the knn function of the package class. Reading the R documentation, I've found out that there

Re: [R] for loop not working

2015-10-13 Thread mnw
Thank you for your comments. I guess what I need is to store movements and then loop through that matrix. Here is a snippet of the 'movement' code. I am trying to fill the matrix below with data from this loop. Is there something else my matrix (bottom) requires? Or its positioning? Many thanks,

Re: [R] 3D matrix columns messed up _ looking for your help

2015-10-13 Thread Jim Lemon
Hi Kristi, This is a bit hard to follow, but I'll try. As you are replacing the numeric values of the intermediate table with the character values of the factor A$date, it looks to me as though the answer is as it should be. 2 -> ENL, 3 -> ENM. I suspect that the solution is not difficult, but I ca

Re: [R] using the apply() family to evaluate nested functions with common arguments

2015-10-13 Thread Maram SAlem
I'm sorry, I forgot to mention that in the last step I have to use the choose() function for all the elements of x not only x[1],x[2], and x[3]. Any suggestions or recommendation for some reference or a package that can help me sort that out. Thanks. Maram Salem On 11 October 2015 at 18:52, Mar