All Questions
3,928 questions
6
votes
1
answer
409
views
Wildly different answers replicating a GEE model from SPSS
Describe the bug
I'm attempting to replicate a GEE model in statsmodels from a published paper that used SPSS (https://pubmed.ncbi.nlm.nih.gov/33279717/). I am getting very different answers for what ...
6
votes
2
answers
21k
views
Binary or Multinomial Logistic Regression?
I wish to analyze the following :
Predictor Variable (IV): Satisfaction of sexual needs as important (4 items scale and respond based on 4 point likert scale. Sum up to get the item score.)
Response ...
6
votes
2
answers
454
views
Multilevel regression model for nested data using "multilevel" and "lme4" R packages?
Hypothesis 1
Independent variable = Thermal inertia (TI)
Dependent variable = Depth diameter ratio (d/D), Radii variation (RV) and Rim irregularity (RI).
Hypothesis 2
Independent variable = Thermal ...
6
votes
2
answers
1k
views
What is the Joint Density Function of a Three-Level Mixed-Effects Model?
This is a follow-up question to a question I posted earlier. Obviously, maximum-likelihood estimation of mixed-effects models requires the joint density function.
Let us assume a two-level mixed ...
6
votes
1
answer
33k
views
Hierarchical or TwoStep cluster analysis for binary data?
(This question is an edited version of a question I previously posted which one user recommended would benefit from more focus).
I have 2000 questionnaires from respondents which ask 33 different ...
6
votes
2
answers
638
views
Correlation between vegetation and erosion
I am wanting to find the correlation between two variables: a categorical vegetation variable (with categories 1-4) and a continuous variable of erosion (with values anywhere between 0.2 and 9). I ...
6
votes
1
answer
12k
views
Inter-rater reliability using Intra-class correlation with ratings for multiple objects on multiple properties
Initial Question:
I am trying to calculate inter-rater reliability. Previous researchers in this area have used intraclass correlation. SPSS has options for two-way random, mixed random and one-way ...
6
votes
2
answers
1k
views
Justifying the Need for Mixed Effects Models (aka. LME, MLM, etc.)
Firstly, I am not an expert in using multi-level modelling (MLM), and I have read this and this questions, however, my question is slightly different in the sense that method 2 below is not mentioned.
...
6
votes
1
answer
3k
views
Distinct results between "emmeans" and "multcomp" - package in multi level model
I have calculated a multi-level model with a biomarker as dependent variable (which was measured three time), a 5-level factor variable called „module“ as predictor (which is an intervention including ...
6
votes
2
answers
1k
views
Conditional intraclass correlation (ICC) from a linear mixed model as evidence for test-retest reliability?
In my experiment with two conditions (between-subjects design), participants completed a single-item scale three times: (1) before the experimental manipulation, (2) after the experimental ...
6
votes
1
answer
748
views
Using Bayesian model diagrams to present both model description and results (posteriors)?
The model diagrams in "Doing Bayesian Data Analysis", John Kruschke creates diagrams like this:
To represent The following BUGS/JAGS code:
He discusses this representation in his related blog post, ...
6
votes
2
answers
812
views
Possible identifiability issue in hierarchical model
I'm trying to fit some data using a hierarchical normal model
$y_i \sim N(\theta_i,\sigma^2)$
$\theta_i \sim N(\mu, \sigma_\theta^2)$
$(\mu,\sigma^2,\sigma_\theta^2) \sim diffuse$
I fit this ...
6
votes
2
answers
899
views
Deriving mathematical model of pLSA
After knowing how LSA works, I went on continue reading on pLSA but couldn't really make sense of the mathematical formula. This is what I get from wikipedia (other academic papers/tutorial show ...
6
votes
3
answers
2k
views
Strange outcomes in binary logistic regression in SPSS
I did a binary logistic regression with SPSS 23 and I found some strange outcomes. This is for NOACprev until No_Prev_treatment, the last 6 variables. First of all they have very high outcomes for B, ...
6
votes
4
answers
35k
views
Intraclass correlation coefficient interpretation
I'm having a look at the intraclass correlation coefficient in SPSS.
Data: 17 participants rated two lists of 9 & 7 items from 0 to 5 (0 being unimportant and 5 being very important).
All ...
6
votes
2
answers
9k
views
How to do meta-regression in SPSS?
I am trying to manage a meta-regression in SPSS17 using the effect size as the dependent variable. I want to explore if my independent variables affects the effect size. Some small practical questions:...
6
votes
1
answer
28k
views
Should I treat these ordinal IVs as covariates or factors, in a regression?
I have a survey problem where the dependent variable (ordinal) is in Likert-type scale (i.e. 1 to 5 from most satisfied to most dissatisfied) and two sets of independent variables. One set has 7 IVs (...
6
votes
1
answer
2k
views
Inconsistency in mixed-effects model estimation results (Stata and SPSS)
For my thesis there's a big chance that I will need some sort of mixed-effects specification. I have some (non-syntax) experience with SPSS but feel that it won't suffice for my analysis. I have very ...
6
votes
2
answers
761
views
Longitudinal relationship between chocolate consumption and happiness: repeated measures ANOVA?
I am using SPSS and having some trouble with a research question which is analogous to the hypothetical question:
Is there a longitudinal Relationship between Happiness and Chocolate Consumption?
...
6
votes
1
answer
127
views
Is there a way to forecast by subgroup without forecasting each subgroup separately?
I am trying to find an appropriate model to forecast the number of applications received at the end of a recruitment cycle based on previous recruitment cycles and the number of applications received ...
6
votes
2
answers
1k
views
Why can't we use top-down methods in forecasting grouped time series?
As I asked in here I was trying to forecast grouped time series with two grouping variables and I find some limitation of hierarchical forecasting methods. In particular, using hts package from R, we ...
6
votes
1
answer
3k
views
Difference between a random slope/intercept model and an ANCOVA with an interaction?
I attempted to run an ANCOVA with one binary predictor, one continuous outcome, and one continuous covariate. I found that there was heterogeneity of regression slopes and thus I concluded that an ...
6
votes
1
answer
5k
views
How to assess multilevel model assumptions using residual plots
I am not clear as to how to assess if a multilevel model fit using lmer satisfies the assumptions of normality and homoscedasticity?
I have used the following <...
6
votes
2
answers
3k
views
Correct specification of longitudinal model in lme4
I am trying to fit a multilevel longitudinal model and i have a question regarding how to specify it.
The data consist of about 8k observations collected from about 3k individuals at four time points. ...
6
votes
1
answer
444
views
Expected Fisher information isn't positive definite for truncated normal with heteroskedasticity
This question is about having a non-positive-definite expected Fisher information in a normal model in which observations have different dispersions.
Consider this simple normal model:
$$Y_i \sim N(\...
6
votes
1
answer
307
views
How to interpret my profile likelihood plots (three-level meta-analysis)?
I conducted a three-level meta-analysis to test effects of non-driving related tasks on take-over quality in highly automated driving. When it now comes to interpretation of the results, I struggle to ...
6
votes
1
answer
2k
views
Model not singular but doesn't converge what could be the reason (lme4 in R)
I'm following up on this great answer regarding running Principal Component Analysis (PCA) to uncover the reason behind lack of convergence and/or singularity for Mixed-Effects models.
My model below ...
6
votes
1
answer
3k
views
How can I make the correct time-series analysis for my data?
I conducted a research with 90 participants who were exposed to a film stimuli for 60 minutes and their physiological responses were measured continuously. I would like to test the causal relationship ...
6
votes
1
answer
7k
views
ICC in a multi level model with two random effects
My understanding is that intraclass correlation gives you an idea of how much variance your level two factor can explain in overall variance of the dependent variable. It is supposed to give an ...
6
votes
1
answer
4k
views
Acceptable values for the intraclass correlation coefficient (empty model)
I'm using xtmixed in Stata to test a Hierarchical Linear Model.
My problem is that variance at level 2 is about 4% of the total variance. So most of the variance is ...
6
votes
4
answers
6k
views
What regression analysis should I perform on my data and why?
I am a law student researching which factors influence the CSR (corporate social responsibility, GSE_RAW) behavior of companies. As my studies didn't offer any ...
6
votes
1
answer
178
views
Role of regression model fit in causal analysis
When analysing causal questions, we use DAGs that give us covariates needed for modelling. But another time we assess model fit to get the best prediction. These two approaches have different purposes ...
6
votes
1
answer
77
views
Doubt regarding mixed modeling format
Say, I have a dataset that looks at how many times my 5 babies chases a cat around the house . I'm trying to estimate 'y' which is the number of times the cat runs one complete round around the house ...
6
votes
2
answers
854
views
How do I interpret the coefficients of a mixed effects multilevel logistic regression differently from regular logistic model?
I am trying to wrap my head around mixed effects multilevel logistic regression. Have a look at my variables:
y: Popularity (0 = Not popular, 1 = Popular)
x1: Extraversion (Continuous)
x2: Teacher ...
6
votes
1
answer
3k
views
Weighted cases in a cluster analysis for cases in SPSS
I am conducting a cluster analysis (of cases) for a database which has weight attributed to the individual cases to ensure that it mirrors the general population in terms of sociodemographic ...
6
votes
1
answer
4k
views
Two-level hierarchical model using time-series cross sectional data?
A question from someone who is fairly new to hierarchical modeling, and I'm looking for the best approach within R, preferably with the package lme4, MCMCpack, or rjags using a BUGS document. I'm ...
6
votes
3
answers
11k
views
Creating ROC curve for multi-level logistic regression model in R
I used the functions from this link for creating ROC curve for logistic regression model. Since the object produced by glmer in ...
6
votes
3
answers
2k
views
How to control for market return in an (SPSS) OLS?
Please consider the following panel dataset:
...
6
votes
1
answer
173
views
Getting an average measurement based on two raters for cases where data is missing for one rater
Context:
I'm investigating behaviour in a clinical study involving children. I had both parents and teachers completing questionnaires to inform an understanding of the same underlying constructs, ...
6
votes
2
answers
96
views
Appropriate regression model when variables were measured on different subjects?
Imagine the following study design: We have propagated $N$ plants. Each plant produces many seeds. From each plant, $M$ of the seeds are used to measure community composition of bacteria living within ...
6
votes
1
answer
400
views
Proper syntax for coding group level variables in mixed effect model using GLMER
I am running a multilevel log odds regression on a dataset that is at the individual level. Each individual belongs to a district, and there are both individual level variables (which vary by ...
6
votes
1
answer
375
views
Increasing multicollinearity in multilevel/hierarchical modeling?
I have a linear model with response variable $\textbf{y}$ and explanatory variable matrix $\textbf{X}$ for which coefficients $\textbf{b}$ are physically meaningful and worth estimating:
\begin{...
6
votes
1
answer
2k
views
Is it possible to use LASSO regression with multi-levlel data?
I have real-time monitoring data where participants report on a variety of variables four times per day for a month. Is it possible to use LASSO regression (e.g,. glmnet r package) with this data? I'm ...
6
votes
3
answers
2k
views
Whether to use a hierarchical linear model
I've been reading through Gelmans book: Data Analysis Using Regression and Multilevel/Hierarchical Models trying to learn more about how to implement hierarchal models. I have a dataset that I think ...
6
votes
1
answer
1k
views
Wald vs. LR $\chi^2$ tests in SPSS generalized linear models
I'm using the generalized linear models function in SPSS with a normal distribution and identity link function. If I choose the likelihood ratio $\chi^2$ statistic, I get the same results as the ...
6
votes
2
answers
5k
views
False discovery rate & q-values: how are q-values to be interpreted when rank of p-values is altered?
I am trying to wrap my head around False Discovery Rate, and its associated q-value; I am new to this technique, but it seems quite promising for my needs.
One sticking point I keep coming across and ...
6
votes
3
answers
477
views
Estimating correlated parameters with multi-level model
I would like to estimate a multi level model in Stata or R (using lmer) where the first level coefficients are the same for all observations, but the coefficients within observation are correlated.
...
6
votes
1
answer
2k
views
Plotting (multilevel) multiple regression [closed]
Lets say I have some data like this:
...
6
votes
1
answer
2k
views
How many levels in multilevel modeling is too many?
This is pretty general, but what are the pros and cons of including additional levels in multilevel model (linear mixed model)?
I have a data containing information on multilevel administrative ...
6
votes
1
answer
10k
views
Mixed model in SPSS with random effect and repeated measures
I am working on analyzing a dataset that involves repeated measures data. The data was previously analyzed by a colleague using custom code written in C++, but I have expanded the dataset and am ...