[討論] 程式碼

[mangogogo (mangogo)]

我之前就一直有一個想法就是

大家可以一起分享使用過不錯的code

降大家也不用再大費周章花同樣的時間找

所以版主置個底吧

大家要分享就推文

有需要再來推文找找看

說不定就有自己要的

當然這需要大家的共享~~

ps.

R code package 可以去http://www.r-project.org/的Search 中有一個google bar

打你要的關鍵字 即可得到介紹code的pdf

如果是你要的 去使用R下載 並安裝即可

--
※ 發信站: 批踢踢實業坊(ptt.cc)
◆ From: 140.113.114.192
→ mangogogo:par(mfrow=c(a,b)) 圖檔存成a*b圖
→ mangogogo:read.table('C:/', header =FALSE, sep ='')  讀檔
→ mangogogo:hist(X)  + lines(density(X)) 即可畫出histo+density
→ mangogogo:library("car") + qq.plot() 可以畫出有邊界的Q-Q plot갠
→ mangogogo:抓"nortest"這個package 裡面有很多常態性檢定的code
→ mangogogo:C=numeric(30)   給C提供size30的大小
→ mangogogo:矩陣和矩陣相乘 用 "%*%" EX A%*%B ;反矩陣 solve()
→ mangogogo:library("ellipse") 這個package可劃出多維C.I.的橢圓
※ 編輯: mangogogo       來自: 140.113.114.192      (11/20 14:18)
→ mangogogo:library("outlier") 這個package有可以找outlier的方法
推 josephw:我每次都有分享喔  沒有暗K喔!!!
→ mangogogo:library("lawstat") 有隨機性檢定的code
→ mangogogo:對呀 文廷都會跟我說 我都跟他學的 哈哈
推 mangogogo:library("rgl") 可以畫3D圖
推 mangogogo:存圖到外面jpeg(file="c:/")+plot()+dev.off()
推 josephw:屁蛋! 我都跟黃翁賢學的
推 aflilfesy:我都跟許仲竹學的
推 Oyin:hist(X,prob=TRUE)  畫出Density
推 Oyin:再加上 lines(density(X,bw=0.1)) 畫出來比較好看
推 mangogogo:hist(r,prob=TRUE,main = paste("Histogram of n=" , n)
→ mangogogo:可以隨著for自動改TITLE
推 mangogogo:text(x label,y label,"....") 在圖中座標(x,y)處加註解
→ mangogogo:plot的裡面加上一個code xlim=c(-1,1) 表圖x範圍由-1~1
推 mangogogo:segments(x0, y0, x1, y1) 可以在圖中加線段^^
推 mangogogo:BOX-COX T =>library("car")>box.cox.powers()
推 mangogogo:package"rrcov"=> screeplot() <----畫screeplot
推 mangogogo:eigen() <-可算出eigenvalue&eigenvector
推 mangogogo:factanal()<-算出factor loading~
→ mangogogo:apply(X,1,mean)<----1表算各row  2表算各col
推 mangogogo:diag("vector")=>  diagonal matrix
→ mangogogo:round(,digits=n) 四捨五入到小數點第n位
推 mangogogo:二維陣列(YA)~Z=matrix(0,a,b)  給一個(a*b)的二維陣列
推 mangogogo:which(x, arr.ind = )
推 mangogogo:save(...,file= ".Rdata")+load(".Rdata")=>存成R檔
推 mangogogo:Tinn-R:http://0rz.tw/e72kD
推 mangogogo:control charts library("qcc")=>qcc(data,type)
推 mangogogo:help(par) => can find all options for plotting
推 mangogogo:help.search("apply")=>可以找到apply系列的相關指令^^
推 mangogogo:一維積分integrate() 二維積分library("adapt")+adapt()
推 shyfang:有一個很多data的package,KMsurv, 載入後用data()就可以抓
推 josephw:Sys.time()   可以抓系統時間
推 mangogogo:proc.time()<-most useful for "timing”the evaluation
推 mangogogo:package "combinat"=>library("combinat") nCm(n,n)
推 josephw:mass裡  fraction可以把小數變成分數
推 josephw:fractions
推 josephw:x先任意給個陣列或位  data.entry(x)  會跳出視窗輸入資料
推 josephw:table(cut(y,pretty(y))) y是你的資料  他會將y做次數分配
→ josephw:如果要像列聯表那樣兩個
→ josephw:table(cut(y,pretty(y)),cut(x,pretty(x))) x,y是你的data
推 ilanglenn:matrix(scan("資料位址"),#row,#col)  <=j大師教學
推 mangogogo:http://www.math.unm.edu/~bedrick/PIBBS/Rsurv.pdf
推 Y0SHIKI:http://0rz.tw/142F4   Tinn-R_1.19.2.3
推 Oyin:http://www.math.ntu.edu.tw/~hchen/Prediction/notes/R-prog
→ Oyin:/R-programming.ppt#291,31,apply
推 mangogogo:legend(X座,Y座,c('名稱1','名稱2'),lty=c(樣1,樣2))
推 mangogogo:mod => "%%"   ex. 7mod3=1  => "7%%3"
→ mangogogo:iris[sort(iris[,4],index.return = TRUE)$ix,]
→ mangogogo:上面就是 EXCEL般的排序 (擴大範圍)
※ 編輯: mangogogo       來自: 140.113.253.144      (01/16 01:45)
→ mangogogo:求mle optimize(f,interval=c(0,1),maximum=T)
→ mangogogo:找多圍的最佳解用  optim
→ mangogogo:bcs:packages:: RODBC/foreign/xlsReadWrite 皆能滿足你
推 aflilfesy:圖存pdf檔 pdf(file="ex.pdf"), width = 4, height = 4)
→ aflilfesy:# [your program]
→ aflilfesy:dev.off()
→ mangogogo:http://cged.genes.nig.ac.jp/RGM2/index.php?clear=all
推 ilanglenn:嗡嗡提供double integrate
→ ilanglenn:http://tolstoy.newcastle.edu.au/R/help/04/10/5951.ht
→ ilanglenn:ml
推 ilanglenn:變數迴圈命名:
→ ilanglenn:X=matrix(rnorm(15),ncol=3)
→ ilanglenn:for(i in 1:5) assign(paste("b",i, sep=""), X[i,])