*********************************************** * This file contains the SPSS syntax file and * the LIMDEP command file that Nishizawa prepared * for his JPSA99 paper. The top half is the * SPSS file and the second hald is the LIMDEP file. * * This file is in copyright by Yoshitaka Nishizawa * 25 Sept 1999 ***********************************************. ************************************************ NOTE: var name convention Starting section 9, I will use prefixs: P for parties C for candidates D for SMDs, and R for pR districts ************************************************. *********************************************** * Section 1 Reading the data ***********************************************. set length=none. GET FILE='C:\My Documents\........\jeds9698j2.sav'. EXECUTE . *********************************************** * Section 2a Grouping by vote pattern: dis and pr. * ***********************************************. *freq b2s1px b2s5. *選挙区、比例区の投票政党のコードの確認. *crosstabs table=b2s1px by b2s5. *crosstabs table=b2s1x by b2s5. *選挙区、比例区の一方のみ棄権について確認. recode b2s1px (1=1)(2 thru 11=0)(else=99) into DISTldp. recode b2s5 (1=1)(2 thru 7=0)(else=99) into PRldp. missing values DISTldp PRldp (99). *crosstabs table=b2s1px by distldp /table=b2s5 by PRldp. *自民投票に関する2つのダミーの作成 選挙区で自民:1, 非自民:0 DISTldp 比例区で自民:1, 非自民:0 PRdp. compute vtype4=-999. if (PRldp eq 1 and DISTldp eq 1) vtype4=1. if (PRldp eq 0 and DISTldp eq 0) vtype4=2. if (PRldp eq 1 and DISTldp eq 0) vtype4=3. if (PRldp eq 0 and DISTldp eq 1) vtype4=4. if (b2s5 eq 8 and DISTldp eq 1) vtype4=4. if (b2s5 eq 8 and DISTldp eq 0) vtype4=2. if (b2s1x eq 1 and PRldp eq 1) vtype4=3. if (b2s1x eq 1 and PRldp eq 0) vtype4=2. if (b2s1x eq 1 and b2s5 eq 8) vtype4=-999. missing values vtype4 (-999). execute. *freq var=vtype4. *crosstabs tables=PRldp by vtype4 /tables=DISTldp by vtype4. * DISTldp * 1 0 * * 1 (1) (3) * PRldp * 0 (4) (2). *(1)から(4)のタイプまで に整理、vote typeという 意味で vtype4とする. *なお、いずれか一方のみ棄権(30+-ケース)は (4),(3)に含めている . *********************************************** * Section 2b Dummy vars for party vote ***********************************************. compute Dldp=0. compute Dnfp=0. compute Ddpj=0. compute Dsdp=0. compute Djcp=0. compute Pldp=0. compute Pnfp=0. compute Pdpj=0. compute Psdp=0. compute Pjcp=0. if (b2s1px eq 1) Dldp=1. if (b2s1px eq 2) Dnfp=1. if (b2s1px eq 3) Ddpj=1. if (b2s1px eq 4) Dsdp=1. if (b2s1px eq 5) Djcp=1. if (b2s5 eq 1) Pldp=1. if (b2s5 eq 2) Pnfp=1. if (b2s5 eq 3) Pdpj=1. if (b2s5 eq 4) Psdp=1. if (b2s5 eq 5) Pjcp=1. *********************************************** * Section 3 Some IVs ***********************************************. *freq var a6 a6s1 a6s2 a9s2 a11 a11s1 a11s2ldp a14 a15ldp b36a to b36d. recode a6 (9,10=-999)(else=copy). recode a6s1 (3,4,9=-999)(else=copy). recode a6s2 (9,10,99=-999)(else=copy). recode a9s2 (6,7,9=-999)(else=copy). recode a11 (3,4=-999)(else=copy). recode a11s1 (8,9=-999)(else=copy). recode a11s2ldp (9=-999)(else=copy). recode a14 a15ldp (11,12=-999)(else=copy). missing values a6 a6s1 a6s2 a9s2 a11 a11s1 a11s2ldp a14 a15ldp(-999). execute. compute pidldp=-999. if (a6 eq 1 and a6s1 eq 1)pidldp=1. if (a6 eq 1 and a6s1 eq 2)pidldp=.75. if (a6 eq 8) pidldp=.5. if (a6 eq 1 and a6s1 eq 2)pidldp=.75. if ((a6 ge 2 and a6 le 7) and a6s1 eq 2)pidldp=.25. if ((a6 ge 2 and a6 le 7) and a6s1 eq 1)pidldp=0. missing values pidldp (-999). execute. *freq var=pidldp. *pidLDP: 自民支持の強さ(1〜0). compute pidstrgx=0. if (a6s1 eq 1) pidstrgx=1. if (a6s1 ge 2 and a6s1 le 4) pidstrgx=.67. if (a6s2 ge 1 and a6s2 le 7) pidstrgx=.33. *pidSTRGx: 政党支持の強さ(1〜0). recode a9s2 (1=1)(2=.75)(3=.5)(4=.25)(5=0)(else=-999) into probWIN. recode a9s2 (1=1)(2=.75)(3=.5)(4=.25)(5=0)(else=.5) into probWINx. missing values probWIN (-999). *probWIN: 意中候補者の当選確率(1〜0). *probWINx: 意中候補者の当選確率(1〜0)NA,DK 救済版. compute coalit=-999. if (a11 eq 1 and a11s1 eq 1) coalit=1. if (a11 eq 2 and a11s2ldp eq 1) coalit=.66. if (a11 eq 2 and a11s1 eq 0) coalit=.33. if (a11 eq 1 and a11s1 ne 1) coalit=0. compute coalitx=.5. if (a11 eq 1 and a11s1 eq 1) coalitx=1. if (a11 eq 2 and a11s2ldp eq 1) coalitx=.66. if (a11 eq 2 and a11s1 eq 0) coalitx=.33. if (a11 eq 1 and a11s1 ne 1) coalitx=0. missing values probWIN coalit (-999). *coalit: 自民単独政権から非自民単独政権(1〜0). *coalitx: 自民単独政権から非自民単独政権(1〜0)NA,DK 救済版. compute ideology=abs(a14 - a15ldp). execute. *ideology: 対自民・イデオロギー距離. *********************************************** * Section 4 Calculating factor scores for A2 an A5 serise ***********************************************. missing values a2s1ldp a2s2ldp a2s3ldp a2s4ldp a2s5ldp a2s6ldp a2s7ldp a2s8ldp a2s9ldp a2s10ldp a2s11ldp a2s12ldp a2s13ldp a2s14ldp a2s15ldp a2s1nfp a2s2nfp a2s3nfp a2s4nfp a2s5nfp a2s6nfp a2s7nfp a2s8nfp a2s9nfp a2s10nfp a2s11nfp a2s12nfp a2s13nfp a2s14nfp a2s15nfp a2s1dpj a2s2dpj a2s3dpj a2s4dpj a2s5dpj a2s6dpj a2s7dpj a2s8dpj a2s9dpj a2s10dpj a2s11dpj a2s12dpj a2s13dpj a2s14dpj a2s15dpj a2s1sdp a2s2sdp a2s3sdp a2s4sdp a2s5sdp a2s6sdp a2s7sdp a2s8sdp a2s9sdp a2s10sdp a2s11sdp a2s12sdp a2s13sdp a2s14sdp a2s15sdp a2s1jcp a2s2jcp a2s3jcp a2s4jcp a2s5jcp a2s6jcp a2s7jcp a2s8jcp a2s9jcp a2s10jcp a2s11jcp a2s12jcp a2s13jcp a2s14jcp a2s15jcp (9). FACTOR /VARIABLES a2s1ldp a2s2ldp a2s3ldp a2s4ldp a2s5ldp a2s6ldp a2s7ldp a2s8ldp a2s9ldp a2s10ldp a2s11ldp a2s12ldp a2s13ldp a2s14ldp a2s15ldp /MISSING LISTWISE /ANALYSIS a2s1ldp a2s2ldp a2s3ldp a2s4ldp a2s5ldp a2s6ldp a2s7ldp a2s8ldp a2s9ldp a2s10ldp a2s11ldp a2s12ldp a2s13ldp a2s14ldp a2s15ldp /PRINT INITIAL EXTRACTION /CRITERIA FACTORS(2) ITERATE(25) /EXTRACTION PC /ROTATION NOROTATE /SAVE REG(all) /METHOD=CORRELATION . FACTOR /VARIABLES a2s1nfp a2s2nfp a2s3nfp a2s4nfp a2s5nfp a2s6nfp a2s7nfp a2s8nfp a2s9nfp a2s10nfp a2s11nfp a2s12nfp a2s13nfp a2s14nfp a2s15nfp /MISSING LISTWISE /ANALYSIS a2s1nfp a2s2nfp a2s3nfp a2s4nfp a2s5nfp a2s6nfp a2s7nfp a2s8nfp a2s9nfp a2s10nfp a2s11nfp a2s12nfp a2s13nfp a2s14nfp a2s15nfp /PRINT INITIAL EXTRACTION /CRITERIA FACTORS(2) ITERATE(25) /EXTRACTION PC /ROTATION NOROTATE /SAVE REG(all) /METHOD=CORRELATION . FACTOR /VARIABLES a2s1dpj a2s2dpj a2s3dpj a2s4dpj a2s5dpj a2s6dpj a2s7dpj a2s8dpj a2s9dpj a2s10dpj a2s11dpj a2s12dpj a2s13dpj a2s14dpj a2s15dpj /MISSING LISTWISE /ANALYSIS a2s1dpj a2s2dpj a2s3dpj a2s4dpj a2s5dpj a2s6dpj a2s7dpj a2s8dpj a2s9dpj a2s10dpj a2s11dpj a2s12dpj a2s13dpj a2s14dpj a2s15dpj /PRINT INITIAL EXTRACTION /CRITERIA FACTORS(2) ITERATE(25) /EXTRACTION PC /ROTATION NOROTATE /SAVE REG(all) /METHOD=CORRELATION . FACTOR /VARIABLES a2s1sdp a2s2sdp a2s3sdp a2s4sdp a2s5sdp a2s6sdp a2s7sdp a2s8sdp a2s9sdp a2s10sdp a2s11sdp a2s12sdp a2s13sdp a2s14sdp a2s15sdp /MISSING LISTWISE /ANALYSIS a2s1sdp a2s2sdp a2s3sdp a2s4sdp a2s5sdp a2s6sdp a2s7sdp a2s8sdp a2s9sdp a2s10sdp a2s11sdp a2s12sdp a2s13sdp a2s14sdp a2s15sdp /PRINT INITIAL EXTRACTION /CRITERIA FACTORS(2) ITERATE(25) /EXTRACTION PC /ROTATION NOROTATE /SAVE REG(all) /METHOD=CORRELATION . FACTOR /VARIABLES a2s1jcp a2s2jcp a2s3jcp a2s4jcp a2s5jcp a2s6jcp a2s7jcp a2s8jcp a2s9jcp a2s10jcp a2s11jcp a2s12jcp a2s13jcp a2s14jcp a2s15jcp /MISSING LISTWISE /ANALYSIS a2s1jcp a2s2jcp a2s3jcp a2s4jcp a2s5jcp a2s6jcp a2s7jcp a2s8jcp a2s9jcp a2s10jcp a2s11jcp a2s12jcp a2s13jcp a2s14jcp a2s15jcp /PRINT INITIAL EXTRACTION /CRITERIA FACTORS(2) ITERATE(25) /EXTRACTION PC /ROTATION NOROTATE /SAVE REG(all) /METHOD=CORRELATION . FACTOR /VARIABLES a5s1ldp a5s2ldp a5s3ldp a5s4ldp a5s5ldp a5s6ldp a5s7ldp a5s8ldp a5s9ldp a5s10ldp a5s11ldp a5s12ldp a5s13ldp a5s14ldp a5s15ldp /MISSING LISTWISE /ANALYSIS a5s1ldp a5s2ldp a5s3ldp a5s4ldp a5s5ldp a5s6ldp a5s7ldp a5s8ldp a5s9ldp a5s10ldp a5s11ldp a5s12ldp a5s13ldp a5s14ldp a5s15ldp /PRINT INITIAL EXTRACTION /CRITERIA FACTORS(2) ITERATE(25) /EXTRACTION PC /ROTATION NOROTATE /SAVE REG(all) /METHOD=CORRELATION . FACTOR /VARIABLES a5s1nfp a5s2nfp a5s3nfp a5s4nfp a5s5nfp a5s6nfp a5s7nfp a5s8nfp a5s9nfp a5s10nfp a5s11nfp a5s12nfp a5s13nfp a5s14nfp a5s15nfp /MISSING LISTWISE /ANALYSIS a5s1nfp a5s2nfp a5s3nfp a5s4nfp a5s5nfp a5s6nfp a5s7nfp a5s8nfp a5s9nfp a5s10nfp a5s11nfp a5s12nfp a5s13nfp a5s14nfp a5s15nfp /PRINT INITIAL EXTRACTION /CRITERIA FACTORS(2) ITERATE(25) /EXTRACTION PC /ROTATION NOROTATE /SAVE REG(all) /METHOD=CORRELATION . FACTOR /VARIABLES a5s1dpj a5s2dpj a5s3dpj a5s4dpj a5s5dpj a5s6dpj a5s7dpj a5s8dpj a5s9dpj a5s10dpj a5s11dpj a5s12dpj a5s13dpj a5s14dpj a5s15dpj /MISSING LISTWISE /ANALYSIS a5s1dpj a5s2dpj a5s3dpj a5s4dpj a5s5dpj a5s6dpj a5s7dpj a5s8dpj a5s9dpj a5s10dpj a5s11dpj a5s12dpj a5s13dpj a5s14dpj a5s15dpj /PRINT INITIAL EXTRACTION /CRITERIA FACTORS(2) ITERATE(25) /EXTRACTION PC /ROTATION NOROTATE /SAVE REG(all) /METHOD=CORRELATION . FACTOR /VARIABLES a5s1sdp a5s2sdp a5s3sdp a5s4sdp a5s5sdp a5s6sdp a5s7sdp a5s8sdp a5s9sdp a5s10sdp a5s11sdp a5s12sdp a5s13sdp a5s14sdp a5s15sdp /MISSING LISTWISE /ANALYSIS a5s1sdp a5s2sdp a5s3sdp a5s4sdp a5s5sdp a5s6sdp a5s7sdp a5s8sdp a5s9sdp a5s10sdp a5s11sdp a5s12sdp a5s13sdp a5s14sdp a5s15sdp /PRINT INITIAL EXTRACTION /CRITERIA FACTORS(2) ITERATE(25) /EXTRACTION PC /ROTATION NOROTATE /SAVE REG(all) /METHOD=CORRELATION . FACTOR /VARIABLES a5s1jcp a5s2jcp a5s3jcp a5s4jcp a5s5jcp a5s6jcp a5s7jcp a5s8jcp a5s9jcp a5s10jcp a5s11jcp a5s12jcp a5s13jcp a5s14jcp a5s15jcp /MISSING LISTWISE /ANALYSIS a5s1jcp a5s2jcp a5s3jcp a5s4jcp a5s5jcp a5s6jcp a5s7jcp a5s8jcp a5s9jcp a5s10jcp a5s11jcp a5s12jcp a5s13jcp a5s14jcp a5s15jcp /PRINT INITIAL EXTRACTION /CRITERIA FACTORS(2) ITERATE(25) /EXTRACTION PC /ROTATION NOROTATE /SAVE REG(all) /METHOD=CORRELATION . rename var (fac1_1 fac2_1 fac1_2 fac2_2 fac1_3 fac2_3 fac1_4 fac2_4 fac1_5 fac2_5 fac1_6 fac2_6 fac1_7 fac2_7 fac1_8 fac2_8 fac1_9 fac2_9 fac1_10 fac2_10 = cldp1 cldp2 cnfp1 cnfp2 cdpj1 cdpj2 csdp1 csdp2 cjcp1 cjcp2 pldp1 pldp2 pnfp1 pnfp2 pdpj1 pdpj2 psdp1 psdp2 pjcp1 pjcp2 ). compute coth1=max(cnfp1,cdpj1,csdp1,cjcp1). if (sysmis(coth1)) coth1=0. compute coth2=min(cnfp2,cdpj2,csdp2,cjcp2). if (sysmis(coth2)) coth2=0. compute poth1=max(pnfp1,pdpj1,psdp1,pjcp1). if (sysmis(poth1)) poth1=0. compute poth2=min(pnfp2,pdpj2,psdp2,pjcp2). if (sysmis(poth2)) poth2=0. *各政党への投票と候補者・政党の評価 についての因子の相関を見たところ、 第一因子についは正の相関、第二因子に                   ついては負の相関があるので、 非自民について、第一は最大値、第二は                   最小値を取ることにした。 また、欠損値については、ゼロを当てた。. *********************************************** * Section 5 Calculating election competitiveness このセクションを実行するには、選挙区データの 追加モジュールが必要。SSJDAからのjeds96j1,savには 含まれない変数があるので注意。 ***********************************************. if (b2s1 eq 1) b2s1cv=unqc1v. *freq var b2s1. compute b2s1Cv=-999. if (b2s1 eq 1) b2s1cv=unqc1v. if (b2s1 eq 2) b2s1cv=unqc2v. if (b2s1 eq 3) b2s1cv=unqc3v. if (b2s1 eq 4) b2s1cv=unqc4v. if (b2s1 eq 5) b2s1cv=unqc5v. if (b2s1 eq 6) b2s1cv=unqc6v. if (b2s1 eq 7) b2s1cv=unqc7v. missing values b2s1cv (-999). compute myCwin=b2s1Cv/max(unqc1v,unqc2v,unqc3v, unqc4v,unqc5v,unqc6v,unqc7v)*100. *list var= b2s1 unqc1v unqc2v unqc3v unqc4v unqc5v unqc6v unqc7v myCwin /cases= from 1 to 50. *execute. *b2s1Cv 投票候補者の得票数 myCwin 投票候補者の得票数/当選者得票数. *********************************************** * Section 6 Considering the role of issue evaluation ***********************************************. missing values a14 a15ldp (). *freq a14 a15ldp to a15jcp a16 a17ldp to a17jcp a18 a19 a20ldp to a20jcp a21 a22x1 to a22x20 a22s1 a22s2px a22s3 a22s4 a23s1px a23s2. missing values a14 a15ldp (-999). missing values a15nfp to a15jcp a16 a17ldp to a17jcp a18 a20ldp to a20jcp (11,12). missing values a19 a21 (5,6). ***** Proximility Model. *イデオロギー距離(iss1Pxxx). do repeat d=iss1Pldp iss1Pnfp iss1Pdpj iss1Psdp iss1Pjcp /p=a15ldp to a15jcp. compute d=(p - a14)*(p - a14). if (sysmis(d)) d=100. end repeat. compute iss1Poth=min(iss1Pnfp, iss1Pdpj, iss1Psdp, iss1Pjcp). *改革・安定-距離(iss2Pxxx). do repeat d=iss2Pldp iss2Pnfp iss2Pdpj iss2Psdp iss2Pjcp /p=a17ldp to a17jcp. compute d=(p - a16)*(p - a16). if (sysmis(d)) d=100. end repeat. compute iss2Poth=min(iss2Pnfp, iss2Pdpj, iss2Psdp, iss2Pjcp). *行政改革-距離(iss3Pxxx). do repeat d=iss3Pldp iss3Pnfp iss3Pdpj iss3Psdp iss3Pjcp /p=a20ldp to a20jcp. compute d=(p - a18)*(p - a18). if (sysmis(d)) d=100. end repeat. compute iss3Poth=min(iss3Pnfp, iss3Pdpj, iss3Psdp, iss3Pjcp). execute. *2争点:イデオロギー/改革・安定-距離(iss4Pxxx). *the square of the Euclidean distance the Euclidean distanceは、 「距離の二乗の合計」について、それの 平方根を取ったもの。そして、それを また二乗するので、結局は、「距離の 二乗の合計」と 同じとなる。 (Merill & Grofman 1999, 21) do repeat d=iss4Pldp iss4Pnfp iss4Pdpj iss4Psdp iss4Pjcp /p1=a15ldp to a15jcp /p2=a17ldp to a17jcp. compute d=(p1 - a14)*(p1 - a14) + (p2 - a16)*(p2 - a16). if (sysmis(d)) d=200. end repeat. compute iss4Poth=min(iss4Pnfp, iss4Pdpj, iss4Psdp, iss4Pjcp). *2争点:イデオロギー/行政改革-距離(iss5Pxxx). do repeat d=iss5Pldp iss5Pnfp iss5Pdpj iss5Psdp iss5Pjcp /p1=a15ldp to a15jcp /p2=a20ldp to a20jcp. compute d=(p1 - a14)*(p1 - a14) + (p2 - a18)*(p2 - a18). if (sysmis(d)) d=200. end repeat. compute iss5Poth=min(iss5Pnfp, iss5Pdpj, iss5Psdp, iss5Pjcp). *2争点:改革・安定/行政改革-距離(iss6Pxxx). do repeat d=iss6Pldp iss6Pnfp iss6Pdpj iss6Psdp iss6Pjcp /p1=a17ldp to a17jcp /p2=a20ldp to a20jcp. compute d=(p1 - a16)*(p1 - a16) + (p2 - a18)*(p2 - a18). if (sysmis(d)) d=200. end repeat. compute iss6Poth=min(iss6Pnfp, iss6Pdpj, iss6Psdp, iss6Pjcp). *3争点:イデオロギー/改革・安定/行政改革-距離(iss7Pxxx). do repeat d=iss7Pldp iss7Pnfp iss7Pdpj iss7Psdp iss7Pjcp /p1=a15ldp to a15jcp /p2=a17ldp to a17jcp /p3=a20ldp to a20jcp. compute d=(p1 - a14)*(p1 - a14) + (p2 - a16)*(p2 - a16) + (p3 - a18)*(p3 - a18). if (sysmis(d)) d=300. end repeat. compute iss7Poth=min(iss7Pnfp, iss7Pdpj, iss7Psdp, iss7Pjcp). ***** Directional Model. *イデオロギー距離(iss1Dxxx). do repeat d=iss1Dldp iss1Dnfp iss1Ddpj iss1Dsdp iss1Djcp /p=a15ldp to a15jcp. compute pp=p-5. compute d=pp*(a14 - 5). if (sysmis(d)) d=-25. end repeat. compute iss1Doth=max(iss1Dnfp, iss1Ddpj, iss1Dsdp, iss1Djcp). *改革・安定-距離(iss2Dxxx). do repeat d=iss2Dldp iss2Dnfp iss2Ddpj iss2Dsdp iss2Djcp /p=a17ldp to a17jcp. compute pp=p-5. compute d=pp*(a16 - 5). if (sysmis(d)) d=-25. end repeat. compute iss2Doth=max(iss2Dnfp, iss2Ddpj, iss2Dsdp, iss2Djcp). *行政改革-距離(iss3Dxxx). do repeat d=iss3Dldp iss3Dnfp iss3Ddpj iss3Dsdp iss3Djcp /p=a20ldp to a20jcp. compute pp=p-5. compute d=pp*(a18 - 5). if (sysmis(d)) d=-25. end repeat. compute iss3Doth=max(iss3Dnfp, iss3Ddpj, iss3Dsdp, iss3Djcp). *2争点:イデオロギー/改革・安定-距離(iss4Dxxx). do repeat d=iss4Dldp iss4Dnfp iss4Ddpj iss4Dsdp iss4Djcp /p1=a15ldp to a15jcp /p2=a17ldp to a17jcp. compute pp1=p1-5. compute pp2=p2-5. compute d=pp1*(a14 - 5) + pp2*(a16 - 5) . if (sysmis(d)) d=-50. end repeat. compute iss4Doth=max(iss4Dnfp, iss4Ddpj, iss4Dsdp, iss4Djcp). *2争点:イデオロギー/行政改革-距離(iss5Dxxx). do repeat d=iss5Dldp iss5Dnfp iss5Ddpj iss5Dsdp iss5Djcp /p1=a15ldp to a15jcp /p2=a20ldp to a20jcp. compute pp1=p1-5. compute pp2=p2-5. compute d=pp1*(a14 - 5) + pp2*(a18 - 5) . if (sysmis(d)) d=-50. end repeat. compute iss5Doth=max(iss5Dnfp, iss5Ddpj, iss5Dsdp, iss5Djcp). *2争点:改革・安定/行政改革-距離(iss6Dxxx). do repeat d=iss6Dldp iss6Dnfp iss6Ddpj iss6Dsdp iss6Djcp /p1=a17ldp to a17jcp /p2=a20ldp to a20jcp. compute pp1=p1-5. compute pp2=p2-5. compute d=pp1*(a16 - 5) + pp2*(a18 - 5) . if (sysmis(d)) d=-50. end repeat. compute iss6Doth=max(iss6Dnfp, iss6Ddpj, iss6Dsdp, iss6Djcp). *3争点:イデオロギー/改革・安定/行政改革-距離(iss7Dxxx). do repeat d=iss7Dldp iss7Dnfp iss7Ddpj iss7Dsdp iss7Djcp /p1=a15ldp to a15jcp /p2=a17ldp to a17jcp /p3=a20ldp to a20jcp. compute pp1=p1-5. compute pp2=p2-5. compute pp3=p3-5. compute d=pp1*(a14 - 5) + pp2*(a16 - 5) + pp3*(a18 - 5) . if (sysmis(d)) d=-75. end repeat. compute iss7Doth=max(iss7Dnfp, iss7Ddpj, iss7Dsdp, iss7Djcp). *取り組んでくれそうな候補者・政党. compute cNisLDP=0. compute cNisOTH=0. compute pNisLDP=0. compute pNisOTH=0. compute cLisLDP=0. compute cLisOTH=0. compute pLisLDP=0. compute pLisOTH=0. if (a22s2px eq 1) cNisLDP=1. if (a22s2px ge 2 and a22s2px le 11) cNisOTH=1. if (a22s3 eq 1) pNisLDP=1. if (a22s3 ge 2 and a22s3 le 7) pNisOTH=1. if (a23s1px eq 1) cLisLDP=1. if (a23s1px ge 2 and a23s1px le 11) cLisOTH=1. if (a23s2 eq 1) pLisLDP=1. if (a23s2 ge 2 and a23s2 le 7) pLisOTH=1. *list var= a22s2px cNisLDP cNisOTH a22s3 pNisLDP pNisOTH /cases= from 1 to 50. *list var= a23s1px cLisLDP cLisOTH a23s2 pLisLDP pLisOTH /cases= from 1 to 50. execute. *********************************************** Section 7: Other independent variables ***********************************************. *民主主義に対する満足度. recode a24 (1=1)(2=.75)(3=.25)(4=0)(5, 6=.5) into satdemo. *生活に対する満足度. recode a28 (1=1)(2=.75)(3,6,7=.5)(4=.25)(5=1) into satlife. *景気動向評価. recode a28 (1=1)(2=.75)(3,6,7=.5)(4=.25)(5=1) into satecon. *一般的内閣評価. recode a21 (1=1)(2=.75)(3, 6, 7=.5)(4=.25)(5=0) into hashimot. *********************************************** Section 8: mobilization variables ***********************************************. *crosstabs tables b8s2p1x by b8s2ldp. *compute ldp=0. *if (b8s2p1x eq 1 or b8s2p2x eq 1 or b8s2p3x eq 1) ldp=1. *crosstabs tables ldp by b8s2ldp. *compute unknown=0. *if (ldp eq 1 and b8s2ldp eq 0) unknown=1. *if (ldp eq 0 and b8s2ldp eq 1) unknown=1. *temporary. *select if (unknown eq 1). *list var= caseid ldp b8s2ldp b8s2c1 b8s2p1x b8s2p2x b8s2p3x. *temporary. *select if (unknown eq 1). *list var= caseid ldp b8s2ldp to b8s2oth. *他者からの投票依頼. compute askldp=0. compute askoth=0. do repeat p=b8s2p1x b8s2p2x b8s2p2x. if (p eq 1) askldp=askldp + 1. if (p ge 2 and p le 11) askoth=askoth + 1. end repeat. *政党・候補者からの動員. compute cmobldp=0. compute cmoboth=0. do repeat b= b9ap1x b9ap2x b9ap3x b9bp1x b9bp2x b9bp3x b9cp1x b9cp2x b9cp3x. if (b eq 1) cmobldp=cmobldp + 1. if (b ge 2 and b le 11) cmoboth=cmoboth + 1. end repeat. *freq var b9aldp b9bldp b9cldp b9anfp b9adpj b9asdp b9ajcp b9anph b9aoth b9bnfp b9bdpj b9bsdp b9bjcp b9bnph b9both b9cnfp b9cdpj b9csdp b9cjcp b9cnph b9coth. recode b9aldp b9bldp b9cldp b9anfp b9adpj b9asdp b9ajcp b9anph b9aoth b9bnfp b9bdpj b9bsdp b9bjcp b9bnph b9both b9cnfp b9cdpj b9csdp b9cjcp b9cnph b9coth (9=0). compute pmobldp=b9aldp + b9bldp + b9cldp. compute pmoboth= b9anfp + b9adpj + b9asdp + b9ajcp + b9anph + b9aoth + b9bnfp + b9bdpj + b9bsdp + b9bjcp + b9bnph + b9both + b9cnfp + b9cdpj + b9csdp + b9cjcp + b9cnph + b9coth. *freq var pmobldp pmoboth. *********************************************** Section 9: demografic variables ***********************************************. compute age = f2age. compute jimae = 0. if (f10 eq 2 or f10 eq 3 or f11s1 eq 2 or f11s1 eq 3) jimae=1. recode f1 (1=1)(2=0)into male. recode f7 (9,10=-999)(else=copy) into income. recode f6 (1, 8=6)(2, 9=9)(3, 10, 11=12) (4, 5, 12, 13=14)(6, 14=16)(7=18)into educ. recode f10s7 (1=1)(else=0) into union. recode b36b b36d b36e (1=1)(else=0) into memberR memberF memberC. *********************************************** * Section 10 write file ***********************************************. do repeat v= vtype4 pidSTRGx probWINx coalitx cldp1 cldp2 coth1 coth2 pldp1 pldp2 poth1 poth2 myCwin iss1Pldp iss1Poth iss2Pldp iss2Poth iss3Pldp iss3Poth iss4Pldp iss4Poth iss5Pldp iss5Poth iss6Pldp iss6Poth iss7Pldp iss7Poth iss1Dldp iss1Doth iss2Dldp iss2Doth iss3Dldp iss3Doth iss4Dldp iss4Doth iss5Dldp iss5Doth iss6Dldp iss6Doth iss7Dldp iss7Doth cNisLDP cNisOTH pNisLDP pNisOTH cLisLDP cLisOTH pLisLDP pLisOTH satdemo satlife satecon hashimot askldp askoth cmobldp cmoboth pmobldp pmoboth age jimae male income educ union memberR memberF memberC. if (missing(v)) v=-999. end repeat. missing values vtype4 pidSTRGx probWINx coalitx cldp1 cldp2 coth1 coth2 pldp1 pldp2 poth1 poth2 myCwin iss1Pldp iss1Poth iss2Pldp iss2Poth iss3Pldp iss3Poth iss4Pldp iss4Poth iss5Pldp iss5Poth iss6Pldp iss6Poth iss7Pldp iss7Poth iss1Dldp iss1Doth iss2Dldp iss2Doth iss3Dldp iss3Doth iss4Dldp iss4Doth iss5Dldp iss5Doth iss6Dldp iss6Doth iss7Dldp iss7Doth cNisLDP cNisOTH pNisLDP pNisOTH cLisLDP cLisOTH pLisLDP pLisOTH satdemo satlife satecon hashimot askldp askoth cmobldp cmoboth pmobldp pmoboth age jimae male income educ union memberR memberF memberC (-999). write OUTFILE='C:\My Documents\.......\jpsa99pub.dat' / caseid vtype4 pidSTRGx probWINx coalitx cldp1 cldp2 coth1 coth2 pldp1 pldp2 poth1 poth2 myCwin iss1Pldp iss1Poth iss2Pldp iss2Poth iss3Pldp iss3Poth iss4Pldp iss4Poth iss5Pldp iss5Poth iss6Pldp iss6Poth iss7Pldp iss7Poth iss1Dldp iss1Doth iss2Dldp iss2Doth iss3Dldp iss3Doth iss4Dldp iss4Doth iss5Dldp iss5Doth iss6Dldp iss6Doth iss7Dldp iss7Doth cNisLDP cNisOTH pNisLDP pNisOTH cLisLDP cLisOTH pLisLDP pLisOTH satdemo satlife satecon hashimot askldp askoth cmobldp cmoboth pmobldp pmoboth age jimae male income educ union memberR memberF memberC (69f12.5). execute. /* ========================================================= LIMDEP program begins here ========================================================= */ /* Best to set font to @MS Gothic, 11 point. */ /* ***************************************** Section 1. Prepareing the panel data ******************************************** */ rows; 4000 $ read; nvar=69; nobs=1535; names=caseid, vtype4, pidSTRGx, probWINx, coalitx, cldp1, cldp2, coth1, coth2, pldp1, pldp2, poth1, poth2, myCwin, iss1Pldp, iss1Poth, iss2Pldp, iss2Poth, iss3Pldp, iss3Poth, iss4Pldp, iss4Poth, iss5Pldp, iss5Poth, iss6Pldp, iss6Poth, iss7Pldp, iss7Poth, iss1Dldp, iss1Doth, iss2Dldp, iss2Doth, iss3Dldp, iss3Doth, iss4Dldp, iss4Doth, iss5Dldp, iss5Doth, iss6Dldp, iss6Doth, iss7Dldp, iss7Doth, cNisLDP, cNisOTH, pNisLDP, pNisOTH, cLisLDP, cLisOTH, pLisLDP, pLisOTH, satdemo, satlife, satecon, hashimot, askldp, askoth, cmobldp, cmoboth, pmobldp, pmoboth, age, jimae, male, income, educ, union, memberR, memberF, memberC; file=jpsa99pub.dat; format=(69f12.5)$ skip create; LDP_ldp=(vtype4=1) ; OTH_oth=(vtype4=2) ; LDP_oth=(vtype4=3) ; OTH_ldp=(vtype4=4)$ create; screen=0$ create; if (vtype4 = 1) screen=1$ create; if (vtype4 = 2) screen=1$ create; if (vtype4 = 3) screen=1$ create; if (vtype4 = 4) screen=1$ reject; screen=0$ /* LDP_ldp ==> PR - Dist */ write; caseid, LDP_ldp, pldp1, cldp1, pldp2, cldp2, cLDP1, pldp1, cLDP2, pldp2, iss1Pldp, iss1Pldp, iss2Pldp, iss2Pldp, iss3Pldp, iss3Pldp, iss4Pldp, iss4Pldp, iss5Pldp, iss5Pldp, iss6Pldp, iss6Pldp, iss7Pldp, iss7Pldp, iss1Dldp, iss1Dldp, iss2Dldp, iss2Dldp, iss3Dldp, iss3Dldp, iss4Dldp, iss4Dldp, iss5Dldp, iss5Dldp, iss6Dldp, iss6Dldp, iss7Dldp, iss7Dldp, pNisLDP, cNisLDP, pLisLDP, cLisLDP, pmobldp, cmobldp, pidSTRGx, probWINx, coalitx, myCwin, satdemo, satlife, satecon, hashimot, askldp, askoth, age, jimae, male, income, educ, union, memberR, memberF, memberC, caseid, OTH_oth, poth1, coth1, poth2, coth2, cOTH1, poth1, cOTH2, poth2, iss1Poth, iss1Poth, iss2Poth, iss2Poth, iss3Poth, iss3Poth, iss4Poth, iss4Poth, iss5Poth, iss5Poth, iss6Poth, iss6Poth, iss7Poth, iss7Poth, iss1Doth, iss1Doth, iss2Doth, iss2Doth, iss3Doth, iss3Doth, iss4Doth, iss4Doth, iss5Doth, iss5Doth, iss6Doth, iss6Doth, iss7Doth, iss7Doth, pNisOTH, cNisOTH, pLisOTH, cLisOTH, pmoboth, cmoboth, pidSTRGx, probWINx, coalitx, myCwin, satdemo, satlife, satecon, hashimot, askldp, askoth, age, jimae, male, income, educ, union, memberR, memberF, memberC, caseid, LDP_oth, pldp1, coth1, pldp2, coth2, cLDP1, poth1, cLDP2, poth2, iss1Pldp, iss1Poth, iss2Pldp, iss2Poth, iss3Pldp, iss3Poth, iss4Pldp, iss4Poth, iss5Pldp, iss5Poth, iss6Pldp, iss6Poth, iss7Pldp, iss7Poth, iss1Dldp, iss1Doth, iss2Dldp, iss2Doth, iss3Dldp, iss3Doth, iss4Dldp, iss4Doth, iss5Dldp, iss5Doth, iss6Dldp, iss6Doth, iss7Dldp, iss7Doth, pNisLDP, cNisOTH, pLisLDP, cLisOTH, pmobldp, cmoboth, pidSTRGx, probWINx, coalitx, myCwin, satdemo, satlife, satecon, hashimot, askldp, askoth, age, jimae, male, income, educ, union, memberR, memberF, memberC, caseid, OTH_ldp, poth1, cldp1, poth2, cldp2, cOTH1, pldp1, cOTH2, pldp2, iss1Poth, iss1Pldp, iss2Poth, iss2Pldp, iss3Poth, iss3Pldp, iss4Poth, iss4Pldp, iss5Poth, iss5Pldp, iss6Poth, iss6Pldp, iss7Poth, iss7Pldp, iss1Doth, iss1Dldp, iss2Doth, iss2Dldp, iss3Doth, iss3Dldp, iss4Doth, iss4Dldp, iss5Doth, iss5Dldp, iss6Doth, iss6Dldp, iss7Doth, iss7Dldp, pNisOTH, cNisLDP, pLisOTH, cLisLDP, pmoboth, cmobldp, pidSTRGx, probWINx, coalitx, myCwin, satdemo, satlife, satecon, hashimot, askldp, askoth, age, jimae, male, income, educ, union, memberR, memberF, memberC; file=pub.dat; format=(63f12.5)$ /* Need to exit file or RESET the dataset, here to re-read the file just created. */ /* ***************************************** Section 2. Reading the panel data ******************************************** */ rows;4000$ read; nvar=63; nobs=3804 ; names=caseid, vote, Rp1, Dc1, Rp2, Dc2, Rc1, Dp1, Rc2, Dp2, Riss1P, Diss1P, Riss2P, Diss2P, Riss3P, Diss3P, Riss4P, Diss4P, Riss5P, Diss5P, Riss6P, Diss6P, Riss7P, Diss7P, Riss1D, Diss1D, Riss2D, Diss2D, Riss3D, Diss3D, Riss4D, Diss4D, Riss5D, Diss5D, Riss6D, Diss6D, Riss7D, Diss7D, RpNis, DcNis, RpLis, DcLis, Rpmob, Dcmob, pidSTRGx, probWINx, coalitx, myCwin, satdemo, satlife, satecon, hashimot, askldp, askoth, age, jimae, male, income, educ, union, memberR, memberF, memberC; file=pub.dat; format=(63f12.5)$ create; bestbet=.35$ create; Riss6M35=(2*(1-bestbet)*Riss6D) - (bestbet*Riss6P)$ create; Diss6M35=(2*(1-bestbet)*Diss6D) - (bestbet*Diss6P)$ /* Analysis begins here */ /* The tree structure assumed below is following: PR level (first choice): LDP OTH _____|____ _____|_____ | | | | SG level (second choice): LDPldp LDPoth OTHldp OTHoth vtype4 (1) (3) (4) (2) */ /* Constants in one of the branch equations seem to work fine, just as in the equation below. */ /* ***************************************** Section 3. Unified spacial models ****************************************** */ nlog; lhs=vote; start=logit; choices=LDPldp, OTHoth, LDPoth, OTHldp; tree=election[LDP(LDPldp, LDPoth), OTH(OTHldp, OTHoth)]; model: u(ldpldp, othoth, ldpoth, othldp)= + bDc1*Dc1 + bDc2*Dc2 + bDiss6P*Diss6P + bDiss6D*Diss6D/ u(ldp, oth)= + bRp1*Rp1 + bRp2*Rp2 + bRiss6P*Riss6P + bRiss6D*Riss6D$ create; bet=.2$ create; Riss6MG2=(2*(1-bet)*Riss6D) - (bet*Riss6P)$ create; Diss6MG2=(2*(1-bet)*Diss6D) - (bet*Diss6P)$ create; bet=.3$ create; Riss6MG3=(2*(1-bet)*Riss6D) - (bet*Riss6P)$ create; Diss6MG3=(2*(1-bet)*Diss6D) - (bet*Diss6P)$ create; bet=.4$ create; Riss6MG4=(2*(1-bet)*Riss6D) - (bet*Riss6P)$ create; Diss6MG4=(2*(1-bet)*Diss6D) - (bet*Diss6P)$ create; bet=.5$ create; Riss6MG5=(2*(1-bet)*Riss6D) - (bet*Riss6P)$ create; Diss6MG5=(2*(1-bet)*Diss6D) - (bet*Diss6P)$ create; bet=.6$ create; Riss6MG6=(2*(1-bet)*Riss6D) - (bet*Riss6P)$ create; Diss6MG6=(2*(1-bet)*Diss6D) - (bet*Diss6P)$ create; bet=.7$ create; Riss6MG7=(2*(1-bet)*Riss6D) - (bet*Riss6P)$ create; Diss6MG7=(2*(1-bet)*Diss6D) - (bet*Diss6P)$ create; bet=.8$ create; Riss6MG8=(2*(1-bet)*Riss6D) - (bet*Riss6P)$ create; Diss6MG8=(2*(1-bet)*Diss6D) - (bet*Diss6P)$ create; bet=.9$ create; Riss6MG9=(2*(1-bet)*Riss6D) - (bet*Riss6P)$ create; Diss6MG9=(2*(1-bet)*Diss6D) - (bet*Diss6P)$ nlog; lhs=vote; start=logit; choices=LDPldp, OTHoth, LDPoth, OTHldp; tree=election[LDP(LDPldp, LDPoth), OTH(OTHldp, OTHoth)]; model: u(ldpldp, othoth, ldpoth, othldp)= + bDc1*Dc1 + bDc2*Dc2 + bDiss6MG2*Diss6MG2/ u(ldp, oth)= + bRp1*Rp1 + bRp2*Rp2 + bRiss6MG2*Riss6MG2$ /* Adj. R2: .37983 */ nlog; lhs=vote; start=logit; choices=LDPldp, OTHoth, LDPoth, OTHldp; tree=election[LDP(LDPldp, LDPoth), OTH(OTHldp, OTHoth)]; model: u(ldpldp, othoth, ldpoth, othldp)= + bDc1*Dc1 + bDc2*Dc2 + bDiss6MG3*Diss6MG3/ u(ldp, oth)= + bRp1*Rp1 + bRp2*Rp2 + bRiss6MG3*Riss6MG3$ /* Adj. R2: .37988 */ nlog; lhs=vote; start=logit; choices=LDPldp, OTHoth, LDPoth, OTHldp; tree=election[LDP(LDPldp, LDPoth), OTH(OTHldp, OTHoth)]; model: u(ldpldp, othoth, ldpoth, othldp)= + bDc1*Dc1 + bDc2*Dc2 + bDiss6MG4*Diss6MG4/ u(ldp, oth)= + bRp1*Rp1 + bRp2*Rp2 + bRiss6MG4*Riss6MG4$ /* Adj. R2: .37987 */ nlog; lhs=vote; start=logit; choices=LDPldp, OTHoth, LDPoth, OTHldp; tree=election[LDP(LDPldp, LDPoth), OTH(OTHldp, OTHoth)]; model: u(ldpldp, othoth, ldpoth, othldp)= + bDc1*Dc1 + bDc2*Dc2 + bDiss6MG5*Diss6MG5/ u(ldp, oth)= + bRp1*Rp1 + bRp2*Rp2 + bRiss6MG5*Riss6MG5$ /* Adj. R2: .37981 */ nlog; lhs=vote; start=logit; choices=LDPldp, OTHoth, LDPoth, OTHldp; tree=election[LDP(LDPldp, LDPoth), OTH(OTHldp, OTHoth)]; model: u(ldpldp, othoth, ldpoth, othldp)= + bDc1*Dc1 + bDc2*Dc2 + bDiss6MG6*Diss6MG6/ u(ldp, oth)= + bRp1*Rp1 + bRp2*Rp2 + bRiss6MG6*Riss6MG6$ /* Adj. R2: .37970 */ nlog; lhs=vote; start=logit; choices=LDPldp, OTHoth, LDPoth, OTHldp; tree=election[LDP(LDPldp, LDPoth), OTH(OTHldp, OTHoth)]; model: u(ldpldp, othoth, ldpoth, othldp)= + bDc1*Dc1 + bDc2*Dc2 + bDiss6MG7*Diss6MG7/ u(ldp, oth)= + bRp1*Rp1 + bRp2*Rp2 + bRiss6MG7*Riss6MG7$ /* Adj. R2: .37955 */ nlog; lhs=vote; start=logit; choices=LDPldp, OTHoth, LDPoth, OTHldp; tree=election[LDP(LDPldp, LDPoth), OTH(OTHldp, OTHoth)]; model: u(ldpldp, othoth, ldpoth, othldp)= + bDc1*Dc1 + bDc2*Dc2 + bDiss6MG8*Diss6MG8/ u(ldp, oth)= + bRp1*Rp1 + bRp2*Rp2 + bRiss6MG8*Riss6MG8$ /* Adj. R2: .37937 */ nlog; lhs=vote; start=logit; choices=LDPldp, OTHoth, LDPoth, OTHldp; tree=election[LDP(LDPldp, LDPoth), OTH(OTHldp, OTHoth)]; model: u(ldpldp, othoth, ldpoth, othldp)= + bDc1*Dc1 + bDc2*Dc2 + bDiss6MG9*Diss6MG9/ u(ldp, oth)= + bRp1*Rp1 + bRp2*Rp2 + bRiss6MG9*Riss6MG9$ /* Adj. R2: .37916 */ /* ************************ testing for beta between 2 and 4 ********************************** */ create; bet=.22$ create; Riss6M22=(2*(1-bet)*Riss6D) - (bet*Riss6P)$ create; Diss6M22=(2*(1-bet)*Diss6D) - (bet*Diss6P)$ create; bet=.24$ create; Riss6M24=(2*(1-bet)*Riss6D) - (bet*Riss6P)$ create; Diss6M24=(2*(1-bet)*Diss6D) - (bet*Diss6P)$ create; bet=.26$ create; Riss6M26=(2*(1-bet)*Riss6D) - (bet*Riss6P)$ create; Diss6M26=(2*(1-bet)*Diss6D) - (bet*Diss6P)$ create; bet=.28$ create; Riss6M28=(2*(1-bet)*Riss6D) - (bet*Riss6P)$ create; Diss6M28=(2*(1-bet)*Diss6D) - (bet*Diss6P)$ create; bet=.32$ create; Riss6M32=(2*(1-bet)*Riss6D) - (bet*Riss6P)$ create; Diss6M32=(2*(1-bet)*Diss6D) - (bet*Diss6P)$ create; bet=.34$ create; Riss6M34=(2*(1-bet)*Riss6D) - (bet*Riss6P)$ create; Diss6M34=(2*(1-bet)*Diss6D) - (bet*Diss6P)$ create; bet=.36$ create; Riss6M36=(2*(1-bet)*Riss6D) - (bet*Riss6P)$ create; Diss6M36=(2*(1-bet)*Diss6D) - (bet*Diss6P)$ create; bet=.38$ create; Riss6M38=(2*(1-bet)*Riss6D) - (bet*Riss6P)$ create; Diss6M38=(2*(1-bet)*Diss6D) - (bet*Diss6P)$ nlog; lhs=vote; start=logit; choices=LDPldp, OTHoth, LDPoth, OTHldp; tree=election[LDP(LDPldp, LDPoth), OTH(OTHldp, OTHoth)]; model: u(ldpldp, othoth, ldpoth, othldp)= + bDc1*Dc1 + bDc2*Dc2 + bDis6M22*Diss6M22/ u(ldp, oth)= + bRp1*Rp1 + bRp2*Rp2 + bRis6M22*Riss6M22$ /* Adj. R2: .37984 */ nlog; lhs=vote; start=logit; choices=LDPldp, OTHoth, LDPoth, OTHldp; tree=election[LDP(LDPldp, LDPoth), OTH(OTHldp, OTHoth)]; model: u(ldpldp, othoth, ldpoth, othldp)= + bDc1*Dc1 + bDc2*Dc2 + bDis6M24*Diss6M24/ u(ldp, oth)= + bRp1*Rp1 + bRp2*Rp2 + bRis6M24*Riss6M24$ /* Adj. R2: .37985 */ nlog; lhs=vote; start=logit; choices=LDPldp, OTHoth, LDPoth, OTHldp; tree=election[LDP(LDPldp, LDPoth), OTH(OTHldp, OTHoth)]; model: u(ldpldp, othoth, ldpoth, othldp)= + bDc1*Dc1 + bDc2*Dc2 + bDis6M26*Diss6M26/ u(ldp, oth)= + bRp1*Rp1 + bRp2*Rp2 + bRis6M26*Riss6M26$ /* Adj. R2: .37986 */ nlog; lhs=vote; start=logit; choices=LDPldp, OTHoth, LDPoth, OTHldp; tree=election[LDP(LDPldp, LDPoth), OTH(OTHldp, OTHoth)]; model: u(ldpldp, othoth, ldpoth, othldp)= + bDc1*Dc1 + bDc2*Dc2 + bDis6M28*Diss6M28/ u(ldp, oth)= + bRp1*Rp1 + bRp2*Rp2 + bRis6M28*Riss6M28$ /* Adj. R2: .37987 */ nlog; lhs=vote; start=logit; choices=LDPldp, OTHoth, LDPoth, OTHldp; tree=election[LDP(LDPldp, LDPoth), OTH(OTHldp, OTHoth)]; model: u(ldpldp, othoth, ldpoth, othldp)= + bDc1*Dc1 + bDc2*Dc2 + bDis6M32*Diss6M32/ u(ldp, oth)= + bRp1*Rp1 + bRp2*Rp2 + bRis6M32*Riss6M32$ /* Adj. R2: .37988 */ nlog; lhs=vote; start=logit; choices=LDPldp, OTHoth, LDPoth, OTHldp; tree=election[LDP(LDPldp, LDPoth), OTH(OTHldp, OTHoth)]; model: u(ldpldp, othoth, ldpoth, othldp)= + bDc1*Dc1 + bDc2*Dc2 + bDis6M34*Diss6M34/ u(ldp, oth)= + bRp1*Rp1 + bRp2*Rp2 + bRis6M34*Riss6M34$ /* Adj. R2: .37988 */ nlog; lhs=vote; start=logit; choices=LDPldp, OTHoth, LDPoth, OTHldp; tree=election[LDP(LDPldp, LDPoth), OTH(OTHldp, OTHoth)]; model: u(ldpldp, othoth, ldpoth, othldp)= + bDc1*Dc1 + bDc2*Dc2 + bDis6M36*Diss6M36/ u(ldp, oth)= + bRp1*Rp1 + bRp2*Rp2 + bRis6M36*Riss6M36$ /* Adj. R2: .37988 */ nlog; lhs=vote; start=logit; choices=LDPldp, OTHoth, LDPoth, OTHldp; tree=election[LDP(LDPldp, LDPoth), OTH(OTHldp, OTHoth)]; model: u(ldpldp, othoth, ldpoth, othldp)= + bDc1*Dc1 + bDc2*Dc2 + bDis6M38*Diss6M38/ u(ldp, oth)= + bRp1*Rp1 + bRp2*Rp2 + bRis6M38*Riss6M38$ /* Adj. R2: .37988 beta at .3-.39, adj r2s level off. Thus let us use beta=.35 */ /* ********** /* ********** The final model adopted for the paper /* ********** */ nlog; lhs=vote; start=logit; choices=LDPldp, OTHoth, LDPoth, OTHldp; tree=election[LDP(LDPldp, LDPoth), OTH(OTHldp, OTHoth)]; model: u(ldpldp)= B1 + b1Dc1*Dc1 + b1DcNis*DcNis + b1Dp1*Dp1/ u(othoth)= B2 + b2Dc1*Dc1 + b2DcNis*DcNis + b2Dp1*Dp1/ u(ldpoth)= b3Dc1*Dc1 + b3DcNis*DcNis + b3Dp1*Dp1/ u(othldp)= b4Dc1*Dc1 + b4DcNis*DcNis + b4Dp1*Dp1/ u(ldp)= B5 + b5Rp1*Rp1 + b5Ris6M3*Riss6M35 + b5RpNis*RpNis + b5Rpmob*Rpmob/ u(oth)= b6Rp1*Rp1 + b6Ris6M3*Riss6M35 + b6RpNis*RpNis + b6Rpmob*Rpmob$ /* ********** Adj. R2: .43306 /* ********** */ /* ***************************************** Section 4. Introduction of PID-strength into the model. *************************************************** */ create; base=.0$ create; Rp1W0=Rp1*(base + (1-base)*pidSTRGx)$ create; Dp1W0=Dp1*(base + (1-base)*pidSTRGx)$ create; RissW0=Riss6M35*(base + (1-base)*(1-pidSTRGx))$ nlog; lhs=vote; start=logit; choices=LDPldp, OTHoth, LDPoth, OTHldp; tree=election[LDP(LDPldp, LDPoth), OTH(OTHldp, OTHoth)]; model: u(ldpldp)= B1 + b1Dc1*Dc1 + b1DcNis*DcNis + b1Dp1*Dp1W0/ u(othoth)= B2 + b2Dc1*Dc1 + b2DcNis*DcNis + b2Dp1*Dp1W0/ u(ldpoth)= b3Dc1*Dc1 + b3DcNis*DcNis + b3Dp1*Dp1W0/ u(othldp)= b4Dc1*Dc1 + b4DcNis*DcNis + b4Dp1*Dp1W0/ u(ldp)= B5 + b5Rp1*Rp1W0 + b5Ris6M3*RissW0 + b5RpNis*RpNis + b5Rpmob*Rpmob/ u(oth)= b6Rp1*Rp1W0 + b6Ris6M3*RissW0 + b6RpNis*RpNis + b6Rpmob*Rpmob$ /* Adj. R2: .42258. */ create; base=.1$ create; Rp1W0=Rp1*(base + (1-base)*pidSTRGx)$ create; Dp1W0=Dp1*(base + (1-base)*pidSTRGx)$ create; RissW0=Riss6M35*(base + (1-base)*(1-pidSTRGx))$ nlog; lhs=vote; start=logit; choices=LDPldp, OTHoth, LDPoth, OTHldp; tree=election[LDP(LDPldp, LDPoth), OTH(OTHldp, OTHoth)]; model: u(ldpldp)= B1 + b1Dc1*Dc1 + b1DcNis*DcNis + b1Dp1*Dp1W0/ u(othoth)= B2 + b2Dc1*Dc1 + b2DcNis*DcNis + b2Dp1*Dp1W0/ u(ldpoth)= b3Dc1*Dc1 + b3DcNis*DcNis + b3Dp1*Dp1W0/ u(othldp)= b4Dc1*Dc1 + b4DcNis*DcNis + b4Dp1*Dp1W0/ u(ldp)= B5 + b5Rp1*Rp1W0 + b5Ris6M3*RissW0 + b5RpNis*RpNis + b5Rpmob*Rpmob/ u(oth)= b6Rp1*Rp1W0 + b6Ris6M3*RissW0 + b6RpNis*RpNis + b6Rpmob*Rpmob$ /* Adj. R2: .42492. */ create; base=.2$ create; Rp1W0=Rp1*(base + (1-base)*pidSTRGx)$ create; Dp1W0=Dp1*(base + (1-base)*pidSTRGx)$ create; RissW0=Riss6M35*(base + (1-base)*(1-pidSTRGx))$ nlog; lhs=vote; start=logit; choices=LDPldp, OTHoth, LDPoth, OTHldp; tree=election[LDP(LDPldp, LDPoth), OTH(OTHldp, OTHoth)]; model: u(ldpldp)= B1 + b1Dc1*Dc1 + b1DcNis*DcNis + b1Dp1*Dp1W0/ u(othoth)= B2 + b2Dc1*Dc1 + b2DcNis*DcNis + b2Dp1*Dp1W0/ u(ldpoth)= b3Dc1*Dc1 + b3DcNis*DcNis + b3Dp1*Dp1W0/ u(othldp)= b4Dc1*Dc1 + b4DcNis*DcNis + b4Dp1*Dp1W0/ u(ldp)= B5 + b5Rp1*Rp1W0 + b5Ris6M3*RissW0 + b5RpNis*RpNis + b5Rpmob*Rpmob/ u(oth)= b6Rp1*Rp1W0 + b6Ris6M3*RissW0 + b6RpNis*RpNis + b6Rpmob*Rpmob$ /* Adj. R2: .42685. */ create; base=.3$ create; Rp1W0=Rp1*(base + (1-base)*pidSTRGx)$ create; Dp1W0=Dp1*(base + (1-base)*pidSTRGx)$ create; RissW0=Riss6M35*(base + (1-base)*(1-pidSTRGx))$ nlog; lhs=vote; start=logit; choices=LDPldp, OTHoth, LDPoth, OTHldp; tree=election[LDP(LDPldp, LDPoth), OTH(OTHldp, OTHoth)]; model: u(ldpldp)= B1 + b1Dc1*Dc1 + b1DcNis*DcNis + b1Dp1*Dp1W0/ u(othoth)= B2 + b2Dc1*Dc1 + b2DcNis*DcNis + b2Dp1*Dp1W0/ u(ldpoth)= b3Dc1*Dc1 + b3DcNis*DcNis + b3Dp1*Dp1W0/ u(othldp)= b4Dc1*Dc1 + b4DcNis*DcNis + b4Dp1*Dp1W0/ u(ldp)= B5 + b5Rp1*Rp1W0 + b5Ris6M3*RissW0 + b5RpNis*RpNis + b5Rpmob*Rpmob/ u(oth)= b6Rp1*Rp1W0 + b6Ris6M3*RissW0 + b6RpNis*RpNis + b6Rpmob*Rpmob$ /* Adj. R2: .42842. */ create; base=.4$ create; Rp1W0=Rp1*(base + (1-base)*pidSTRGx)$ create; Dp1W0=Dp1*(base + (1-base)*pidSTRGx)$ create; RissW0=Riss6M35*(base + (1-base)*(1-pidSTRGx))$ nlog; lhs=vote; start=logit; choices=LDPldp, OTHoth, LDPoth, OTHldp; tree=election[LDP(LDPldp, LDPoth), OTH(OTHldp, OTHoth)]; model: u(ldpldp)= B1 + b1Dc1*Dc1 + b1DcNis*DcNis + b1Dp1*Dp1W0/ u(othoth)= B2 + b2Dc1*Dc1 + b2DcNis*DcNis + b2Dp1*Dp1W0/ u(ldpoth)= b3Dc1*Dc1 + b3DcNis*DcNis + b3Dp1*Dp1W0/ u(othldp)= b4Dc1*Dc1 + b4DcNis*DcNis + b4Dp1*Dp1W0/ u(ldp)= B5 + b5Rp1*Rp1W0 + b5Ris6M3*RissW0 + b5RpNis*RpNis + b5Rpmob*Rpmob/ u(oth)= b6Rp1*Rp1W0 + b6Ris6M3*RissW0 + b6RpNis*RpNis + b6Rpmob*Rpmob$ /* Adj. R2: .42966. */ create; base=.5$ create; Rp1W0=Rp1*(base + (1-base)*pidSTRGx)$ create; Dp1W0=Dp1*(base + (1-base)*pidSTRGx)$ create; RissW0=Riss6M35*(base + (1-base)*(1-pidSTRGx))$ nlog; lhs=vote; start=logit; choices=LDPldp, OTHoth, LDPoth, OTHldp; tree=election[LDP(LDPldp, LDPoth), OTH(OTHldp, OTHoth)]; model: u(ldpldp)= B1 + b1Dc1*Dc1 + b1DcNis*DcNis + b1Dp1*Dp1W0/ u(othoth)= B2 + b2Dc1*Dc1 + b2DcNis*DcNis + b2Dp1*Dp1W0/ u(ldpoth)= b3Dc1*Dc1 + b3DcNis*DcNis + b3Dp1*Dp1W0/ u(othldp)= b4Dc1*Dc1 + b4DcNis*DcNis + b4Dp1*Dp1W0/ u(ldp)= B5 + b5Rp1*Rp1W0 + b5Ris6M3*RissW0 + b5RpNis*RpNis + b5Rpmob*Rpmob/ u(oth)= b6Rp1*Rp1W0 + b6Ris6M3*RissW0 + b6RpNis*RpNis + b6Rpmob*Rpmob$ /* Adj. R2: .43064. */ create; base=.6$ create; Rp1W0=Rp1*(base + (1-base)*pidSTRGx)$ create; Dp1W0=Dp1*(base + (1-base)*pidSTRGx)$ create; RissW0=Riss6M35*(base + (1-base)*(1-pidSTRGx))$ nlog; lhs=vote; start=logit; choices=LDPldp, OTHoth, LDPoth, OTHldp; tree=election[LDP(LDPldp, LDPoth), OTH(OTHldp, OTHoth)]; model: u(ldpldp)= B1 + b1Dc1*Dc1 + b1DcNis*DcNis + b1Dp1*Dp1W0/ u(othoth)= B2 + b2Dc1*Dc1 + b2DcNis*DcNis + b2Dp1*Dp1W0/ u(ldpoth)= b3Dc1*Dc1 + b3DcNis*DcNis + b3Dp1*Dp1W0/ u(othldp)= b4Dc1*Dc1 + b4DcNis*DcNis + b4Dp1*Dp1W0/ u(ldp)= B5 + b5Rp1*Rp1W0 + b5Ris6M3*RissW0 + b5RpNis*RpNis + b5Rpmob*Rpmob/ u(oth)= b6Rp1*Rp1W0 + b6Ris6M3*RissW0 + b6RpNis*RpNis + b6Rpmob*Rpmob$ /* Adj. R2: .43141. */ create; base=.7$ create; Rp1W0=Rp1*(base + (1-base)*pidSTRGx)$ create; Dp1W0=Dp1*(base + (1-base)*pidSTRGx)$ create; RissW0=Riss6M35*(base + (1-base)*(1-pidSTRGx))$ nlog; lhs=vote; start=logit; choices=LDPldp, OTHoth, LDPoth, OTHldp; tree=election[LDP(LDPldp, LDPoth), OTH(OTHldp, OTHoth)]; model: u(ldpldp)= B1 + b1Dc1*Dc1 + b1DcNis*DcNis + b1Dp1*Dp1W0/ u(othoth)= B2 + b2Dc1*Dc1 + b2DcNis*DcNis + b2Dp1*Dp1W0/ u(ldpoth)= b3Dc1*Dc1 + b3DcNis*DcNis + b3Dp1*Dp1W0/ u(othldp)= b4Dc1*Dc1 + b4DcNis*DcNis + b4Dp1*Dp1W0/ u(ldp)= B5 + b5Rp1*Rp1W0 + b5Ris6M3*RissW0 + b5RpNis*RpNis + b5Rpmob*Rpmob/ u(oth)= b6Rp1*Rp1W0 + b6Ris6M3*RissW0 + b6RpNis*RpNis + b6Rpmob*Rpmob$ /* Adj. R2: .43200. */ create; base=.8$ create; Rp1W0=Rp1*(base + (1-base)*pidSTRGx)$ create; Dp1W0=Dp1*(base + (1-base)*pidSTRGx)$ create; RissW0=Riss6M35*(base + (1-base)*(1-pidSTRGx))$ nlog; lhs=vote; start=logit; choices=LDPldp, OTHoth, LDPoth, OTHldp; tree=election[LDP(LDPldp, LDPoth), OTH(OTHldp, OTHoth)]; model: u(ldpldp)= B1 + b1Dc1*Dc1 + b1DcNis*DcNis + b1Dp1*Dp1W0/ u(othoth)= B2 + b2Dc1*Dc1 + b2DcNis*DcNis + b2Dp1*Dp1W0/ u(ldpoth)= b3Dc1*Dc1 + b3DcNis*DcNis + b3Dp1*Dp1W0/ u(othldp)= b4Dc1*Dc1 + b4DcNis*DcNis + b4Dp1*Dp1W0/ u(ldp)= B5 + b5Rp1*Rp1W0 + b5Ris6M3*RissW0 + b5RpNis*RpNis + b5Rpmob*Rpmob/ u(oth)= b6Rp1*Rp1W0 + b6Ris6M3*RissW0 + b6RpNis*RpNis + b6Rpmob*Rpmob$ /* Adj. R2: .43245. */ create; base=.9$ create; Rp1W0=Rp1*(base + (1-base)*pidSTRGx)$ create; Dp1W0=Dp1*(base + (1-base)*pidSTRGx)$ create; RissW0=Riss6M35*(base + (1-base)*(1-pidSTRGx))$ nlog; lhs=vote; start=logit; choices=LDPldp, OTHoth, LDPoth, OTHldp; tree=election[LDP(LDPldp, LDPoth), OTH(OTHldp, OTHoth)]; model: u(ldpldp)= B1 + b1Dc1*Dc1 + b1DcNis*DcNis + b1Dp1*Dp1W0/ u(othoth)= B2 + b2Dc1*Dc1 + b2DcNis*DcNis + b2Dp1*Dp1W0/ u(ldpoth)= b3Dc1*Dc1 + b3DcNis*DcNis + b3Dp1*Dp1W0/ u(othldp)= b4Dc1*Dc1 + b4DcNis*DcNis + b4Dp1*Dp1W0/ u(ldp)= B5 + b5Rp1*Rp1W0 + b5Ris6M3*RissW0 + b5RpNis*RpNis + b5Rpmob*Rpmob/ u(oth)= b6Rp1*Rp1W0 + b6Ris6M3*RissW0 + b6RpNis*RpNis + b6Rpmob*Rpmob$ /* Adj. R2: .43280. */ /* ***************************************** Section 5. Post simulation estimation B. Estimating change in prob(LDPoth)= exp(B'Xij)/sum[exp(B'Xij)] ******************************************** */ /* Effect of IV on the 4 choices. The first set is a template. ******************************************** */ dstat; Rhs= Dc1, DcNis, Dp1, Rp1, Riss6M35, RpNis, Rpmob$ Descriptive Statistics All results based on nonmissing observations. Variable Mean Std.Dev. Minimum Maximum Cases ------------------------------------------------------------------------------- DC1 .386750383 1.28359216 -.751230000 9.13725000 3652 DCNIS .168243954 .374131959 .000000000 1.00000000 3804 DP1 .591278113 1.30362859 -1.35509000 8.58405000 3730 RP1 .591278113 1.30362859 -1.35509000 8.58405000 3730 RISS6M35 -49.2693481 64.7875042 -135.000000 65.0000000 3804 RPNIS .280757098 .449428134 .000000000 1.00000000 3804 RPMOB .804416404 1.13353558 .000000000 7.00000000 3804 /* /* *** IV (max: xxx, min: nnn ***************************** */ create; IVmax=xxx$ create; IVmin=nnn$ create; sumEbXx= exp( 2.330112598 + DC1 * .6899220787 + DCNIS * .9310144323 + DP1 * .1103967457) +exp( .8519967372 + DC1 * .3021260969 + DCNIS * .9303025610 + DP1 * .2760329276) +exp(DC1 * .4379134483 + DCNIS * 1.708758736 + DP1 * .4864035791) +exp(DC1 * .7886561438 + DCNIS * 1.108511324 + DP1 * .3407188575)$ create; pLlx= exp( 2.330112598 + DC1 * .6899220787 + DCNIS * .9310144323 + DP1 * .1103967457)/sumEbXx$ create; pOox= exp( .8368348233 + DC1 * .2983179912 + DCNIS * .9303025610 + DP1 * .2760329276)/sumEbXx$ create; pLox= exp(DC1 * .4379134483 + DCNIS * 1.708758736 + DP1 * .4864035791)/sumEbXx$ create; pOlx= exp(DC1 * .7886561438 + DCNIS * 1.108511324 + DP1 * .3407188575)/sumEbXx$ create; sumEbXn= exp( 2.330112598 + DC1 * .6899220787 + DCNIS * .9310144323 + DP1 * .1103967457) +exp( .8368348233 + DC1 * .2983179912 + DCNIS * .9303025610 + DP1 * .2760329276) +exp(DC1 * .4379134483 + DCNIS * 1.708758736 + DP1 * .4864035791) +exp(DC1 * .7886561438 + DCNIS * 1.108511324 + DP1 * .3407188575)$ create; pLln= exp( 2.330112598 + DC1 * .6899220787 + DCNIS * .9310144323 + DP1 * .1103967457)/sumEbXn$ create; pOon= exp( .8368348233 + DC1 * .2983179912 + DCNIS * .9303025610 + DP1 * .2760329276)/sumEbXn$ create; pLon= exp(DC1 * .4379134483 + DCNIS * 1.708758736 + DP1 * .4864035791)/sumEbXn$ create; pOln= exp(DC1 * .7886561438 + DCNIS * 1.108511324 + DP1 * .3407188575)/sumEbXn$ dstat; Rhs= pLlx, pLln, pOox, pOon, pLox, pLon, pOlx, pOln$ /* ***1. Dc1 (max: 9.14, min: -.751 ***************************** */ create; IVmax=9.14$ create; IVmin=-.751$ create; sumEbXx= exp( 2.330112598 + IVmax * .6899220787 + DCNIS * .9310144323 + DP1 * .1103967457) +exp( .8368348233 + IVmax * .2983179912 + DCNIS * .9303025610 + DP1 * .2760329276) +exp(IVmax * .4379134483 + DCNIS * 1.708758736 + DP1 * .4864035791) +exp(IVmax * .7886561438 + DCNIS * 1.108511324 + DP1 * .3407188575)$ create; pLlx= exp( 2.330112598 + IVmax * .6899220787 + DCNIS * .9310144323 + DP1 * .1103967457)/sumEbXx$ create; pOox= exp( .8368348233 + IVmax * .2983179912 + DCNIS * .9303025610 + DP1 * .2760329276)/sumEbXx$ create; pLox= exp(IVmax * .4379134483 + DCNIS * 1.708758736 + DP1 * .4864035791)/sumEbXx$ create; pOlx= exp(IVmax * .7886561438 + DCNIS * 1.108511324 + DP1 * .3407188575)/sumEbXx$ create; sumEbXn= exp( 2.330112598 + IVmin * .6899220787 + DCNIS * .9310144323 + DP1 * .1103967457) +exp( .8368348233 + IVmin * .2983179912 + DCNIS * .9303025610 + DP1 * .2760329276) +exp(IVmin * .4379134483 + DCNIS * 1.708758736 + DP1 * .4864035791) +exp(IVmin * .7886561438 + DCNIS * 1.108511324 + DP1 * .3407188575)$ create; pLln= exp( 2.330112598 + IVmin * .6899220787 + DCNIS * .9310144323 + DP1 * .1103967457)/sumEbXn$ create; pOon= exp( .8368348233 + IVmin * .2983179912 + DCNIS * .9303025610 + DP1 * .2760329276)/sumEbXn$ create; pLon= exp(IVmin * .4379134483 + DCNIS * 1.708758736 + DP1 * .4864035791)/sumEbXn$ create; pOln= exp(IVmin * .7886561438 + DCNIS * 1.108511324 + DP1 * .3407188575)/sumEbXn$ dstat; Rhs= pLlx, pLln, pOox, pOon, pLox, pLon, pOlx, pOln$ /* ***2. DcNis (max: 1, min: 0 ***************************** */ create; IVmax=1$ create; IVmin=0$ create; sumEbXx= exp( 2.330112598 + DC1 * .6899220787 + IVmax * .9310144323 + DP1 * .1103967457) +exp( .8368348233 + DC1 * .2983179912 + IVmax * .9303025610 + DP1 * .2760329276) +exp(DC1 * .4379134483 + IVmax * 1.708758736 + DP1 * .4864035791) +exp(DC1 * .7886561438 + IVmax * 1.108511324 + DP1 * .3407188575)$ create; pLlx= exp( 2.330112598 + DC1 * .6899220787 + IVmax * .9310144323 + DP1 * .1103967457)/sumEbXx$ create; pOox= exp( .8368348233 + DC1 * .2983179912 + IVmax * .9303025610 + DP1 * .2760329276)/sumEbXx$ create; pLox= exp(DC1 * .4379134483 + IVmax * 1.708758736 + DP1 * .4864035791)/sumEbXx$ create; pOlx= exp(DC1 * .7886561438 + IVmax * 1.108511324 + DP1 * .3407188575)/sumEbXx$ create; sumEbXn= exp( 2.330112598 + DC1 * .6899220787 + IVmin * .9310144323 + DP1 * .1103967457) +exp( .8368348233 + DC1 * .2983179912 + IVmin * .9303025610 + DP1 * .2760329276) +exp(DC1 * .4379134483 + IVmin * 1.708758736 + DP1 * .4864035791) +exp(DC1 * .7886561438 + IVmin * 1.108511324 + DP1 * .3407188575)$ create; pLln= exp( 2.330112598 + DC1 * .6899220787 + IVmin * .9310144323 + DP1 * .1103967457)/sumEbXn$ create; pOon= exp( .8368348233 + DC1 * .2983179912 + IVmin * .9303025610 + DP1 * .2760329276)/sumEbXn$ create; pLon= exp(DC1 * .4379134483 + IVmin * 1.708758736 + DP1 * .4864035791)/sumEbXn$ create; pOln= exp(DC1 * .7886561438 + IVmin * 1.108511324 + DP1 * .3407188575)/sumEbXn$ dstat; Rhs= pLlx, pLln, pOox, pOon, pLox, pLon, pOlx, pOln$ /* ***3. Dp1 (max: 8.58, min: -1.36 ***************************** */ create; IVmax=8.58$ create; IVmin=-1.36$ create; sumEbXx= exp( 2.330112598 + DC1 * .6899220787 + DCNIS * .9310144323 + IVmax * .1103967457) +exp( .8368348233 + DC1 * .2983179912 + DCNIS * .9303025610 + IVmax * .2760329276) +exp(DC1 * .4379134483 + DCNIS * 1.708758736 + IVmax * .4864035791) +exp(DC1 * .7886561438 + DCNIS * 1.108511324 + IVmax * .3407188575)$ create; pLlx= exp( 2.330112598 + DC1 * .6899220787 + DCNIS * .9310144323 + IVmax * .1103967457)/sumEbXx$ create; pOox= exp( .8368348233 + DC1 * .2983179912 + DCNIS * .9303025610 + IVmax * .2760329276)/sumEbXx$ create; pLox= exp(DC1 * .4379134483 + DCNIS * 1.708758736 + IVmax * .4864035791)/sumEbXx$ create; pOlx= exp(DC1 * .7886561438 + DCNIS * 1.108511324 + IVmax * .3407188575)/sumEbXx$ create; sumEbXn= exp( 2.330112598 + DC1 * .6899220787 + DCNIS * .9310144323 + IVmin * .1103967457) +exp( .8368348233 + DC1 * .2983179912 + DCNIS * .9303025610 + IVmin * .2760329276) +exp(DC1 * .4379134483 + DCNIS * 1.708758736 + DP1 * .4864035791) +exp(DC1 * .7886561438 + DCNIS * 1.108511324 + IVmin * .3407188575)$ create; pLln= exp( 2.330112598 + DC1 * .6899220787 + DCNIS * .9310144323 + IVmin * .1103967457)/sumEbXn$ create; pOon= exp( .8368348233 + DC1 * .2983179912 + DCNIS * .9303025610 + IVmin * .2760329276)/sumEbXn$ create; pLon= exp(DC1 * .4379134483 + DCNIS * 1.708758736 + IVmin * .4864035791)/sumEbXn$ create; pOln= exp(DC1 * .7886561438 + DCNIS * 1.108511324 + IVmin * .3407188575)/sumEbXn$ dstat; Rhs= pLlx, pLln, pOox, pOon, pLox, pLon, pOlx, pOln$ /* ************************* For PR equations ************************** */ /* ***4. Rp1 (max: 8.58, min: -1.36 ***************************** */ create; IVmax=8.58$ create; IVmin=-1.36$ create; sumEbXx= exp( -1.839883441 + IVmax * .7227771089 + RISs6M35 * .007107381046 + RPNIS * 1.044896618 + RPMOB * .2721177448) +exp( IVmax * .3862251208 + RISs6M35 * .007856606591 + RPNIS * .7155853178 + RPMOB * .2563995165)$ create; pLx= exp( -1.839883441 + IVmax * .7227771089 + RISs6M35 * .007107381046 + RPNIS * 1.044896618 + RPMOB * .2721177448)/sumEbXx$ create; pOx= +exp( IVmax * .3862251208 + RISs6M35 * .007856606591 + RPNIS * .7155853178 + RPMOB * .2563995165)/sumEbXx$ create; sumEbXn= exp( -1.839883441 + IVmin * .7227771089 + RISs6M35 * .007107381046 + RPNIS * 1.044896618 + RPMOB * .2721177448) +exp( IVmin * .3862251208 + RISs6M35 * .007856606591 + RPNIS * .7155853178 + RPMOB * .2563995165)$ create; pLn= exp( -1.839883441 + IVmin * .7227771089 + RISs6M35 * .007107381046 + RPNIS * 1.044896618 + RPMOB * .2721177448)/sumEbXn$ create; pOn= +exp( IVmin * .3862251208 + RISs6M35 * .007856606591 + RPNIS * .7155853178 + RPMOB * .2563995165)/sumEbXn$ dstat; Rhs= pLx, pLn, pOx, pOn$ /* ***5. Iss (max: 65, min: -135 ***************************** */ create; IVmax=65$ create; IVmin=-135$ create; sumEbXx= exp( -1.839883441 + RP1 * .7227771089 + IVmax * .007107381046 + RPNIS * 1.044896618 + RPMOB * .2721177448) +exp( RP1 * .3862251208 + IVmax * .007856606591 + RPNIS * .7155853178 + RPMOB * .2563995165)$ create; pLx= exp( -1.839883441 + RP1 * .7227771089 + IVmax * .007107381046 + RPNIS * 1.044896618 + RPMOB * .2721177448)/sumEbXx$ create; pOx= +exp( RP1 * .3862251208 + IVmax * .007856606591 + RPNIS * .7155853178 + RPMOB * .2563995165)/sumEbXx$ create; sumEbXn= exp( -1.839883441 + RP1 * .7227771089 + IVmin * .007107381046 + RPNIS * 1.044896618 + RPMOB * .2721177448) +exp( RP1 * .3862251208 + IVmin * .007856606591 + RPNIS * .7155853178 + RPMOB * .2563995165)$ create; pLn= exp( -1.839883441 + RP1 * .7227771089 + IVmin * .007107381046 + RPNIS * 1.044896618 + RPMOB * .2721177448)/sumEbXn$ create; pOn= +exp( RP1 * .3862251208 + IVmin * .007856606591 + RPNIS * .7155853178 + RPMOB * .2563995165)/sumEbXn$ dstat; Rhs= pLx, pLn, pOx, pOn$ /* ***6. Rpnis (max: 1, min: 0 ***************************** */ create; IVmax=1$ create; IVmin=0$ create; sumEbXx= exp( -1.839883441 + RP1 * .7227771089 + RISs6M35 * .007107381046 + IVmax * 1.044896618 + RPMOB * .2721177448) +exp( RP1 * .3862251208 + RISs6M35 * .007856606591 + IVmax * .7155853178 + RPMOB * .2563995165)$ create; pLx= exp( -1.839883441 + RP1 * .7227771089 + RISs6M35 * .007107381046 + IVmax * 1.044896618 + RPMOB * .2721177448)/sumEbXx$ create; pOx= +exp( RP1 * .3862251208 + RISs6M35 * .007856606591 + IVmax * .7155853178 + RPMOB * .2563995165)/sumEbXx$ create; sumEbXn= exp( -1.839883441 + RP1 * .7227771089 + RISs6M35 * .007107381046 + IVmin * 1.044896618 + RPMOB * .2721177448) +exp( RP1 * .3862251208 + RISs6M35 * .007856606591 + IVmin * .7155853178 + RPMOB * .2563995165)$ create; pLn= exp( -1.839883441 + RP1 * .7227771089 + RISs6M35 * .007107381046 + IVmin * 1.044896618 + RPMOB * .2721177448)/sumEbXn$ create; pOn= +exp( RP1 * .3862251208 + RISs6M35 * .007856606591 + IVmin * .7155853178 + RPMOB * .2563995165)/sumEbXn$ dstat; Rhs= pLx, pLn, pOx, pOn$ /* ***7. Rp1 (max: 7, min: 0 ***************************** */ create; IVmax=7$ create; IVmin=0$ create; sumEbXx= exp( -1.839883441 + RP1 * .7227771089 + RISs6M35 * .007107381046 + RPNIS * 1.044896618 + IVmax * .2721177448) +exp( RP1 * .3862251208 + RISs6M35 * .007856606591 + RPNIS * .7155853178 + IVmax * .2563995165)$ create; pLx= exp( -1.839883441 + RP1 * .7227771089 + RISs6M35 * .007107381046 + RPNIS * 1.044896618 + IVmax * .2721177448)/sumEbXx$ create; pOx= +exp( RP1 * .3862251208 + RISs6M35 * .007856606591 + RPNIS * .7155853178 + IVmax * .2563995165)/sumEbXx$ create; sumEbXn= exp( -1.839883441 + RP1 * .7227771089 + RISs6M35 * .007107381046 + RPNIS * 1.044896618 + IVmin * .2721177448) +exp( RP1 * .3862251208 + RISs6M35 * .007856606591 + RPNIS * .7155853178 + IVmin * .2563995165)$ create; pLn= exp( -1.839883441 + RP1 * .7227771089 + RISs6M35 * .007107381046 + RPNIS * 1.044896618 + IVmin * .2721177448)/sumEbXn$ create; pOn= +exp( RP1 * .3862251208 + RISs6M35 * .007856606591 + RPNIS * .7155853178 + IVmin * .2563995165)/sumEbXn$ dstat; Rhs= pLx, pLn, pOx, pOn$