カイ二乗値の等しいテーブルを作る

• 予定パッケージ名
  • sphere
• 関数名
  • makeTblSameChiNM
• タイトル

makeTblSameChiNM 

• 説明

観測テーブルを与え、それと周辺度数を同じくするテーブルで、指定のカイ二乗値をとるテーブルを指定個作成する

• 使用例

makeTblSameChiNM (O,k,N)

• ソース

makeTblSameChiNM <-
function (O = matrix(c(10, 20, 30, 40, 50, 60), nrow = 2, byrow = TRUE), 
    k = 1, N = 1000) 
{
    m1 <- apply(O, 1, sum)
    m2 <- apply(O, 2, sum)
    ms <- sum(O)
    E <- outer(m1, m2, FUN = "*")/ms
    D <- O - E
    K <- sum(D^2/E)
    A <- matrix(runif(length(O) * N), nrow = N)
    A <- t(apply(A, 1, FUN = "standardizeNMfromLine", nrow = length(O[, 
        1]), byrow = TRUE))
    Kp <- sqrt(apply(t(t(A^2)/c(t(E))), 1, FUN = "sum"))
    t(t(A/Kp * sqrt(K) * k) + c(t(E)))
}

• Rdファイル

\name{makeTblSameChiNM }
\alias{makeTblSameChiNM }
%- Also NEED an '\alias' for EACH other topic documented here.
\title{
makeTblSameChiNM 
}
\description{
観測テーブルを与え、それと周辺度数を同じくするテーブルで、指定のカイ二乗値をとるテーブルを指定個作成する
}
\usage{
makeTblSameChiNM (O,k,N)
}
%- maybe also 'usage' for other objects documented here.
\details{
%%  ~~ If necessary, more details than the description above ~~
}
\value{
(N-1)x(M-1)ベクトル
}
\references{
%% ~put references to the literature/web site here ~
}
\author{
%%  ~~who you are~~
}
\note{
%%  ~~further notes~~
}

%% ~Make other sections like Warning with \section{Warning }{....} ~

\seealso{
}
\examples{
N<-3
M<-4
O<-matrix(sample(10:50,N*M,replace=TRUE),N,M)
tableExpAndMarginals(O)

A<-makeTblSameChiNM(O)
plot3d(A[,1],A[,2],A[,3])
plot(as.data.frame(A))
}