F_MC_powerCircadian_norm: F-statistics-based Monte-Carlo Circadian Power Calculation...

In weiiizong/powerCircadian: An R package for circadian data power calculation

Description

Monte-Carlo circadian power calculation based on F-statistics assuming independent Normal errors

Usage

F_MC_powerCircadian_norm(
  B_MC = 10000,
  n,
  A,
  sigma,
  phi = 0,
  period = 24,
  C = 10,
  cts = NULL,
  alpha = 0.05
)

Arguments

B_MC	Numer of Monte-Carlo simulations. Default is 10000.
n	Sample size.
A	Amplitude of the sine curve: A * sin(2π/period * (phi + cts)) + C.
sigma	\sigma in the independent Normal error term N(0,σ).
phi	Phase of the sine curve. Default is 0.
period	Period of the sine curve. Default is 24.
C	Offset of the sine curve. Default is 10.
cts	Circadian time design vector.
alpha	Type I error control. Default is 0.05.

Details

Calculate power of circaidan data from Monte-Carlo simulated data assuming independent Normal errors.

Value

A vector of empirical type I error and Monte-Carlo power based on F-statistics.

Author(s)

Wei Zong, Zhiguang Huo

Examples

B_MC = 10000
n = 100
cts = seq(0,24,length.out = n+1)[-1]
A = 1
sigma = 1
phi = 0
C = 10
F_MC_powerCircadian_norm(B_MC, n, A, sigma, phi=0, period = 24, cts=cts, alpha = 0.05)

weiiizong/powerCircadian documentation built on Dec. 23, 2021, 5:09 p.m.