The response variable is level of free fatty acid in a sample of prepubescent boys. The explanatory variables are age (in months), weight (in lbs), and skin fold thickness.
Hettmansperger, T.P. and McKean J.W. (2011), Robust Nonparametric Statistical Methods, 2nd ed., New York: Chapman-Hall.
Examples
data(ffa)
summary(rfit(ffa~age+weight+skin,data=ffa)) #using the default (Wilcoxon scores)
summary(rfit(ffa~age+weight+skin,data=ffa,scores=bentscores1))
Results
R version 3.3.1 (2016-06-21) -- "Bug in Your Hair"
Copyright (C) 2016 The R Foundation for Statistical Computing
Platform: x86_64-pc-linux-gnu (64-bit)
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(Rfit)
Loading required package: quantreg
Loading required package: SparseM
Attaching package: 'SparseM'
The following object is masked from 'package:base':
backsolve
> png(filename="/home/ddbj/snapshot/RGM3/R_CC/result/Rfit/ffa.Rd_%03d_medium.png", width=480, height=480)
> ### Name: ffa
> ### Title: Free Fatty Acid Data
> ### Aliases: ffa
> ### Keywords: datasets
>
> ### ** Examples
>
> data(ffa)
> summary(rfit(ffa~age+weight+skin,data=ffa)) #using the default (Wilcoxon scores)
Call:
rfit.default(formula = ffa ~ age + weight + skin, data = ffa)
Coefficients:
Estimate Std. Error t.value p.value
(Intercept) 1.4899763 0.2694023 5.5307 2.714e-06 ***
age -0.0011321 0.0026360 -0.4295 0.6700657
weight -0.0153345 0.0038481 -3.9849 0.0003048 ***
skin 0.2743114 0.1342777 2.0429 0.0482380 *
---
Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
Multiple R-squared (Robust): 0.3756868
Reduction in Dispersion Test: 7.42171 p-value: 0.00052
> summary(rfit(ffa~age+weight+skin,data=ffa,scores=bentscores1))
Call:
rfit.default(formula = ffa ~ age + weight + skin, data = ffa,
scores = bentscores1)
Coefficients:
Estimate Std. Error t.value p.value
(Intercept) 1.35964045 0.18873747 7.2039 1.519e-08 ***
age -0.00048203 0.00178357 -0.2703 0.7884625
weight -0.01539475 0.00260369 -5.9127 8.244e-07 ***
skin 0.35616269 0.09085428 3.9202 0.0003686 ***
---
Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
Multiple R-squared (Robust): 0.4758888
Reduction in Dispersion Test: 11.19857 p-value: 2e-05
>
>
>
>
>
> dev.off()
null device
1
>