美国社区调查，SAS EG 误差范围代码答案

【问题标题】：American Community Survey, SAS EG code for margin of error美国社区调查，SAS EG 误差范围代码
【发布时间】：2015-10-06 04:37:08
【问题描述】：

人口普查局提供了计算美国社区调查误差幅度的数学公式，但不包括其 SAS 代码。该公式位于此处文档的第 24 页：http://www2.census.gov/programs-surveys/acs/tech_docs/accuracy/ACS_Accuracy_of_Data_2014.pdf

有人知道误差范围的 SAS 代码吗？它必须包含所有 80 个 pwgtp。

【问题讨论】：

标签： sas census

【解决方案1】：

这里是相关代码。它使用 90% 的置信区间，因为这是人口普查局在 American FactFinder 上公布的误差范围所使用的。您可以在“1.64537”所在的开头更改置信区间。

/* Margin of Error 90% confidence code*/

        1.64537*(SQRT(.05*(SUM((SUM(t1.pwgtp1)-(SUM(t1.PWGTP)))**2,(SUM(t1.pwgtp2)-(SUM(t1.PWGTP)))**2,(SUM(t1.pwgtp3)-(SUM(t1.PWGTP)))**2,(SUM(t1.pwgtp4)-(SUM(t1.PWGTP)))**2,(SUM(t1.pwgtp5)-(SUM(t1.PWGTP)))**2
            ,(SUM(t1.pwgtp6)-(SUM(t1.PWGTP)))**2,(SUM(t1.pwgtp7)-(SUM(t1.PWGTP)))**2,(SUM(t1.pwgtp8)-(SUM(t1.PWGTP)))**2,(SUM(t1.pwgtp9)-(SUM(t1.PWGTP)))**2,(SUM(t1.pwgtp10)-(SUM(t1.PWGTP)))**2
            ,(SUM(t1.pwgtp11)-(SUM(t1.PWGTP)))**2,(SUM(t1.pwgtp12)-(SUM(t1.PWGTP)))**2,(SUM(t1.pwgtp13)-(SUM(t1.PWGTP)))**2,(SUM(t1.pwgtp14)-(SUM(t1.PWGTP)))**2,(SUM(t1.pwgtp15)-(SUM(t1.PWGTP)))**2
            ,(SUM(t1.pwgtp16)-(SUM(t1.PWGTP)))**2,(SUM(t1.pwgtp17)-(SUM(t1.PWGTP)))**2,(SUM(t1.pwgtp18)-(SUM(t1.PWGTP)))**2,(SUM(t1.pwgtp19)-(SUM(t1.PWGTP)))**2,(SUM(t1.pwgtp20)-(SUM(t1.PWGTP)))**2
        ,(SUM(t1.pwgtp21)-(SUM(t1.PWGTP)))**2,(SUM(t1.pwgtp22)-(SUM(t1.PWGTP)))**2,(SUM(t1.pwgtp23)-(SUM(t1.PWGTP)))**2,(SUM(t1.pwgtp24)-(SUM(t1.PWGTP)))**2,(SUM(t1.pwgtp25)-(SUM(t1.PWGTP)))**2
        ,(SUM(t1.pwgtp26)-(SUM(t1.PWGTP)))**2,(SUM(t1.pwgtp27)-(SUM(t1.PWGTP)))**2,(SUM(t1.pwgtp28)-(SUM(t1.PWGTP)))**2,(SUM(t1.pwgtp29)-(SUM(t1.PWGTP)))**2,(SUM(t1.pwgtp30)-(SUM(t1.PWGTP)))**2
            ,(SUM(t1.pwgtp31)-(SUM(t1.PWGTP)))**2,(SUM(t1.pwgtp32)-(SUM(t1.PWGTP)))**2,(SUM(t1.pwgtp33)-(SUM(t1.PWGTP)))**2,(SUM(t1.pwgtp34)-(SUM(t1.PWGTP)))**2,(SUM(t1.pwgtp35)-(SUM(t1.PWGTP)))**2
        ,(SUM(t1.pwgtp36)-(SUM(t1.PWGTP)))**2,(SUM(t1.pwgtp37)-(SUM(t1.PWGTP)))**2,(SUM(t1.pwgtp38)-(SUM(t1.PWGTP)))**2,(SUM(t1.pwgtp39)-(SUM(t1.PWGTP)))**2,(SUM(t1.pwgtp40)-(SUM(t1.PWGTP)))**2
        ,(SUM(t1.pwgtp41)-(SUM(t1.PWGTP)))**2,(SUM(t1.pwgtp42)-(SUM(t1.PWGTP)))**2,(SUM(t1.pwgtp43)-(SUM(t1.PWGTP)))**2,(SUM(t1.pwgtp44)-(SUM(t1.PWGTP)))**2,(SUM(t1.pwgtp45)-(SUM(t1.PWGTP)))**2
        ,(SUM(t1.pwgtp46)-(SUM(t1.PWGTP)))**2,(SUM(t1.pwgtp47)-(SUM(t1.PWGTP)))**2,(SUM(t1.pwgtp48)-(SUM(t1.PWGTP)))**2,(SUM(t1.pwgtp49)-(SUM(t1.PWGTP)))**2,(SUM(t1.pwgtp50)-(SUM(t1.PWGTP)))**2
        ,(SUM(t1.pwgtp51)-(SUM(t1.PWGTP)))**2,(SUM(t1.pwgtp52)-(SUM(t1.PWGTP)))**2,(SUM(t1.pwgtp53)-(SUM(t1.PWGTP)))**2,(SUM(t1.pwgtp54)-(SUM(t1.PWGTP)))**2,(SUM(t1.pwgtp55)-(SUM(t1.PWGTP)))**2
        ,(SUM(t1.pwgtp56)-(SUM(t1.PWGTP)))**2,(SUM(t1.pwgtp57)-(SUM(t1.PWGTP)))**2,(SUM(t1.pwgtp58)-(SUM(t1.PWGTP)))**2,(SUM(t1.pwgtp59)-(SUM(t1.PWGTP)))**2,(SUM(t1.pwgtp60)-(SUM(t1.PWGTP)))**2
        ,(SUM(t1.pwgtp61)-(SUM(t1.PWGTP)))**2,(SUM(t1.pwgtp62)-(SUM(t1.PWGTP)))**2,(SUM(t1.pwgtp63)-(SUM(t1.PWGTP)))**2,(SUM(t1.pwgtp64)-(SUM(t1.PWGTP)))**2,(SUM(t1.pwgtp65)-(SUM(t1.PWGTP)))**2
        ,(SUM(t1.pwgtp66)-(SUM(t1.PWGTP)))**2,(SUM(t1.pwgtp67)-(SUM(t1.PWGTP)))**2,(SUM(t1.pwgtp68)-(SUM(t1.PWGTP)))**2,(SUM(t1.pwgtp69)-(SUM(t1.PWGTP)))**2,(SUM(t1.pwgtp70)-(SUM(t1.PWGTP)))**2
        ,(SUM(t1.pwgtp71)-(SUM(t1.PWGTP)))**2,(SUM(t1.pwgtp72)-(SUM(t1.PWGTP)))**2,(SUM(t1.pwgtp73)-(SUM(t1.PWGTP)))**2,(SUM(t1.pwgtp74)-(SUM(t1.PWGTP)))**2,(SUM(t1.pwgtp75)-(SUM(t1.PWGTP)))**2
        ,(SUM(t1.pwgtp76)-(SUM(t1.PWGTP)))**2,(SUM(t1.pwgtp77)-(SUM(t1.PWGTP)))**2,(SUM(t1.pwgtp78)-(SUM(t1.PWGTP)))**2,(SUM(t1.pwgtp79)-(SUM(t1.PWGTP)))**2,(SUM(t1.pwgtp80)-(SUM(t1.PWGTP)))**2)))) 
            AS Margin_of_Error,

        /* Plus_Minus */

        (1.64537*(SQRT(.05*(SUM((SUM(t1.pwgtp1)-(SUM(t1.PWGTP)))**2,(SUM(t1.pwgtp2)-(SUM(t1.PWGTP)))**2,(SUM(t1.pwgtp3)-(SUM(t1.PWGTP)))**2,(SUM(t1.pwgtp4)-(SUM(t1.PWGTP)))**2,(SUM(t1.pwgtp5)-(SUM(t1.PWGTP)))**2
            ,(SUM(t1.pwgtp6)-(SUM(t1.PWGTP)))**2,(SUM(t1.pwgtp7)-(SUM(t1.PWGTP)))**2,(SUM(t1.pwgtp8)-(SUM(t1.PWGTP)))**2,(SUM(t1.pwgtp9)-(SUM(t1.PWGTP)))**2,(SUM(t1.pwgtp10)-(SUM(t1.PWGTP)))**2
            ,(SUM(t1.pwgtp11)-(SUM(t1.PWGTP)))**2,(SUM(t1.pwgtp12)-(SUM(t1.PWGTP)))**2,(SUM(t1.pwgtp13)-(SUM(t1.PWGTP)))**2,(SUM(t1.pwgtp14)-(SUM(t1.PWGTP)))**2,(SUM(t1.pwgtp15)-(SUM(t1.PWGTP)))**2
            ,(SUM(t1.pwgtp16)-(SUM(t1.PWGTP)))**2,(SUM(t1.pwgtp17)-(SUM(t1.PWGTP)))**2,(SUM(t1.pwgtp18)-(SUM(t1.PWGTP)))**2,(SUM(t1.pwgtp19)-(SUM(t1.PWGTP)))**2,(SUM(t1.pwgtp20)-(SUM(t1.PWGTP)))**2
        ,(SUM(t1.pwgtp21)-(SUM(t1.PWGTP)))**2,(SUM(t1.pwgtp22)-(SUM(t1.PWGTP)))**2,(SUM(t1.pwgtp23)-(SUM(t1.PWGTP)))**2,(SUM(t1.pwgtp24)-(SUM(t1.PWGTP)))**2,(SUM(t1.pwgtp25)-(SUM(t1.PWGTP)))**2
        ,(SUM(t1.pwgtp26)-(SUM(t1.PWGTP)))**2,(SUM(t1.pwgtp27)-(SUM(t1.PWGTP)))**2,(SUM(t1.pwgtp28)-(SUM(t1.PWGTP)))**2,(SUM(t1.pwgtp29)-(SUM(t1.PWGTP)))**2,(SUM(t1.pwgtp30)-(SUM(t1.PWGTP)))**2
            ,(SUM(t1.pwgtp31)-(SUM(t1.PWGTP)))**2,(SUM(t1.pwgtp32)-(SUM(t1.PWGTP)))**2,(SUM(t1.pwgtp33)-(SUM(t1.PWGTP)))**2,(SUM(t1.pwgtp34)-(SUM(t1.PWGTP)))**2,(SUM(t1.pwgtp35)-(SUM(t1.PWGTP)))**2
        ,(SUM(t1.pwgtp36)-(SUM(t1.PWGTP)))**2,(SUM(t1.pwgtp37)-(SUM(t1.PWGTP)))**2,(SUM(t1.pwgtp38)-(SUM(t1.PWGTP)))**2,(SUM(t1.pwgtp39)-(SUM(t1.PWGTP)))**2,(SUM(t1.pwgtp40)-(SUM(t1.PWGTP)))**2
        ,(SUM(t1.pwgtp41)-(SUM(t1.PWGTP)))**2,(SUM(t1.pwgtp42)-(SUM(t1.PWGTP)))**2,(SUM(t1.pwgtp43)-(SUM(t1.PWGTP)))**2,(SUM(t1.pwgtp44)-(SUM(t1.PWGTP)))**2,(SUM(t1.pwgtp45)-(SUM(t1.PWGTP)))**2
        ,(SUM(t1.pwgtp46)-(SUM(t1.PWGTP)))**2,(SUM(t1.pwgtp47)-(SUM(t1.PWGTP)))**2,(SUM(t1.pwgtp48)-(SUM(t1.PWGTP)))**2,(SUM(t1.pwgtp49)-(SUM(t1.PWGTP)))**2,(SUM(t1.pwgtp50)-(SUM(t1.PWGTP)))**2
        ,(SUM(t1.pwgtp51)-(SUM(t1.PWGTP)))**2,(SUM(t1.pwgtp52)-(SUM(t1.PWGTP)))**2,(SUM(t1.pwgtp53)-(SUM(t1.PWGTP)))**2,(SUM(t1.pwgtp54)-(SUM(t1.PWGTP)))**2,(SUM(t1.pwgtp55)-(SUM(t1.PWGTP)))**2
        ,(SUM(t1.pwgtp56)-(SUM(t1.PWGTP)))**2,(SUM(t1.pwgtp57)-(SUM(t1.PWGTP)))**2,(SUM(t1.pwgtp58)-(SUM(t1.PWGTP)))**2,(SUM(t1.pwgtp59)-(SUM(t1.PWGTP)))**2,(SUM(t1.pwgtp60)-(SUM(t1.PWGTP)))**2
        ,(SUM(t1.pwgtp61)-(SUM(t1.PWGTP)))**2,(SUM(t1.pwgtp62)-(SUM(t1.PWGTP)))**2,(SUM(t1.pwgtp63)-(SUM(t1.PWGTP)))**2,(SUM(t1.pwgtp64)-(SUM(t1.PWGTP)))**2,(SUM(t1.pwgtp65)-(SUM(t1.PWGTP)))**2
        ,(SUM(t1.pwgtp66)-(SUM(t1.PWGTP)))**2,(SUM(t1.pwgtp67)-(SUM(t1.PWGTP)))**2,(SUM(t1.pwgtp68)-(SUM(t1.PWGTP)))**2,(SUM(t1.pwgtp69)-(SUM(t1.PWGTP)))**2,(SUM(t1.pwgtp70)-(SUM(t1.PWGTP)))**2
        ,(SUM(t1.pwgtp71)-(SUM(t1.PWGTP)))**2,(SUM(t1.pwgtp72)-(SUM(t1.PWGTP)))**2,(SUM(t1.pwgtp73)-(SUM(t1.PWGTP)))**2,(SUM(t1.pwgtp74)-(SUM(t1.PWGTP)))**2,(SUM(t1.pwgtp75)-(SUM(t1.PWGTP)))**2
        ,(SUM(t1.pwgtp76)-(SUM(t1.PWGTP)))**2,(SUM(t1.pwgtp77)-(SUM(t1.PWGTP)))**2,(SUM(t1.pwgtp78)-(SUM(t1.PWGTP)))**2,(SUM(t1.pwgtp79)-(SUM(t1.PWGTP)))**2,(SUM(t1.pwgtp80)-(SUM(t1.PWGTP)))**2)))))/ (SUM(t1.PWGTP))
            FORMAT=PERCENT6.1 AS Plus_Minus_Percent

/* End of Margin of Error Code */

要查看这适合查询的位置，这里是嵌入了错误代码边际的查询的完整代码。

/* Example program */


PROC SQL;
   CREATE TABLE WORK.QUERY_FOR_COMBINEDACS2013_SAS7BD(label="QUERY_FOR_combinedacs2013.sas7bdat") AS 
   SELECT /* SUM_of_PWGTP */
            (SUM(t1.PWGTP)) FORMAT=Z5. AS SUM_of_PWGTP, 
          t1.SCHL,
        1.64537*(SQRT(.05*(SUM((SUM(t1.pwgtp1)-(SUM(t1.PWGTP)))**2,(SUM(t1.pwgtp2)-(SUM(t1.PWGTP)))**2,(SUM(t1.pwgtp3)-(SUM(t1.PWGTP)))**2,(SUM(t1.pwgtp4)-(SUM(t1.PWGTP)))**2,(SUM(t1.pwgtp5)-(SUM(t1.PWGTP)))**2
            ,(SUM(t1.pwgtp6)-(SUM(t1.PWGTP)))**2,(SUM(t1.pwgtp7)-(SUM(t1.PWGTP)))**2,(SUM(t1.pwgtp8)-(SUM(t1.PWGTP)))**2,(SUM(t1.pwgtp9)-(SUM(t1.PWGTP)))**2,(SUM(t1.pwgtp10)-(SUM(t1.PWGTP)))**2
            ,(SUM(t1.pwgtp11)-(SUM(t1.PWGTP)))**2,(SUM(t1.pwgtp12)-(SUM(t1.PWGTP)))**2,(SUM(t1.pwgtp13)-(SUM(t1.PWGTP)))**2,(SUM(t1.pwgtp14)-(SUM(t1.PWGTP)))**2,(SUM(t1.pwgtp15)-(SUM(t1.PWGTP)))**2
            ,(SUM(t1.pwgtp16)-(SUM(t1.PWGTP)))**2,(SUM(t1.pwgtp17)-(SUM(t1.PWGTP)))**2,(SUM(t1.pwgtp18)-(SUM(t1.PWGTP)))**2,(SUM(t1.pwgtp19)-(SUM(t1.PWGTP)))**2,(SUM(t1.pwgtp20)-(SUM(t1.PWGTP)))**2
        ,(SUM(t1.pwgtp21)-(SUM(t1.PWGTP)))**2,(SUM(t1.pwgtp22)-(SUM(t1.PWGTP)))**2,(SUM(t1.pwgtp23)-(SUM(t1.PWGTP)))**2,(SUM(t1.pwgtp24)-(SUM(t1.PWGTP)))**2,(SUM(t1.pwgtp25)-(SUM(t1.PWGTP)))**2
        ,(SUM(t1.pwgtp26)-(SUM(t1.PWGTP)))**2,(SUM(t1.pwgtp27)-(SUM(t1.PWGTP)))**2,(SUM(t1.pwgtp28)-(SUM(t1.PWGTP)))**2,(SUM(t1.pwgtp29)-(SUM(t1.PWGTP)))**2,(SUM(t1.pwgtp30)-(SUM(t1.PWGTP)))**2
            ,(SUM(t1.pwgtp31)-(SUM(t1.PWGTP)))**2,(SUM(t1.pwgtp32)-(SUM(t1.PWGTP)))**2,(SUM(t1.pwgtp33)-(SUM(t1.PWGTP)))**2,(SUM(t1.pwgtp34)-(SUM(t1.PWGTP)))**2,(SUM(t1.pwgtp35)-(SUM(t1.PWGTP)))**2
        ,(SUM(t1.pwgtp36)-(SUM(t1.PWGTP)))**2,(SUM(t1.pwgtp37)-(SUM(t1.PWGTP)))**2,(SUM(t1.pwgtp38)-(SUM(t1.PWGTP)))**2,(SUM(t1.pwgtp39)-(SUM(t1.PWGTP)))**2,(SUM(t1.pwgtp40)-(SUM(t1.PWGTP)))**2
        ,(SUM(t1.pwgtp41)-(SUM(t1.PWGTP)))**2,(SUM(t1.pwgtp42)-(SUM(t1.PWGTP)))**2,(SUM(t1.pwgtp43)-(SUM(t1.PWGTP)))**2,(SUM(t1.pwgtp44)-(SUM(t1.PWGTP)))**2,(SUM(t1.pwgtp45)-(SUM(t1.PWGTP)))**2
        ,(SUM(t1.pwgtp46)-(SUM(t1.PWGTP)))**2,(SUM(t1.pwgtp47)-(SUM(t1.PWGTP)))**2,(SUM(t1.pwgtp48)-(SUM(t1.PWGTP)))**2,(SUM(t1.pwgtp49)-(SUM(t1.PWGTP)))**2,(SUM(t1.pwgtp50)-(SUM(t1.PWGTP)))**2
        ,(SUM(t1.pwgtp51)-(SUM(t1.PWGTP)))**2,(SUM(t1.pwgtp52)-(SUM(t1.PWGTP)))**2,(SUM(t1.pwgtp53)-(SUM(t1.PWGTP)))**2,(SUM(t1.pwgtp54)-(SUM(t1.PWGTP)))**2,(SUM(t1.pwgtp55)-(SUM(t1.PWGTP)))**2
        ,(SUM(t1.pwgtp56)-(SUM(t1.PWGTP)))**2,(SUM(t1.pwgtp57)-(SUM(t1.PWGTP)))**2,(SUM(t1.pwgtp58)-(SUM(t1.PWGTP)))**2,(SUM(t1.pwgtp59)-(SUM(t1.PWGTP)))**2,(SUM(t1.pwgtp60)-(SUM(t1.PWGTP)))**2
        ,(SUM(t1.pwgtp61)-(SUM(t1.PWGTP)))**2,(SUM(t1.pwgtp62)-(SUM(t1.PWGTP)))**2,(SUM(t1.pwgtp63)-(SUM(t1.PWGTP)))**2,(SUM(t1.pwgtp64)-(SUM(t1.PWGTP)))**2,(SUM(t1.pwgtp65)-(SUM(t1.PWGTP)))**2
        ,(SUM(t1.pwgtp66)-(SUM(t1.PWGTP)))**2,(SUM(t1.pwgtp67)-(SUM(t1.PWGTP)))**2,(SUM(t1.pwgtp68)-(SUM(t1.PWGTP)))**2,(SUM(t1.pwgtp69)-(SUM(t1.PWGTP)))**2,(SUM(t1.pwgtp70)-(SUM(t1.PWGTP)))**2
        ,(SUM(t1.pwgtp71)-(SUM(t1.PWGTP)))**2,(SUM(t1.pwgtp72)-(SUM(t1.PWGTP)))**2,(SUM(t1.pwgtp73)-(SUM(t1.PWGTP)))**2,(SUM(t1.pwgtp74)-(SUM(t1.PWGTP)))**2,(SUM(t1.pwgtp75)-(SUM(t1.PWGTP)))**2
        ,(SUM(t1.pwgtp76)-(SUM(t1.PWGTP)))**2,(SUM(t1.pwgtp77)-(SUM(t1.PWGTP)))**2,(SUM(t1.pwgtp78)-(SUM(t1.PWGTP)))**2,(SUM(t1.pwgtp79)-(SUM(t1.PWGTP)))**2,(SUM(t1.pwgtp80)-(SUM(t1.PWGTP)))**2)))) 
            AS Margin_of_Error,

        /* Plus_Minus */

        (1.64537*(SQRT(.05*(SUM((SUM(t1.pwgtp1)-(SUM(t1.PWGTP)))**2,(SUM(t1.pwgtp2)-(SUM(t1.PWGTP)))**2,(SUM(t1.pwgtp3)-(SUM(t1.PWGTP)))**2,(SUM(t1.pwgtp4)-(SUM(t1.PWGTP)))**2,(SUM(t1.pwgtp5)-(SUM(t1.PWGTP)))**2
            ,(SUM(t1.pwgtp6)-(SUM(t1.PWGTP)))**2,(SUM(t1.pwgtp7)-(SUM(t1.PWGTP)))**2,(SUM(t1.pwgtp8)-(SUM(t1.PWGTP)))**2,(SUM(t1.pwgtp9)-(SUM(t1.PWGTP)))**2,(SUM(t1.pwgtp10)-(SUM(t1.PWGTP)))**2
            ,(SUM(t1.pwgtp11)-(SUM(t1.PWGTP)))**2,(SUM(t1.pwgtp12)-(SUM(t1.PWGTP)))**2,(SUM(t1.pwgtp13)-(SUM(t1.PWGTP)))**2,(SUM(t1.pwgtp14)-(SUM(t1.PWGTP)))**2,(SUM(t1.pwgtp15)-(SUM(t1.PWGTP)))**2
            ,(SUM(t1.pwgtp16)-(SUM(t1.PWGTP)))**2,(SUM(t1.pwgtp17)-(SUM(t1.PWGTP)))**2,(SUM(t1.pwgtp18)-(SUM(t1.PWGTP)))**2,(SUM(t1.pwgtp19)-(SUM(t1.PWGTP)))**2,(SUM(t1.pwgtp20)-(SUM(t1.PWGTP)))**2
        ,(SUM(t1.pwgtp21)-(SUM(t1.PWGTP)))**2,(SUM(t1.pwgtp22)-(SUM(t1.PWGTP)))**2,(SUM(t1.pwgtp23)-(SUM(t1.PWGTP)))**2,(SUM(t1.pwgtp24)-(SUM(t1.PWGTP)))**2,(SUM(t1.pwgtp25)-(SUM(t1.PWGTP)))**2
        ,(SUM(t1.pwgtp26)-(SUM(t1.PWGTP)))**2,(SUM(t1.pwgtp27)-(SUM(t1.PWGTP)))**2,(SUM(t1.pwgtp28)-(SUM(t1.PWGTP)))**2,(SUM(t1.pwgtp29)-(SUM(t1.PWGTP)))**2,(SUM(t1.pwgtp30)-(SUM(t1.PWGTP)))**2
            ,(SUM(t1.pwgtp31)-(SUM(t1.PWGTP)))**2,(SUM(t1.pwgtp32)-(SUM(t1.PWGTP)))**2,(SUM(t1.pwgtp33)-(SUM(t1.PWGTP)))**2,(SUM(t1.pwgtp34)-(SUM(t1.PWGTP)))**2,(SUM(t1.pwgtp35)-(SUM(t1.PWGTP)))**2
        ,(SUM(t1.pwgtp36)-(SUM(t1.PWGTP)))**2,(SUM(t1.pwgtp37)-(SUM(t1.PWGTP)))**2,(SUM(t1.pwgtp38)-(SUM(t1.PWGTP)))**2,(SUM(t1.pwgtp39)-(SUM(t1.PWGTP)))**2,(SUM(t1.pwgtp40)-(SUM(t1.PWGTP)))**2
        ,(SUM(t1.pwgtp41)-(SUM(t1.PWGTP)))**2,(SUM(t1.pwgtp42)-(SUM(t1.PWGTP)))**2,(SUM(t1.pwgtp43)-(SUM(t1.PWGTP)))**2,(SUM(t1.pwgtp44)-(SUM(t1.PWGTP)))**2,(SUM(t1.pwgtp45)-(SUM(t1.PWGTP)))**2
        ,(SUM(t1.pwgtp46)-(SUM(t1.PWGTP)))**2,(SUM(t1.pwgtp47)-(SUM(t1.PWGTP)))**2,(SUM(t1.pwgtp48)-(SUM(t1.PWGTP)))**2,(SUM(t1.pwgtp49)-(SUM(t1.PWGTP)))**2,(SUM(t1.pwgtp50)-(SUM(t1.PWGTP)))**2
        ,(SUM(t1.pwgtp51)-(SUM(t1.PWGTP)))**2,(SUM(t1.pwgtp52)-(SUM(t1.PWGTP)))**2,(SUM(t1.pwgtp53)-(SUM(t1.PWGTP)))**2,(SUM(t1.pwgtp54)-(SUM(t1.PWGTP)))**2,(SUM(t1.pwgtp55)-(SUM(t1.PWGTP)))**2
        ,(SUM(t1.pwgtp56)-(SUM(t1.PWGTP)))**2,(SUM(t1.pwgtp57)-(SUM(t1.PWGTP)))**2,(SUM(t1.pwgtp58)-(SUM(t1.PWGTP)))**2,(SUM(t1.pwgtp59)-(SUM(t1.PWGTP)))**2,(SUM(t1.pwgtp60)-(SUM(t1.PWGTP)))**2
        ,(SUM(t1.pwgtp61)-(SUM(t1.PWGTP)))**2,(SUM(t1.pwgtp62)-(SUM(t1.PWGTP)))**2,(SUM(t1.pwgtp63)-(SUM(t1.PWGTP)))**2,(SUM(t1.pwgtp64)-(SUM(t1.PWGTP)))**2,(SUM(t1.pwgtp65)-(SUM(t1.PWGTP)))**2
        ,(SUM(t1.pwgtp66)-(SUM(t1.PWGTP)))**2,(SUM(t1.pwgtp67)-(SUM(t1.PWGTP)))**2,(SUM(t1.pwgtp68)-(SUM(t1.PWGTP)))**2,(SUM(t1.pwgtp69)-(SUM(t1.PWGTP)))**2,(SUM(t1.pwgtp70)-(SUM(t1.PWGTP)))**2
        ,(SUM(t1.pwgtp71)-(SUM(t1.PWGTP)))**2,(SUM(t1.pwgtp72)-(SUM(t1.PWGTP)))**2,(SUM(t1.pwgtp73)-(SUM(t1.PWGTP)))**2,(SUM(t1.pwgtp74)-(SUM(t1.PWGTP)))**2,(SUM(t1.pwgtp75)-(SUM(t1.PWGTP)))**2
        ,(SUM(t1.pwgtp76)-(SUM(t1.PWGTP)))**2,(SUM(t1.pwgtp77)-(SUM(t1.PWGTP)))**2,(SUM(t1.pwgtp78)-(SUM(t1.PWGTP)))**2,(SUM(t1.pwgtp79)-(SUM(t1.PWGTP)))**2,(SUM(t1.pwgtp80)-(SUM(t1.PWGTP)))**2)))))/ (SUM(t1.PWGTP))
            FORMAT=PERCENT6.1 AS Plus_Minus_Percent
      FROM EC100005.combinedacs2013 t1
      GROUP BY t1.SCHL;
QUIT;

/* End example program */

【讨论】：

这可能是正确的，但它相当笨拙。在使用数组的数据步骤中执行此操作会更整洁。