Re: [R] survey package -- doesn't appear to match svy

2014-10-28 Thread Anthony Damico
could you provide a minimal reproducible example? perhaps use ?dput. in general the survey package matches all other languages http://journal.r-project.org/archive/2009-2/RJournal_2009-2_Damico.pdf here's an example of a minimal reproducible example that does match

Re: [R] Trace of product of matrices

2014-10-28 Thread Martin Maechler
peter dalgaard pda...@gmail.com on Sun, 19 Oct 2014 21:26:39 +0200 writes: On 19 Oct 2014, at 19:00 , Spencer Graves spencer.gra...@structuremonitoring.com wrote: On 10/19/2014 8:42 AM, peter dalgaard wrote: On 19 Oct 2014, at 16:43 , Wagner Bonat wbo...@gmail.com

[R] how to structure data to use a cengaussian/cenpoisson distribution in MCMCglmm

2014-10-28 Thread David Villegas Ríos
Hi. I was wondering how data has to be structured in the database to be able to run a mixed-model with the MCMCglmm funcion using censored family distributions, like cengaussian or cenpoisson. I know my response variable should have two columns, but can anyone provide an example of this? In my

Re: [R] Trace of product of matrices

2014-10-28 Thread Berend Hasselman
On 28-10-2014, at 11:05, Martin Maechler maech...@stat.math.ethz.ch wrote: …... Thank you, Peter, and Spencer. For a few years now, I have had in my TODO file for the Matrix package: ** TODO tr(A %*% B) {and even tr(A %*% B %*% C) ...} are also needed frequently in some computations

[R] Putting R script in a function.

2014-10-28 Thread Frederic Ntirenganya
Hi All, I wrote this script to calculate the water balance using the following formula: Water_Balance = Water_Balance yesterday + Rainfall - Evaporation. The code works well and I want to put into a function. conditions: Water_Balance0 is equal to 0. Water_Balance100 is equal

[R] Nonparametric Estimation of Tail Dependence Coefficients

2014-10-28 Thread Gildas Mazo
Dear all, in which package can I find an implementation of the nonparametric estimation of tail dependence coefficients: lambda_L = lim_{u\to 0} P[F_1(X_1)u|F_2(X_2)u] lambda_U = lim_{u\to 1} P[F_1(X_1)u|F_2(X_2)u], where (X_1,X_2) has marginal distribution functions F_1 and F_2? (The

[R] breakpoints

2014-10-28 Thread Erasmo Papagni
Hello, Using the package strucchange to implement Bai-Perron methods I found that sometimes intervals between 2 break dates are exactly equal to the minimum that I set in the command breakpoints. I think this is a problem and a signal of possible different breaks with a trimming value

[R] spectral coherence non-overlapping windows: confidence intervals?

2014-10-28 Thread Martin Ivanov
Dear R users, The confidence intervals for the squared coherence and phase, plotted by plot.spec.coherence() and plot.spec.phase, respectively, use the formulae 10.14 and 10.11 from Bloomfield, P. (1976) Fourier Analysis of Time Series: An Introduction. Wiley: gg - 2/x$df se - sqrt(gg/2)  z -

Re: [R] breakpoints

2014-10-28 Thread Achim Zeileis
On Tue, 28 Oct 2014, Erasmo Papagni wrote: Hello, Using the package strucchange to implement Bai-Perron methods I found that sometimes intervals between 2 break dates are exactly equal to the minimum that I set in the command breakpoints. I think this is a problem and a signal of possible

Re: [R] Trace of product of matrices

2014-10-28 Thread Martin Maechler
Berend Hasselman b...@xs4all.nl on Tue, 28 Oct 2014 11:38:29 +0100 writes: On 28-10-2014, at 11:05, Martin Maechler maech...@stat.math.ethz.ch wrote: …... Thank you, Peter, and Spencer. For a few years now, I have had in my TODO file for the Matrix package:

Re: [R] plot hclust object

2014-10-28 Thread Martin Maechler
Greg Snow 538...@gmail.com on Mon, 27 Oct 2014 12:33:18 -0600 writes: I don't know of any tools that automate this process. For small sample sizes it may be easiest to just do this by hand, for large sample sizes that plot will probably be to complicated to make sense

[R] merge coefficients from a glmlist of models

2014-10-28 Thread Michael Friendly
In the vcdExtra package, I have a function glmlist to collect a set of glm() models as a glmlist object, and other functions that generate fit such a collection of models. This is my working example, fitting a set of models to the Donner data # install.packages(vcdExtra,

[R] Adding labels to ColSums

2014-10-28 Thread Muhuri, Pradip (SAMHSA/CBHSQ)
Hello, I was trying to add labels to the colSums of the integers variable corresponding to a factor. Below are the warning message and the reproducible code. How would I tweak the code to replace the NA with the Total in the output? Your advice toward resolving the issue would be greatly

Re: [R] merge coefficients from a glmlist of models

2014-10-28 Thread John Fox
Hi Michael, How about this? coef.glmlist - function(object, result=c(list, matrix, data.frame), ...){ result - match.arg(result) coefs - lapply(object, coef) if (result == list) return(coefs) coef.names - unique(unlist(lapply(coefs, names))) n.mods - length(object)

Re: [R] plot hclust object

2014-10-28 Thread Greg Snow
Thanks Martin, It is always great to learn that I don't need to reinvent the wheel (especially when I learn that before reinventing). Do you know if there are any help pages that point to cophenetic (see also or other sections). Maybe it is just the way that my brain is wired (along with being

Re: [R] merge coefficients from a glmlist of models

2014-10-28 Thread Michael Friendly
On 10/28/2014 12:13 PM, John Fox wrote: Hi Michael, How about this? That's perfect! And more general than I had hoped for. You probably used mindreader() :-) -Michael coef.glmlist - function(object, result=c(list, matrix, data.frame), ...){ result - match.arg(result) coefs -

Re: [R] Putting R script in a function.

2014-10-28 Thread Adams, Jean
Frederic, You can simplify your code somewhat using vectorization and the cumsum() function instead of a for loop. Since you focus in on a single year, perhaps you could do the calculations for all the years, and store the result in an array, instead. Then you could grab any day*month matrix

[R] R and Newest version of Java

2014-10-28 Thread Antonio Paredes
Hell All, Last night I updated Java to it newest version and this morning when I got to my office some of the R packages, that I am using in a current project, are not loading at all. For example Loading required package: XLConnectJars Error : .onLoad failed in loadNamespace() for 'rJava',

Re: [R] plot hclust object

2014-10-28 Thread Martin Maechler
Greg Snow 538...@gmail.com on Tue, 28 Oct 2014 10:31:27 -0600 writes: Thanks Martin, It is always great to learn that I don't need to reinvent the wheel (especially when I learn that before reinventing). Do you know if there are any help pages that point to cophenetic (see

Re: [R] R and Newest version of Java

2014-10-28 Thread Prof Brian Ripley
On 28/10/2014 16:47, Antonio Paredes wrote: Hell All, Last night I updated Java to it newest version and this morning when I got to my office some of the R packages, that I am using in a current project, are not loading at all. For example Loading required package: XLConnectJars Error :

Re: [R] plot hclust object

2014-10-28 Thread Greg Snow
I would suggest links to cophenetic on the help pages for dendrogram and possibly plot.hclust and related functions. I was not complaining. I always enjoy learning new things, it is a testament to the breadth and depth of R that even after more than 25 years using S and R, that I can still be

Re: [R] plot hclust object

2014-10-28 Thread Bert Gunter
Copacetic cophenetics are a way To better see much genetics. ;-) -- Bert Bert Gunter Genentech Nonclinical Biostatistics (650) 467-7374 Data is not information. Information is not knowledge. And knowledge is certainly not wisdom. Clifford Stoll On Tue, Oct 28, 2014 at 1:28 PM, Greg Snow

[R] big datasets for R

2014-10-28 Thread Tim Hoolihan
Another source of large datasets is the Public Data Sets on AWS http://aws.amazon.com/public-data-sets/ Tim Hoolihan @thoolihan http://linkedin.com/in/timhoolihan On Oct 28, 2014, at 7:00 AM, r-help-requ...@r-project.org wrote: -- Message: 2 Date: Mon, 27 Oct

[R] weighting histograms

2014-10-28 Thread Chris Jackson-Jordan
I am attempting to weight a histogram with another variable in the spreadsheet. does anyone know an easy way of doing this. I have been attempting to use weighted.hist from the plotrix package but I can't figure out how to use the other variable as the weight. I need to weight the Teton_Co_V

[R-es] Encuesta

2014-10-28 Thread web
  -- Puede desuscribirse de este boletín siguiendo este enlace: http://r-es.org/tiki-newsletters.php?unsubscribe=df02c5ccbde821d5c05c50c46c705946 [[alternative HTML version deleted]] ___ R-help-es mailing list R-help-es@r-project.org

[R-es] duda (URGENTE) R

2014-10-28 Thread Nicolás Corral Gómez
Hola quisiera saber de que manera puedo generar numeros aleatorios enteros entre 0 y 4, para completar una matriz de un automata celular Muchas gracias Juan Nicol�s Corral G�mezEstudiante Carrera de Ecolog�aPontificia Universidad Javeriana

Re: [R-es] duda (URGENTE) R

2014-10-28 Thread Isidro Hidalgo
Si quieres 100: sample(0:4, 100, replace = TRUE) Si quieres más, cambia la cantidad. Un saludo, Isidro Hidalgo Arellano Observatorio Regional de Empleo Consejería de Empleo y Economía http://www.jccm.es -Mensaje original- De: r-help-es-boun...@r-project.org

[R-es] predict en objetos creados por kknn

2014-10-28 Thread Raúl Vaquerizo
Compañeros, me siento superado por este problema. A ver si me podéis ayudar. No me funciona (o no entiendo) predict en objetos creados por kknn. Os pongo un ejemplo sencillo que trata de predecir un cuadrado en una nube de puntos en el plano. Muy sencillo viendo el código: #Datos #Se trata de una

Re: [R-es] predict en objetos creados por kknn

2014-10-28 Thread Jose Luis Cañadas Reche
Hola Raúl. No he visto tu código en profundidad, pero puede ser que en el último gráfico estás llamando a g, que es ggplot creado para los datos de entrenamiento y le pintes la predicción de los datos de test. Deberías pintar los datos de test y su predicción. Quizá esto te valga. g2 -