ASI.AMMI: AMMI Stability Index

Description Usage Arguments Details Value References See Also Examples

View source: R/ASI.AMMI.R

Description

ASI.AMMI computes the AMMI Stability Index (ASI) (Jambhulkar et al., 2014; Jambhulkar et al., 2015; Jambhulkar et al., 2017) considering the first two interaction principal components (IPCs) in the AMMI model. Using ASI, the Simultaneous Selection Index for Yield and Stability (SSI) is also calculated according to the argument ssi.method.

Usage

1
ASI.AMMI(model, ssi.method = c("farshadfar", "rao"), a = 1)

Arguments

model

The AMMI model (An object of class AMMI generated by AMMI).

ssi.method

The method for the computation of simultaneous selection index. Either "farshadfar" or "rao" (See SSI).

a

The ratio of the weights given to the stability components for computation of SSI when method = "rao" (See SSI).

Details

The AMMI Stability Index (ASI) is computed as follows:

ASI = √[[PC21 × θ21] + [PC22 × θ22]]

Where, PC1 and PC2 are the scores of 1st and 2nd IPCs respectively; and θ1 and θ2 are percentage sum of squares explained by the 1st and 2nd principal component interaction effect respectively.

Value

A data frame with the following columns:

ASI

The ASI values.

SSI

The computed values of simultaneous selection index for yield and stability.

rASI

The ranks of ASI values.

rY

The ranks of the mean yield of genotypes.

means

The mean yield of the genotypes.

The names of the genotypes are indicated as the row names of the data frame.

References

\insertRef

jambhulkar_ammi_2014ammistability

\insertRef

jambhulkar_genotype_2015ammistability

\insertRef

jambhulkar_stability_2017ammistability

See Also

AMMI, SSI

Examples

 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
26
library(agricolae)
data(plrv)

# AMMI model
model <- with(plrv, AMMI(Locality, Genotype, Rep, Yield, console = FALSE))

# ANOVA
model$ANOVA

# IPC F test
model$analysis

# Mean yield and IPC scores
model$biplot

# G*E matrix (deviations from mean)
array(model$genXenv, dim(model$genXenv), dimnames(model$genXenv))

# With default ssi.method (farshadfar)
ASI.AMMI(model)

# With  ssi.method = "rao"
ASI.AMMI(model, ssi.method = "rao")

# Changing the ratio of weights for Rao's SSI
ASI.AMMI(model, ssi.method = "rao", a = 0.43)

ajaygpb/AMMIStbp documentation built on Jan. 27, 2020, 8:53 a.m.