20
votes
0answers
751 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: ...
11
votes
0answers
1k views

Inverting the Fourier Transform for a Fisher distribution

The characteristic function of Fisher $\mathcal{F}(1,\alpha)$ distribution is: $$C(t)=\frac{\Gamma \left(\frac{\alpha +1}{2}\right) U\left(\frac{1}{2},1-\frac{\alpha }{2},-i t \alpha \right)}{\Gamma ...
10
votes
0answers
228 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 ...
10
votes
0answers
127 views

Testing certain contrasts: Is this provably a hard problem, or not?

I posted this to mathoverflow and no one's answering: Scheffé's method for identifying statistically significant contrasts is widely known. A contrast among the means $\mu_i$, $i=1,\ldots,r$ of $r$ ...
10
votes
0answers
1k views

Random forests for multivariate regression

I have a multi-output regression problem with $d_x$ input features and $d_y$ outputs. The outputs have a complex, non-linear correlation structure. I'd like to use random forests to do the ...
9
votes
0answers
207 views

AIC/BIC: how many parameters does a permutation count for?

Let's say I have a model selection problem and I am trying to use AIC or BIC to evaluate the models. This is straightforward for models that have some number $k$ of real-valued parameters. However, ...
9
votes
0answers
344 views

What are the options in proportional hazard regression model when Schoenfeld residuals are not good?

I am doing a Cox proportional hazards regression in R using coxph, which includes many variables. The Martingale residuals look great, and the Schoenfeld residuals ...
9
votes
0answers
826 views

Intraclass Correlation Coefficients (ICC) with Multiple Variables

Suppose I have measured some variable in siblings, which are nested within families. The data structure looks like this: family sibling value ------ ------- ----- 1 1 y_11 1 2 ...
8
votes
0answers
242 views

Advanced regression modeling examples

I'm looking for an advanced linear regression case study illustrating the steps required to model complex, multiple non-linear relationships using GLM or OLS. It is surprisingly difficult to find ...
8
votes
0answers
349 views

Can the Burnham-Anderson book on multimodel inference be recommended?

As motivated by the recent change of the default model selection statistic in the R's forecast package from AIC to AICc, I am curious whether the latter is indeed applicable wherever the former is. ...
8
votes
0answers
483 views

Clustered standard errors vs. multilevel modeling?

I've skimmed through several books (Raudenbush & Bryk, Snijders & Bosker, Gelman & Hill, etc.) and several articles (Gelman, Jusko, Primo & Jacobsmeier, etc.), and I still haven't ...
8
votes
0answers
152 views

Penalized spline confidence intervals based on cluster-sandwich VCV

This is my first post here, but I've benefited a lot from this forum's results popping up in google search results. I've been teaching myself semi-parametric regression using penalized splines. ...
8
votes
0answers
303 views

Physical/pictoral interpretation of higher-order moments

I'm preparing a presentation about parallel statistics. I plan to illustrate the formulas for distributed computation of the mean and variance with examples involving center of gravity and moment of ...
7
votes
0answers
150 views

Blind source separation of convex mixture?

Suppose I have $n$ independent sources, $X_1, X_2, ..., X_n$ and I observe $m$ convex mixtures: \begin{align} Y_1 &= a_{11}X_1 + a_{12}X_2 + \cdots + a_{1n}X_n\\ ...&\\ Y_m &= a_{m1}X_1 + ...
7
votes
0answers
74 views

Are data handling errors already 'priced in' to statistical analysis?

Ok, fair warning--this is a philosophical question that involves no numbers. I've been thinking a lot about how errors creep into data sets over time and how that should be treated by analysts--or if ...
7
votes
0answers
120 views

Why does the scikit-learn bootstrap function resample the test set?

When using bootstrapping for model evaluation, I always thought the out-of-bag samples were directly used as a test set. However, this appears not to be the case for the scikit-learn bootstrap ...
7
votes
0answers
231 views

Model selection with Firth logistic regression

In a small data set ($n\sim100$ ) that I am working with, several variables give me perfect prediction/separation. I thus use Firth logistic regression to deal with the issue. If I select the best ...
7
votes
0answers
416 views

What is the difference between conditional and unconditional quantile regression?

The conditional quantile regression estimator by Koenker and Basset (1978) for the $\tau^{th}$ quantile is defined as $$\widehat{\beta}_{QR} = \min_{b} \sum^{n}_{i=1} \rho_\tau (y_i - X'_i b_\tau)$$ ...
7
votes
0answers
2k views

What does the anova() command do with a lmer model object?

Hopefully this is a question that someone here can answer for me on the nature of decomposing sums of squares from a mixed-effects model fit with lmer (from the ...
7
votes
0answers
163 views

Bound for Arithmetic Harmonic mean inequality for matrices?

NOTE: This question has originally been posted in MSE, but it did not generate any interest. It was first posted there, because the question itself is a pure matrix-algebra question. Nevertheless, ...
7
votes
0answers
218 views

p-value adjustment for Local Moran's I statistic (LISA)

I'm working with some exploratory spatial analysis in R using spdep package. I came across an option to adjust p-values of local indicators of spatial association (LISA) calculated using ...
7
votes
0answers
287 views

Unique (?) idea for forecasting sales

I'm working on developing a model to predict total sales of a product. I have about a year and a half of bookings data, so I could do a standard time series analysis. However, I also have a lot of ...
7
votes
0answers
142 views

Analyze a football match: similar players with DBSCAN and similar trajectories with TRACLUS

I'm trying to analyze a dataset that originates from sensors located near players' shoes in a match (http://www.orgs.ttu.edu/debs2013/index.php?goto=cfchallengedetails). I decided to look at ...
7
votes
0answers
259 views

Luce choice axiom, question about conditional probability

I'm reading Luce (1959). Then I found this statement: When a person chooses among alternatives, very often their responses appear to be governed by probabilities that are conditioned on the ...
7
votes
0answers
118 views

Identification of parameters problem

I always struggle to get the true essence of identification in econometrics. I know that we state that a parameter (say $\hat{\theta}$) can be identified if by simply looking at its (joint) ...
7
votes
0answers
348 views

Can these data be aggregated into a proportion for a binomial glm?

We asked 60 people to list as many restaurant franchises in Atlanta as they could. The overall list included over 70 restaurants, but we eliminated those that were mentioned by fewer than 10% of the ...
7
votes
0answers
219 views

Can I use optimally scaled variables for a factor analysis to account for rotation? If I can then how?

I have discussed this issue several times in this site, but I am asking it again for a final justification from the experts of our community. I wanted to extract four factors (I should call dimensions ...
7
votes
0answers
944 views

Multiple imputation questions for multiple regression in SPSS

I am currently running a multiple regression model using imputed data and have a few questions. Background: Using SPSS 18. My data appears to be MAR. Listwise deletion of cases leaves me with only ...
7
votes
0answers
124 views

Bayesian, MDL or ML interpretation of cross-validation?

Is there any known Bayesian, ML or MDL interpretation of cross-validation? Can I interpret cross validation as performing the right update on a specifically crafted prior?
7
votes
0answers
235 views

How can we bound the probability that a random variable is maximal?

Suppose we have $N$ independent random variables $X_1$, $\ldots$, $X_n$ with finite means $\mu_1 \leq \ldots \leq \mu_N$ and variances $\sigma_1^2$, $\ldots$, $\sigma_N^2$. I am looking for ...
7
votes
0answers
253 views

Distribution of inverse Wishart to a power?

In a related question, I had asked about the norm induced by an inverse Wishart matrix. I am interested in generalizing that result somewhat. Let $A\sim\mathcal{W}_p\left(I,n\right)$, a Wishart matrix ...
7
votes
0answers
77 views

Estimation and functional space

In the first chapter of the book Algebraic Geometry and Statistical Learning Theory which talks about the convergence of estimations in different functional space, it mentions that the Bayesian ...
7
votes
0answers
1k views

How to analyze longitudinal count data: accounting for temporal autocorrelation in GLMM?

Hello statistical gurus and R programming wizards, I am interested in modeling animal captures as a function of environmental conditions and day of the year. As part of another study, I have counts ...
7
votes
0answers
460 views

What data and statistics skills are currently in high demand and where are they in high demand?

I have a job doing data analysis in finance. My current job is such that I don't have much exposure to things happening in the rest of my industry or other industries. I have a fair amount of ...
7
votes
0answers
956 views

How to use G Power 3 to calculate statistical power in mixed design ANOVA with unequal group sample sizes

In G power 3, ANOVA repeated measures within-between interaction: Only the total sample size is reported assuming equal sample size for the two groups. My questions are: How would it work if the ...
7
votes
0answers
388 views

Marginal distribution of the diagonal of an inverse Wishart distributed matrix

Suppose $X\sim InvWishart(\nu, \Sigma_0)$. I'm interested in the marginal distribution of the diagonal elements $diag(X) = (x_{11}, \dots, x_{pp})$. There are a few simple results on the distribution ...
7
votes
0answers
524 views

How can I assess GEE/logistic model fit when covariates have some missing data?

I have fit two generalized estimating equation (GEE) models to my data: 1) Model 1: Outcome is longitudinal Yes/No variable (A) (year 1,2,3,4,5) with longitudinal continuous predictor (B) for years ...
7
votes
0answers
200 views

Fisher information in a hierarchical model

Given the following hierarchical model, $$ X \sim {\mathcal N}(\mu,1), $$ and, $$ \mu \sim {\rm Laplace}(0, c) $$ where $\mathcal{N}(\cdot,\cdot)$ is a normal distribution. Is there a way to get an ...
6
votes
0answers
112 views

Overdispersion and modeling alternatives in Poisson random effect models with offsets

I have run into a number of practical questions when modeling count data from experimental research using a within-subject experiment. I briefly describe the experiment, data, and what I have done so ...
6
votes
0answers
69 views

Prediction interval based on CV

In the text books and youtube lectures I learned a lot about iterative models such as boosting, but I never saw anything about deriving a prediction interval. Cross validation is used for the ...
6
votes
0answers
66 views

What distribution results in adding two Pareto distributions

I'm wondering what distribution results in adding two (or more) type-one Pareto distributions of the form $x^{-\alpha}$. Experimentally, it looks like a two-mode power-law, asymptotic to the ...
6
votes
0answers
140 views

Applying a variance-stabilizing transform to a fitted function (rather than data)

Outline I'm working with data corrupted by a mixed Poisson-Gaussian noise model (for example with images gathered in astronomy or electron microscopy), and have been using the generalized Anscombe ...
6
votes
0answers
132 views

Confidence bands for QQ line

This question doesn't specifically pertain to R, but I chose to use R to illustrate it. Consider the code for producing ...
6
votes
0answers
54 views

Thesaurus for statistics and machine learning terms

Does there exist any reference thesaurus for statistics and machine learning terms? I know that Wikipedia articles often contain synonyms, but I would like to have a mere thesaurus that I could go ...
6
votes
0answers
113 views

Why doesn't correlation of residuals matter when testing for normality?

When $Y = AX + \varepsilon$ (i.e., $Y$ comes from linear regression model), $$\varepsilon \sim \mathcal{N}(0, \sigma^2 I) \hspace{1em} \Rightarrow \hspace{1em} \hat{e} = (I - H) Y \sim \mathcal{N}(0, ...
6
votes
0answers
356 views

Instrumental variables and mixed/multilevel models

I want to estimate a growth model to model the growth trajectories of individuals $j$ over multiple time points $t$ by applying a standard mixed/mutilevel model (also known as random coefficient ...
6
votes
0answers
101 views

Standard errors from flexsurvreg

I'm using the flexsurv package (in R) to fit an exponential distribution to the veteran ...
6
votes
0answers
177 views

Detecting outliers using 95% PI around a natural spline fit

Please see the picture below: I wanted to mark the points that are not consistent with their adjacent points as outlier. What I did was to fit a natural spline fit to 1000 observations (the purple ...
6
votes
0answers
231 views

When an analytical Jacobian is available, is it better to approximate the Hessian by $J^TJ$, or by finite differences of the Jacobian?

Let's say I'm computing some model parameters my minimizing the sum squared residuals, and I'm assuming my errors are Gaussian. My model produces analytical derivatives, so the optimizer does not need ...
6
votes
0answers
114 views

Using percentiles as predictors - good idea?

I am thinking about a problem which is to predict log(spend) of a customer using linear regression. I am considering what features to use as input and wondering if it would be OK to use the ...

15 30 50 per page