Error in BEEFPRODACT[ENDDAY[j], j] : subscript out of bounds

[Chris1992]

Hi I am using LiGAPS beef model to analyze my data. The model is long and coded in R. The following error is preventing the whole code from running. May you assist me, I keep getting Error in BEEFPRODACT[ENDDAY[j], j] : subscript out of bounds

[nirgrahamuk]

j is your subscript
A number indicating which element from a set of elements is being pulled out. Out of bounds means trying to pull out an element that is not there.
Like if I tell you have a box with 7 apples, and you ask for the 8th apple, your request is out of bounds.

[Chris1992]

So how do i deal with that, in the code J is coming from here
for (j in 1:jmax){ # Loop for individual animals starts here (j = jth animal)
And the jmax code is the one below.
if(SCALE== 1) SEX <- SEX_ANIMAL else SEX <- c(1,SEX_CALVES)
jmax <- MAXCALFNR + 1 # Number of animals in the simulation at the animal level or the herd
# level
MAXCALFNR is 8 meaning eventully my jmax value is 9. So where do i edit now and how

[nirgrahamuk]

So then the question becomes whether ENDDAY has 9 things and whether beefproduct has 9 in its 2nd dimension, and if all the values of ENDDAY from 1 to 9 are themselves present in the first dimension of beefproduct

[Chris1992]

here is the ENDDAY code with 9 things
ENDDAY = c(1,1,1,1,1,1,1,1,1)

Another related code hereunder

Day the animal is slaughtered (days)

# Is initially 9999, but is replaced by the correct number of days at slaughter
ENDDAY[j] <- min(SLAUGHTERDAYACT[1:i,j]) 
  
# Beef production at slaughter (kg per animal)
if(ENDDAY[j] < 9999) BEEFPROD[j] <- BEEFPRODACT[ENDDAY[j],j] 
# Beef production (kg beef per head per year)
if(ENDDAY[j] < 9999) BEEFPRODYEAR[j] <- BEEFPRODACT[ENDDAY[j],j]/(TIME[ENDDAY[j],j]/365) 
# Live weight production at slaughter (kg per animal)
if(ENDDAY[j] < 9999) LWPROD[j] <- LWPRODACT[ENDDAY[j],j] 
# Live weight production (kg live weight per head per year)
if(ENDDAY[j] < 9999) LWPRODYEAR[j] <- LWPRODACT[ENDDAY[j],j]/(TIME[ENDDAY[j],j]/365) 

The model is very long with many codes. Would you be available to assist me via an online short meeting where i can share my screen, please contact me via direct message.