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

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] unsuscribe

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] Fitting a curve to weibull distribution in R using nls

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

Re: [R] unsuscribe

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

Re: [R] heritability parameter

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

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

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

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

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

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

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,

[R] k-nearest neighbour classification

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

Re: [R] for loop not working

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] Windows release R 3.2.2 : internet.dll module

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

Re: [R] k-nearest neighbour classification

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

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

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

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

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] 3D matrix columns messed up _ looking for your help

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

Re: [R] Parsing XML File

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

Re: [R] R help

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

[R] quantile regression: warning message

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

Re: [R] quantile regression: warning message

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-6678

Re: [R] contrasts among simple effects

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

Re: [R] R help

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()

[R] contrasts among simple effects

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

[R] R help

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")

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

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

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

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

[R] NaNs produced running icc()

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

Re: [R] k-nearest neighbour classification

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

[R] heritability parameter

Is there in any library in R to calculate gcv, pcv and heritability parameters in R? -- *Dr R N SarmaProfessorDepartment of Plant Breeding and

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

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] quantile regression: warning message

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] contrasts among simple effects - 2

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

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

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

Re: [R] k-nearest neighbour classification

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