This message is in MIME format. The first part should be readable text, while the remaining parts are likely unreadable without MIME-aware tools. Send mail to mime@docserver.cac.washington.edu for more info. --1465796871-645734172-974620647=:18687 Content-Type: TEXT/PLAIN; charset=US-ASCII Hi all, I've discovered and (hopefully) fixed a bug in the stars function I wrote some time ago. A corrected function is attached. Please let me know if you have any questions or comments. Tom Thomas S. Dye, Ph.D. http://www.tsdye.com Home: 812A 19th Avenue, Honolulu, Hawaii 96816. (808) 739-1367 or 387-9352. Work: International Archaeological Research Institute, Inc., 2081 Young St., Honolulu, Hawaii 96826. Voice (808) 946-2548; Fax 943-0716. --1465796871-645734172-974620647=:18687 Content-Type: TEXT/PLAIN; charset=US-ASCII; name="new-stars.r" Content-Transfer-Encoding: BASE64 Content-ID: <Pine.BSI.4.10.10011182157270.18687@malasada.lava.net> Content-Description: bug fix Content-Disposition: attachment; filename="new-stars.r" ZnVuY3Rpb24gKHgsIGZ1bGwgPSBUUlVFLCBzY2FsZSA9IFRSVUUsIHJhZGl1 cyA9IFRSVUUsIGxhYmVscyA9IGRpbW5hbWVzKHgpW1sxXV0sIA0KICAgICAg ICAgICAgbG9jYXRpb25zID0gTlVMTCwgeGxpbWl0ID0gTlVMTCwgeWxpbWl0 ID0gTlVMTCwgbGVuID0gMSwgDQogICAgICAgICAgICBjb2xvcnMgPSBOVUxM LCBrZXkubG9jID0gTlVMTCwga2V5LmxhYmVscyA9IE5VTEwsIGRyYXcuc2Vn bWVudHMgPSBGQUxTRSwgDQogICAgICAgICAgICBkcmF3LmF4ZXMgPSBGQUxT RSwgLi4uKSANCnsNCiAgaWYgKGlzLmRhdGEuZnJhbWUoeCkpIA0KICAgIHgg PC0gYXMubWF0cml4KHgpDQogIGVsc2UgaWYgKCFpcy5tYXRyaXgoeCkpIA0K ICAgIHN0b3AoInggbXVzdCBiZSBhIG1hdHJpeCBvciBhIGRhdGEgZnJhbWUi KQ0KICBpZiAoIWlzLm51bWVyaWMoeCkpIA0KICAgIHN0b3AoImRhdGEgaW4g eCBtdXN0IGJlIG51bWVyaWMiKQ0KICBuLmxvYyA8LSBucm93KHgpDQogIG4u c2VnIDwtIG5jb2woeCkNCiAgZGVnIDwtIHBpLzE4MA0KICBzZWcuY29sb3Jz IDwtIGlmICghaXMubnVsbChjb2xvcnMpKSANCiAgICBjb2xvcnMNCiAgZWxz ZSAxOm4uc2VnDQogIGlmIChpcy5udWxsKGxvY2F0aW9ucykpIHsNCiAgICBt YXQuZGltIDwtIGNlaWxpbmcoc3FydChuLmxvYykpDQogICAgdGVtcC5sb2Mu MSA8LSByZXAoeCA9IDIuMSAqIDE6bWF0LmRpbSwgdGltZXMgPSBtYXQuZGlt LCANCiAgICAgICAgICAgICAgICAgICAgICBsZW5ndGggPSBuLmxvYykNCiAg ICB0ZW1wLmxvYy4yIDwtIHJlcCh4ID0gMi4xICogbWF0LmRpbToxLCByZXAo eCA9IG1hdC5kaW0sIA0KICAgICAgICAgICAgICAgICAgICAgICAgdGltZXMg PSBtYXQuZGltKSwgbGVuZ3RoID0gbi5sb2MpDQogICAgbG9jIDwtIG1hdHJp eChkYXRhID0gYyh0ZW1wLmxvYy4xLCB0ZW1wLmxvYy4yKSwgbmNvbCA9IDIp DQogIH0NCiAgZWxzZSB7DQogICAgaWYgKCFpcy5tYXRyaXgobG9jYXRpb25z KSB8fCBuY29sKGxvY2F0aW9ucykgIT0gMikgDQogICAgICBzdG9wKCJsb2Nh dGlvbnMgbXVzdCBiZSBhIDItY29sdW1uIG1hdHJpeC4iKQ0KICAgIGxvYyA8 LSAuQWxpYXMobG9jYXRpb25zKQ0KICB9DQogIGlmIChuLmxvYyAhPSBucm93 KGxvYykpIA0KICAgIHN0b3AoIm51bWJlciBvZiByb3dzIG9mIGxvY2F0aW9u cyBhbmQgeCBtdXN0IGJlIGVxdWFsLiIpDQogIGFuZ2xlcyA8LSBpZiAoZnVs bCkgDQogICAgc2VxKDAsIDIgKiBwaSwgbGVuZ3RoID0gbi5zZWcgKyAxKVst KG4uc2VnICsgMSldDQogIGVsc2UgaWYgKGRyYXcuc2VnbWVudHMpIA0KICAg IHNlcSgwLCBwaSwgbGVuZ3RoID0gbi5zZWcgKyAxKVstKG4uc2VnICsgMSld DQogIGVsc2Ugc2VxKDAsIHBpLCBsZW5ndGggPSBuLnNlZykNCiAgaWYgKGxl bmd0aChhbmdsZXMpICE9IG4uc2VnKSANCiAgICBzdG9wKCJsZW5ndGgoYW5n bGVzKSBtdXN0IGJlIHRoZSBzYW1lIGFzIG5jb2woeCkiKQ0KICB4W2lzLm5h KHgpXSA8LSAwDQogIGlmIChzY2FsZSkgDQogICAgeCA8LSBzd2VlcCh4LCAy LCBhcHBseSh4LCAyLCBtYXgpLCBGVU4gPSAiLyIpDQojIENvbHVtbnMgb2Yg MHMgd2lsbCBwdXQgTkFzIGluIHgsIG5leHQgbGluZSBnZXRzIHJpZCBvZiB0 aGVtDQogIHhbaXMubmEoeCldIDwtIDANCiAgeCA8LSB4ICogbGVuDQogIHRl bXAueGxpbSA8LSBpZiAoaXMubnVsbCh4bGltaXQpKSANCiAgICByYW5nZShs b2NbLCAxXSArIG1heCh4KSwgbG9jWywgMV0gLSBtYXgoeCkpDQogIGVsc2Ug eGxpbWl0DQogIHRlbXAueWxpbSA8LSBpZiAoaXMubnVsbCh5bGltaXQpKSAN CiAgICByYW5nZShsb2NbLCAyXSArIG1heCh4KSwgbG9jWywgMl0gLSBtYXgo eCkpDQogIGVsc2UgeWxpbWl0DQogIG9wYXIgPC0gcGFyKG5vLnJlYWRvbmx5 ID0gVFJVRSkNCiAgb24uZXhpdChwYXIob3BhcikpDQogIHBsb3QoMCwgdHlw ZSA9ICJuIiwgLi4uLCB4bGltID0gdGVtcC54bGltLCB5bGltID0gdGVtcC55 bGltLCANCiAgICAgICB4bGFiID0gIiIsIHlsYWIgPSAiIiwgYXNwID0gMSwg YXhlcyA9IGRyYXcuYXhlcykNCiAgaWYgKGRyYXcuc2VnbWVudHMpIHsNCiAg ICBmb3IgKGkgaW4gMTpuLmxvYykgew0KICAgICAgcG9seS54IDwtIE5BDQog ICAgICBwb2x5LnkgPC0gTkENCiAgICAgIHN0YXJ0LnguY29vcmQgPC0geFtp LCBdICogY29zKGFuZ2xlcykgKyBsb2NbaSwgMV0NCiAgICAgIHN0YXJ0Lnku Y29vcmQgPC0geFtpLCBdICogc2luKGFuZ2xlcykgKyBsb2NbaSwgMl0NCiAg ICAgIGZvciAoaiBpbiAxOm4uc2VnKSB7DQogICAgICAgIHBvbHkueCA8LSBj KHBvbHkueCwgbG9jW2ksIDFdLCBzdGFydC54LmNvb3JkW2pdKQ0KICAgICAg ICBwb2x5LnkgPC0gYyhwb2x5LnksIGxvY1tpLCAyXSwgc3RhcnQueS5jb29y ZFtqXSkNCiAgICAgICAgbmV4dC5hbmdsZSA8LSBpZiAoaiA8IG4uc2VnKSAN CiAgICAgICAgICBhbmdsZXNbaiArIDFdDQogICAgICAgIGVsc2UgKGlmIChm dWxsKSANCiAgICAgICAgICAgICAgMzYwDQogICAgICAgIGVsc2UgMTgwKSAq IGRlZw0KICAgICAgICBrIDwtIHNlcShmcm9tID0gYW5nbGVzW2pdLCB0byA9 IG5leHQuYW5nbGUsIGJ5ID0gZGVnKQ0KICAgICAgICBwb2x5LnggPC0gYyhw b2x5LngsIHhbaSwgal0gKiBjb3MoaykgKyBsb2NbaSwgDQogICAgICAgICAg ICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAxXSwg TkEpDQogICAgICAgIHBvbHkueSA8LSBjKHBvbHkueSwgeFtpLCBqXSAqIHNp bihrKSArIGxvY1tpLCANCiAgICAgICAgICAgICAgICAgICAgICAgICAgICAg ICAgICAgICAgICAgICAgICAgICAgIDJdLCBOQSkNCiAgICAgIH0NCiAgICAg IHBhcihsd2QgPSAwLjI1KQ0KICAgICAgcG9seWdvbihwb2x5LngsIHBvbHku eSwgY29sID0gc2VnLmNvbG9ycykNCiAgICAgIHBhcihsd2QgPSAxKQ0KICAg ICAgaWYgKCFpcy5udWxsKGxhYmVscykpIA0KICAgICAgICB0ZXh0KGxvY1tp LCAxXSwgbG9jW2ksIDJdIC0gaWYgKGZ1bGwpIA0KICAgICAgICAgICAgIG1h eCh4KQ0KICAgICAgICBlbHNlIDAuMSAqIG1heCh4KSwgbGFiZWxzW2ldLCBj ZXggPSAwLjUsIGFkaiA9IGMoMC41LCANCiAgICAgICAgICAgICAgICAgICAg ICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgIDEpKQ0KICAgIH0NCiAg fQ0KICBlbHNlIHsNCiAgICBmb3IgKGkgaW4gMTpuLmxvYykgew0KICAgICAg dGVtcC54LmNvb3JkIDwtIHhbaSwgXSAqIGNvcyhhbmdsZXMpICsgbG9jW2ks IDFdDQogICAgICB0ZW1wLnkuY29vcmQgPC0geFtpLCBdICogc2luKGFuZ2xl cykgKyBsb2NbaSwgMl0NCiAgICAgIGlmIChyYWRpdXMpIHsNCiAgICAgICAg cGFyKGx3ZCA9IDAuMjUpDQogICAgICAgIHNlZ21lbnRzKHJlcChsb2NbaSwg MV0sIG4uc2VnKSwgcmVwKGxvY1tpLCAyXSwgDQogICAgICAgICAgICAgICAg ICAgICAgICAgICAgICAgICAgICAgICAgICAgIG4uc2VnKSwgdGVtcC54LmNv b3JkLCB0ZW1wLnkuY29vcmQpDQogICAgICAgIHBhcihsd2QgPSAxKQ0KICAg ICAgfQ0KICAgICAgbGluZXMoYyh0ZW1wLnguY29vcmQsIHRlbXAueC5jb29y ZFsxXSksIGModGVtcC55LmNvb3JkLCANCiAgICAgICAgICAgICAgICAgICAg ICAgICAgICAgICAgICAgICAgICAgICAgICAgIHRlbXAueS5jb29yZFsxXSks IGx3ZCA9IDAuMjUpDQogICAgICBpZiAoIWlzLm51bGwobGFiZWxzKSkgDQog ICAgICAgIHRleHQobG9jW2ksIDFdLCBsb2NbaSwgMl0gLSBpZiAoZnVsbCkg DQogICAgICAgICAgICAgbWF4KHgpDQogICAgICAgIGVsc2UgMC4xICogbWF4 KHgpLCBsYWJlbHNbaV0sIGNleCA9IDAuNSwgYWRqID0gYygwLjUsIA0KICAg ICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAg ICAgMSkpDQogICAgfQ0KICB9DQogIGlmICghaXMubnVsbChrZXkubG9jKSkg ew0KICAgIGlmIChkcmF3LnNlZ21lbnRzKSB7DQogICAgICBrZXkueCA8LSBO QQ0KICAgICAga2V5LnkgPC0gTkENCiAgICAgIGtleS54LmNvb3JkIDwtIGNv cyhhbmdsZXMpICogbGVuICsga2V5LmxvY1sxXQ0KICAgICAga2V5LnkuY29v cmQgPC0gc2luKGFuZ2xlcykgKiBsZW4gKyBrZXkubG9jWzJdDQogICAgICBm b3IgKGogaW4gMTpuLnNlZykgew0KICAgICAgICBrZXkueCA8LSBjKGtleS54 LCBrZXkubG9jWzFdLCBrZXkueC5jb29yZFtqXSkNCiAgICAgICAga2V5Lnkg PC0gYyhrZXkueSwga2V5LmxvY1syXSwga2V5LnkuY29vcmRbal0pDQogICAg ICAgIGsgPC0gYW5nbGVzW2pdICsgZGVnDQogICAgICAgIG5leHQuYW5nbGUg PC0gaWYgKGogPCBuLnNlZykgDQogICAgICAgICAgYW5nbGVzW2ogKyAxXQ0K ICAgICAgICBlbHNlIChpZiAoZnVsbCkgDQogICAgICAgICAgICAgIDM2MA0K ICAgICAgICBlbHNlIDE4MCkgKiBkZWcNCiAgICAgICAgd2hpbGUgKGsgPCBu ZXh0LmFuZ2xlKSB7DQogICAgICAgICAga2V5LnggPC0gYyhrZXkueCwgbGVu ICogY29zKGspICsga2V5LmxvY1sxXSkNCiAgICAgICAgICBrZXkueSA8LSBj KGtleS55LCBsZW4gKiBzaW4oaykgKyBrZXkubG9jWzJdKQ0KICAgICAgICAg IGsgPC0gayArIGRlZw0KICAgICAgICB9DQogICAgICAgIGtleS54IDwtIGMo a2V5LngsIGxlbiAqIGNvcyhuZXh0LmFuZ2xlKSArIGtleS5sb2NbMV0sIA0K ICAgICAgICAgICAgICAgICAgIE5BKQ0KICAgICAgICBrZXkueSA8LSBjKGtl eS55LCBsZW4gKiBzaW4obmV4dC5hbmdsZSkgKyBrZXkubG9jWzJdLCANCiAg ICAgICAgICAgICAgICAgICBOQSkNCiAgICAgIH0NCiAgICAgIHBhcihsd2Qg PSAwLjI1KQ0KICAgICAgcG9seWdvbihrZXkueCwga2V5LnksIGNvbCA9IHNl Zy5jb2xvcnMpDQogICAgICBwYXIobHdkID0gMSkNCiAgICB9DQogICAgZWxz ZSB7DQogICAgICB0ZW1wLnguY29vcmQgPC0gY29zKGFuZ2xlcykgKiBsZW4g KyBrZXkubG9jWzFdDQogICAgICB0ZW1wLnkuY29vcmQgPC0gc2luKGFuZ2xl cykgKiBsZW4gKyBrZXkubG9jWzJdDQogICAgICBwYXIobHdkID0gMC4yNSkN CiAgICAgIGlmIChyYWRpdXMpIA0KICAgICAgICBzZWdtZW50cyhyZXAoa2V5 LmxvY1sxXSwgbi5zZWcpLCByZXAoa2V5LmxvY1syXSwgDQogICAgICAgICAg ICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICBuLnNlZyksIHRl bXAueC5jb29yZCwgdGVtcC55LmNvb3JkKQ0KICAgICAgbGluZXMoYyh0ZW1w LnguY29vcmQsIHRlbXAueC5jb29yZFsxXSksIGModGVtcC55LmNvb3JkLCAN CiAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAg ICAgIHRlbXAueS5jb29yZFsxXSkpDQogICAgICBwYXIobHdkID0gMSkNCiAg ICB9DQogICAgaWYgKGlzLm51bGwoa2V5LmxhYmVscykpIA0KICAgICAga2V5 LmxhYmVscyA8LSBkaW1uYW1lcyh4KVtbMl1dDQogICAgbGFiLmFuZ2wgPC0g YW5nbGVzICsgaWYgKGRyYXcuc2VnbWVudHMpIA0KICAgICAgKGFuZ2xlc1sy XSAtIGFuZ2xlc1sxXSkvMg0KICAgIGVsc2UgMA0KICAgIGxhYmVsLnguY29v cmQgPC0gY29zKGxhYi5hbmdsKSAqIGxlbiAqIDEuMSArIGtleS5sb2NbMV0N CiAgICBsYWJlbC55LmNvb3JkIDwtIHNpbihsYWIuYW5nbCkgKiBsZW4gKiAx LjEgKyBrZXkubG9jWzJdDQogICAgZm9yIChrIGluIDE6bi5zZWcpIHsNCiAg ICAgIHRleHQuYWRqIDwtIGlmIChsYWIuYW5nbFtrXSA8ICg5MCAqIGRlZykg fHwgbGFiLmFuZ2xba10gPiANCiAgICAgICAgICAgICAgICAgICAgICAoMjcw ICogZGVnKSkgDQogICAgICAgIDANCiAgICAgIGVsc2UgaWYgKGxhYi5hbmds W2tdID4gKDkwICogZGVnKSAmJiBsYWIuYW5nbFtrXSA8IA0KICAgICAgICAg ICAgICAgKDI3MCAqIGRlZykpIA0KICAgICAgICAxDQogICAgICBlbHNlIDAu NQ0KICAgICAgaWYgKGxhYi5hbmdsW2tdIDw9ICg5MCAqIGRlZykpIA0KICAg ICAgICB0ZXh0LmFkaiA8LSBjKHRleHQuYWRqLCAwLjUgKiAoMSAtIGxhYi5h bmdsW2tdLyg5MCAqIA0KICAgICAgICAgICAgICAgICAgICAgICAgICAgICAg ICAgICAgICAgICAgICAgICAgICAgICAgICBkZWcpKSkNCiAgICAgIGVsc2Ug aWYgKGxhYi5hbmdsW2tdID4gKDkwICogZGVnKSAmIGxhYi5hbmdsW2tdIDw9 IA0KICAgICAgICAgICAgICAgKDI3MCAqIGRlZykpIA0KICAgICAgICB0ZXh0 LmFkaiA8LSBjKHRleHQuYWRqLCAobGFiLmFuZ2xba10gLSAoOTAgKiANCiAg ICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAg IGRlZykpLygxODAgKiBkZWcpKQ0KICAgICAgZWxzZSBpZiAobGFiLmFuZ2xb a10gPiAoMjcwICogZGVnKSkgDQogICAgICAgIHRleHQuYWRqIDwtIGModGV4 dC5hZGosIDEgLSAoMC41ICogKGxhYi5hbmdsW2tdIC0gDQogICAgICAgICAg ICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICgyNzAgKiBkZWcp KS8oOTAgKiBkZWcpKSkNCiAgICAgIHRleHQuZGVmYXVsdCh4ID0gbGFiZWwu eC5jb29yZFtrXSwgeSA9IGxhYmVsLnkuY29vcmRba10sIA0KICAgICAgICAg ICAgICAgICAgIGxhYmVscyA9IGtleS5sYWJlbHNba10sIGNleCA9IDAuNSwg YWRqID0gdGV4dC5hZGopDQogICAgfQ0KICB9DQogIGludmlzaWJsZSgpDQp9 DQo= --1465796871-645734172-974620647=:18687-- -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-devel mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html Send "info", "help", or "[un]subscribe" (in the "body", not the subject !) To: r-devel-request@stat.math.ethz.ch _._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._
Bug in stars function (PR#738)
2 messages · Tom Dye, Peter Dalgaard
Please don't use MIME attachments when reporting bugs. The resending mechanism renders them undecodable when seen on r-devel. tdye@lava.net writes:
This message is in MIME format. The first part should be readable text, while the remaining parts are likely unreadable without MIME-aware tools. Send mail to mime@docserver.cac.washington.edu for more info. --1465796871-645734172-974620647=:18687 Content-Type: TEXT/PLAIN; charset=US-ASCII Hi all, I've discovered and (hopefully) fixed a bug in the stars function I wrote some time ago. A corrected function is attached. Please let me know if you have any questions or comments. Tom Thomas S. Dye, Ph.D. http://www.tsdye.com Home: 812A 19th Avenue, Honolulu, Hawaii 96816. (808) 739-1367 or 387-9352. Work: International Archaeological Research Institute, Inc., 2081 Young St., Honolulu, Hawaii 96826. Voice (808) 946-2548; Fax 943-0716. --1465796871-645734172-974620647=:18687 Content-Type: TEXT/PLAIN; charset=US-ASCII; name="new-stars.r" Content-Transfer-Encoding: BASE64 Content-ID: <Pine.BSI.4.10.10011182157270.18687@malasada.lava.net> Content-Description: bug fix Content-Disposition: attachment; filename="new-stars.r" ZnVuY3Rpb24gKHgsIGZ1bGwgPSBUUlVFLCBzY2FsZSA9IFRSVUUsIHJhZGl1 cyA9IFRSVUUsIGxhYmVscyA9IGRpbW5hbWVzKHgpW1sxXV0sIA0KICAgICAg ICAgICAgbG9jYXRpb25zID0gTlVMTCwgeGxpbWl0ID0gTlVMTCwgeWxpbWl0 ID0gTlVMTCwgbGVuID0gMSwgDQogICAgICAgICAgICBjb2xvcnMgPSBOVUxM
O__ ---- Peter Dalgaard Blegdamsvej 3 c/ /'_ --- Dept. of Biostatistics 2200 Cph. N (*) \(*) -- University of Copenhagen Denmark Ph: (+45) 35327918 ~~~~~~~~~~ - (p.dalgaard@biostat.ku.dk) FAX: (+45) 35327907 -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-devel mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html Send "info", "help", or "[un]subscribe" (in the "body", not the subject !) To: r-devel-request@stat.math.ethz.ch _._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._