能不能幫忙看一下code 哪里有問(wèn)題，run不出結(jié)果來(lái)

是 否 +2 論壇幣

k人 參與回答

經(jīng)管之家送您一份

應(yīng)屆畢業(yè)生專屬福利!

求職就業(yè)群

趙安豆老師微信：zhaoandou666

經(jīng)管之家聯(lián)合CDA

送您一個(gè)全額獎(jiǎng)學(xué)金名額~ !

立即領(lǐng)取

感謝您參與論壇問(wèn)題回答

經(jīng)管之家送您兩個(gè)論壇幣！

+2 論壇幣

title ’Choice of Chocolate Candies’;

data chocs;

input Subj c Dark Soft Nuts @@;

Set = 1;

datalines;

1 2 0 0 0 1 2 0 0 1 1 2 0 1 0 1 2 0 1 1

1 1 1 0 0 1 2 1 0 1 1 2 1 1 0 1 2 1 1 1

2 2 0 0 0 2 2 0 0 1 2 2 0 1 0 2 2 0 1 1

2 2 1 0 0 2 1 1 0 1 2 2 1 1 0 2 2 1 1 1

3 2 0 0 0 3 2 0 0 1 3 2 0 1 0 3 2 0 1 1

3 2 1 0 0 3 2 1 0 1 3 1 1 1 0 3 2 1 1 1

4 2 0 0 0 4 2 0 0 1 4 2 0 1 0 4 2 0 1 1

4 1 1 0 0 4 2 1 0 1 4 2 1 1 0 4 2 1 1 1

5 2 0 0 0 5 1 0 0 1 5 2 0 1 0 5 2 0 1 1

5 2 1 0 0 5 2 1 0 1 5 2 1 1 0 5 2 1 1 1

6 2 0 0 0 6 2 0 0 1 6 2 0 1 0 6 2 0 1 1

6 2 1 0 0 6 1 1 0 1 6 2 1 1 0 6 2 1 1 1

7 2 0 0 0 7 1 0 0 1 7 2 0 1 0 7 2 0 1 1

7 2 1 0 0 7 2 1 0 1 7 2 1 1 0 7 2 1 1 1

8 2 0 0 0 8 2 0 0 1 8 2 0 1 0 8 2 0 1 1

8 2 1 0 0 8 1 1 0 1 8 2 1 1 0 8 2 1 1 1

9 2 0 0 0 9 2 0 0 1 9 2 0 1 0 9 2 0 1 1

9 2 1 0 0 9 1 1 0 1 9 2 1 1 0 9 2 1 1 1

10 2 0 0 0 10 2 0 0 1 10 2 0 1 0 10 2 0 1 1

10 2 1 0 0 10 1 1 0 1 10 2 1 1 0 10 2 1 1 1

;

proc print data=chocs noobs;

where subj <= 2;

var subj set c dark soft nuts;

run;

* Alternative Form of Data Entry;

data combos; /* Read the design matrix. */

input Dark Soft Nuts;

datalines;

0 0 0

0 0 1

0 1 0

0 1 1

1 0 0

1 0 1

1 1 0

1 1 1

;

data chocs; /* Create the data set. */

input Choice @@; drop choice; /* Read the chosen combo num. */

Subj = _n_; Set = 1; /* Store subj, choice set num. */

do i = 1 to 8; /* Loop over alternatives. */

c = 2 - (i eq choice); /* Designate chosen alt. */

set combos point=i; /* Read design matrix. */

output; /* Output the results. */

end;

datalines;

5 6 7 5 2 6 2 6 6 6

;

proc phreg data=chocs outest=betas;

strata subj set;

model c*c(2) = dark soft nuts / ties=breslow;

label dark = ’Dark Chocolate’ soft = ’Soft Center’

nuts = ’With Nuts’;

run;

data chocs2;

set chocs;

Milk = 1 - dark; Chewy = 1 - Soft; NoNuts = 1 - nuts;

label dark = ’Dark Chocolate’ milk = ’Milk Chocolate’

soft = ’Soft Center’ chewy = ’Chewy Center’

nuts = ’With Nuts’ nonuts = ’No Nuts’;

run;

proc phreg data=chocs2;

strata subj set;

model c*c(2) = dark milk soft chewy nuts nonuts / ties=breslow;

run;

* Estimate the probability that each alternative will be chosen; （在這之前的應(yīng)該問(wèn)題不是很大，就是在求probability中出不來(lái)結(jié)果）

data p;

retain sum 0;

set combos end=eof;

* On the first pass through the DATA step (_n_ is the pass

number), get the regression coefficients in B1-B3.

Note that they are automatically retained so that they

can be used in all passes through the DATA step.;

if _n_ = 1 then

set betas(rename=(dark=b1 soft=b2 nuts=b3));

keep dark soft nuts p;

array x[3] dark soft nuts;

array b[3] b1-b3;

* For each combination, create x * b;

p = 0;

do j = 1 to 3;

p = p + x[j] * b[j];

end;

* Exponentiate x * b and sum them up;

p = exp(p);

sum = sum + p;

* Output sum exp(x * b) in the macro variable ’&sum’;

if eof then call symput(’sum’,put(sum,best12.));

run;

proc format;

value df 1 = ’Dark’ 0 = ’Milk’;

value sf 1 = ’Soft’ 0 = ’Chewy’;

value nf 1 = ’Nuts’ 0 = ’No Nuts’;

run;

* Divide each exp(x * b) by sum exp(x * b);

data p;

set p;

p = p / (&sum);

format dark df. soft sf. nuts nf.;

run;

proc sort;

by descending p;

run;

proc print;

run;

紅字的部分是有問(wèn)題的部分，怎樣修改呢？ 謝謝

[此貼子已經(jīng)被作者于2008-12-4 14:05:30編輯過(guò)]

掃碼加我 拉你入群

請(qǐng)注明：姓名-公司-職位

以便審核進(jìn)群資格，未注明則拒絕

分享0 收藏0 回帖

關(guān)鍵詞：code ODE COD Run 哪里有 結(jié)果 幫忙 code Run

p = p / (∑); 應(yīng)該是p=p/(&sum)      &sum貼不出來(lái)

[此貼子已經(jīng)被作者于2008-12-4 14:08:59編輯過(guò)]

樓主發(fā)完題目就跑了 這種春天播種秋天收獲的模式需要改進(jìn)

以下是引用sushe1527在2008-12-4 14:44:00的發(fā)言：
樓主發(fā)完題目就跑了 這種春天播種秋天收獲的模式需要改進(jìn)

那應(yīng)該是哪種模式呢？

以下是引用chasedreamtrue在2008-12-4 22:50:00的發(fā)言：

哈哈 你來(lái)得滿早 

[此貼子已經(jīng)被作者于2008-12-5 1:13:17編輯過(guò)]

以下是引用sushe1527在2008-12-4 22:51:00的發(fā)言：

哈哈 你來(lái)得滿早

起個(gè)大早看到有人回帖以為問(wèn)題解決了呢，原來(lái)不沾邊

以下是引用chasedreamtrue在2008-12-4 22:58:00的發(fā)言：

有個(gè) label 標(biāo)錯(cuò)了 sum沒(méi)有賦值  也迷糊了 是不是這個(gè)��？

	Dark	Soft	Nuts	p
1	1	0	1	0.5039995573
2	1	0	0	0.2159998103
3	0	0	1	0.1259998894
4	1	1	1	0.0560004427
5	0	0	0	0.0539999526
6	1	1	0	0.0240001897
7	0	1	1	0.0140001107
8	0	1	0	0.0060000474

以下是引用sushe1527在2008-12-4 23:48:00的發(fā)言：

有個(gè) label 標(biāo)錯(cuò)了 sum沒(méi)有賦值  也迷糊了 是不是這個(gè)啊？

	Dark	Soft	Nuts	p
1	1	0	1	0.5039995573
2	1	0	0	0.2159998103
3	0	0	1	0.1259998894
4	1	1	1	0.0560004427
5	0	0	0	0.0539999526
6	1	1	0	0.0240001897
7	0	1	1	0.0140001107
8	0	1	0	0.0060000474

是呀，但是怎莫改那個(gè)code，能過(guò)run出這個(gè)結(jié)果呢？

有個(gè)sum = sum + p，先求P ，然后sum 試試  

別的一些格式錯(cuò)誤不致命  呼呼去了 我這可不是大早

以下是引用sushe1527在2008-12-5 1:47:00的發(fā)言：

有個(gè)sum = sum + p，先求P ，然后sum 試試  

別的一些格式錯(cuò)誤不致命  呼呼去了 我這可不是大早

 code 是p = p / (&sum) ，但是一貼出來(lái)就是變了，在SAS里總是說(shuō)這個(gè)有問(wèn)題（&與sum相乘，如果兩個(gè)放到一起就成了∑這個(gè)符號(hào)了，其實(shí)code是&與sum相乘，不知道說(shuō)清楚沒(méi)有，反正這個(gè)貼出來(lái)的與我發(fā)的不一樣