【正文】 2) 2( 2) 3( 3) 3( 3) 17( 4) 17( 4) 60( 5) 60( 5) Frequency Table Percents Percents Value Count Cell Cum Value Count Cell Cum 1 1 17 1 2 1 60 1 3 1 原始的 5個(gè)數(shù)目清楚 , 還是這 50多個(gè)數(shù)目清楚 ? 這些輸出都有意義嗎 ? 實(shí)際上 … ?對(duì)于一個(gè)統(tǒng)計(jì)數(shù)據(jù)，如果選擇“傻瓜”軟件適合該數(shù)據(jù)的選項(xiàng)的全部輸出，那么 … ?其中可能包含（很多，甚至多數(shù)）沒(méi)有意義的結(jié)果 （比如SPSS中的對(duì)數(shù)線(xiàn)性模型、非參數(shù)檢驗(yàn)中的多重相關(guān)樣本的檢驗(yàn)等等，各種回歸中的適合正態(tài)樣本的檢驗(yàn)） 而且 … ?這些眾多的輸出，給了那些諸如檢驗(yàn)中的“少數(shù)服從多數(shù)”的“創(chuàng)新者”以似是而非的“論據(jù)”。 proc univariate freq normal。 cards。 data test。 ?因此，就應(yīng)該特別留神，明白自己是在干什么。 ?這就使得使用者大傷腦筋。 ?可能統(tǒng)計(jì)學(xué)家缺乏商業(yè)頭腦。 ?但統(tǒng)計(jì)軟件則不那么負(fù)責(zé)。 統(tǒng)計(jì)軟件的問(wèn)題 ?諸如法律和醫(yī)學(xué)的軟件都有不少警告，不時(shí)提醒你去咨詢(xún)專(zhuān)家。 統(tǒng)計(jì)軟件 ?是否傻瓜式的統(tǒng)計(jì)軟件使用可以代替統(tǒng)計(jì)課程了？ ?當(dāng)然不是。 統(tǒng)計(jì)軟件 ?統(tǒng)計(jì)軟件的發(fā)展，也使得統(tǒng)計(jì)從統(tǒng)計(jì)學(xué)家的圈內(nèi)游戲變成了大眾的游戲。 ?計(jì)算機(jī)的使用 ， 也從過(guò)去必須學(xué)會(huì)計(jì)算機(jī)語(yǔ)言到只需要 “ 傻瓜式 ”地點(diǎn)擊鼠標(biāo) 。 ?而現(xiàn)在統(tǒng)計(jì)仍然是進(jìn)行數(shù)字計(jì)算最多的用戶(hù) 。 ?最初的計(jì)算機(jī) 僅僅是為科學(xué)計(jì)算而建造的 。統(tǒng)計(jì)軟件和R語(yǔ)言 裝了 R沒(méi)有？ 一個(gè)廣泛接受的統(tǒng)計(jì)定義為： ?統(tǒng)計(jì) 是用以收集 數(shù)據(jù) 、分析 數(shù)據(jù) 和由 數(shù)據(jù) 得出結(jié)論的一組概念、原則和方法 . 這個(gè)定義決定了統(tǒng)計(jì)的命運(yùn)： ?和數(shù)學(xué)及音樂(lè)不同 , 統(tǒng)計(jì)不能欣賞自己 , 它不為實(shí)際服務(wù)就 沒(méi)有存在必要 ?統(tǒng)計(jì)必須 為各個(gè)領(lǐng)域服務(wù) ?統(tǒng)計(jì)必須和 數(shù)據(jù) 打交道 ?因此 ,統(tǒng)計(jì)必須和 計(jì)算機(jī) 結(jié)合 統(tǒng)計(jì)和計(jì)算機(jī) ?現(xiàn)代生活已離不開(kāi)計(jì)算機(jī)了 。 ?但最早使用計(jì)算機(jī)的是統(tǒng)計(jì) 。 大型計(jì)算機(jī)的最早一批用戶(hù)就包含統(tǒng)計(jì) 。 統(tǒng)計(jì)和計(jì)算機(jī) ?計(jì)算機(jī)現(xiàn)在早已脫離了僅有計(jì)算功能的單一模式 ， 而成為 百姓生活的一部分 。 結(jié)果也從單純的數(shù)字輸出到包括漂亮的表格和圖形的各種形式 。 ?只要輸入數(shù)據(jù)，點(diǎn)幾下鼠標(biāo)，做一些選項(xiàng)，馬上就得到令人驚嘆的 漂亮結(jié)果了 。 ?數(shù)據(jù)的整理和識(shí)別，方法的選用，計(jì)算機(jī)輸出結(jié)果的理解都不象使用傻瓜相機(jī)那樣簡(jiǎn)單可靠。 ?這是注意飯碗的律師和大夫的高明之處。只要數(shù)據(jù)格式無(wú)誤、方法不矛盾而且不用零作為除數(shù)就一定給你結(jié)果，而且 沒(méi)有任何警告 。 統(tǒng)計(jì)軟件的問(wèn)題 ?另外，統(tǒng)計(jì)軟件 輸出的結(jié)果太多 ； ?即使是同樣的方法，不同軟件輸出的內(nèi)容還不一樣；有時(shí)同樣的內(nèi)容名稱(chēng)也不一樣。即使統(tǒng)計(jì)學(xué)家也不一定能解釋所有的輸出。 不要在得到一堆毫無(wú)意義的垃圾之后還沾沾自喜。 input x。 1 2 3 17 60 run。 run。 ?許多“傻瓜”軟件的輸出找不到確切數(shù)學(xué)背景。 統(tǒng)計(jì)軟件的種類(lèi) ?統(tǒng)計(jì)軟件的種類(lèi)很多。 ?還有些是專(zhuān)門(mén)的軟件，只處理某一類(lèi)統(tǒng)計(jì)問(wèn)題。 統(tǒng)計(jì)軟件的種類(lèi) ?SPSS： 這是一個(gè)很受歡迎的統(tǒng)計(jì)軟件； ?它容易操作 ， ?輸出漂亮 ， ?功能齊全 ， ?價(jià)格合理 。 統(tǒng)計(jì)軟件的種類(lèi) ? Excel： 嚴(yán)格說(shuō)來(lái)不是統(tǒng)計(jì)軟件 ，但作為數(shù)據(jù)表格軟件 ， 必然有一定統(tǒng)計(jì)計(jì)算功能 。 ? 對(duì)于簡(jiǎn)單分析 ， Excel還算方便 ，但隨著問(wèn)題的深入 ， 就不那么 “ 傻瓜 ” ， 需要很麻煩地使用函數(shù) ， 甚至根本沒(méi)有相應(yīng)的方法了 。 ?價(jià)格不菲 ,每年交費(fèi) .即使贈(zèng)送 ,條件苛刻 ? SAS公司傲慢無(wú)禮 ?盡管現(xiàn)在已經(jīng)盡量 “ 傻瓜化 ” ， 仍然需要一定的訓(xùn)練才可以進(jìn)入 。 統(tǒng)計(jì)軟件的種類(lèi) ? Splus：這是 統(tǒng)計(jì)學(xué)家 喜愛(ài)的軟件 。 也可以自己加入算法 . ? 強(qiáng)大而又方便的編程功能 ， 使得研究人員可以編制自己的程序來(lái)實(shí)現(xiàn)自己的理論和方法 。 但主要以其 方便的編程 為顧客所青睞 。 ? 但是對(duì)于不會(huì)編程者 ,不那么 “ 傻瓜 ” 統(tǒng)計(jì)軟件的種類(lèi) ? R軟件： 免費(fèi)的 ,志愿者管理的軟件 。 也可以自己加入自己算法的軟件包 . ? 從網(wǎng)上可以不斷更新和增加有關(guān)的軟件包 、 數(shù)據(jù)和程序 。 ? 可以通過(guò)網(wǎng)上服務(wù)把程序變成機(jī)器語(yǔ)言 ? 這是 發(fā)展最快 的軟件 ， 受到世界上統(tǒng)計(jì)師生的歡迎 。 ? 對(duì)于一般非統(tǒng)計(jì)工作者來(lái)說(shuō) ， 主要問(wèn)題是它沒(méi)有“ 傻瓜化 ” 。 ? Statistica： 也是功能強(qiáng)大而齊全的 “ 傻瓜化 ” 的軟件 ， 在我國(guó)用的也不如 SAS與 SPSS那么普遍 。 ? Systat: 已經(jīng)被 SPSS購(gòu)買(mǎi)而消失 ? GAUSS： 這是一個(gè)很好用的統(tǒng)計(jì)軟件 ，許多搞經(jīng)濟(jì)的喜歡它 。 目前在我國(guó)使用的人不多 。 編程類(lèi)似于 S和 R。 ? FORTRAN： 這是應(yīng)用于各個(gè)領(lǐng)域的歷史很長(zhǎng)的非常優(yōu)秀的編程軟件 ， 功能強(qiáng)大 ，有大量數(shù)學(xué)軟件包及一定的統(tǒng)計(jì)軟件包 。 但需要編程和編譯 。 ? C/C++: 是許多軟件 ,諸如 SAS的基礎(chǔ) 。只能輸出軟件規(guī)定的輸出 . ?無(wú)法在中間 插入任何算法 . ?無(wú)法實(shí)現(xiàn)軟件所 沒(méi)有的 計(jì)算 .無(wú)法實(shí)現(xiàn)任何方法或計(jì)算方面的創(chuàng)新 . ?是輸入輸出皆有限制的 黑盒子 . ?用語(yǔ)句的任何計(jì)算 (即使 1+1=2) 都需類(lèi)似八股文的 “ 花架子 ” . 毫無(wú)例外，沒(méi)有任何新的重要統(tǒng)計(jì)方法是利用諸如 SAS或 SPSS等傻瓜軟件產(chǎn)生的。 使用什么軟件編程呢？ ? SAS、 SPSS、 Minitab等也有程序，但很難學(xué)，不好記，也沒(méi)有任何普遍性和規(guī)律性。容易舉一反三。特別是 Splus軟件、 R軟件和Gauss軟件。T實(shí)驗(yàn)室開(kāi)發(fā) . ? R 項(xiàng)目由 Auckland 大學(xué)統(tǒng)計(jì)系的 Robert Gentleman和 Ross Ihaka于 1995年開(kāi)始的 . ?它很快得到廣泛用戶(hù)的歡迎 . 目前它是由 R核心發(fā)展團(tuán)隊(duì)維持 。 Jones (1995) ? lattice Lattice Graphics Interface ? tools Tools for Package Development ? utils The R Utils Package Packages （繼續(xù)） ? MASSMain Package of Venables and Ripley39。 ?希望你是下一個(gè)加盟這些packages的作者之一。, |, ! 還有 ? 輸入輸出數(shù)據(jù) : scan, , dump, save, load, write, ? letters, LETTERS ? list, matrix, array, cbind, rbind, merge ? sort, order, , rev, stack, unstack , reshape 序列和向量 ? z=seq(1,10,length=100) ? z=seq(1,10, len=100) ? z=seq(10,1,1) ? z=10:1 ? x=rep(3,1:3) ? x=rep(3:5,1:3) ? x ? [1] 3 4 4 5 5 5 ? x=rep(c(1,10),c(4,5)) ? w=c(1,3,x,z)。dnorm(,2,1)。rnorm(10,0,1)rnorm(10) ? t分布 :pt(,1)。qt(.7,1)。z=1:3。names(x)=LETTERS[1:3] ? x ? A B C ? 1 2 3 ? append(x,runif(3),after=2) ? A B C ? 向量賦值 ? z=1:5 ? z[7]=8。 z ? [1] 1 NA 2 NA 3 ? rnorm(10)[c(2,5)] ? z[c(1,3)] 去掉第 3元素 . ? z[(length(z)4):length(z)] 最后五個(gè)元素 . 向量的大小次序 ? z=sample(1:100,10)。x [,1] [,2] [,3] [,4] [,5] [1,] 1 5 9 13 17 [2,] 2 6 10 14 18 [3,] 3 7 11 15 19 [4,] 4 8 12 16 20 ? x=matrix(1:20,4,5,byrow=T)。 ncol(x)。x[,3]=1,1]。amp。all(x!=0)。(1:10)[x0] ? x=sample(1:7,5,rep=T)。x [,1] [,2] [,3] [1,] [2,] [3,] ? t(x) [,1] [,2] [,3] [1,] [2,] [3,] ? solve(x) solve(a,b)可以解 ax=b方程 [,1] [,2] [,3] [1,] [2,] [3,] 警告 :計(jì)算機(jī)中的 0是什么 ? ? x%*%solve(x) [,1] [,2] [,3] [1,] +00 [2,] +00 [3,] +00 ?要用線(xiàn)性代數(shù)的知識(shí)來(lái)判斷諸如有多少非零特征根等問(wèn)題 .假定 v是特征 值 組 成 的 向 量 , 不 能 用 諸 如sum(v!=0) 等方法來(lái)判斷非零特征根的數(shù)目 ! Matrix amp。 x [,1] [,2] [,3] [,4] [,5] [1,] [2,] [3,] [4,] ? (x) ? [1] TRUE ? x[1,2] ? x[1,] ? x[,2] ? dim(x)得到維數(shù) (4,5) Array ? x=array(runif(24),c(4,3,2)) ? (x) 可由 dim(x)得到維數(shù) (4,3,2) ? [1] FALSE ? x , , 1 [,1] [,2] [,3] [1,] [2,] [3,] [4,] , , 2 [,1] [,2] [,3] [1,] [2,] [3,] [4,] Array的子集 ? x=array(1:24,c(4,3,2)) ? x[c(1,3),] , , 1 [,1] [,2] [,3] [1,] 1 5 9 [2,] 3 7 11 , , 2 [,1] [,2] [,3] [1,] 13 17 21 [2,] 15 19 23 矩陣乘法及行列運(yùn)算 ? x=matrix(1:30,5,6)。z[[2]]。z$T$a