R petrol N. L. Prater 的炼油厂数据

说明

具有四个协变量的石油精炼过程的产量。原油似乎仅来自 10 个不同的样品。

Prater (1956) 最初使用这些数据来建立原油炼制汽油过程中产量的估计方程。

用法

petrol

格式

变量如下

No

原油样品识别标签。 (因子。)

SG

比重，API 度。 (样本内恒定。)

VP

蒸气压，单位为磅每平方英寸。 (样本内恒定。)

V10

原油的波动性； ASTM 10% 点。 (样本内恒定。)

EP

汽油的所需挥发性。 (终点。样本内有所不同。)

Y

产量占原油的百分比。

例子

library(nlme)
Petrol <- petrol
Petrol[, 2:5] <- scale(as.matrix(Petrol[, 2:5]), scale = FALSE)
pet3.lme <- lme(Y ~ SG + VP + V10 + EP,
                random = ~ 1 | No, data = Petrol)
pet3.lme <- update(pet3.lme, method = "ML")
pet4.lme <- update(pet3.lme, fixed. = Y ~ V10 + EP)
anova(pet4.lme, pet3.lme)

来源

N. H. Prater (1956) 估计原油的汽油产量。石油精炼厂 35, 236-238。

该数据集也在 D. J. Hand、F. Daly、K. McConway、D. Lunn 和 E. Ostrowski (eds) (1994) A Handbook of Small Data Sets 中给出。查普曼和霍尔。

参考

Venables, W. N. and Ripley, B. D. (2002) Modern Applied Statistics with S. Fourth edition. Springer.