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
http://www.ats.ucla.edu/stat
> peter dalgaard
> on Sun, 19 Oct 2014 21:26:39 +0200 writes:
>> On 19 Oct 2014, at 19:00 , Spencer Graves
wrote:
>>
>> On 10/19/2014 8:42 AM, peter dalgaard wrote:
On 19 Oct 2014, at 16:43 , Wagner Bonat wrote:
Dear,
I ha
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
On 28-10-2014, at 11:05, Martin Maechler 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 {conditional normal
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_Balance<0 is equal to 0.
Water_Balance>100 is equal
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],
where (X_1,X_2) has marginal distribution functions F_1 and F_2?
(The nonparametric estimators typically require to choose a
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 low
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 <
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
> Berend Hasselman
> on Tue, 28 Oct 2014 11:38:29 +0100 writes:
> On 28-10-2014, at 11:05, Martin Maechler
wrote:
>> …...
>> Thank you, Peter, and Spencer.
>>
>> For a few years now, I have had in my TODO file for the Matrix
>> package:
>>
>> ** TO
> 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 mak
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", repos="http://R-F
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
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(objec
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 a
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)
coe
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 fro
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', detail
> 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 co
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 : .onLoad
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 ple
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
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 Oc
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 colum
24 matches
Mail list logo
