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.