Skip to main content
Stack Overflow
  1. About
  3. For Teams
Filter by
Sorted by
Tagged with
0 votes
0 answers
61 views

I am modelling an health outcome, percentage days alive and out of hospital (%DAOH), with gamlss and family BEINF (zero-one-inflated beta). This outcome has point mass at 0 and 1, and clusters near ...
2 votes
0 answers
120 views

I'm trying to fit a log-logistic distribution (truncated) with GAMLSS and can't find how to do it. In the documentation, it says that for Generalized Beta 2 (GB2) : Setting ν = 1 and τ = 1 in (15.20) ...
0 votes
1 answer
95 views

I am a seasoned R programmer trying to help my friend to visualize percentile growth curves using gamlss. I have been successfully produce the curve for each group. Now the problem is she asked me to ...
0 votes
0 answers
73 views

I am using the R package GAMLSS to fit a GAM to illustrate the seasonality (by month) and trend (week_t number in a year) in a time series data (mydat) Here is the code I ran: fit_lss <- gamlss(...
Trang Nguyen's user avatar
0 votes
2 answers
537 views

I have a GAMLSS trained model which I have saved using saveRDS() in .rda format. For eg, I trained model as: gamlss_model<- gamlss(res~pb(x)+pb(y), family=BCTo, data = test) When I load the above ...
1 vote
0 answers
625 views

I’ll preface this with I am NOT savvy with statistics, so I appreciate any help, especially in simple terms. I hope to understand and not just get the "correct" answer. I also admit I am most ...
Clara's user avatar
  • 43
1 vote
1 answer
237 views

I have fitted a model with m0 <- lms(weight,age, data=df) After, I use the model to predict percentiles for same ages: newages <- c(1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15) centiles....
Juan's user avatar
  • 452
0 votes
1 answer
259 views

I'm trying to write functions that wrap basic gamlss modeling (because I'm preparing a higher-level, application-specific package). But I've run into issues wrapping predict. The example from the ...
0 votes
1 answer
49 views

I am trying to use qbcpeo to produce lower limit than normal for a gamless bcpe object, however, some of the sigma values (sigma link is log) yields negative value so when I attempt to use qbcpeo I ...
0 votes
0 answers
300 views

My shiny app works fine in the console, but when I clean the R environment and runApp(), I get Error: C stack usage 15924528 is too close to the limit. One related post suggests to type 'ulimit -s ...
Meer's user avatar
  • 1
0 votes
0 answers
55 views

I'm looking to run the fitDist function on columns in a dataframe using sapply and save the output to a dataframe, using something like the following: data.frame(sapply(df, fitDist)) Run singularly, ...
0 votes
1 answer
166 views

I am trying to fit a mixed effects multinomial model for my 4 level categorical outcome variable (y). I am using the gamlss function from gamlss package. The call looks like: library(mixcat) data(&...
1 vote
0 answers
74 views

Although a long time user of R, I have little experience to write complex programs in R. During one of my first tries, I experience a unexpected behavior of my code. Namely, if I try to run a function ...
0 votes
2 answers
37 views

my dataset can be found here: https://raw.githubusercontent.com/yuliaUU/test/main/test.csv library(gamlss) library(tidyverse) data_final<- read_csv("https://raw.githubusercontent.com/yuliaUU/...
yuliaUU's user avatar
  • 1,743
0 votes
1 answer
52 views

I'm trying to fit centile curves for a censored gamma distribution. The censoring goes well I think, but when trying to make the percentile curves, I get the following result: I tried the following ...

15 30 50 per page
1
2 3 4 5 6

AltStyle によって変換されたページ (->オリジナル) /