254 questions
-1 votes
0 answers
17 views
Is MANCOVA the right model to use in the experiment? [closed]
I'm conducting an experiment and need to choose a reasonable statistical model. Study design: 3 groups (each group will be presented with varying instructions to the same task). This is the IV. 2 DVs ...
Advice
0 votes
1 replies
30 views
Check assumptions of mixed effects cox model
I am running a mixed effects cox model using coxme in R, and I cannot find any information on how to check model assumptions for this type of model. I have seen the older chats, which have recommended ...
0 votes
1 answer
81 views
Testing equality of multivariate coefficients with linearHypothesis()
I would like to conduct a F-test (Wilks' Lambda test, ideally) to test the equality of the slope coefficients for a single independent variable on two dependent variables. Using the mtcars dataset, ...
0 votes
1 answer
92 views
Stepdown analysis of multivariate models in R
Roy - Bargmann stepdown tests are commonly recommended for multivariate linear models (MANOVA, MMRA) fit with lm() as: lm(cbind(y1, y2, y3, ...) ~ x1 + x2 + x3 + ...) But I can't find any ...
1 vote
1 answer
57 views
is it posible to fit a multivariate linear model using lmtree?
I want to fit a lmtree where the response is multivariable. You know that you can fit a Multivariate Multiple Regression using lm function as follows: lm(cbind(Y1, Y2) ~ X1+X2+X3, data = my_df) so I ...
1 vote
1 answer
68 views
Need help using "reframe" to run multivariate regression across multiple data sets (changing model from one year to multiple years)
I have a simple multivariate regression model testing the impact of two variables on another over one year (1971). I have a large data set that includes these data over a number of years and am being ...
1 vote
1 answer
160 views
Optimizing variable applied to discrete data in order to minimize error function
I'm trying to optimize a function based on discrete data sets I got in the lab. Spectrum_1 and Spectrum_2 are experimental data sets with length N. Those two arrays contain values taken for each ...
0 votes
0 answers
42 views
Creating a vector of in which Mardia's b1p can be saved
i want to create a vector that have 'q0' rows and 1 column with N repetitions. My minimum working example is: library(Matrix) library(psych) library(MASS) set.seed(10) N0 <- 2 n0 <- 5 p0 <- ...
0 votes
1 answer
610 views
Multivariate Linear Mixed Model
I want to make a multivariate linear mixed model with the variables slope_reactivity, peak_reactivity, and slope_recovery as dependent variables and normalized_CTseverity, gender, age_years and ...
2 votes
1 answer
49 views
Want to create a vector in which multivariate skewness of the four multivariate normal populations can be stored
Minimum Working example: library(Matrix) library(psych) library(MASS) set.seed(10) N0 <- 1 n0 <- 5 p0 <- 2 q0 <- 4 n <- n0*q0 m2 <- matrix(c(0, 0), p0, 1) s2 <- matrix(c(1, 0, 0,...
0 votes
1 answer
76 views
How to save Mardia's coefficient of skewness in a vector when the data is divided into 4 parts
Minimum Working Example: ############################################################################## set.seed(10) ############################################################################## ...
0 votes
1 answer
182 views
Coefplot.gllvm shows different plots per categorical factor level rather than each factor
I'm not well versed in R but am using the GLLVM package to assess the influence of season and water mass on plankton distribution - when I try to make the coefficient plots using coeffplot.gllvm() I ...
0 votes
1 answer
74 views
Stepwise Multivariate Linear Regression in RStudio Stalling
My data has 6 independent variables and 5 dependent variables, and I am running stepwise multivariate linear regressions on them. My code removes any independent variables with p-values > 0.05 in ...
1 vote
1 answer
232 views
Multivariate regression tree with "mvpart" (in R) and plots for each leaf of the tree visualization
I want to perform an MVRT to understand how some abiotic factors affect the cover percentage of different benthic species. So, I used the 'mvpart' package in R version 3.2.5, as it is not possible to ...
1 vote
1 answer
72 views
mvabund fourth corner model using categorical environmental levels
Within mvabund::traitglm() is it appropriate to use a categorical R matrix rather than continuous environmental variables? I am interested in the effect of the interaction between MHW period (before, ...