funs.base.class.spec | R Documentation |

Two functions: `base.class.spec.ct`

and `base.class.spec`

.

Both functions are objects of class `"classhtest"`

but with different arguments (see the parameter list below).
Each one performs class specific segregation tests due to Dixon for *k ≥ 2* classes. That is,
each one performs hypothesis tests of deviations of
entries in each row of NNCT from the expected values under RL or CSR for each row.
Recall that row labels in the NNCT are base class labels.
The test for each row *i* is based on the chi-squared approximation of the corresponding quadratic form
and are due to \insertCitedixon:NNCTEco2002;textualnnspat.

Each function yields the test statistic, *p*-value and `df`

for each base class *i*, description of the
alternative with the corresponding null values (i.e. expected values) for the row *i*, estimates for the entries in row *i*
for *i=1,…,k*. The functions also provide names of the test statistics, the method and the data set used.

The null hypothesis for each row is that the corresponding *N_{ij}* entries in row *i* are equal to their
expected values under RL or CSR.

See also (\insertCitedixon:NNCTEco2002,ceyhan:stat-neer-class2009;textualnnspat) and the references therein.

base.class.spec.ct(ct, covN) base.class.spec(dat, lab, ...)

`ct` |
A nearest neighbor contingency table, used in |

`covN` |
The |

`dat` |
The data set in one or higher dimensions, each row corresponds to a data point,
used in |

`lab` |
The |

`...` |
are for further arguments, such as |

A `list`

with the elements

`type` |
Type of the class-specific test, which is |

`statistic` |
The |

`stat.names` |
Name of the test statistics |

`p.value` |
The |

`df` |
Degrees of freedom for the chi-squared test, which is |

`estimate` |
Estimates of the parameters, NNCT, i.e., matrix of the observed |

`null.value` |
Matrix of hypothesized null values for the parameters which are expected values of
the |

`null.name` |
Name of the null values |

`method` |
Description of the hypothesis test |

`ct.name` |
Name of the contingency table, |

`data.name` |
Name of the data set, |

Elvan Ceyhan

`NN.class.spec.ct`

, `NN.class.spec`

, `class.spec.ct`

and `class.spec`

n<-20 #or try sample(1:20,1) Y<-matrix(runif(3*n),ncol=3) ipd<-ipd.mat(Y) cls<-sample(1:2,n,replace = TRUE) #or try cls<-rep(1:2,c(10,10)) ct<-nnct(ipd,cls) W<-Wmat(ipd) Qv<-Qvec(W)$q Rv<-Rval(W) varN<-var.nnct(ct,Qv,Rv) covN<-cov.nnct(ct,varN,Qv,Rv) base.class.spec(Y,cls) base.class.spec.ct(ct,covN) base.class.spec(Y,cls,method="max") #cls as a factor na<-floor(n/2); nb<-n-na fcls<-rep(c("a","b"),c(na,nb)) ct<-nnct(ipd,fcls) base.class.spec(Y,fcls) base.class.spec.ct(ct,covN) ############# n<-40 Y<-matrix(runif(3*n),ncol=3) ipd<-ipd.mat(Y) cls<-sample(1:4,n,replace = TRUE) #or try cls<-rep(1:2,c(10,10)) ct<-nnct(ipd,cls) W<-Wmat(ipd) Qv<-Qvec(W)$q Rv<-Rval(W) varN<-var.nnct(ct,Qv,Rv) covN<-cov.nnct(ct,varN,Qv,Rv) base.class.spec(Y,cls) base.class.spec.ct(ct,covN)

Embedding an R snippet on your website

Add the following code to your website.

For more information on customizing the embed code, read Embedding Snippets.