Search Results
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 |
Generalized additive models for location, scale and shape (GAMLSS).
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 …
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? …
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 …
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) + …
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) …
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) …
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 …
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? …
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 …