Location via proxy:   [ UP ]  
[Report a bug]   [Manage cookies]                
Scribd
Upload
14 views

Source Code

Source Codes

Uploaded by

Teguh Prasetyo
Copyright
© Attribution Non-Commercial (BY-NC)
Available Formats
Download as TXT, PDF, TXT or read online on Scribd
14 views

Source Code

Source Codes

Uploaded by

Teguh Prasetyo
Copyright
© Attribution Non-Commercial (BY-NC)
Available Formats
Download as TXT, PDF, TXT or read online on Scribd
You are on page 1/ 2

varian<-function(x)

{
a=length(x)
b=sum(x)
c=b/a
d=x-c
e=d^2
f=sum(e)
g=a-1
h=f/g
print(h)
}
uji_t=function(x,mu0,arah)
{
df=length(x)-1
T=abs((mean(x)-mu0)/(sd(x))/sqrt(length(x)))
if (arah==0) P=2*(1-pt(abs(T),df))
else if (arah==-1) P=pt(T,df)
else P=pt(T,df,lower.tail=FALSE)
cat ("Uji t Sampel Tunggal","\n")
cat ("H0:mu =",mu0,"\n")
if (arah==0) cat("H1:mu tidak=",mu0,"\n")
if (arah==-1) cat("H1:mu <",mu0,"\n")
if (arah==1) cat("H1:mu >",mu0,"\n")
cat ("Mean =",mean(x),",Stdev =",sd(x),"n =",length(x),"\n")
cat ("T =",T,",df =",df,",P value =",P,"\n")
}
atau
uji_t=function(x,mu0,arah)
{
df=length(x)-1
T=abs((mean(x)-mu0)/(sd(x))/sqrt(length(x)))
if (arah==0) P=2*(1-pt(abs(T),df))
else if (arah==-1) P=pt(T,df)
else P=pt(T,df,upper.tail=TRUE)
cat ("Uji t Sampel Tunggal","\n")
cat ("H0:mu =",mu0,"\n")
if (arah==0) cat("H1:mu tidak=",mu0,"\n")
if (arah==-1) cat("H1:mu <",mu0,"\n")
if (arah==1) cat("H1:mu >",mu0,"\n")
cat ("Mean =",mean(x),",Stdev =",sd(x),"n =",length(x),"\n")
cat ("T =",T,",df =",df,",P value =",P,"\n")
}

regresi_berganda=function(x,y)
{
a=t(x)
b=a%*%x
c=solve(b)
d=a%*%y
e=c%*%d
cat ("===========================================","\n")
cat ("Pendugaan Parameter Regresi Linier Berganda","\n")

cat ("===========================================","\n")
cat ("Beta =",e,"\n")
cat ("===========================================","\n")
}

uji_t_independen_var_sama=function(x1,x2,mu0,arah)
{
df=length(x1)+length(x2)-2
s=(((length(x1)-1)*var(x1))+((length(x2)-1)*var(x2)))/(((length(x1))-1)+((length
(x2))-1))
T=abs((mean(x1)-mean(x2))-(mu0))/sqrt(s*(1/length(x1))*(length(x2)))
if (arah==0) P=2*(1-pt(abs(T),df))
else if (arah==-1) P=pt(T,df)
else P=pt(T,df,lower.tail=FALSE)
cat ("=============================================","\n")
cat ("Uji t 2 Sampel Independen Ragam Dianggap Sama","\n")
cat ("=============================================","\n")
cat ("H0:mu a - mu b =",mu0,"\n")
if (arah==0) cat("H1:mu a - mu b tidak=",mu0,"\n")
if (arah==-1) cat("H1:mu a - mu b <",mu0,"\n")
if (arah==1) cat("H1:mu a - mu b >",mu0,"\n")
cat ("Mean x1=",mean(x1),",Mean x2=",mean(x2),",Stdev x1 =",sd(x1),",Stdev x2 ="
,sd(x2),",n =",length(x1)+length(x2),"\n")
cat ("Statistik t =",T,",df =",df,",P value =",P,"\n")
cat ("=============================================","\n")
}

You might also like