# Tag Info

1 vote

### Interpretation of binomial regression in R

[...] why is the data for smkcigst NA, even though the dataset has no NAs for it and instead has 1 for all observations? The answer lies in your question: smkcigst ...

### How to test linear hypotheses with inequality in r using the function glht for a model in log-form?

I've happened to find the answer by searching for a solution on the internet. For all with the same problem: Just put \ and ` together before and after your variable name. ...

### T-distribution parameters with QRM package

par.ests referes to the parameter estimates. par.ses refers to the parameter standard errors. It is maybe not so obvious, since ...
1 vote

### How to test linear hypotheses with inequality in r using the function glht for a model in log-form?

Classical hypothesis test are a bit funny for this kind of thing. Testing a null hypothesis that consists of an inequality is not all that different from testing the corresponding equation --- the ...

### Are these effects missing from my glm output because of a possible dummy trap?

With only 5 values of Lineage you are getting yourself into trouble trying to treat it as a random effect. That typically leads to singularity problems with the ...
1 vote

### How to perform lasso on a wide matrix?

The error message indeed suggests that you don't have enough RAM. However, this answer suggests that glmnet should be able to handle this amount of data, and suggests h2o for situations where glmnet ...

### Maximum Likelihood with Categorical Variables - Does this Change Anything?

On the one hand, there are features, which are attributes of your data, stuff that you measure, and on the other hand, there are the parameters to your model. Both can be continuous or discrete. If ...
Accepted

### Compute two proportion tests (with three variables?)

Comment. One more thing: In R, the procedure prop.test is essentially the same as a chi-squared test of homogeneity, using your $3 \times 2$ table. Except for ...

### Maximum Likelihood with Categorical Variables - Does this Change Anything?

The algorithm used is exactly the same: although the features are not continuous in your example, the coefficients in a regression model still are. We are optimizing the coefficients in the model, the ...

### Discrepancy in degrees of freedom from R svyglm vs glm

You don't say anything about your survey design. The residual df for a svyglm object is the design df plus one, minus the number of parameters estimated. From the ...
1 vote
Accepted

### Doubt regarding inverse CDF/quantile function or qpois in R

First, here is a 4-place table of PDF and CDF values for $\mathsf{Pois}(\lambda = 9.29),$ x = 0:5 pdf = round(dpois(x, 9.29), 5) cdf = round(ppois(x, 9.29), 5) (...
1 vote
Accepted

1 vote

### I have a group that contains zero data, and I want to know whether it is considered normally or not normally distributed.?

Comment in Answer format--to allow showing output from R: Clearly, there is a difference between Groups 1 and 2. You could take the position (1) that the 0's in Group 2 are real data, but simply too ...

### Ideal Settings for Longitudinal Models?

There is a simple counterexample to your suggestion. Imagine a case where some units are affected by large measurement error, while others have only small measurement error. It is more statistically ...

### R: Anova and Linear Regression

ANOVA is a specific regression model with categorical regressors. It only models means of different groups i.e. intercepts. There is no slope parameter. Regression is a general methodology to estimate ...

### QQ plot result doesn't correspond to normality test

Just on first look, this distribution looks very short tailed, as you can see it looks kind of like this simulation with a uniform distribution

### QQ plot result doesn't correspond to normality test

Your Q-Q plot does not show that the data is normally distributed. In fact, it shows that the distribution diverges from Gaussian (values lower than -1 and higher than 1.5 diverge from the diagonal ...
1 vote

### Why do fixed effects in a logistic regression model differ depending on the presence of a random slope?

Fixed effects change You do get that the fixed effects change when you add a random effect. With your example both the intercept and slope change. Below is an example of the situation. black line: ...
1 vote

### How to tell if data has homgenous variance?

For an ANOVA model like this, you have 1 mean value estimate for each of 3 groups. The residuals are the differences of each individual observation from its corresponding group mean. Roughly, the ...
Accepted

### How to determine how much variation was due to the differences between a group?

Analysis of variance (ANOVA) provides the answer to your specific question. Here, the proper outcome to examine is Time in the ...
1 vote
Accepted

### How to fit parametric distribution to baseline hazard in cox PH model

Several thoughts: First, even if you could generate a smooth baseline survival curve beyond 130 months, it would be unreliable. That's implicit in the large confidence bands around the survival curve ...

### Time series model in production - Re-train on the fly as as batch process?

What you can do is basically train the model for example once a day (depending on the interval), save it and append its forecasts to a copy of the array you're storing your data in so all of it could ...
Accepted

### Missing Prediction Intervals by R package autoTS

Not all methods of producing point forecasts are probabilistic in nature, in which case there may really be no prediction interval to report. This package is focused on point forecasts and picks ...

### Interpreting coefficients of ordinal independent variables in logistic regression in R

I understand that the influence is likely to be insignificant, but which P value do I use? What can I say about the other coefficients? The model turns the single categorical/ordinal variable ...

### Given a specific value for a variable, how do you find the predicted value of a fixed effects multivariate regression?

You have to use I(N^2) instead of N^2, because otherwise lm() interprets ...

### XGBoost interpretation of the plot in R

XGBoost creates an ensemble of binary trees, which are built "on top of each other". I.e. you first try to find a binary tree that best approximates your data. Then you build further trees ...
Accepted

### Non-positive definite matrix problem for desired correlation structure

You have stumbled into the fact that you cannot simply make a correlation matrix by assembling individually valid pairwise correlations. There are many questions on the site related to this, have a ...

### Cox PH controlling for multiple events

It would be simpler to think of this as a multi-state survival model, as explained in the vignette, and study both Zones together. The id values keep track directly ...

### Fit an ARMAX model in R

As of 2022, Hyndman is using fable::ARIMA, but in his otherwise excellent time series regression guide (https://otexts.com/fpp3/) is still only showing how to do a linear regression with arima errors. ...

### Interpreting interaction effects for categorical reference group in regression

@EdM makes valid points - read those first. Just for reference, you could get the effects of interest along with their confidence intervals using the emmeans package. For example: ...

### Interpreting interaction effects for categorical reference group in regression

First, with the default R treatment coding of your categorical predictors, the individual coefficients for things like Story Vision are their associations with ...

### Anova shows a main effect but not the t-test

I do not have your data but you are testing two different things. In the ANOVA model you are testing for the effect of room, time and their interaction in anxiety, in an analysis that uses the ...
1 vote

### Valid to compare variable importance ranks across RF with different responses?

I suggest thinking about this problem like an experiment. What is being manipulated/changed and what is being controlled? Thinking about the set of models like this, you might be able to see where ...
1 vote

### SE Interpretation

The standard error it's a metric that show the deviation that your estimated coefficient could have from the estimated value (mean), that is, the error. It's a number in the same units as the ...
1 vote

### pointwise envelopes not including Theoretical line Foxall J

Please provide a minimal working example. Pooling of objects (such as simulation envelopes) is only justified when the individual objects (envelopes) were created under identical conditions. In your ...

Top 50 recent answers are included