Er der nogen der ved det og jeg mangler Bartlett test. Nogen der overhovedet kan SAS herinde? Det er et program til databehandling af statistiske data.
Næ, men google kan bruges til mange ting; søg på bartlett test sas , jeg fik f.eks. følgende %macro Bartlett(data=_last_,var=,gruppe=,bygruppe=_ingen,tabel=ja);
/****************************************************************/ /* */ /* Laver Bartletts test for varianshomogenitet for variablen */ /* "var" over grupperne givet ved klassifikationsvariablen */ /* "gruppe" for hvert niveau af "bygruppe"-variablen. */ /* */ /* Grupper med ens observationer eller med et antal ikke- */ /* manglende observationer mindre end 2 ekskluderes fra testet */ /* (men relevant statistik beregnes alligevel). */ /* */ /* Af hensyn til plads i udskriften tilrådes det, at "gruppe"- */ /* variablen enten er numerisk eller højst består af 2 tegn. */ /* */ /* Datasettet "_bartlet" indeholder næsten alle de udskrevne */ /* værdier. */ /* */ /* Reference: */ /* Preben Blaesild & Jqrgen Granfeldt: */ /* Noter til Biostatistik og Geostatistik, */ /* 1994 */ /* */ /* Jørn Attermann 19-01-94, 14-04-94, 08-06-94 */ /* */ /* Modificeret af Asger Roer Pedersen 21-02-96. */ /* */ /****************************************************************/
data _mact1; set &data; _ingen=1; run;
proc sort data=_mact1; by &bygruppe &gruppe; run;
proc summary data=_mact1; by &bygruppe &gruppe; var &var; output out=_mact2 n=_ni sum=_Si USS=_USSi mean=_meani; run;
data _bartlet; set _mact2; by &bygruppe; retain _n _S _USS _CT _SSD _f _s_invfi _s_flns _k . _ekskl 0; if first.&bygruppe then do; _n=.; _S=.; _USS=.; _CT=.; _SSD=.; _f=.; _s_invfi=.; _s_flns=.; _k=.; _ekskl=0; end; if (_ni>0) then _CTi=_Si**2/_ni; else _CTi=.; _SSDi=_USSi-_CTi; _fi=_ni-1; if ((not (_fi<1)) and (_SSDi>0)) then do; _vari=_SSDi/_fi; _s_invfi+1/_fi; _s_flns+_fi*log(_vari); _k+1; end; else do; _ekskl+1; _fi=0; _vari=.; end; _n+_ni; _S+_Si; _USS+_USSi; _CT+_CTi; _SSD+_SSDi; _f+_fi; if last.&bygruppe then do; _var=_SSD/_f; if (_k>1) then do; _C=1+1/3/(_k-1)*(_s_invfi-1/_f); _2lnQ=_f*log(_var)-_s_flns; _Ba=_2lnQ/_C; _testss=1-probchi(_Ba,_k-1); end; end;
/* put / @33 'Macro "BARTLETT"'; put / @33 "Dataset: %upcase(&data)"; put @25 "Responsvariabel: %upcase(&var)"; put @26 "Gruppevariabel: %upcase(&gruppe)"; */ /* put @24 "Bygruppevariabel: %upcase(&bygruppe)"; */ /* put // @33 "&bygruppe = " @44 _bychar $28.; put ' ____________________________________________________'; %end; */
%if &tabel=ja %then %do; file print; if (_n_=1) then do;
put / @10 'Macro "BARTLETT"' " Dataset: %upcase(&data)"; put @14 "Responsvariabel : %upcase(&var) Gruppevariabel : %upcase(&gruppe)";
Tilladte BB-code-tags: [b]fed[/b] [i]kursiv[/i] [u]understreget[/u] Web- og emailadresser omdannes automatisk til links. Der sættes "nofollow" på alle links.