# R/is.f.R In DnE: Distribution and Equation

#### Documented in is.f

is.f <-
function(x,m,a=10,k1=NULL,k2=NULL)
{
re=1;
for(i in 1:length(x))
if(x[i]<0)
re=-1;
p=rep(0,m+2);
y=rep(0,m+2);
q=0;
if(re==-1)
{
return(data.frame("state"=-1,"pvalue"=1));
}
else
{
x1=mean(x);sum=0;
for(i in 1:length(x))
{
sum=sum+x[i]^2;
}
x2=sum/length(x);
if(is.null(k1)&&is.null(k2))
{
k2=2*x1/(x1-1);
k1=2*x1^2*(k2-2)/((k2-4)*(x2-x1^2)-2*x1^2);
df=m-1;
}
else if(is.null(k1))
{
k1=2*x1^2*(k2-2)/((k2-4)*(x2-x1^2)-2*x1^2);
df=m;
}
else if(is.null(k2))
{
k2=2*x1/(x1-1);
df=m;
}
else
{
df=m+1;
}
if(k1>0&&k2>4)
{
di=max(x)-min(x);
for(i in 1:m)
{
p[i]=pf(min(x)+di*i/m,k1,k2)-pf(min(x)+di*(i-1)/m,k1,k2);
if(p[i]==0)
{
break;
}
for(j in 1:length(x))
if(x[j]>(min(x)+di*(i-1)/m) && x[j]<=(min(x)+di*i/m))
y[i]=y[i]+1;
q=q+(y[i]-(length(x)*p[i]))^2/(length(x)*p[i]);
}
p[m+1]=pf(Inf,k1,k2)-pf(max(x),k1,k2);
p[m+2]=pf(min(x),k1,k2);
y[m+1]=length(which(x==min(x)));
if(p[m+1]!=0)
q=q+(y[m+1]-(length(x)*p[m+1]))^2/(length(x)*p[m+1]);
if(p[m+2]!=0)
q=q+(y[m+2]-(length(x)*p[m+2]))^2/(length(x)*p[m+2]);
q0=qchisq(1-a,df);
pvalue=pchisq(q,df);
if(q<=q0)
{
return(data.frame("qchisq"=q,"pvalue"=pvalue));
}
else
{
return(data.frame("state"=-1,"pvalue"=1));
}
}
else
{
return(data.frame("state"=-1,"pvalue"=1));
}
}
}

## Try the DnE package in your browser

Any scripts or data that you put into this service are public.

DnE documentation built on May 2, 2019, 4:03 p.m.