빈도의 선 그래프와 SAS의 여러 질병에 걸린 발병 연령

Question

질병의 발병 연령과 관련된 환자 데이터 세트가 있으며 각 질병의 발병 연령 빈도를 다른 줄로 계획하고 싶습니다. x 축은 발병 연령이고 y 축은 빈도이며, 각 선은 다른 질병을 나타냅니다. 나이가 0 인 환자는 그 병이 없다는 것을 나타냅니다. 이렇게하려면 SAS에 어떤 코드가 있어야합니까? 정말 고마워!

HoHTAge HoGDAge AddDAge CelDAge 
0 0 32 0 
0 0 0 0 
12 0 23 0 
0 20 0 0 
25 0 0 0 
0 0 0 0 
32 0 0 0 
0 0 0 0 
0 0 0 0 
0 0 0 0 
0 0 0 0 
0 0 0 35 
45 0 0 0 
0 0 0 0 
0 0 0 0 
43 0 0 0 
0 23 0 0 
0 18 0 0 
0 0 0 0 
0 0 0 0 
0 0 0 0 
0 12 0 0 
30 26 0 0 
0 40 46 0 
0 0 0 30 
57 0 0 0 
0 0 0 0 

Answer 1

정확하게 문제를 이해하면 100 % 확신 할 수 없지만 해결 방법을 제공하려고 시도했습니다.

이것은 아마도 복잡한 해결책 일 것입니다. 나는 많은 짧은/쉬운 해결책이 있다고 생각합니다. 나는 각 질병에 대한 freqs를 계산, 하나 개의 데이터 세트로 병합하고 gplot로 그릴이이 그래프로 이어질 것입니다 샘플 데이터와

data x; 
input HoHTAge HoGDAge AddDAge CelDAge; 
datalines; 
0 0 32 0 
0 0 0 0 
12 0 23 0 
0 20 0 0 
25 0 0 0 
0 0 0 0 
32 0 0 0 
0 0 0 0 
0 0 0 0 
0 0 0 0 
0 0 0 0 
0 0 0 35 
45 0 0 0 
0 0 0 0 
0 0 0 0 
43 0 0 0 
0 23 0 0 
0 18 0 0 
0 0 0 0 
0 0 0 0 
0 0 0 0 
0 12 0 0 
30 26 0 0 
0 40 46 0 
0 0 0 30 
57 0 0 0 
0 0 0 0 
; 
run; 

proc freq data=x noprint ; 
tables HoHTAge/out=a; 

run; 
proc freq data=x noprint ; 
tables HoGDAge/out=b; 

run; 
proc freq data=x noprint ; 
tables AddDAge/out=c; 

run; 
proc freq data=x noprint ; 
tables CelDAge/out=d; 

run; 

data res (drop =percent count); 
merge a (in=a rename=(HoHTAge=age)) b (in=b rename=(HoGDAge=age)) c (in=c rename=(AddDAge=age)) d(in=d rename=(CelDAge=age)); 
by age; 
*if age=0 then count=0; /*if you want to exclude age 0*/ 
if a then HoHTAge=count; else HoHTAge=0; 
if b then HoGDAge=count; else HoGDAge=0; 
if c then AddDAge=count; else AddDAge=0; 
if d then CelDAge=count; else CelDAge=0; 

ruN; 
/* Set the graphics environment */                          
goptions reset=all cback=white border htext=10pt htitle=12pt; 

axis1 label=("age");               
axis2 label=("Count");                      

symbol1 interpol=join color=R height=14pt font='Arial' ;               
symbol2 interpol=join color=B height=14pt font='Arial';              
symbol3 interpol=join color=O height=14pt font='Arial' ;               
symbol4 interpol=join color=BL height=14pt font='Arial' ;                                                          
legend1 repeat=1 label=none frame;                     

proc gplot data=res;                            
    plot (HoHTAge HoGDAge AddDAge CelDAge)*age/ overlay legend=legend1 haxis=axis1 vaxis=axis2;                   
run;                                  

, 나는이 더 잘 보일 것 실제 데이터로 생각 때문에 지금 우리 내가 이해로까지 당신처럼 출력을 원

으로 간단한 대안은 당신이 PROC 주파수의 dotplot를 사용할 수 있지만 다음 그래프 만 점을 구분 한 각 질병에 대해 한 번 더 한 후 연령이 없습니다 긴 해결책 :

ods graphics on; 
proc freq data=x ; 
tables HoHTAge HoGDAge AddDAge CelDAge/plots=freqplot(type=dot ORIENT = VERTICAL); 
run;