I used the codes below trying to fit a model at .95 c.i. of size relationship but it fails.
library(quantreg) range(len$length) range(len$preyl) x<-rq(len$preyl~len$length,tau=0.95) plot(x, type="b", xlab="length (cm)",ylab="preyl (cm)",ylim=c(5,35), xlim=c(10,100), data=len) Cheers, Nuru
length weight mat sex preyl 87 4929.377934 3 1 21 87 4119.583815 3 1 25 82 3271.873879 3 1 31 77 3404.157313 3 1 22 66 1597.807053 2 0 17 60 1597.807053 2 0 15 62 1679.924649 2 0 19 66 2133.177309 2 1 24 66 1096.918041 2 1 17 65 1764.823773 2 0 18 64 1943.154908 2 0 28 69 2440.95732 2 0 13 69 2440.95732 2 0 12 63 1852.551494 2 1 11 63 1852.551494 2 1 12 64 1943.154908 2 0 26 68 2335.268203 2 0 27 58 1441.728286 2 0 22 62 1764.823773 4 0 28 64 1943.154908 2 0 18 56 1296.211412 2 0 26 73 2895.773935 2 0 29 98 7072.406331 3 1 31 68 2335.268203 2 1 18 59 1518.423939 2 0 17 79 3679.38373 2 1 31 88 5103.183732 2 1 24 78 3539.978123 3 1 12 81 3969.139107 2 1 24 82 4119.583815 3 1 9 82 4119.583815 3 1 14 64 1943.154908 2 1 17 76 3271.873879 2 0 17 76 3271.873879 2 0 12 63 1852.551494 2 0 13 63 1852.551494 2 0 12 48 812.2564608 1 1 12 28 158.4721295 1 1 12 50 919.27785 1 1 13 50 919.27785 1 1 13 50 919.27785 1 1 12 46 713.9231888 1 0 12 46 713.9231888 1 1 18 32 237.5661048 1 1 12 50 919.27785 1 1 12 50 919.27785 1 1 12 50 919.27785 1 1 13 39 432.7904124 1 0 13 48 812.2564608 2 0 12 61 1350 2 1 11.8 63.6 1500 2 0 25.9 63.4 1755 2 1 24.2 62.9 1920 3 0 27.5 56.4 1380 2 1 25.2 61.3 1570 2 0 16.5 61.3 1570 2 0 13.4 62 1550 2 1 11.5 52 1480 2 1 21.4 63 1790 2 1 28 63 1790 2 1 15.5 62.3 500 1 1 22 93 6550 4 1 33 93 6550 4 1 17.9
______________________________________________ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.
