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 ...
509 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: ...
413 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 ...
657 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 ...
282 views

### How much smaller can $p$ values from ANOVA's $F$-test be vs. those from multiple $t$-tests on the same data?

Intro: Having noted the attention received today by this question, "Can ANOVA be significant when none of the pairwise t-tests is?," I thought I might be able to reframe it in an interesting way that ...
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 ...
609 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 ...
100 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 ...
168 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 ...
673 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 ...
208 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 ...
### 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 ...