R Under development (unstable) (2025-04-26 r88181 ucrt) -- "Unsuffered Consequences" Copyright (C) 2025 The R Foundation for Statistical Computing Platform: x86_64-w64-mingw32/x64 R is free software and comes with ABSOLUTELY NO WARRANTY. You are welcome to redistribute it under certain conditions. Type 'license()' or 'licence()' for distribution details. R is a collaborative project with many contributors. Type 'contributors()' for more information and 'citation()' on how to cite R or R packages in publications. Type 'demo()' for some demos, 'help()' for on-line help, or 'help.start()' for an HTML browser interface to help. Type 'q()' to quit R. > > > library(msme) Loading required package: MASS Loading required package: lattice > > data(ufc) > ufc <- na.omit(ufc) > > ufc.g.reg <- ml_g(height.m ~ dbh.cm, data = ufc) > > ufc.g.lm <- lm(height.m ~ dbh.cm, data = ufc) > > ufc.g.reg Call: ml_g(formula = height.m ~ dbh.cm, data = ufc) No coefficients > ufc.g.lm Call: lm(formula = height.m ~ dbh.cm, data = ufc) Coefficients: (Intercept) dbh.cm 12.6757 0.3126 > > summary(ufc.g.reg) Call: ml_g(formula = height.m ~ dbh.cm, data = ufc) Residuals: Min 1Q Median 3Q Max -33.525 -2.863 0.132 2.851 13.320 Coefficients: Estimate SE Z LCL UCL (Intercept) 12.67708 0.56271 22.52853 11.57416 13.780 dbh.cm 0.31256 0.01384 22.57707 0.28543 0.340 Log Sigma 1.59521 0.03577 44.60167 1.52511 1.665 Residual standard error: 4.929 on 391 degrees of freedom > summary(ufc.g.lm) Call: lm(formula = height.m ~ dbh.cm, data = ufc) Residuals: Min 1Q Median 3Q Max -33.526 -2.862 0.132 2.851 13.321 Coefficients: Estimate Std. Error t value Pr(>|t|) (Intercept) 12.67570 0.56406 22.47 <2e-16 *** dbh.cm 0.31259 0.01388 22.52 <2e-16 *** --- Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1 Residual standard error: 4.941 on 389 degrees of freedom Multiple R-squared: 0.566, Adjusted R-squared: 0.5649 F-statistic: 507.4 on 1 and 389 DF, p-value: < 2.2e-16 > > coef(ufc.g.reg) NULL > > logLik(ufc.g.reg) 'log Lik.' NaN (df=) > logLik(ufc.g.lm) 'log Lik.' -1178.469 (df=3) > > AIC(ufc.g.reg) numeric(0) > AIC(ufc.g.lm) [1] 2362.938 > > > > > > proc.time() user system elapsed 0.28 0.12 0.39