For循環錯誤 - 下標出界

Question

的我一直在尋找了幾個小時的互聯網試圖解決以下錯誤：

Error in Dataset[i, Year] : subscript out of bounds 

下面是我的代碼產生錯誤的部分：

for(i in (2*YF):1){ 
    if(Dataset[i,Year] < 0){ 
     Dataset[i,Total_Births] <- Dataset[i,Male_Births] + Dataset[i,Female_Births] 
    }else{ 
     Dataset[i,Total_Births] <- with(Dataset, sum(Dataset[Female_Births > (i-AEB) & Female_Births <= (i-ABB),Female_Births])) 
     Dataset[i,Male_Births] <- MBR * Dataset[i,Total_Births] 
     Dataset[i,Female_Births] <- FBR * Dataset[i,Total_Births] 
    } 
} 

'Year'開始於500，結束於-500跳過0.

我需要從下往上填充列。

任何幫助將不勝感激。

以下是完整代碼：

#DECLARE PARAMETERS 
YF <- 500  #Years Ago Flood Occurred 
FBR <- .5  #Historical Female Birth Rate 
MBR <- (1-FBR) #Historical Male Birth Rate 
ABF <- 12  #Average Births per Female 
AL <- 60  #Average Lifespan 
NF <- 4  #Number of Females at Debark 
NM <- 4  #Number of Males at Debark 
NAD <- 40  #Noah's Age at Debark 
NWA <- 40  #Emzara's (Noah's Wife) Age at Debark 
SA <- 35  #Shem's (Noah's Son) Age at Debark 
SWA <- 35  #Sedeqetelebab's (Shem's Wife) Age at Debark 
HA <- 30  #Ham (Noah's Son) Age at Debark 
HWA <- 30  #Ne'elatama'uk's (Ham's Wife) Age at Debark 
JA <- 25  #Japheth (Noah's Son) Age at Debark 
JWA <- 25  #Adataneses's (Japheth's Wife) Age at Debark 
CP <- 20  #Current Population of World Today 
ABB <- 18  #Age Begin Births 
AEB <- 30  #Age End Births 


###CREATE MATRIX - YEAR### 
YearA <- YF:1 
YearB <- -1:-YF 

Year <- c(YearA,YearB) 

###POPULATE BIRTHS DATA### 
BF_Male_Births <- rep.int(0,YF) 
BF_Female_Births <- rep.int(0,YF) 

BF_Male_Births[NAD] <- 1 
BF_Female_Births[NWA] <- 1 
BF_Male_Births[SA] <- 1 
BF_Female_Births[SWA] <- 1 
BF_Male_Births[HA] <- 1 
BF_Female_Births[HWA] <- 1 
BF_Male_Births[JA] <- 1 
BF_Female_Births[JWA] <- 1 

BF_Total_Births <- rep.int(0,YF) 

AF_Male_Births <- rep.int(0,YF) 
AF_Female_Births <- rep.int(0,YF) 
AF_Total_Births <- c(rep.int(4,9),rep.int(3,3),rep.int(0,(YF-12))) #make parameters for 9,3,12 

Male_Births <- c(AF_Male_Births,BF_Male_Births) 
Female_Births <- c(AF_Female_Births,BF_Female_Births) 
Total_Births <- c(AF_Total_Births,BF_Total_Births) 

Total_Births <- rep.int(0,YF) 

Births <- cbind(Male_Births,Female_Births,Total_Births) 
Dataset <- cbind(Year,Births) 

for(i in (2*YF):1){ 
    if(Dataset[i,Year] < 0){ 
    Dataset[i,Total_Births] <- Dataset[i,Male_Births] + Dataset[i,Female_Births] 
    }else{ 
    Dataset[i,Total_Births] <- with(Dataset, sum(Dataset[Female_Births > (i-AEB) & Female_Births <= (i-ABB),Female_Births])) 
    Dataset[i,Male_Births] <- MBR * Dataset[i,Total_Births] 
    Dataset[i,Female_Births] <- FBR * Dataset[i,Total_Births] 
    } 
} 

Answer 1

當你這樣做Dataset[i,Year]，既i和Year被視爲indicies。由於year是數字，所以子集操作符假定您希望返回這些列值。但是，Year是一個值爲500：-500的向量，並且只有4列，這就是爲什麼您會遇到越界錯誤。這有點令人困惑，因爲你在矩陣外部定義了一個Year矢量，在矩陣內部有一個Year列。

如果使用字符向量進行索引，它將匹配行名稱和列名稱。所以你認爲你的意思是

Dataset[i,"Year"] 

所有其他地方，你索引Dataset。現在，沒有引號，它使用你在矩陣之外定義的值，這不是你想要的值。

另外，您不能在矩陣中使用with()。 with()應該與data.frames或環境或列表一起使用，並且Dataset不是這些中的任何一個。這更接近於正確的語法

for(i in (2*YF):1){ 
    if(Dataset[i,"Year"] < 0){ 
    Dataset[i,"Total_Births"] <- Dataset[i,"Male_Births"] + Dataset[i,"Female_Births"] 
    }else{ 
    Dataset[i,"Total_Births"] <- sum(Dataset[Dataset[, "Female_Births"] > (i-AEB) & Dataset[,"Female_Births"] <= (i-ABB),"Female_Births"]) 
    Dataset[i,"Male_Births"] <- MBR * Dataset[i,"Total_Births"] 
    Dataset[i,"Female_Births"] <- FBR * Dataset[i,"Total_Births"] 
    } 
} 

但似乎你可能仍然存在一些邏輯錯誤。