2014-08-01から1ヶ月間の記事一覧
##dplyr##magrittr http://d.hatena.ne.jp/teramonagi/20140814/1408024760http://d.hatena.ne.jp/teramonagi/20140403/1396529027http://rpubs.com/teramonagi/1112http://meme.biology.tohoku.ac.jp/students/iwasaki/rstats/tidyr.htmlhttp://meme.biolog…
「第5回デスマコロシアム」問題のトーナメント結果発表です！──優勝者は…！ #デスマコロシアム｜CodeIQ MAGAZINE 標準出力に、余計な記号やスペースなどなしに 12fizz4buzzfizz78fizzbuzz11fizz1314FIZZBUZZ1617fizz19buzzfizz2223fizzbuzz26fizz2829FIZZBUZ…
単なる覚書 ####比率の差の信頼区間 多分Miettinen and Nurminenlibrary(PropCIs)diffscoreci(10, 100, 20, 150, 0.95) ###正確な比率の差の信頼区間 時間がかかるので注意#library(ExactCIdiff)#BinomCI(100,150,10,20,0.95,CItype = "Two.sided", precisio…
厚生労働省の所にデータが有ったので、やり直す。 ##男女別生命表#データフレーム作成 厚生労働省H22国勢調査による完全生命表dat <- data.frame(numeric( 226 ))dat $ x = c( 0,1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20,21,22,23,24,25,26,27,2…
mean survial time; AUC of survival curve r-statstics-fan Tuesday, August 12, 2014 1年生存率と5年生存率 - 驚異のアニヲタ社会復帰への道 に触発されて、５年生存率ではなく平均寿命やMSTを計算したくなった。 前回の記事を使って、簡単にコピペでデー…
1年生存率と5年生存率 - 驚異のアニヲタ社会復帰への道 に触発されて、平均寿命とか計算したくなった。 平均寿命といえば 生存曲線下面積RMST(Restricted mean survival time)～比例ハザード性が仮定できない生存曲線での代替手法 - r-statistics-fanの日記 …
多倍長精度計算パッケージ Rmpfr - 裏 RjpWiki #追記2014.9.20 よく見ると自分の関数は小数点以下も出力されていて読みにくい。 裏 RjpWikiさんの所で、より良い方法の記事が出ました。 <a href="http://blog.goo.ne.jp/r-de-r/e/0dee8bc64bfb50b39e796e207d50969e" data-mce-href="http://blog.goo.ne.jp/r-de-r/e/0dee8bc64bfb50b39e796e207d50969e">多倍長精度計算結果を見やすく出力 - 裏 RjpWiki</a> 多倍長精度計算…
□＝１～９（重複しない） ちょ、難しいんだけど。 唯一解なのか分からない。総当りで確認する。 library(e1071)f1 <- function(){s <- permutations(9)ans <- s[,1] / (10*s[,2] + s[,3]) + s[,4] / (10*s[,5] + s[,6]) + s[,7] / (10*s[,8] + s[,9])return(…
生存曲線下面積RMST(Restricted mean survival time)というのを聞いた。 論文の多くは田舎病院では入手できなかったが、下記は読めた。 Royston, P. & Parmar, M.K., 2013. Restricted mean survival time: an alternative to the hazard ratio for the desi…
