Tương quan và hồi quy với SAS

Bạn đang xem bản rút gọn của tài liệu. Xem và tải ngay bản đầy đủ của tài liệu tại đây (73.03 KB, 5 trang )

(1)<div class='page_container' data-page=1>

Ph

ầ

n 3

T

ươ

ng quan và H

ồ

i quy v

ớ

i ph

ầ

n m

ề

m SAS

3.1.

T

ươ

ng quan và h

ồ

i quy

ðể

tính h

ệ

s

ố

t

ươ

ng quan và xây d

ự

ng ph

ươ

ng trình h

ồ

i quy, s

ố

li

ệ

u ln ph

ả

i t

ạ

o thành t

ừ

ng

c

ặ

p và

ñượ

c nh

ậ

p vào t

ừ

ng c

ộ

t

ñố

i t

ừ

ng ch

ỉ

tiêu.

Ví d

ụ

11: Ti

ế

n hành cân kh

ố

i l

ượ

ng (P),

đ

o

đườ

ng kính l

ớ

n (D) và

đườ

ng kính bé (d) c

ủ

a 22

qu

ả

tr

ứ

ng gà. S

ố

li

ệ

u thu

đượ

c trình bày

ở

b

ả

ng d

ướ

i

ñ

ây.

P (gram)

66,80 60,10 71,20 61,60 61,20 59,00 67,90 59,00 51,50 62,60 64,20

D (mm)

58,37 54,95 60,58 56,73 57,36 53,26 57,07 58,17 52,28 55,62 56,82

d (mm)

45,12 44,35 45,56 44,34 43,57 44,86 46,27 42,82 41,91 44,95 44,79

P (gram)

71,20 54,20 54,50 69,10 55,90 66,00 68,00 62,00 56,70 67,00 53,80

D (mm)

61,15 54,24 54,99 60,99 54,41 58,19 59,93 56,80 55,66 58,49 52,44

d (mm)

46,00 42,58 42,32 44,85 42,62 45,69 45,50 44,20 42,41 45,56 43,38

3.1.1.

H

ệ

s

ố

t

ươ

ng quan

Gi

ả

thi

ế

t

ñố

i v

ớ

i ki

ể

m

đị

nh hai phía H

: ρ = 0 và

đố

i thi

ế

t H

: ρ ≠ 0, trong

đ

ó ρ là t

ươ

ng quan

gi

ữ

a 2 bi

ế

n nghiên c

ứ

u.

SAS CODE

DATA

SAS12;

INPUT

KL DL DN;

CARDS

;

</div>
(2)<div class='page_container' data-page=2>

PROC

CORR

;

VAR

KL DL DN;

RUN

;

K

ế

t qu

ả

t

ừ

SAS:

The CORR Procedure

3 Variables: KL DL DN

Simple Statistics

Variable N Mean Std Dev Sum Minimum Maximum
KL 22 61.97727 5.94426 1364 51.50000 71.20000
DL 22 56.75000 2.61551 1249 52.28000 61.15000
DN 22 44.25682 1.34160 973.65000 41.91000 46.27000

Pearson Correlation Coefficients, N = 22
Prob > |r| under H0: Rho=0

KL DL DN
KL 1.00000 0.89667 0.90527
<.0001 <.0001
DL 0.89667 1.00000 0.64802
<.0001 0.0011
DN 0.90527 0.64802 1.00000
<.0001 0.0011

H

ệ

s

ố

t

ươ

ng quan m

ẫ

u gi

ữ

a kh

ố

i l

ượ

ng và

đườ

ng kính l

ớ

n là 0,89667; kh

ố

i l

ượ

ng và

đườ

ng kính bé là 0,90527;

đườ

ng kính l

ớ

n và

đườ

ng kính bé là 0,64802. Xác su

ấ

t

ñố

i v

ớ

i

t

ừ

ng h

ệ

s

ố

t

ươ

ng quan (p-values)

đề

u bé h

ơ

n 0,05 (α) vì v

ậ

y k

ế

t lu

ậ

n t

ồ

n t

ạ

i m

ố

i quan h

ệ

gi

ữ

a

các ch

ỉ

tiêu này.

3.1.2.

Ph

ươ

ng trình h

ồ

i quy tuy

ế

n tính

Có th

ể

xây d

ự

ng h

ồ

i quy

ñơ

n bi

ế

n y = a + bx ho

ặ

c

ñ

a bi

ế

n y = a + b

x

+ b

x

+...+b

x

. V

ớ

i

ví d

ụ

11, ta có th

ể

xây d

ự

ng ph

ươ

ng trình h

ồ

i tuy

ế

n tính

đơ

n bi

ế

n quy

ướ

c tính kh

ố

i l

ượ

ng

tr

ứ

ng thơng qua

đườ

ng kính l

ớ

n/

đườ

ng kính bé ho

ặ

c

đ

a bi

ế

n thơng qua

đườ

ng kính l

ớ

n và

đườ

ng kính bé.

V

ớ

i SAS CODE sau

ñ

ây ta có th

ể

xây d

ự

ng ph

ươ

ng trình h

ồ

i quy tuy

ế

n tính v

ớ

i bi

ế

n ph

ụ

thu

ộ

c Y là kh

ố

i l

ượ

ng tr

ứ

ng và bi

ế

n

độ

c l

ậ

p X là

đườ

ng kính l

ớ

n:

SAS CODE

PROC

REG

;

MODEL

KL = DL;

</div>
(3)<div class='page_container' data-page=3>

K

ế

t qu

ả

t

ừ

SAS

The REG Procedure
Model: MODEL1
Dependent Variable: KL

Analysis of Variance

Sum of Mean

Source DF Squares Square F Value Pr > F
Model 1 596.59551 596.59551 82.05 <.0001
Error 20 145.42313 7.27116

Corrected Total 21 742.01864

Root MSE 2.69651 R-Square 0.8040
Dependent Mean 61.97727 Adj R-Sq 0.7942
Coeff Var 4.35080

Parameter Estimates
Parameter Standard

Variable DF Estimate Error t Value Pr > |t|
Intercept 1 -53.67124 12.78032 -4.20 0.0004
DL 1 2.03786 0.22498 9.06 <.0001

Ph

ươ

ng trình h

ồ

i quy v

ớ

i bi

ế

n ph

ụ

thu

ộ

c - kh

ố

i l

ượ

ng (Y) và bi

ế

n

ñộ

c l

ậ

p -

ñườ

ng kính l

ớ

n

(X) t

ừ

ph

ầ

n m

ề

m SAS thu

đượ

c nh

ư

Y = -53,7 + 2,04X.

Ph

ầ

n

Parameter Estimates

, cho ta k

ế

t qu

ả

ki

ể

m

ñị

nh các h

ệ

s

ố

c

ủ

a ph

ươ

ng trình h

ồ

i quy. V

ớ

i

xác su

ấ

t P = 0,0004 và P < 0,0001 ta có th

ể

k

ế

t lu

ậ

n các h

ệ

s

ố

trong ph

ươ

ng trình h

ồ

i quy

khác 0 (P < 0,05). H

ệ

s

ố

xác

đị

nh c

ủ

a ph

ươ

ng trình R² = 80,40%, hi

ệ

u ch

ỉ

nh R² = 79,4%.

Ta c

ũ

ng có th

ể

xây d

ự

ng h

ồ

i quy v

ớ

i 2 bi

ế

n

độ

c l

ậ

p là

đườ

ng kính l

ớ

n và

đườ

ng kính bé v

ớ

i

SAS CODE sau:

SAS CODE

PROC

REG

;

MODEL

KL = DL DN;

RUN

;

K

ế

t qu

ả

t

ừ

SAS:

</div>
(4)<div class='page_container' data-page=4>

Sum of Mean

Source DF Squares Square F Value Pr > F
Model 2 731.05126 365.52563 633.24 <.0001
Error 19 10.96738 0.57723

Corrected Total 21 742.01864

Root MSE 0.75976 R-Square 0.9852
Dependent Mean 61.97727 Adj R-Sq 0.9837
Coeff Var 1.22586

Parameter Estimates
Parameter Standard

Variable DF Estimate Error t Value Pr > |t|
Intercept 1 -116.55512 5.47204 -21.30 <.0001
DL 1 1.21473 0.08323 14.60 <.0001
DN 1 2.47638 0.16226 15.26 <.0001

Ta có k

ế

t qu

ả

hồn tồn t

ươ

ng t

ự

nh

ư

vi

ệ

c xây d

ự

ng ph

ươ

ng trình h

ồ

i quy

đơ

n gi

ả

n.

Y = - 116,55512 + 1.21473X

+ 2.47638X

Trong

đ

ó:

Y = kh

ố

i l

ượ

ng tr

ứ

ng,

X

=

đườ

ng kính l

ớ

n và

X

=

đườ

ng kính bé

ð

i

ề

u khác bi

ệ

t trong tr

ườ

ng h

ợ

p này là h

ệ

s

ố

xác

ñị

nh R² = 98,5% l

ớ

n h

ơ

n nhi

ề

u so v

ớ

i

tr

ườ

ng h

ợ

p

đơ

n bi

ế

n R² = 80,4%.

V

ớ

i tr

ườ

ng h

ợ

p ph

ươ

ng trình

ñượ

c xây d

ự

ng t

ừ

nhi

ề

u bi

ế

n

ñộ

c l

ậ

p, c

ầ

n ph

ả

i xác

đị

nh rõ vai

trị c

ủ

a t

ừ

ng bi

ế

n. Có nh

ữ

ng bi

ế

n

độ

c l

ậ

p khi

đư

a vào ph

ươ

ng trình

đ

óng góp khơng

đ

áng k

ể

vào mơ hình vì v

ậ

y c

ầ

n lo

ạ

i b

ỏ

khi xây d

ự

ng ph

ươ

ng trình.

ðể

có th

ể

tìm

đượ

c ph

ươ

ng trình

h

ồ

i quy t

ố

t nh

ấ

t c

ầ

n có thêm các l

ệ

nh ph

ụ

khác.

SAS CODE

PROC

REG

;

MODEL

KL = DL DN /

SSE

CP

AIC

SELECTION

=

CP

;

RUN

;

K

ế

t qu

ả

t

ừ

SAS :

The REG Procedure
Model: MODEL1
Dependent Variable: KL
C(p) Selection Method
Number in

Model C(p) R-Square AIC SSE Variables in Model
2 3.0000 0.9852 -9.3146 10.96738 DL DN

</div>
(5)<div class='page_container' data-page=5>

TÀI LI

Ệ

U THAM KH

Ả

O

Ti

ế

ng Vi

ệ

t

1. Nguy

ễ

n M

ạ

nh

ðứ

c, Giáo trình tin h

ọ

c

ứ

ng d

ụ

ng (dùng cho sinh viên các ngành kinh

t

ế

, nông – lâm nghi

ệ

p), Nxb. Nơng nghi

ệ

p, 2000.

2. Nguy

ễ

n

ð

ình Hi

ề

n (ch

ủ

biên),

ðỗ

ðứ

c L

ự

c, Giáo trình thi

ế

t k

ế

thí nghi

ệ

m (dùng cho

các ngành Thú y, Ch

ă

n nuôi

−

Thú y và Nuôi tr

ồ

ng thu

ỷ

s

ả

n), Nxb. Nông nghi

ệ

p,

2007.

3. ðỗ

ðứ

c L

ự

c, Bài gi

ả

ng ph

ươ

ng pháp thí nghi

ệ

m và x

ử

lý d

ữ

li

ệ

u v

ớ

i ph

ầ

n m

ề

m SAS,

Khoa Ch

ă

n nuôi và Nuôi tr

ồ

ng thu

ỷ

s

ả

n, Tr

ườ

ng

ðạ

i h

ọ

c Nông nghi

ệ

p Hà N

ộ

i, 2008.

4. Lê V

ă

n Ti

ế

n (1999), Lý thuy

ế

t xác su

ấ

t và th

ố

ng kê tốn h

ọ

c, Nxb. Nơng nghi

ệ

p,

1999.

5. Nguy

ễ

n H

ả

i Thanh (ch

ủ

biên),

ðỗ

Th

ị

M

ơ

và các tác gi

ả

khác, Tin h

ọ

c

ứ

ng d

ụ

ng

trong ngành nông nghi

ệ

p (Ch

ươ

ng I), Nxb. Khoa h

ọ

c và K

ỹ

thu

ậ

t, 2005.

Ti

ế

ng Anh

6. Kaps M. and Lamberson W.R. (2004). Biostatistics for animal science. CABI

Publishing

7. Kenneth N. Berk, Patrick Carey, Data Analysis with Microsoft Excel, Duxbury Press,

1998.

8. Peter Thomson and Louise Helby, Short Course Manual on Development of

Research and Teaching Skills in Experimental Design and Statistical Analysis – Part

A: Basic Applied Statistícs using Minitab, Uni. of Sydney, 2001.

9. Ramon C. Littell, Rudolf J. Freund, Philip C. Spector (1991). SAS System for Linear

Models, third edition (1991). SAS Institute Inc.

10. Richard A. Johnson and Dean W. Wichern (1992), Applied Multivariate Statistical

Analysis, Prentice Hall, New Jersay, 1992.

11. Robert D. Manson, Douglas A. Lind, Statistical Technique in Business and

Economics, Irwin Inc., 1990.

</div>

Tương quan và hồi quy với SAS

<b>Ph</b>

ầ

<b>n 3 </b>

<b>T</b>

ươ

<b>ng quan và H</b>

ồ

<b>i quy v</b>

ớ

<b>i ph</b>

ầ

<b>n m</b>

ề

<b>m SAS </b>

<b>3.1.</b>

<b>T</b>

ươ

<b>ng quan và h</b>

ồ

<b>i quy </b>

ðể

tính h

ệ

s

ố

t

ươ

ng quan và xây d

ự

ng ph

ươ

ng trình h

ồ

i quy, s

ố

li

ệ

u ln ph

ả

i t

ạ

o thành t

ừ

ng

c

ặ

p và

ñượ

c nh

ậ

p vào t

ừ

ng c

ộ

t

ñố

i t

ừ

ng ch

ỉ

tiêu.

<b>Ví d</b>

ụ

<b> 11: Ti</b>

ế

n hành cân kh

ố

i l

ượ

ng (P),

đ

o

đườ

ng kính l

ớ

n (D) và

đườ

ng kính bé (d) c

ủ