Skip to main content
Search type Search syntax
Tags [tag]
Exact "words here"
Author user:1234
user:me (yours)
Score score:3 (3+)
score:0 (none)
Answers answers:3 (3+)
answers:0 (none)
isaccepted:yes
hasaccepted:no
inquestion:1234
Views views:250
Code code:"if (foo != bar)"
Sections title:apples
body:"apples oranges"
URL url:"*.example.com"
Saves in:saves
Status closed:yes
duplicate:no
migrated:no
wiki:no
Types is:question
is:answer
Exclude -[tag]
-apples
For more details on advanced search visit our help page
Results tagged with
Search options not deleted user 364552

Generalized additive models for location, scale and shape (GAMLSS).

3 votes
1 answer
224 views

How to compare centiles from different models?

I am comparing the centiles from different GAMLSS models. Which model has a better performance. Why? For model 1, the number of cases below 0.4 centile is 0.5. Is this possible? …
stats's user avatar
  • 89
2 votes
1 answer
486 views

GAMLSS vs VGAM for percentile curves (Growth chart)

. -- GAMLSS (Generalized Additive Model for Location, Scale and Shape) -- VGAM (Vector generalized additive model) Which one of these two techniques and their corresponding R packages can estimate more …
stats's user avatar
  • 89
1 vote
0 answers
92 views

A question about applying GAMLSS models

So I used the following function in R m4 <- gamlss(y ~ log(Ph1_Groesse) + pb(log(Ph1_Alter_2)), sigma.fo =~pb(log(Ph1_Alter_2)), nu.fix=T, nu.start=1, family = BCCGo(mu.link = "log"), … A4<- gamlss(y ~ log(Ph1_Groesse) + pb(log(Ph1_Alter_2)), sigma.fo =~log(Ph1_Alter_2), nu.fix=T, nu.start=1, family = BCCGo(mu.link = "log"), data=DAT1.F) …
stats's user avatar
  • 89
1 vote
1 answer
128 views

How to find the appropriate model to apply (GAMLSS)/Approving statistical thinking

BCCGo(mu.link = "log"), data=df) m3 <- gamlss(TLC ~ log(Height) + pb(Age), sigma.fo =~pb(Age), nu.fo =~pb(Age), family = BCCGo(mu.link = "log"), data=df) m1.1 <- gamlss(TLC ~ Height + pb(Age), sigma.fo … =~pb(Age), nu.fo =~1, family = BCCGo(mu.link = "log"), data=df) m2.1 <- gamlss(TLC ~ Height + pb(log(Age)), sigma.fo =~pb(log(Age)), nu.fo =~1, family = BCCGo(mu.link = "log"), data=df) m3.1 <- gamlss
stats's user avatar
  • 89
1 vote
1 answer
118 views

Choosing whether to eliminate or keep a predictor in a GAMLSS model

I need to calculate the centile curve for y using a GAMLSS model with age and height as predictors. The plots below depict the relationship between log(y) and each of the independent variables. … m8 <- gamlss(y ~ log(Ph1_Groesse) + pb(log(Ph1_Alter_2)), sigma.fo =~pb(log(Ph1_Alter_2)), nu.fo =~1, family = BCCGo(mu.link = "log"), data=DAT1.F) …
stats's user avatar
  • 89
0 votes
1 answer
94 views

Proportion/Ratio response variable

In order to calcualte the curves with LMS method from gamlss package we use the following code (forexample) library(gamlss) A<-lungFunction m1 <- gamlss(slf ~ log(height) + pb(log(age)), sigma.fo =~pb …
stats's user avatar
  • 89
0 votes
1 answer
184 views

Find the more appropriate fit by using term.plot in GAMLSS

I'm comparing several GAMLSS models and was wondering how I might use the term.plot to identify the best model/fit. What type of information can we get from this graph? …
stats's user avatar
  • 89
0 votes
2 answers
103 views

qBCCGo function in GAMLSS package (quantile function)

I am working with GAMLSS technique and I use LMS (with three parameters of mean,variation and skewness) method. … LLN= exp(log(mu)+log(1-1.645*nu*sigma)/nu) In GAMLSS package, there is a function qBCCGo(p, mu = 1, sigma = 0.1, nu = 1, lower.tail = TRUE, log.p = FALSE) to calculate the lower limit of normal for Box-Cox …
stats's user avatar
  • 89
0 votes
0 answers
58 views

Need some help to identify the appropriate distribution

Thanks so much for your time and help in advance. m4 <- gamlss(y ~ log(Ph1_Groesse) + pb(log(Ph1_Alter_2)), sigma.fo =~pb(log(Ph1_Alter_2)), nu.fo=~log(Ph1_Alter_2), family = BCCGo(mu.link = "log"), data … =DAT1.F) m17 <- gamlss(y ~ log(Ph1_Groesse) + pb(Ph1_Alter_2), sigma.fo =~pb(Ph1_Alter_2), nu.fo =~pb(Ph1_Alter_2), family = BCCGo(mu.link = "log"), data=DAT1.F) m24 <- gamlss(y ~ log(Ph1_Groesse) + …
stats's user avatar
  • 89