filename mydata 'se97prct.txt'; data da1; infile mydata; input #1 AREA 1-2 IDP 3-4 KESIGETU 6 KESINITI 7-8 Q1SEX 9 Q1NEN 10-11 Q1GETU 12-13 Q2 15 Q3 16 Q4 17 Q5AIMP 18 Q5BIMP 19 Q5CIMP 20 Q5DIMP 21 Q5EIMP 22 Q5FIMP 23 Q5GIMP 24 Q5ASTS 25 Q5BSTS 26 Q5CSTS 27 Q5DSTS 28 Q5ESTS 29 Q5FSTS 30 Q5GSTS 31 Q6A 33 Q6B 34 Q6C 35 Q6D 36 Q6E 37 Q6F 38 Q6G 39 Q6H 40 Q6I 41 Q6J 42 Q6K 43 Q6L 44 Q6M 45 Q6N 46 Q7A 47 Q7B 48 Q7C 49 Q8A 51 Q8B 52 Q8C 53 Q8D 54 Q8E 55 Q8F 56 Q8G 57 Q9A 58 Q9B 59 Q9C 60 Q9D 61 Q10 62 Q11 63-64 #2 Q12 6 Q14S01 11 Q14S02 12 Q14S03 13 Q14S04 14 Q14S05 15 Q14S06 16 Q14S07 17 Q14S08 18 Q14S09 19 Q14S10 20 Q14S11 21 Q14S12 22 Q14S13 23 Q14S14 24 Q14S15 25 Q14S16 26 Q15 27-28 Q16 30 Q17ASEX 31 Q17AAGE 32-33 Q17AEDU 34 Q17AEMPL 35 Q17ACOMP 36 Q17AJOB 37 Q17BSEX 38 Q17BAGE 39-40 Q17BEDU 41 Q17BEMPL 42 Q17BCOMP 43 Q17BJOB 44 Q17CSEX 45 Q17CAGE 46-47 Q17CEDU 48 Q17CEMPL 49 Q17CCOMP 50 Q17CJOB 51 Q17DSEX 52 Q17DAGE 53-54 Q17DEDU 55 Q17DEMPL 56 Q17DCOMP 57 Q17DJOB 58 Q18AB 59 Q18AC 60 Q18AD 61 Q18BC 62 Q18BD 63 Q18CD 64 Q19A 66 Q19BS1 67 Q19BS2 68 Q19CS1 69 Q19CS2 70 Q19DS1 71 Q19DS2 72 OPINION 73; /***** SYAKAI ISIKI NI KANSURU SENDAI SIMIN TYOUSA BY MURASE *****/ /***** BUNSEKI YOU SAS PROGRAM Version1.0 19980114 *****/ /***** *****/ /***** 「社会意識に関する仙台市民調査」 *****/ /***** 調査主体   東北大学文学部行動科学研究室(担当 村瀬洋一)*****/ /***** 調査実施費用 1996年度文部省科学研究費補助金 *****/ /***** (特別研究員奨励費:研究代表者 村瀬洋一) *****/ /***** *****/ /***** 母集団    仙台市青葉区と泉区の20歳以上70歳未満の男女 *****/ /***** サンプル数  1557人(農村部補充用サンプル57人を含む) *****/ /***** 有効回収数  1093人 回収率70.2% *****/ /***** 回収期間   1997年03月01日から04月30日 *****/ /***** 回収法    郵送調査法 *****/ /***** 標本抽出法  無作為抽出法(確率比例2段抽出法) *****/ /***** 1行目は、データファイルとしてSE97PRCT.TXT *****/ /*****              を使うことを意味している *****/ /*****  このデータは練習用のもので本物のデータではない *****/ /***** ここ以下が欠損値処理            *****/ /*****  無回答は欠損値となっている。       *****/ /***** SASは欠損値をピリオドで表す。      *****/ /***** データ中には、欠損値として9または99が入っているので *****/ /*****  これらをピリオドに置き換える       *****/ ARRAY MISS{65} Q2--Q9D Q14S01--Q14S16 Q16 Q19A Q19BS1 Q19BS2 ; DO J=1 TO 65; IF MISS{J}=9 THEN MISS{J}=.; END; ARRAY MI2{2} Q11 Q15 ; DO K=1 TO 2; IF MI2{K}=99 THEN MI2{K}=.; IF MI2{K}=98 THEN MI2{K}=.; END; ARRAY MI3{4} Q10 Q12 Q19DS1 Q19DS2 ; DO L=1 TO 4; IF MI3{L}=0 THEN MI3{L}=.; IF MI3{L}=0 THEN MI3{L}=.; END; /***** *****/ /***** データの加工 カテゴリー合併や新変数作成 *****/ /***** *****/ /***** 新変数として N3 を設定し、           *****/ /***** Q3の 4段階の値を、2カテゴリーに合併        *****/ N3=0; IF (Q3=1 OR Q3=2) THEN N3=1; IF (Q3=3 OR Q3=4) THEN N3=2; IF (Q2=1 and Q3=1) THEN N2 =1; else N2 =0; /***** NENREI TO NENREI KATEGORI SAKUSEI *****/ AGE =71-Q1NEN; IF Q1GETU<5 THEN AGE=AGE+1; /***** 新変数として NENDAI を設定し、                *****/ /***** 20〜70までの51段階の年齢を、2から6までの5カテゴリーに合併 *****/ select; when (19 < age < 30) nendai=2; when (29 < age < 40) nendai=3; when (39 < age < 50) nendai=4; when (49 < age < 60) nendai=5; otherwise nendai=6; end; /***** GAKUREKI WO KYOUIKU NENSUU HE HENKAN *****/ if q16=1 then edu=6; if q16=2 then edu=9; if q16=3 then edu=12; if q16=4 then edu=15; if q16=5 then edu=16; if q16=6 then edu=18; if q16=7 then edu=.; if q16=9 then edu=.; /***** GAKUREKI KATEGORI SAKUSEI *****/ if q16=1 or q16=2 then educt=1; if q16=3 or q16=4 then educt=2; if q16=5 or q16=6 then educt=3; if q16=7 or q16=9 then educt=.; if q16=1 or q16=2 then educt2=1; if q16=3 or q16=4 then educt2=2; if q16=5 then educt2=3; if q16=6 then educt2=4; if q16=7 or q16=9 then educt2=.; /***** *****/ /***** 新変数の作成 *****/ /***** *****/ /***** 問6の値1234を4321に変換 *****/ OTOKOSO=5-Q6A; YARITA =5-Q6B; DENTO =5-Q6C; SEX = Q1SEX-1; /***** 問15の世帯年収カテゴリーを年収何万円の単位に変換 *****/ FAMINCOM = Q15; IF Q15 = 1 THEN FAMINCOM = 50; IF Q15 = 2 THEN FAMINCOM = 150; IF Q15 = 3 THEN FAMINCOM = 250; IF Q15 = 4 THEN FAMINCOM = 350; IF Q15 = 5 THEN FAMINCOM = 450; IF Q15 = 6 THEN FAMINCOM = 550; IF Q15 = 7 THEN FAMINCOM = 650; IF Q15 = 8 THEN FAMINCOM = 750; IF Q15 = 9 THEN FAMINCOM = 850; IF Q15 =10 THEN FAMINCOM = 950; IF Q15 =11 THEN FAMINCOM = 1100; IF Q15 =12 THEN FAMINCOM = 1300; IF Q15 =13 THEN FAMINCOM = 1500; IF Q15 =14 THEN FAMINCOM = 1700; /***** 財産保有数の変数を作成 *****/ property = Q14S01 + Q14S02 + Q14S03 + Q14S04 + Q14S05 + Q14S06 + Q14S07 + Q14S08 + Q14S09 + Q14S10 + Q14S11 + Q14S12 + Q14S13 + Q14S14 + Q14S15; /***** 以下は行幅のオプション どこに書いてもよい  *****/ options ls=78; /***** ★★★ この下が proc step  ★★★        *****/ /***** sas プログラムは data step と proc step に分かれる *****/ /***** 以下が分析のための procステップの例 *****/ /***** *****/ /***** procステップはいくつあっても良い   *****/ /***** はじめの3つは度数分布表(単純集計とクロス集計)*****/ /***** 4つめは基本統計量を出力 *****/ /***** *****/ /***** プログラムの最後に必ずrun;をつける(実行せよ、の意味)*****/ /***** 単純集計   *****/ PROC FREQ; tables Q2--Q4 Q7A Q7B*Q7C; /***** クロス集計  *****/ PROC FREQ; tables Q1SEX*Q3 /nopercent nocol; PROC FREQ; tables Q1SEX*opinion /nopercent nocol; /***** 基本統計量 *****/ proc means; var q6a--q6c; by Q1SEX; /***** 以下は回帰分析のプログラム *****/ /*****   stbオプションを付けると、標準偏回帰係数を出力する *****/ proc reg data=da1; model q3 = q6a q6b q6c/stb; /***** 以下は相関係数のプログラム。 *****/ /*****   回帰分析の際は、説明変数間の相関が高いと *****/ /*****   多重共線性が出るのでよく注意すること *****/ proc corr; var age q6d--q6i; run;