| Type: | Package |
| Title: | Neutrosophic Analysis of Completely Randomized Designs and Randomized Complete Block Designs |
| Version: | 0.0.1 |
| Maintainer: | Vinaykumar L.N. <vinaymandya123@gmail.com> |
| Description: | Provides neutrosophic statistical methods for interval-valued data from completely randomized and randomized complete block designs. Methods include neutrosophic analysis of variance, analysis of covariance, multivariate analysis of variance, pooled analysis of variance, Levene's test, and Aitken transformation. When the lower and upper bounds are equal (crisp data), the methods reduce to their corresponding classical statistical analyses. The basic concept of neutrosophic statistics is based on Smarandache (2014) https://fs.unm.edu/NeutrosophicStatistics.pdf, while the statistical analysis procedures implemented in this package are newly developed. |
| License: | GPL (≥ 3) |
| Encoding: | UTF-8 |
| Imports: | MASS, stats |
| RoxygenNote: | 7.3.3 |
| NeedsCompilation: | no |
| Packaged: | 2026-08-05 04:16:59 UTC; admin |
| Author: | Neethu R.S. [aut, ctb], Boyina Devi Priyanka [aut, ctb], Cini Varghese [aut, ctb], Susheel Kumar Sarkar [aut, ctb], Mohd Harun [aut, ctb], Vinaykumar L.N. [aut, cre], Anindita Datta [aut, ctb] |
| Repository: | CRAN |
| Date/Publication: | 2026-08-09 08:40:13 UTC |
Neutrosophic Analysis of Covariance for Completely Randomized Design
Description
Performs Neutrosophic Analysis of Covariance (NANCOVA) for interval-valued response and covariate data from a Completely Randomized Design (CRD).
Usage
CRDnsANCOVA(
Lower_y,
Upper_y,
Lower_z,
Upper_z,
design,
alpha = 0.05,
verbose = FALSE
)
Arguments
Lower_y |
Numeric matrix of lower bounds of the response variable. |
Upper_y |
Numeric matrix of upper bounds of the response variable. |
Lower_z |
Numeric matrix of lower bounds of the covariate. |
Upper_z |
Numeric matrix of upper bounds of the covariate. |
design |
Numeric matrix representing the CRD treatment layout. |
alpha |
Significance level for interval-based LSD comparisons. Default is 0.05. |
verbose |
Logical. If |
Value
A list containing the Neutrosophic ANCOVA table, interval LSD comparisons (if treatment effects are significant), and the interval LSD value.
Examples
Lower_y <- matrix(c(
20,25,30,35,
21,26,31,36,
19,24,29,34,
20,25,30,35,
21,26,31,36
), nrow = 5, byrow = TRUE)
Upper_y <- matrix(c(
22,27,32,37,
23,28,33,38,
21,26,31,36,
22,27,32,37,
23,28,33,38
), nrow = 5, byrow = TRUE)
Lower_z <- matrix(c(
8,9,10,11,
9,10,11,12,
8,9,10,11,
9,10,11,12,
8,9,10,11
), nrow = 5, byrow = TRUE)
Upper_z <- matrix(c(
9,10,11,12,
10,11,12,13,
9,10,11,12,
10,11,12,13,
9,10,11,12
), nrow = 5, byrow = TRUE)
design <- matrix(c(
1,2,3,4,
1,2,3,4,
1,2,3,4,
1,2,3,4,
1,2,3,4
), nrow = 5, byrow = TRUE)
CRDnsANCOVA(Lower_y, Upper_y, Lower_z, Upper_z, design)
Neutrosophic Analysis of Variance for Completely Randomized Design
Description
Performs neutrosophic analysis of variance (NANOVA) for a completely randomized design using interval-valued observations.
Usage
CRDnsANOVA(Lower_y, Upper_y, design, alpha = 0.05, verbose = FALSE)
Arguments
Lower_y |
Matrix of lower bounds of the response variable. |
Upper_y |
Matrix of upper bounds of the response variable. |
design |
Matrix specifying treatment allocation. |
alpha |
Significance level for LSD test. |
verbose |
Logical. If TRUE, prints the analysis. |
Value
A list containing the NANOVA table, treatment means, pairwise comparisons, and LSD interval.
Examples
Lower_y <- matrix(c(
9.5, 19.5, 29.5, 39.5,
10.0, 20.0, 30.0, 40.0,
10.5, 20.5, 30.5, 40.5,
9.8, 19.8, 29.8, 39.8,
10.2, 20.2, 30.2, 40.2
), nrow = 5, byrow = TRUE)
Upper_y <- matrix(c(
10.5, 20.5, 30.5, 40.5,
11.0, 21.0, 31.0, 41.0,
11.5, 21.5, 31.5, 41.5,
10.8, 20.8, 30.8, 40.8,
11.2, 21.2, 31.2, 41.2
), nrow = 5, byrow = TRUE)
design <- matrix(c(
1, 2, 3, 4,
1, 2, 3, 4,
1, 2, 3, 4,
1, 2, 3, 4,
1, 2, 3, 4
), nrow = 5, byrow = TRUE)
CRDnsANOVA(
Lower_y = Lower_y,
Upper_y = Upper_y,
design = design,
alpha = 0.05,
verbose = TRUE
)
Neutrosophic MANOVA for a Completely Randomized Design
Description
Performs neutrosophic multivariate analysis of variance (MANOVA) for interval-valued multivariate observations under a completely randomized design using Wilks' Lambda as the multivariate test criterion.
Usage
CRDnsMANOVA(Lower_y, Upper_y, group, alpha = 0.05, verbose = TRUE)
Arguments
Lower_y |
A numeric matrix containing the lower bounds of the response variables. Rows represent observations and columns represent response variables. |
Upper_y |
A numeric matrix containing the upper bounds of the response
variables. Must have the same dimensions as |
group |
A vector or factor identifying the treatment group for each
observation. Its length must equal |
alpha |
Significance level used to obtain the critical F value.
Defaults to |
verbose |
Logical; if |
Value
A list containing the results of the neutrosophic MANOVA with the following components:
- MANOVA.table
A data frame containing the neutrosophic MANOVA table, including the source of variation, degrees of freedom, Wilks' Lambda interval, approximate F-statistic interval, and decision.
- Wilks.Lambda
A named numeric vector containing the lower and upper bounds of Wilks' Lambda.
- Approx.F
A named numeric vector containing the lower and upper bounds of the approximate F-statistic.
- F.critical
The critical F value at the specified significance level.
- Decision
The hypothesis testing decision ("Significant", "Not Significant", or "Indeterminate").
- treatment_ss
A list containing the lower, upper, determinate, and indeterminate treatment sum of squares matrices.
- total_ss
A list containing the lower, upper, determinate, and indeterminate total sum of squares matrices.
- residual_ss
A list containing the lower, upper, determinate, and indeterminate residual sum of squares matrices.
Examples
YL <- matrix(c(8.27, 2.50, 5.30, 1.41,
8.40, 6.27, 0.00, 3.42),
ncol = 2, byrow = TRUE)
YU <- matrix(c(9.73, 3.50, 6.70, 2.59,
9.60, 7.73, 0.00, 4.58),
ncol = 2, byrow = TRUE)
g <- c(1, 1, 2, 2)
fit <- CRDnsMANOVA(
Lower_y = YL,
Upper_y = YU,
group = g
)
fit$MANOVA.table
fit$Decision
Pooled Neutrosophic ANOVA for a Completely Randomized Design
Description
Performs pooled neutrosophic analysis of variance for a completely randomized design.
Usage
PooledCRDnsANOVA(Lower_y, Upper_y, design, alpha = 0.05, verbose = FALSE)
Arguments
Lower_y |
Numeric matrix of lower response bounds. Columns must be arranged in consecutive year or location blocks, each containing all treatments. |
Upper_y |
Numeric matrix of upper response bounds. |
design |
Numeric treatment-layout matrix with the same dimensions as
|
alpha |
Significance level. Default is |
verbose |
Logical; print results when |
Value
A list containing the initial and final NANOVA tables, interaction test, decision, and pooled MSE when applicable.
Examples
# Three replications, three treatments, and two years.
# Columns 1-3 are Year 1; columns 4-6 are Year 2.
Lower_y <- matrix(c(
20, 25, 30, 21, 26, 31,
22, 27, 32, 23, 28, 33,
21, 26, 31, 22, 27, 32
), nrow = 3, byrow = TRUE)
Upper_y <- matrix(c(
22, 27, 32, 23, 28, 33,
24, 29, 34, 25, 30, 35,
23, 28, 33, 24, 29, 34
), nrow = 3, byrow = TRUE)
# Treatment layout:
# Year 1 = Treatments 1, 2, 3
# Year 2 = Treatments 1, 2, 3
design <- matrix(c(
1, 2, 3, 1, 2, 3,
1, 2, 3, 1, 2, 3,
1, 2, 3, 1, 2, 3
), nrow = 3, byrow = TRUE)
result <- PooledCRDnsANOVA(
Lower_y = Lower_y,
Upper_y = Upper_y,
design = design,
alpha = 0.05,
verbose = TRUE
)
result$NANOVA.table
Pooled Neutrosophic Analysis of Variance for a Randomized Block Design
Description
Performs the pooled neutrosophic ANOVA calculation for interval-valued observations measured in a Randomized Block Design over multiple years or locations or seasons.
Usage
PooledRBDnsANOVA(Lower_y, Upper_y, design, alpha = 0.05, verbose = FALSE)
Arguments
Lower_y |
Numeric matrix of lower response bounds. Columns must be ordered in consecutive year blocks, each containing all treatments. |
Upper_y |
Numeric matrix of upper response bounds, with the same
dimensions as |
design |
Numeric matrix giving the treatment labels for |
alpha |
Significance level for the treatment-by-year interaction test. |
verbose |
Logical; print the neutrosophic ANOVA table when |
Value
An invisible list containing the neutrosophic ANOVA table, interaction F interval and critical value, interaction decision, and the mean-square interval used for subsequent F intervals.
Examples
# Three blocks, three treatments, and two years.
# Columns 1-3 are Year 1; columns 4-6 are Year 2.
# Each block contains every treatment once in each year.
Lower_y <- matrix(c(
20, 25, 30, 22, 28, 23,
27, 32, 21, 33, 24, 29,
31, 22, 26, 25, 30, 34
), nrow = 3, byrow = TRUE)
Upper_y <- matrix(c(
22, 27, 32, 24, 30, 25,
29, 34, 23, 35, 26, 31,
33, 24, 28, 27, 32, 36
), nrow = 3, byrow = TRUE)
# Treatment layout for each block:
# Year 1 occupies columns 1-3.
# Year 2 occupies columns 4-6.
design <- matrix(c(
1, 2, 3, 2, 3, 1,
2, 3, 1, 3, 1, 2,
3, 1, 2, 1, 2, 3
), nrow = 3, byrow = TRUE)
result <- PooledRBDnsANOVA(
Lower_y = Lower_y,
Upper_y = Upper_y,
design = design,
alpha = 0.05,
verbose = TRUE
)
result$NANOVA.table
Neutrosophic Analysis of Covariance for RCBD
Description
Performs neutrosophic analysis of covariance for a Randomized Complete Block Design using interval-valued response and covariate data.
Usage
RCBDnsANCOVA(
Lower_y,
Upper_y,
Lower_z,
Upper_z,
design,
alpha = 0.05,
verbose = FALSE
)
Arguments
Lower_y |
Numeric matrix of lower response values. |
Upper_y |
Numeric matrix of upper response values. |
Lower_z |
Numeric matrix of lower covariate values. |
Upper_z |
Numeric matrix of upper covariate values. |
design |
Numeric matrix representing the RCBD layout. |
alpha |
Significance level for LSD test. Default is |
verbose |
Logical; if |
Value
A list containing the NANCOVA table, LSD interval, and treatment comparisons (if treatment effect is significant).
Examples
Lower_y <- matrix(c(
46.84, 97.69, 39.30, 49.20,
51.52, 83.38, 56.48, 42.15,
44.42, 49.64, 44.21, 35.79,
32.77, 74.30, 55.72, 70.55
), nrow = 4, byrow = TRUE)
Upper_y <- matrix(c(
49.66,101.31,47.70,55.80,
60.98,87.62,60.52,44.85,
52.26,59.36,52.79,44.61,
42.23,82.70,59.28,75.45
), nrow = 4, byrow = TRUE)
Lower_z <- matrix(c(
224.95,245.87,245.19,259.41,
222.06,213.83,253.10,247.51,
255.67,230.55,265.05,246.67,
137.98,214.68,251.49,257.70
), nrow = 4, byrow = TRUE)
Upper_z <- matrix(c(
229.05,250.13,252.81,268.59,
229.94,222.17,258.90,256.49,
262.33,237.45,274.95,249.33,
142.02,219.32,260.51,264.30
), nrow = 4, byrow = TRUE)
design <- matrix(c(
1,2,3,4,
1,2,3,4,
1,2,3,4,
1,2,3,4
), nrow = 4, byrow = TRUE)
result <- RCBDnsANCOVA(
Lower_y,
Upper_y,
Lower_z,
Upper_z,
design,
alpha = 0.05,
verbose = TRUE
)
Neutrosophic Analysis of Variance for Randomized Complete Block Design
Description
Performs Neutrosophic Analysis of Variance (NANOVA) for interval-valued response data from a Randomized Complete Block Design (RCBD).
Usage
RCBDnsANOVA(Lower_y, Upper_y, design, alpha = 0.05, verbose = FALSE)
Arguments
Lower_y |
Numeric matrix of lower bounds of the response variable. |
Upper_y |
Numeric matrix of upper bounds of the response variable. |
design |
Numeric matrix representing the RCBD layout. |
alpha |
Significance level for interval-based LSD test. Default is 0.05. |
verbose |
Logical. If |
Value
A list containing the Neutrosophic ANOVA table, interval-based LSD comparisons (if applicable), and the interval LSD.
Examples
Lower_y <- matrix(c(
120.230,125.488,132.987,127.086,127.672,128.013,
122.594,121.009,123.969,120.358,120.424,122.197,
121.183,130.671,128.794,114.863,122.595,122.073,
127.620,124.532,132.893,125.528,125.850,127.550
), nrow = 4, byrow = TRUE)
Upper_y <- matrix(c(
127.6967536,131.2116955,141.2127373,136.1540904,130.6884772,136.8474149,
129.8264289,130.3314544,133.3113414,126.5063118,128.4362999,130.2714433,
124.5068016,139.3287297,134.1060197,124.2774447,127.2248520,130.3948469,
131.0638721,129.8884785,135.5666716,127.7580663,132.0178679,133.3903886
), nrow = 4, byrow = TRUE)
design <- matrix(c(
1,2,3,4,5,6,
1,2,3,4,5,6,
1,2,3,4,5,6,
1,2,3,4,5,6
), nrow = 4, byrow = TRUE)
RCBDnsANOVA(Lower_y, Upper_y, design)
Aitken Transformation for Neutrosophic Interval Data
Description
Applies an Aitken variance-stabilizing transformation to interval-valued observations using group-specific interval MSE values.
Usage
nsAitkenTransform(Lower_y, Upper_y, group, MSE_lower, MSE_upper)
Arguments
Lower_y |
Numeric matrix of lower interval bounds. |
Upper_y |
Numeric matrix of upper interval bounds. |
group |
Vector identifying the group/year for every row of the input. |
MSE_lower |
Numeric vector of lower MSE bounds, one per group. |
MSE_upper |
Numeric vector of upper MSE bounds, one per group. |
Value
A list containing transformed lower/upper matrices and their determinate and indeterminate components.
Examples
Lower_y <- matrix(
c(8.2, 5.1,
7.9, 5.4,
9.0, 6.2,
8.8, 6.0,
7.5, 4.8,
7.8, 5.0,
8.6, 5.8,
8.4, 5.6),
ncol = 2,
byrow = TRUE
)
Upper_y <- matrix(
c(8.8, 5.7,
8.5, 6.0,
9.6, 6.8,
9.4, 6.6,
8.1, 5.4,
8.4, 5.6,
9.2, 6.4,
9.0, 6.2),
ncol = 2,
byrow = TRUE
)
group <- rep(1:4, each = 2)
transformed <- nsAitkenTransform(
Lower_y = Lower_y,
Upper_y = Upper_y,
group = group,
MSE_lower = c(29.83, 16.30, 21.21, 3.62),
MSE_upper = c(30.02, 16.43, 21.24, 3.89)
)
transformed$weights
Neutrosophic Levene's Test for Interval-Valued Data
Description
Tests homogeneity of variances for interval-valued observations using neutrosophic lower and upper test-statistic bounds.
Usage
nsLeveneTest(Lower_y, Upper_y, design, alpha = 0.05, verbose = TRUE)
Arguments
Lower_y |
Numeric matrix of lower interval bounds. |
Upper_y |
Numeric matrix of upper interval bounds; same dimensions as
|
design |
Treatment-design matrix. Treatment labels in each row must be arranged in consecutive year blocks. |
alpha |
Significance level. Default is |
verbose |
Logical; print the test table when |
Value
An invisible list containing the Levene statistic interval, critical F value, homogeneity decision, and ANOVA table.
Examples
# Two replications, three treatments, and two years.
# Columns are ordered as:
# Year 1 (T1, T2, T3), Year 2 (T1, T2, T3).
Lower_y <- matrix(c(
10, 12, 11, 13, 14, 12,
9, 11, 10, 12, 13, 11
), nrow = 2, byrow = TRUE)
Upper_y <- matrix(c(
12, 14, 13, 15, 16, 14,
11, 13, 12, 14, 15, 13
), nrow = 2, byrow = TRUE)
design <- matrix(c(
1, 2, 3, 1, 2, 3,
1, 2, 3, 1, 2, 3
), nrow = 2, byrow = TRUE)
levene_result <- nsLeveneTest(
Lower_y = Lower_y,
Upper_y = Upper_y,
design = design,
alpha = 0.05,
verbose = TRUE
)
levene_result$Homogeneity