Blog

Latest News and Updates

Rlogo-1.png

(Ellips) 100(1-a)% – software R

Kali ini akan dibahas mengenai cara membuat Daerah (ellips) Kepercayaan 100(1-alpha)% bagi Vektor Nilai Tengah dalam software R.

Agar lebih mudah, buat fungsi yang bisa berlaku umum untuk suatu peubah acak X yang menyebar normal ganda (X merupakan matriks berukuran (n x 2)).

Berikut adalah fungsi yang telah dibuat dengan nama “conf.ell” (nama fungsi bisa disesuaikan). X merupakan matriks berukuran (n x 2) yang terdiri dari peubah x1 dan x2, alpha merupakan taraf nyata yang dipergunakan dengan default 0.05, x1 dan x2 berturut-turut merupakan nama peubah acak x1 dan x2.
[box] conf.ell <- function( X, alpha = 0.05, x1 = “X1”, x2 = “X2” )
{

judul <- paste(“Daerah (ellips) Kepercayaan “, 100*(1-alpha),”% bagi Vektor Nilai Tengah”, sep = “”)

dimnames(X)[[2]] <- c(x1,x2)

n <- nrow(X)

mu <- c(mean(X[,1]), mean(X[,2]))

sigma2 <- var(X)

plot(X, main = judul)

evals <- eigen(sigma2)$values

evecs <- eigen(sigma2)$vectors

a <- seq(0, 2*pi, len=1000000)

c2 <- qf(1-alpha, 2, n-2)*(n-1)*2/(n-2)

c <- sqrt(c2)

ma <-  c * sqrt(evals[1]/n)

mi <- c * sqrt(evals[2]/n)

xT <- ma * cos(a)

yT <- mi * sin(a)

M <- cbind(xT, yT)

sma <- cbind(mu+ma*evecs[,1], mu-ma*evecs[,1])

smi <- cbind(mu+mi*evecs[,2], mu-mi*evecs[,2])

transM <- evecs %*% t(M)

transM <- t(transM)

lines(transM + matrix(rep(mu,50), 1000000, 2, byrow=T), lwd = 2)

lines(t(sma), lwd = 2)

lines(t(smi), lwd = 2)

}[/box]

Contoh :
[box]X <- scan()

72.8

46

59.2

66.7

84.2

50.4

49.6

77.9

63.9

55.1

69.9

68.9

58.4

78.2

63.9

54.6

66.5

71.6

77.2

56.8

X <- matrix(X, 10, 2)

conf.ell(X, 0.1, “Matematika”, “Fisika”)
[/box]

gsbipb(Ellips) 100(1-a)% – software R
Share this post

Join the conversation