[R] order by which the eigenvalues are presented

2015-05-06 Thread Luis Borda de Agua
Thank you, Bill, for your reply. However, I'm afraid I didn't explain myself properly. Imagine you have a 2x2 matrix Then the eigenvalues lambda_1 and lambda_2 are analytically calculated from lambda_1 = (-b+sqrt(delta))/2a lambda_2 = (-b-sqrt(delta))/2a where delta = b^2-4ac If

[R] order by which the eigenvalues are presented

2015-05-05 Thread Luis Borda de Agua
I’m using R 3.1.3 in OS X 10.10.3 (Yosemite) I use the function “eigen” to calculate the eigenvalues of matrix where each element is sampled from a given distribution (normal, beta, etc.). According to the information provided: values a vector containing the p eigenvalues of x, sorted in

Re: [R] keep information on the number of warnings

2014-08-06 Thread Luis Borda de Agua
Thank you very much, David, Luís On 04 Aug 2014, at 17:24, Luis Borda de Agua lbagua.cl...@gmail.com wrote: Dear David Thank you very much for your reply. I’ve only seen it now. I tried length(warnings) and I got a strange result. When I used lw - length(warnings) print(lw) I

[R] keep information on the number of warnings

2014-08-05 Thread Luis Borda de Agua
Dear David, Once again, thank you very much for your email. I believe I understand why I got a different number of warnings”. It seems that when you type length(warnings) you keep the information on the run before the last. I come up with a simple example called testing (you will need to

[R] keep information on the number of warnings

2014-08-04 Thread Luis Borda de Agua
Dear David Thank you very much for your reply. I’ve only seen it now. I tried length(warnings) and I got a strange result. When I used lw - length(warnings) print(lw) I obtained lw=36 however, the number of warnings was 38 according to message to screen: There were 38 warnings (use

[R] keep information on the number of warnings

2014-07-31 Thread Luis Borda de Agua
I’m using R 3.1.0 in OS X 10.9.4 (Mavericks) I’m running a function Y that calls a function X which occasionally generates a warning. Say that I call the function X 1000 times, and out of these 1000 times I get the following message: There were 36 warnings (use warnings() to see them) How

[R] fitting a Bi-Weibull

2012-06-25 Thread Luis Borda de Agua
I want to fit a bi-Weibull distribution, as defined by Evans et al. in their book Statistical distributions: http://www.scribd.com/doc/57260630/164/FIVE-PARAMETER-BI-WEIBULL-DISTRIBUTION I would like to know if there is any package in R that already does it, or any quick procedure to

[R] how to extract information from anova results

2010-07-13 Thread Luis Borda de Agua
Hi, I have used the instruction aov in the following manner: res - aov(qwe ~ asd) when I typed res I get: _ Call: aov(formula = qwe ~ asd) Terms: asd Residuals Sum of Squares 0.0708704 0.5255957 Deg. of Freedom 1 8 Residual standard error:

[R] changing the number of elements in a list

2010-06-09 Thread Luis Borda de Agua
I would like to have a list where each element is a matrix, for example: my.list - list(matrix(0, ncol=3, nrow=3), matrix(0, ncol=3, nrow=3), matrix(0, ncol=3, nrow=3)) The problem is, I would like to be able to change automatically the number of