我想使用nlme
软件包以混合效果的方式拟合幂数法。我能够使用brms
来拟合贝叶斯模型,见下文,但我想要适应频率论者'使用nlme
或其他一些能够适应非线性混合效果模型的包的相同模型的版本。 a
和b
应允许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")