作者mangogogo (mangogo)
看板NCTU-STAT95G
标题[讨论] 程式码
时间Mon Nov 20 14:04:32 2006
我之前就一直有一个想法就是
大家可以一起分享使用过不错的code
降大家也不用再大费周章花同样的时间找
所以版主置个底吧
大家要分享就推文
有需要再来推文找找看
说不定就有自己要的
当然这需要大家的共享~~
ps.
R code package 可以去
http://www.r-project.org/的Search 中有一个google bar
打你要的关键字 即可得到介绍code的pdf
如果是你要的 去使用R下载 并安装即可
--
※ 发信站: 批踢踢实业坊(ptt.cc)
◆ From: 140.113.114.192
1F：→ mangogogo:par(mfrow=c(a,b)) 图档存成a*b图 11/20 14:04
2F：→ mangogogo:read.table('C:/', header =FALSE, sep ='') 读档 11/20 14:06
3F：→ mangogogo:hist(X) + lines(density(X)) 即可画出histo+density 11/20 14:07
4F：→ mangogogo:library("car") + qq.plot() 可以画出有边界的Q-Q plot갠 11/20 14:09
5F：→ mangogogo:抓"nortest"这个package 里面有很多常态性检定的code 11/20 14:11
6F：→ mangogogo:C=numeric(30) 给C提供size30的大小 11/20 14:13
7F：→ mangogogo:矩阵和矩阵相乘 用 "%*%" EX A%*%B ;反矩阵 solve() 11/20 14:14
8F：→ mangogogo:library("ellipse") 这个package可划出多维C.I.的椭圆 11/20 14:16
※ 编辑: mangogogo 来自: 140.113.114.192 (11/20 14:18)
9F：→ mangogogo:library("outlier") 这个package有可以找outlier的方法 11/20 14:19
10F：推 josephw:我每次都有分享喔 没有暗K喔!!! 11/20 14:22
11F：→ mangogogo:library("lawstat") 有随机性检定的code 11/20 14:21
12F：→ mangogogo:对呀 文廷都会跟我说 我都跟他学的 哈哈 11/20 14:23
13F：推 mangogogo:library("rgl") 可以画3D图 11/20 14:26
14F：推 mangogogo:存图到外面jpeg(file="c:/")+plot()+dev.off() 11/20 14:31
15F：推 josephw:屁蛋! 我都跟黄翁贤学的 11/20 15:54
16F：推 aflilfesy:我都跟许仲竹学的 11/20 16:38
17F：推 Oyin:hist(X,prob=TRUE) 画出Density 11/21 10:17
18F：推 Oyin:再加上 lines(density(X,bw=0.1)) 画出来比较好看 11/21 10:22
19F：推 mangogogo:hist(r,prob=TRUE,main = paste("Histogram of n=" , n) 11/21 22:13
20F：→ mangogogo:可以随着for自动改TITLE 11/21 22:13
21F：推 mangogogo:text(x label,y label,"....") 在图中座标(x,y)处加注解 11/27 23:17
22F：→ mangogogo:plot的里面加上一个code xlim=c(-1,1) 表图x范围由-1~1 11/27 23:19
23F：推 mangogogo:segments(x0, y0, x1, y1) 可以在图中加线段^^ 12/04 14:19
24F：推 mangogogo:BOX-COX T =>library("car")>box.cox.powers() 12/07 09:35
25F：推 mangogogo:package"rrcov"=> screeplot() <----画screeplot 12/14 10:00
26F：推 mangogogo:eigen() <-可算出eigenvalue&eigenvector 12/26 15:26
27F：推 mangogogo:factanal()<-算出factor loading~ 12/26 17:36
28F：→ mangogogo:apply(X,1,mean)<----1表算各row 2表算各col 12/29 20:04
29F：推 mangogogo:diag("vector")=> diagonal matrix 01/16 15:02
30F：→ mangogogo:round(,digits=n) 四舍五入到小数点第n位 01/16 15:02
31F：推 mangogogo:二维阵列(YA)~Z=matrix(0,a,b) 给一个(a*b)的二维阵列 01/17 12:30
32F：推 mangogogo:which(x, arr.ind = ) 01/20 10:20
33F：推 mangogogo:save(...,file= ".Rdata")+load(".Rdata")=>存成R档 01/26 12:28
35F：推 mangogogo:control charts library("qcc")=>qcc(data,type) 02/05 16:50
36F：推 mangogogo:help(par) => can find all options for plotting 02/06 09:26
37F：推 mangogogo:help.search("apply")=>可以找到apply系列的相关指令^^ 02/10 21:54
38F：推 mangogogo:一维积分integrate() 二维积分library("adapt")+adapt() 03/06 23:02
39F：推 shyfang:有一个很多data的package,KMsurv, 载入後用data()就可以抓 03/13 15:31
40F：推 josephw:Sys.time() 可以抓系统时间 03/28 00:36
41F：推 mangogogo:proc.time()<-most useful for "timing”the evaluation 03/28 01:36
42F：推 mangogogo:package "combinat"=>library("combinat") nCm(n,n) 04/04 09:12
43F：推 josephw:mass里 fraction可以把小数变成分数 04/13 20:55
44F：推 josephw:fractions 04/13 20:59
45F：推 josephw:x先任意给个阵列或位 data.entry(x) 会跳出视窗输入资料 04/13 21:05
46F：推 josephw:table(cut(y,pretty(y))) y是你的资料 他会将y做次数分配 04/13 21:14
47F：→ josephw:如果要像列联表那样两个 04/13 21:15
48F：→ josephw:table(cut(y,pretty(y)),cut(x,pretty(x))) x,y是你的data 04/13 21:16
49F：推 ilanglenn:matrix(scan("资料位址"),#row,#col) <=j大师教学 04/18 11:37
53F：→ Oyin:/R-programming.ppt#291,31,apply 07/22 12:49
54F：推 mangogogo:legend(X座,Y座,c('名称1','名称2'),lty=c(样1,样2)) 09/04 12:33
55F：推 mangogogo:mod => "%%" ex. 7mod3=1 => "7%%3" 10/31 14:21
56F：→ mangogogo:iris[sort(iris[,4],index.return = TRUE)$ix,] 12/12 17:02
57F：→ mangogogo:上面就是 EXCEL般的排序 (扩大范围) 12/12 17:02
※ 编辑: mangogogo 来自: 140.113.253.144 (01/16 01:45)
58F：→ mangogogo:求mle optimize(f,interval=c(0,1),maximum=T) 01/16 23:35
59F：→ mangogogo:找多围的最佳解用 optim 01/16 23:35
60F：→ mangogogo:bcs:packages:: RODBC/foreign/xlsReadWrite 皆能满足你 02/13 23:11
61F：推 aflilfesy:图存pdf档 pdf(file="ex.pdf"), width = 4, height = 4) 02/20 13:53
62F：→ aflilfesy:# [your program] 02/20 13:54
63F：→ aflilfesy:dev.off() 02/20 13:55
65F：推 ilanglenn:嗡嗡提供double integrate 04/18 16:48
67F：→ ilanglenn:ml 04/18 16:48
68F：推 ilanglenn:变数回圈命名: 10/05 12:35
69F：→ ilanglenn:X=matrix(rnorm(15),ncol=3) 10/05 12:35
70F：→ ilanglenn:for(i in 1:5) assign(paste("b",i, sep=""), X[i,]) 10/05 12:36