1k views

Variance on the sum of predicted values from a mixed effect model on a timeseries

I have a mixed effect model (in fact a generalized additive mixed model) that gives me predictions for a timeseries. To counter the autocorrelation, I use a corCAR1 model, given the fact I have ...
514 views

SVD of correlated matrix should be additive but doesn't appear to be

I'm just trying to replicate a claim made in the following paper, Finding Correlated Biclusters from Gene Expression Data, which is: Proposition 4. If $X_{IJ}=R_{I}C^{T}_{J}$. then we have: ...
414 views

Stability of cross-validation in Bayesian models

I'm fitting a Bayesian HLM in JAGS using k-fold cross-validation (k=5). I'd like to know whether estimates of parameter $\beta$ are stable across all folds. What's the best way to do this? One idea ...
661 views

EM maximum likelihood estimation for Weibull distribution

Note: I am posting a question from a former student of mine unable to post on his own for technical reasons. Given an iid sample $x_1,\ldots,x_n$ from a Weibull distribution with pdf  f(x) = k ...
210 views

When is binomial distribution function above/below its limiting Poisson distribution function?

Let $B(n,p,r)$ denote the binomial distribution function (DF) with parameters $n \in \mathbb N$ and $p \in (0,1)$ evaluated at $r \in \{0,1,\ldots,n\}$: B(n,p,r) = \sum_{i=0}^r ...
612 views

Gaussian Ratio Distribution: Derivatives wrt underlying $\mu$'s and $\sigma^2$s

I'm working with two independent normal distributions $X$ and $Y$, with means $\mu_x$ and $\mu_y$ and variances $\sigma^2_x$ and $\sigma^2_y$. I'm interested in the distribution of their ratio ...
102 views

Training a basic Markov Random Field for classifying pixels in an image

I am attempting to learn how to use Markov Random Fields to segment regions in an image. I do not understand some of the parameters in the MRF or why the expectation maximisation I perform fails to ...
169 views

Estimating variance of center-censored Normal samples

I have normally-distributed processes from which I get small samples (n typically 10-30) that I want to use to estimate variance. But frequently the samples are so close together that we can't ...
680 views

Brant test in R

In testing the parallel regression assumption in ordinal logistic regression I find there are several approaches. I've used both the graphical approach (as detailed in Harrell´s book) and the approach ...
210 views

Can the Mantel test be extended to asymmetric matrices?

The Mantel test is usually applied to symmetric distance/difference matrices. As far as I understand, an assumption of the test is that the measure used to define differences must be at least a ...
3k views

How to estimate variance components with lmer for models with random effects and compare them with lme results

I performed an experiment where I raised different families coming from two different source populations, where each family was split up into a different treatments. After the experiment I measured ...
205 views

Phylogenetic dependent variables: ANOVA?

I understand deriving a covariance matrix from phylogenetic data to make $cov(X,Y) = 0$ for two variables you're making a regression on. But what happens if you have one continuous variable, that ...
500 views

How to test equality of variances with circular data

I am interested in comparing the amount of variability within 8 different samples (each from a different population). I am aware that this can be done by several methods with ratio data: F-test ...
Tail bounds on Euclidean norm for uniform distribution on $\{-n,-(n-1),…,n-1,n\}^d$
What are known upper bounds on how often the Euclidean norm of a uniformly chosen element of $\:\{-n,~-(n-1),~...,~n-1,~n\}^d\:$ will be larger than a given threshold? I'm mainly interested in bounds ...