統計分析

上回這篇Stata: (Kuso) 性伴侶人數研究 帶來不少流量，看來大家都對這檔事滿有興趣的…不過沒人留言要看結果，再加上我快被報告逼瘋了，於是就一拖再拖 @@ Anyway, 跑regression很簡單，接著上次的code，我跑的程式如下，有興趣的自己可以replicate結果：

regress dpartners dmale age dblack dother dincome if (dwidowed!=.)…

人在寫報告的時候，總是會有一堆奇思幻想。關於報告本身的聯想與進展並不多，想的都是一些有的沒的。小的在下為了推廣Stata，並為了娛樂大家，不惜犧牲形象，想出勁爆話題，讓大家學習統計時，印象更加深刻。請不要以有色眼光看待此kuso研究和想出這研究的人…

所用的資料庫仍是Stata: 敘述統計(descriptive statistics)所用的GSS 2000。一開始資料要做很多recode和產生很多虛擬變項(dummy variable)，這方面不熟的就從Stata: 敘述統計(descriptive statistics) 看起。變數後面括號的英文為該變數在GSS 2000裡的variable name。

我們想知道的是：有什麼因素可以預測性伴侶人數(partners)？很直覺的就會想到幾個因素：種族(race)、年紀(age)、收入(income)、性別(sex)、婚姻狀況(marital)。

如果是這樣，就不叫惡搞了 XD 所以小弟想了兩個假設：

當控制了上述變數之後：
1.…

最近發文不是很勤，一方面是在忙著作自己的研究，一方面期末又到了，要開始擠報告了 🙁 。另外還有一件正在忙的事就是：開發程式！開發什麼程式？就像文章標題所說的，是跟表格有關的程式。之前不是介紹過Stata: 輸出regression table到Word或Excel 嗎？秉持著懶還可以更懶的原則，我就突發奇想：既然表格都可以直接匯出，那分析也一併好了！於是乎：ReadTable(名稱暫訂)就這樣產生了。

先來看看畫面。輸入視窗如下：

當然不會要你輸入的。你可以從Stata或Excel將資料貼到此程式。

資料直接貼上。

 

 

 

然後按下分析：

目前只作了OLS regression，未來還會加logistic regression或其它的。至於中文的，嗯啊…，我自己翻譯怎麼都不順。有人要提供翻譯的話非常歡迎。英文到時候我還會再修過，現在權當測試用。

載點在哪？就說給你聞香而已，沒有要給你喝 XD 目前還有一些小問題，比如說coefficient還沒設定要四捨五入到哪一位，p的值會不小心最後多個0，表格目前還沒有column name……。另外，雖然是用java開發的，但因為我是用jdk 1.6，mac的相容性還不足，還沒搞定向下相容的問題。多語界面也還在想，還沒動手做。

第一次用java寫程式，所以速度比較慢。大夥就慢慢等哪時候我有空，完成第一個版本吧！

…

上次說到：仔細看下圖，問題有兩個：

一個是問卷的觀察數(obs)不同，一個是race、sex這些必須轉換成虛擬變量(dummy variable)才有意義。後者上次已解決了，前者則是要利用if指令來處理missing values。

在Stata裡，missing value通常都是以 .…

學統計程式，有個dataset可以玩會比較能夠進入狀況，不然指令看過就忘了。此外，我覺得學習統計程式也是了解統計觀念的一個好方法，因為你可以直接看到結果。當結果與common sense不合，就會再去找找是不是哪做錯了。

找了很久，終於找到一個不大不小剛剛好，可以拿來練手的dataset，這dataset是general social survey 2000年所做的，你可以到此下載(記得選2000年)，或是直接用我上傳到借來的空間。Gss2000的大小為3.61mb，原汁原味，沒經過任何修改。由於gss的網站好像沒看到可否自由散佈的相關資訊，就假設可以好了 XD 真的不行的話，那到時候再找別的dataset。

use http://twtcsl.org/dataset/gss2000.dta

…