Description Usage Arguments Details Value Author(s) References Examples
A profile likelihood based method of estimation and hypothesis testing on the correlation coefficient of bivariate data with different types of cencoring.
1 2 3 4 5 6 7 8 9 | clikcorr(data, lower1, upper1, lower2, upper2, cp = 0.95, dist = "n",
df = 4, sv = NA, nlm = FALSE, ...)
## Default S3 method:
clikcorr(data, lower1, upper1, lower2, upper2, cp = 0.95, dist = "n",
df = 4, sv = NA, nlm = FALSE, ...)
## S3 method for class 'clikcorr'
print(x, ...)
## S3 method for class 'clikcorr'
summary(object, ...)
|
data |
a data frame name. |
lower1 |
the lower bound of the first of the two variables whose correlation coefficient to be calculated. |
upper1 |
the upper bound of the first of the two variables whose correlation coefficient to be calculated. |
lower2 |
the lower bound of the second of the two variables whose correlation coefficient to be calculated. |
upper2 |
the upper bound of the second of the two variables whose correlation coefficient to be calculated. |
cp |
confidence level for the confidence interval. |
dist |
working distribution. By default, |
df |
degree of freedom of the bivariate t-distribution when |
sv |
user specified starting values for the vector of (mean1, mean2, var1, corr, var2). |
nlm |
use |
x |
an object of class |
object |
an object of class |
... |
not used. |
clikcorr
conducts point estimation and hypothesis testing on the correlation coefficient of bivariate data with different types of cencoring.
A list with components:
pairName |
variable names for the input paired data structure in the |
pairData |
a paired data structure in the |
dist |
Normal or t distribution. |
df |
degree of freedom for t distribution. |
coefficients |
maximum likelihood estimate (MLE) of the correlation coefficient. |
Cov |
estimated variance covariance matrix. |
Mean |
estimated means. |
CI |
unsymmetric profile confidence interval for the estimated correlation coefficient. |
P0 |
p-value for likelihood ratio test with null hypothesis says that the true correlation coefficient equals zero. |
logLik |
the value of the log likelihood at MLE. |
Yanming Li, Kerby Shedden, Brenda W. Gillespie and John A. Gillespie.
Yanming Li, Kerby Shedden, Brenda W. Gillespie and John A. Gillespie (2016). Calculating Profile Likelihood Estimates of the Correlation Coefficient in the Presence of Left, Right or Interval Censoring and Missing Data.
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 | data(ND)
logND <- log(ND)
logND1 <- logND[51:90,]
obj <- clikcorr(logND1, "t1_OCDD", "t2_OCDD", "t1_HxCDF_234678", "t2_HxCDF_234678")
## Not run:
clikcorr(logND, "t1_OCDD", "t2_OCDD", "t1_HxCDF_234678", "t2_HxCDF_234678")
clikcorr(logND, "t1_OCDD", "t2_OCDD", "t1_HxCDF_234678", "t2_HxCDF_234678",
nlm=TRUE)
clikcorr(logND, "t1_OCDD", "t2_OCDD", "t1_HxCDF_234678", "t2_HxCDF_234678",
method="BFGS")
clikcorr(logND, "t1_OCDD", "t2_OCDD", "t1_HxCDF_234678", "t2_HxCDF_234678",
sv=c(5,-0.5,0.6,0.5,0.6))
clikcorr(logND, "t1_OCDD", "t2_OCDD", "t1_HxCDF_234678", "t2_HxCDF_234678",
dist="t", df=10, nlm=TRUE)
## End(Not run)
print(obj)
summary(obj)
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.