在nlme中指定幂律的非线性混合效应模型

时间:2018-06-06 14:09:00

标签: r mixed-models non-linear-regression nlme

我想使用nlme软件包以混合效果的方式拟合幂数法。我能够使用brms来拟合贝叶斯模型,见下文,但我想要适应频率论者'使用nlme或其他一些能够适应非线性混合效果模型的包的相同模型的版本。 ab应允许File变化。

这是有效的brms语法:

prior1 <- prior(normal(0, 5), nlpar = "a") + prior(normal(1, 0.5), nlpar = "b")
fit1 <- brm(bf(y ~ a * cd^b , a~(1|File), b ~(1|File), nl = TRUE),data = df,prior = prior1)

数据:

dput(df)
structure(list(File = c("001", "001", "001", "001", "001", "001", 
"001", "001", "001", "001", "001", "001", "001", "001", "001", 
"001", "001", "001", "001", "001", "001", "001", "001", "001", 
"001", "001", "001", "001", "001", "001", "001", "001", "001", 
"001", "001", "001", "001", "001", "001", "001", "002", "002", 
"002", "002", "002", "002", "002", "002", "002", "002", "002", 
"002", "002", "002", "002", "002", "002", "002", "002", "002", 
"002", "002", "002", "002", "002", "002", "002", "002", "002", 
"002", "002", "002", "002", "002", "002", "002", "002", "002", 
"002", "002", "003", "003", "003", "003", "003", "003", "003", 
"003", "003", "003", "003", "003", "003", "003", "003", "003", 
"003", "003", "003", "003", "003", "003", "003", "003", "003", 
"003", "003", "003", "003", "003", "003", "003", "003", "003", 
"003", "003", "003", "003", "003", "003", "004", "004", "004", 
"004", "004", "004", "004", "004", "004", "004", "004", "004", 
"004", "004", "004", "004", "004", "004", "004", "004", "004", 
"004", "004", "004", "004", "004", "004", "004", "004", "004", 
"004", "004", "004", "004", "004", "004", "004", "004", "004", 
"004", "005", "005", "005", "005", "005", "005", "005", "005", 
"005", "005", "005", "005", "005", "005", "005", "005", "005", 
"005", "005", "005", "005", "005", "005", "005", "005", "005", 
"005", "005", "005", "005", "005", "005", "005", "005", "005", 
"005", "005", "005", "005", "005", "006", "006", "006", "006", 
"006", "006", "006", "006", "006", "006", "006", "006", "006", 
"006", "006", "006", "006", "006", "006", "006", "006", "006", 
"006", "006", "006", "006", "006", "006", "006", "006", "006", 
"006", "006", "006", "006", "006", "006", "006", "006", "006", 
"007", "007", "007", "007", "007", "007", "007", "007", "007", 
"007", "007", "007", "007", "007", "007", "007", "007", "007", 
"007", "007", "007", "007", "007", "007", "007", "007", "007", 
"007", "007", "007", "007", "007", "007", "007", "007", "007", 
"007", "007", "007", "007", "008", "008", "008", "008", "008", 
"008", "008", "008", "008", "008", "008", "008", "008", "008", 
"008", "008", "008", "008", "008", "008", "008", "008", "008", 
"008", "008", "008", "008", "008", "008", "008", "008", "008", 
"008", "008", "008", "008", "008", "008", "009", "009", "009", 
"009", "009", "009", "009", "009", "009", "009", "009", "009", 
"009", "009", "009", "009", "009", "009", "009", "009", "009", 
"009", "009", "009", "009", "009", "009", "009", "009", "009", 
"009", "009", "009", "009", "009", "009", "009", "009", "009", 
"009", "010", "010", "010", "010", "010", "010", "010", "010", 
"010", "010", "010", "010", "010", "010", "010", "010", "010", 
"010", "010", "010", "010", "010", "010", "010", "010", "010", 
"010", "010", "010", "010", "010", "010", "010", "010", "010", 
"010", "010", "010", "010", "010", "011", "011", "011", "011", 
"011", "011", "011", "011", "011", "011", "011", "011", "011", 
"011", "011", "011", "011", "011", "011", "011", "011", "011", 
"011", "011", "011", "011", "011", "011", "011", "011", "011", 
"011", "011", "011", "011", "011", "011", "011", "011", "011", 
"012", "012", "012", "012", "012", "012", "012", "012", "012", 
"012", "012", "012", "012", "012", "012", "012", "012", "012", 
"012", "012", "012", "012", "012", "012", "012", "012", "012", 
"012", "012", "012", "012", "012", "012", "012", "012", "012", 
"012", "012", "012", "012", "013", "013", "013", "013", "013", 
"013", "013", "013", "013", "013", "013", "013", "013", "013", 
"013", "013", "013", "013", "013", "013", "013", "013", "013", 
"013", "013", "013", "013", "013", "013", "013", "013", "013", 
"013", "013", "013", "013", "013", "013", "013", "013", "014", 
"014", "014", "014", "014", "014", "014", "014", "014", "014", 
"014", "014", "014", "014", "014", "014", "014", "014", "014", 
"014", "014", "014", "014", "014", "014", "014", "014", "014", 
"014", "014", "014", "014", "014", "014", "014", "014", "014", 
"014", "014", "014", "015", "015", "015", "015", "015", "015", 
"015", "015", "015", "015", "015", "015", "015", "015", "015", 
"015", "015", "015", "015", "015", "015", "015", "015", "015", 
"015", "015", "015", "015", "015", "015", "015", "015", "015", 
"015", "015", "015", "015", "015", "015", "015", "016", "016", 
"016", "016", "016", "016", "016", "016", "016", "016", "016", 
"016", "016", "016", "016", "016", "016", "016", "016", "016", 
"016", "016", "016", "016", "016", "016", "016", "016", "016", 
"016", "016", "016", "016", "016", "016", "016", "016", "016", 
"016", "017", "017", "017", "017", "017", "017", "017", "017", 
"017", "017", "017", "017", "017", "017", "017", "017", "017", 
"017", "017", "017", "017", "017", "017", "017", "017", "017", 
"017", "017", "017", "017", "017", "017", "017", "017", "017", 
"017", "017", "017", "017", "017", "019", "019", "019", "019", 
"019", "019", "019", "019", "019", "019", "019", "019", "019", 
"019", "019", "019", "019", "019", "019", "019", "019", "019", 
"019", "019", "019", "019", "019", "019", "019", "019", "019", 
"019", "019", "019", "019", "019", "019", "019", "019", "019", 
"020", "020", "020", "020", "020", "020", "020", "020", "020", 
"020", "020", "020", "020", "020", "020", "020", "020", "020", 
"020", "020", "020", "020", "020", "020", "020", "020", "020", 
"020", "020", "020", "020", "020", "020", "020", "020", "020", 
"020", "020", "020", "020", "021", "021", "021", "021", "021", 
"021", "021", "021", "021", "021", "021", "021", "021", "021", 
"021", "021", "021", "021", "021", "021", "021", "021", "021", 
"021", "021", "021", "021", "021", "021", "021", "021", "021", 
"021", "021", "021", "021", "021", "021", "021", "021", "022", 
"022", "022", "022", "022", "022", "022", "022", "022", "022", 
"022", "022", "022", "022", "022", "022", "022", "022", "022", 
"022", "022", "022", "022", "022", "022", "022", "022", "022", 
"022", "022", "022", "022", "022", "022", "022", "022", "022", 
"022", "022", "022", "023", "023", "023", "023", "023", "023", 
"023", "023", "023", "023", "023", "023", "023", "023", "023", 
"023", "023", "023", "023", "023", "023", "023", "023", "023", 
"023", "023", "023", "023", "023", "023", "023", "023", "023", 
"023", "023", "023", "023", "023", "023", "023", "024", "024", 
"024", "024", "024", "024", "024", "024", "024", "024", "024", 
"024", "024", "024", "024", "024", "024", "024", "024", "024", 
"024", "024", "024", "024", "024", "024", "024", "024", "024", 
"024", "024", "024", "024", "024", "024", "024", "024", "024", 
"024", "024", "025", "025", "025", "025", "025", "025", "025", 
"025", "025", "025", "025", "025", "025", "025", "025", "025", 
"025", "025", "025", "025", "025", "025", "025", "025", "025", 
"025", "025", "025", "025", "025", "025", "025", "025", "025", 
"025", "025", "025", "025", "025", "025", "026", "026", "026", 
"026", "026", "026", "026", "026", "026", "026", "026", "026", 
"026", "026", "026", "026", "026", "026", "026", "026", "026", 
"026", "026", "026", "026", "026", "026", "026", "026", "026", 
"026", "026", "026", "026", "026", "026", "026", "026", "026", 
"026", "027", "027", "027", "027", "027", "027", "027", "027", 
"027", "027", "027", "027", "027", "027", "027", "027", "027", 
"027", "027", "027", "027", "027", "027", "027", "027", "027", 
"027", "027", "027", "027", "027", "027", "027", "027", "027", 
"027", "027", "027", "027", "027", "028", "028", "028", "028", 
"028", "028", "028", "028", "028", "028", "028", "028", "028", 
"028", "028", "028", "028", "028", "028", "028", "028", "028", 
"028", "028", "028", "028", "028", "028", "028", "028", "028", 
"028", "028", "028", "028", "028", "028", "028", "028", "028", 
"029", "029", "029", "029", "029", "029", "029", "029", "029", 
"029", "029", "029", "029", "029", "029", "029", "029", "029", 
"029", "029", "029", "029", "029", "029", "029", "029", "029", 
"029", "029", "029", "029", "029", "029", "029", "029", "029", 
"029", "029", "029", "029", "031", "031", "031", "031", "031", 
"031", "031", "031", "031", "031", "031", "031", "031", "031", 
"031", "031", "031", "031", "031", "031", "031", "031", "031", 
"031", "031", "031", "031", "031", "031", "031", "031", "031", 
"031", "031", "031", "031", "031", "031", "031", "031", "032", 
"032", "032", "032", "032", "032", "032", "032", "032", "032", 
"032", "032", "032", "032", "032", "032", "032", "032", "032", 
"032", "032", "032", "032", "032", "032", "032", "032", "032", 
"032", "032", "032", "032", "032", "032", "032", "032", "032", 
"032", "032", "032", "033", "033", "033", "033", "033", "033", 
"033", "033", "033", "033", "033", "033", "033", "033", "033", 
"033", "033", "033", "033", "033", "033", "033", "033", "033", 
"033", "033", "033", "033", "033", "033", "033", "033", "033", 
"033", "033", "033", "033", "033", "033", "033", "034", "034", 
"034", "034", "034", "034", "034", "034", "034", "034", "034", 
"034", "034", "034", "034", "034", "034", "034", "034", "034", 
"034", "034", "034", "034", "034", "034", "034", "034", "034", 
"034", "034", "034", "034", "034", "034", "034", "034", "034", 
"034", "034", "035", "035", "035", "035", "035", "035", "035", 
"035", "035", "035", "035", "035", "035", "035", "035", "035", 
"035", "035", "035", "035", "035", "035", "035", "035", "035", 
"035", "035", "035", "035", "035", "035", "035", "035", "035", 
"035", "035", "035", "035", "035", "035", "036", "036", "036", 
"036", "036", "036", "036", "036", "036", "036", "036", "036", 
"036", "036", "036", "036", "036", "036", "036", "036", "036", 
"036", "036", "036", "036", "036", "036", "036", "036", "036", 
"036", "036", "036", "036", "036", "036", "036", "036", "036", 
"036"), cd = c(57.2, 32, 3.2, 57.2, 17.9, 1.8, 32, 1, 32, 1.8, 
1.8, 5.7, 100, 1, 3.2, 3.2, 5.7, 5.7, 57.2, 17.9, 17.9, 17.9, 
100, 1, 1, 32, 3.2, 57.2, 32, 100, 100, 57.2, 1.8, 5.7, 1.8, 
1, 100, 5.7, 17.9, 3.2, 1, 5.7, 17.9, 5.7, 100, 3.2, 1, 32, 17.9, 
57.2, 3.2, 57.2, 17.9, 17.9, 1.8, 3.2, 32, 1.8, 3.2, 3.2, 57.2, 
1.8, 57.2, 5.7, 100, 1, 57.2, 1.8, 100, 1, 1, 17.9, 5.7, 32, 
100, 100, 1.8, 32, 32, 5.7, 32, 5.7, 17.9, 32, 3.2, 57.2, 5.7, 
17.9, 3.2, 5.7, 17.9, 3.2, 1.8, 57.2, 1, 3.2, 100, 100, 32, 1, 
57.2, 57.2, 1.8, 32, 1, 100, 100, 32, 1, 3.2, 1.8, 1.8, 17.9, 
17.9, 57.2, 5.7, 100, 5.7, 1.8, 1, 3.2, 5.7, 32, 57.2, 1, 1.8, 
3.2, 17.9, 100, 1.8, 32, 57.2, 3.2, 100, 5.7, 1.8, 32, 1, 32, 
100, 17.9, 100, 1.8, 1, 5.7, 32, 5.7, 1, 17.9, 57.2, 57.2, 57.2, 
17.9, 17.9, 5.7, 1.8, 1, 100, 3.2, 3.2, 57.2, 57.2, 17.9, 1.8, 
3.2, 3.2, 100, 3.2, 32, 57.2, 5.7, 1, 32, 1.8, 5.7, 17.9, 1, 
5.7, 32, 1.8, 57.2, 32, 1.8, 17.9, 1, 3.2, 100, 1, 100, 100, 
1, 17.9, 3.2, 1.8, 57.2, 100, 32, 5.7, 5.7, 17.9, 100, 57.2, 
3.2, 17.9, 5.7, 3.2, 1.8, 5.7, 1.8, 100, 1, 32, 32, 5.7, 1.8, 
1.8, 100, 17.9, 57.2, 17.9, 1, 3.2, 17.9, 32, 100, 17.9, 5.7, 
5.7, 57.2, 100, 3.2, 1, 1, 57.2, 3.2, 32, 57.2, 1, 1.8, 32, 3.2, 
5.7, 57.2, 1, 3.2, 1.8, 1.8, 100, 32, 57.2, 5.7, 1.8, 1, 57.2, 
5.7, 100, 17.9, 32, 3.2, 1, 17.9, 17.9, 100, 57.2, 1, 32, 3.2, 
32, 100, 1.8, 3.2, 57.2, 5.7, 17.9, 32, 100, 5.7, 1, 1.8, 17.9, 
5.7, 3.2, 17.9, 3.2, 5.7, 32, 17.9, 17.9, 3.2, 100, 5.7, 1.8, 
57.2, 1.8, 1.8, 57.2, 5.7, 1.8, 3.2, 100, 100, 5.7, 32, 57.2, 
32, 1, 57.2, 3.2, 1, 32, 100, 1, 100, 57.2, 17.9, 32, 17.9, 1.8, 
1, 17.9, 32, 3.2, 100, 32, 1.8, 3.2, 5.7, 1.8, 1, 1.8, 3.2, 17.9, 
3.2, 32, 17.9, 100, 57.2, 1.8, 1.8, 5.7, 100, 5.7, 32, 5.7, 17.9, 
1, 5.7, 57.2, 57.2, 100, 57.2, 100, 17.9, 1, 57.2, 1, 32, 3.2, 
1.8, 57.2, 3.2, 1, 1.8, 32, 1.8, 3.2, 5.7, 5.7, 1.8, 100, 32, 
1, 5.7, 100, 5.7, 3.2, 3.2, 17.9, 1, 100, 17.9, 1.8, 100, 17.9, 
57.2, 100, 3.2, 32, 32, 57.2, 32, 1, 57.2, 1, 17.9, 17.9, 57.2, 
5.7, 17.9, 100, 1, 1.8, 1.8, 3.2, 32, 5.7, 100, 57.2, 5.7, 5.7, 
57.2, 17.9, 3.2, 100, 1, 57.2, 17.9, 1.8, 3.2, 1, 17.9, 32, 5.7, 
57.2, 100, 1.8, 5.7, 32, 3.2, 57.2, 17.9, 1, 100, 32, 3.2, 1, 
32, 1.8, 17.9, 100, 5.7, 1.8, 1, 1, 3.2, 57.2, 57.2, 32, 32, 
32, 57.2, 5.7, 100, 100, 32, 1.8, 100, 1.8, 3.2, 5.7, 1.8, 17.9, 
3.2, 17.9, 1, 3.2, 17.9, 17.9, 1, 1, 57.2, 1.8, 57.2, 100, 32, 
5.7, 5.7, 3.2, 32, 17.9, 1.8, 100, 100, 1.8, 3.2, 57.2, 1, 100, 
3.2, 1.8, 5.7, 1, 5.7, 3.2, 17.9, 1.8, 3.2, 1.8, 1, 32, 17.9, 
5.7, 17.9, 57.2, 5.7, 32, 57.2, 17.9, 57.2, 100, 1, 57.2, 32, 
1, 32, 3.2, 100, 5.7, 3.2, 100, 57.2, 1, 17.9, 1, 1, 100, 1.8, 
17.9, 1, 5.7, 3.2, 1.8, 32, 3.2, 5.7, 5.7, 17.9, 3.2, 1.8, 100, 
57.2, 100, 3.2, 32, 5.7, 32, 57.2, 1.8, 57.2, 17.9, 100, 5.7, 
1.8, 57.2, 32, 17.9, 32, 1, 5.7, 5.7, 100, 32, 1, 1, 5.7, 5.7, 
1.8, 1.8, 17.9, 1.8, 57.2, 3.2, 17.9, 5.7, 32, 57.2, 1, 100, 
32, 17.9, 57.2, 17.9, 3.2, 100, 17.9, 32, 100, 57.2, 3.2, 57.2, 
1, 3.2, 1.8, 3.2, 1, 32, 1.8, 100, 100, 1, 32, 32, 57.2, 5.7, 
5.7, 57.2, 17.9, 57.2, 3.2, 1, 1.8, 17.9, 3.2, 3.2, 100, 5.7, 
5.7, 3.2, 57.2, 100, 3.2, 17.9, 100, 1, 1.8, 17.9, 57.2, 1.8, 
32, 100, 1.8, 32, 32, 1, 1.8, 5.7, 17.9, 100, 32, 3.2, 1.8, 17.9, 
3.2, 5.7, 17.9, 57.2, 3.2, 1.8, 32, 100, 57.2, 32, 1.8, 1, 57.2, 
5.7, 1, 1, 17.9, 1, 5.7, 17.9, 1.8, 1.8, 3.2, 100, 100, 5.7, 
32, 57.2, 5.7, 3.2, 1, 100, 32, 57.2, 17.9, 32, 1.8, 3.2, 1, 
1, 17.9, 57.2, 1, 100, 100, 57.2, 1.8, 17.9, 5.7, 1.8, 32, 3.2, 
32, 5.7, 100, 32, 100, 5.7, 1.8, 57.2, 100, 5.7, 17.9, 1, 5.7, 
57.2, 3.2, 57.2, 1, 3.2, 17.9, 32, 3.2, 17.9, 1.8, 3.2, 3.2, 
5.7, 1.8, 17.9, 17.9, 1, 5.7, 1.8, 3.2, 5.7, 32, 57.2, 1, 100, 
5.7, 1.8, 100, 100, 17.9, 1, 1.8, 32, 32, 1, 17.9, 57.2, 57.2, 
32, 3.2, 3.2, 1, 32, 100, 100, 5.7, 57.2, 17.9, 1.8, 57.2, 1.8, 
57.2, 57.2, 100, 32, 1, 1, 5.7, 3.2, 5.7, 32, 32, 1.8, 1, 17.9, 
1.8, 17.9, 57.2, 3.2, 5.7, 57.2, 17.9, 17.9, 1.8, 5.7, 32, 3.2, 
3.2, 57.2, 32, 100, 100, 1, 17.9, 100, 1.8, 5.7, 1, 3.2, 100, 
3.2, 32, 100, 100, 1.8, 5.7, 57.2, 5.7, 1, 32, 32, 57.2, 17.9, 
100, 3.2, 5.7, 5.7, 1.8, 5.7, 57.2, 3.2, 1, 3.2, 1.8, 57.2, 1, 
1, 100, 1, 100, 57.2, 3.2, 32, 32, 1.8, 1.8, 17.9, 17.9, 17.9, 
17.9, 1.8, 17.9, 5.7, 100, 1.8, 3.2, 3.2, 5.7, 57.2, 1.8, 17.9, 
57.2, 57.2, 32, 3.2, 3.2, 1, 17.9, 100, 32, 17.9, 3.2, 57.2, 
1, 32, 5.7, 17.9, 5.7, 1, 1, 1, 5.7, 100, 100, 1.8, 1.8, 57.2, 
100, 32, 32, 100, 17.9, 57.2, 1, 3.2, 17.9, 3.2, 100, 32, 5.7, 
1.8, 1, 57.2, 3.2, 32, 1.8, 1, 1, 57.2, 100, 32, 5.7, 32, 1, 
100, 1.8, 5.7, 57.2, 3.2, 1.8, 100, 5.7, 17.9, 57.2, 3.2, 5.7, 
32, 1.8, 17.9, 17.9, 1, 5.7, 100, 3.2, 1, 1.8, 100, 5.7, 1.8, 
1, 17.9, 17.9, 32, 3.2, 32, 3.2, 5.7, 1.8, 57.2, 100, 57.2, 32, 
32, 3.2, 32, 17.9, 1, 17.9, 100, 17.9, 5.7, 57.2, 1.8, 57.2, 
5.7, 3.2, 100, 1.8, 1, 57.2, 1.8, 32, 17.9, 32, 17.9, 1, 1, 5.7, 
32, 3.2, 57.2, 57.2, 32, 5.7, 100, 100, 57.2, 3.2, 1.8, 3.2, 
1.8, 17.9, 1, 1, 5.7, 100, 5.7, 100, 100, 3.2, 1, 1.8, 17.9, 
57.2, 57.2, 5.7, 17.9, 1.8, 3.2, 32, 57.2, 100, 57.2, 5.7, 1, 
5.7, 1, 17.9, 5.7, 3.2, 17.9, 1.8, 17.9, 32, 17.9, 100, 57.2, 
100, 3.2, 1.8, 32, 1, 3.2, 1, 32, 1.8, 1, 57.2, 32, 100, 57.2, 
3.2, 3.2, 5.7, 1.8, 100, 1.8, 32, 17.9, 5.7, 57.2, 5.7, 5.7, 
57.2, 100, 1.8, 5.7, 100, 100, 57.2, 17.9, 100, 3.2, 32, 57.2, 
32, 3.2, 32, 1.8, 1, 32, 1.8, 17.9, 1, 5.7, 1.8, 3.2, 17.9, 5.7, 
17.9, 1, 1.8, 17.9, 100, 3.2, 32, 3.2, 1, 57.2, 1, 5.7, 17.9, 
3.2, 100, 5.7, 32, 1, 32, 100, 100, 3.2, 17.9, 1.8, 57.2, 17.9, 
5.7, 17.9, 3.2, 3.2, 32, 100, 3.2, 5.7, 57.2, 1, 57.2, 100, 32, 
1.8, 1, 1, 1, 17.9, 1.8, 1.8, 5.7, 1.8, 57.2, 57.2, 32, 1.8, 
3.2, 1, 57.2, 3.2, 100, 17.9, 1.8, 17.9, 100, 3.2, 3.2, 57.2, 
5.7, 17.9, 57.2, 1, 5.7, 32, 1.8, 1, 100, 5.7, 32, 32, 32, 32, 
3.2, 1, 57.2, 1.8, 5.7, 100, 17.9, 1, 57.2, 1.8, 100, 17.9, 5.7, 
1, 57.2, 5.7, 57.2, 32, 1.8, 100, 17.9, 57.2, 17.9, 17.9, 5.7, 
100, 57.2, 3.2, 5.7, 1.8, 32, 1.8, 1, 32, 3.2, 100, 1, 3.2, 57.2, 
32, 100, 1, 5.7, 5.7, 1.8, 3.2, 1.8, 17.9, 32, 3.2, 100, 17.9, 
1, 1, 1.8, 57.2, 32, 57.2, 32, 17.9, 5.7, 1.8, 57.2, 1.8, 17.9, 
1, 100, 1.8, 32, 17.9, 57.2, 3.2, 5.7, 32, 1, 3.2, 57.2, 100, 
1, 100, 3.2, 17.9, 5.7, 1, 3.2, 3.2, 1.8, 17.9, 5.7, 32, 5.7, 
100, 100, 3.2, 100, 3.2, 1.8, 17.9, 1, 1.8, 32, 100, 100, 57.2, 
57.2, 32, 17.9, 1, 3.2, 1, 57.2, 1.8, 17.9, 1, 32, 5.7, 5.7, 
5.7, 32, 5.7, 57.2, 17.9, 3.2, 32, 57.2, 100, 1.8, 1.8, 3.2, 
5.7, 1, 17.9, 100, 17.9, 3.2, 1, 100, 3.2, 17.9, 32, 5.7, 1, 
100, 1.8, 57.2, 5.7, 5.7, 100, 3.2, 17.9, 1.8, 57.2, 1, 100, 
32, 57.2, 57.2, 32, 5.7, 17.9, 32, 1, 1.8, 5.7, 1.8, 3.2, 1, 
17.9, 100, 3.2, 57.2, 1.8, 32, 32, 5.7, 1.8, 100, 32, 100, 5.7, 
17.9, 100, 5.7, 1, 1, 17.9, 17.9, 1, 17.9, 100, 1, 1.8, 3.2, 
17.9, 57.2, 1.8, 57.2, 32, 3.2, 1.8, 32, 3.2, 57.2, 5.7, 3.2, 
3.2, 100, 57.2, 1, 32, 57.2, 1.8, 5.7), y = c(15, 13, 7.69230769230769, 
15, 11, 5, 12, 3.33333333333333, 12, 6.66666666666667, 6.66666666666667, 
9.09090909090909, 25, 5, 7.69230769230769, 7.69230769230769, 
8.33333333333333, 8.33333333333333, 15, 11, 11, 11, 20, 6.66666666666667, 
6.66666666666667, 13, 7.14285714285714, 13, 12, 16, 15, 13, 7.69230769230769, 
8.33333333333333, 6.25, 5, 15, 8.33333333333333, 11, 8.33333333333333, 
8, 9, 6, 7, 20, 5, 3, 15, 12, 18, 3, 14, 11, 10.5, 5, 6, 12, 
5, 6, 6, 20, 5, 20, 8, 20, 5, 17, 4, 20, 5, 4, 11, 9, 12, 18, 
18, 5, 13, 14, 7, 20, 6.66666666666667, 13, 14, 6.25, 20, 8.33333333333333, 
11, 5.88235294117647, 7.69230769230769, 11, 6.66666666666667, 
5.55555555555556, 20, 5, 5.55555555555556, 20, 19, 14, 5.55555555555556, 
14, 13, 6.66666666666667, 14, 5.55555555555556, 18, 16, 11, 5.26315789473684, 
7.69230769230769, 7.69230769230769, 7.14285714285714, 11, 11, 
14, 9.09090909090909, 15, 7.69230769230769, 7.14285714285714, 
6.66666666666667, 7, 7, 13, 15, 2, 3, 4, 11, 20, 2, 12, 16, 5, 
20, 7, 3, 12, 3, 14, 20, 11, 20, 2, 1, 8, 12, 8, 2, 10, 20, 19, 
19, 10, 10, 7, 3, 1, 20, 7, 8, 20, 15, 11, 5, 5, 5, 20, 5, 15, 
20, 6.66666666666667, 5, 12, 6.66666666666667, 8.33333333333333, 
11, 2.85714285714286, 6.66666666666667, 15, 6.66666666666667, 
20, 12, 6.66666666666667, 11, 5, 6.66666666666667, 25, 4, 25, 
20, 2.85714285714286, 12, 6.66666666666667, 5.88235294117647, 
20, 30, 12, 7.69230769230769, 6.66666666666667, 11, 50, 30, 2.5, 
12.5, 3, 1, 0.5, 5, 1, 100, 0.5, 15, 12.5, 11, 4, 2, 40, 5, 30, 
12.5, 1, 7.5, 11, 15, 50, 11, 7.5, 7.5, 30, 40, 5, 1, 1, 30, 
5, 15, 20, 2, 2, 15, 6.66666666666667, 8.33333333333333, 20, 
3.33333333333333, 6.66666666666667, 5, 4, 40, 15, 20, 8.33333333333333, 
5, 3.33333333333333, 20, 8.33333333333333, 50, 11, 15, 6.66666666666667, 
3.33333333333333, 11, 11, 50, 30, 3.33333333333333, 15, 6.66666666666667, 
12, 40, 5, 6.66666666666667, 20, 8.33333333333333, 11, 14, 50, 
6.66666666666667, 2, 2.22222222222222, 11, 5, 3, 11, 8, 9.5, 
13, 11.8, 10.5, 3.5, 20, 9.5, 6, 16, 5, 4.5, 13, 9.25, 6.5, 9.25, 
17, 18, 8, 10.75, 12, 11, 3, 14, 7, 2, 12, 18, 1, 18, 12, 10.5, 
11, 10.1, 5.5, 4, 14, 17, 5.55555555555556, 30, 20, 5.88235294117647, 
6.66666666666667, 8.33333333333333, 5.55555555555556, 4, 5, 6.66666666666667, 
12, 7.14285714285714, 16, 12, 30, 25, 5, 4.54545454545455, 8.33333333333333, 
23, 8.33333333333333, 13, 7.14285714285714, 11, 4.34782608695652, 
7.69230769230769, 17, 16, 22, 20, 25, 11, 4, 17, 4, 13, 6.66666666666667, 
7, 16, 8, 4, 5, 14, 7, 7, 9, 9, 6, 30, 20, 4, 8, 30, 9, 9, 9, 
6, 3, 30, 11, 6, 14, 11, 16, 20, 8, 12, 12, 14, 13, 7, 13, 4, 
11, 11, 15, 7, 14, 100, 2, 2.5, 2.5, 3.33333333333333, 20, 5, 
90, 30, 6.25, 6.25, 20, 12, 5, 40, 1.11111111111111, 40, 11, 
3.33333333333333, 5, 2, 12, 20, 7.14285714285714, 40, 90, 3.33333333333333, 
6.25, 20, 5, 30, 12, 3.33333333333333, 80, 30, 5, 1.66666666666667, 
15, 3.33333333333333, 12, 50, 6, 4, 2.5, 2.5, 6, 20, 18, 15, 
14, 17, 25, 7, 23, 24, 15, 5, 30, 4, 6, 7, 4, 11, 7, 12, 4, 6, 
11, 12, 4, 4, 18, 6, 20, 30, 13, 7, 7, 6, 15, 12, 5, 30, 30, 
6.66666666666667, 7.69230769230769, 13, 3.33333333333333, 25, 
5.55555555555556, 6.66666666666667, 8.33333333333333, 5, 8.33333333333333, 
6.66666666666667, 12, 5.88235294117647, 7.69230769230769, 6.66666666666667, 
5, 15, 12, 9.09090909090909, 11, 20, 9.09090909090909, 13, 18, 
11, 20, 20, 6.66666666666667, 15, 12, 5.55555555555556, 13, 6.66666666666667, 
15, 9.09090909090909, 8, 20, 14, 4, 11, 2, 2, 20, 4, 11, 3, 9, 
9, 9, 12, 8, 9, 9.5, 11, 8, 4, 20, 17, 18, 7, 14, 9, 11, 12, 
7, 13, 10.5, 19, 9.5, 6, 18, 12, 11, 11.5, 2, 5, 6.66666666666667, 
25, 20, 2.5, 2.22222222222222, 6.66666666666667, 8.33333333333333, 
5, 5, 12, 5.88235294117647, 20, 6.66666666666667, 12, 8.33333333333333, 
12, 30, 2, 60, 18, 12, 20, 13, 5, 45, 11, 14, 25, 25, 4, 30, 
1.66666666666667, 2.5, 2, 2.5, 1.42857142857143, 20, 3.33333333333333, 
80, 20, 5, 15, 13, 14, 8, 9, 17, 9, 19, 5, 1, 2, 11, 7, 7, 20, 
9, 9, 7, 19, 20, 6, 6, 20, 2, 2, 3, 18, 1, 13, 17, 1, 6, 13, 
1, 2, 8, 11, 30, 13, 5, 5, 11, 6.66666666666667, 7.69230769230769, 
11, 20, 5.55555555555556, 5, 15, 25, 18, 13, 4, 2.85714285714286, 
20, 6.66666666666667, 3.33333333333333, 3.33333333333333, 11, 
3.84615384615385, 8.33333333333333, 11, 5, 4.16666666666667, 
6.66666666666667, 30, 30, 8.33333333333333, 15, 25, 7.14285714285714, 
5, 3.33333333333333, 30, 23, 28, 18, 15, 6.66666666666667, 7.14285714285714, 
5, 5.88235294117647, 12, 15, 5.88235294117647, 18, 20, 16, 8.33333333333333, 
11, 7.69230769230769, 7.14285714285714, 16, 7.69230769230769, 
12, 8.33333333333333, 18, 14, 17, 8.33333333333333, 7.14285714285714, 
17, 17, 8.33333333333333, 11, 7.69230769230769, 8.33333333333333, 
16, 7.69230769230769, 15, 7.69230769230769, 8.33333333333333, 
11, 13, 8.33333333333333, 11, 8.33333333333333, 5, 7, 9, 3, 12, 
13, 2, 8, 3, 5, 9, 20, 21, 4, 30, 7, 3, 20, 21, 11, 5, 6, 15, 
13, 5, 11, 20, 21, 13, 8, 7, 4, 14, 20, 21, 8, 15, 11, 6, 20, 
5, 20, 15, 20, 12, 6.66666666666667, 6.66666666666667, 9.09090909090909, 
9.09090909090909, 9.09090909090909, 11, 11, 7.69230769230769, 
7.14285714285714, 12, 7.69230769230769, 11, 14, 8.33333333333333, 
7.69230769230769, 14, 11, 11, 6.66666666666667, 8.33333333333333, 
13, 7.69230769230769, 7.69230769230769, 14, 13, 20, 20, 5, 12, 
20, 7.14285714285714, 8.33333333333333, 5, 6.66666666666667, 
20, 7, 12, 30, 30, 6, 7, 13, 7, 1, 8, 7, 8, 9, 13, 6, 7, 8, 3, 
7, 13, 5, 2, 3, 2, 15, 1, 2, 17, 2, 18, 15, 8, 14, 13, 4, 5, 
11, 11, 12, 12, 4, 15, 6.66666666666667, 40, 3.33333333333333, 
5, 5, 6.66666666666667, 30, 3.33333333333333, 15, 30, 30, 20, 
5, 4, 3.33333333333333, 13, 30, 20, 12, 5, 30, 2.85714285714286, 
17, 7.69230769230769, 11, 7.14285714285714, 2.85714285714286, 
2.85714285714286, 2.85714285714286, 6.66666666666667, 30, 30, 
5, 5, 25, 30, 20, 18, 20, 12, 17, 2, 4, 6, 5, 19, 15, 9, 3, 2, 
17, 7, 14, 4, 3, 4, 16, 18, 13, 9, 13, 1, 20, 4, 7, 19, 6, 5, 
20, 9, 11, 17, 7, 8, 13, 6, 11, 12, 5, 8.33333333333333, 40, 
6.66666666666667, 3.33333333333333, 6.66666666666667, 35, 7.69230769230769, 
5, 2.5, 11, 12, 15, 5.88235294117647, 15, 5.26315789473684, 7.14285714285714, 
4, 23, 30, 20, 17, 16, 4.54545454545455, 17, 11, 2, 12, 35, 11, 
7.14285714285714, 25, 3.33333333333333, 30, 5.55555555555556, 
3.84615384615385, 40, 2.63157894736842, 2.32558139534884, 25, 
5, 15, 12, 14, 11, 2, 2, 7, 13, 6, 18, 18, 15, 7, 20, 19, 16, 
7, 5, 6, 5, 11, 2, 3, 8, 20, 8, 20, 20, 7, 5, 6, 11, 19, 19, 
8, 11, 6, 7, 13, 17, 20, 16, 7.69230769230769, 5, 8.33333333333333, 
4.54545454545455, 11, 8.33333333333333, 7.14285714285714, 11, 
5.55555555555556, 11, 12, 11, 24, 18, 18, 6.66666666666667, 5.88235294117647, 
12, 5, 7.14285714285714, 4.76190476190476, 12, 5.26315789473684, 
4.54545454545455, 16, 12, 18, 15, 7.14285714285714, 7.14285714285714, 
9.09090909090909, 7.69230769230769, 20, 6.66666666666667, 12, 
11, 8.33333333333333, 20, 8, 8, 17, 25, 5, 7, 23, 22, 19, 13, 
23, 5, 11, 15, 12, 6, 11, 4, 3, 12, 4, 10.5, 3.2, 9.8, 6, 6, 
10.4, 9.7, 10.7, 2.6, 7, 10.2, 23, 7.6, 11, 9, 2.3, 12.4, 2, 
5, 13, 5.55555555555556, 35, 6.66666666666667, 17, 2.5, 15, 30, 
30, 5, 12, 6.66666666666667, 20, 11, 6.66666666666667, 12, 6.66666666666667, 
6.25, 15, 30, 6.66666666666667, 7.69230769230769, 25, 3.33333333333333, 
25, 30, 15, 5.55555555555556, 4, 3.33333333333333, 3.33333333333333, 
12, 6.66666666666667, 6.66666666666667, 8.33333333333333, 6.66666666666667, 
25, 25, 20, 6.66666666666667, 5, 5, 20, 7.69230769230769, 20, 
11, 5.55555555555556, 11, 20, 7.14285714285714, 8.33333333333333, 
14, 7.14285714285714, 11, 17, 4.54545454545455, 8.33333333333333, 
13, 5.55555555555556, 4, 20, 7.14285714285714, 12, 11, 11, 12, 
7.69230769230769, 4.16666666666667, 15, 5, 8.33333333333333, 
20, 11, 4.34782608695652, 14, 5.55555555555556, 20, 11, 8.33333333333333, 
3, 20, 7, 18, 15, 5, 20, 11, 15, 11, 11, 7, 20, 15, 5, 7, 3, 
15, 3, 2, 15, 7, 20, 1, 6, 17, 15, 20, 3, 7, 7, 5, 6, 5, 11, 
15, 6, 18, 11, 2, 4, 5, 30, 15, 20, 15, 11, 6.66666666666667, 
2.85714285714286, 30, 3.57142857142857, 11, 2.5, 35, 3.33333333333333, 
15, 11, 35, 6.66666666666667, 8.33333333333333, 20, 2.85714285714286, 
3.33333333333333, 40, 50, 2.5, 50, 5, 11, 5, 2.5, 3.33333333333333, 
4, 3.33333333333333, 11, 6.66666666666667, 15, 5, 50, 40, 8, 
30, 7, 6, 11, 5, 6, 13, 20, 20, 16, 15, 12, 10, 3, 6, 4, 20, 
3, 11, 2, 12, 9, 9, 9, 15, 9, 24, 10, 7, 20, 27, 30, 4, 4, 5, 
7, 3, 9, 26, 15, 7.69230769230769, 5, 40, 5, 11, 12, 8.33333333333333, 
6.66666666666667, 15, 6.66666666666667, 18, 8.33333333333333, 
9.09090909090909, 100, 3.33333333333333, 11, 5, 30, 2.5, 30, 
13, 14, 15, 13, 7.14285714285714, 11, 13, 1.25, 3.33333333333333, 
2.5, 5, 6.66666666666667, 2, 11, 100, 7.14285714285714, 40, 5, 
30, 15, 8, 5, 20, 9, 20, 7, 11, 20, 5, 2, 1, 12, 11, 0.5, 11, 
19, 0.5, 1, 4, 11, 18, 1, 17, 14, 8, 5, 13, 6, 18, 7, 8, 7, 20, 
16, 3, 11, 17, 1, 5)), .Names = c("File", "cd", "y"), row.names = c(NA, 
-1357L), class = "data.frame")

0 个答案:

没有答案