用戶定義函數內部和外部的不同結果，爲什麼？

Question

我有一個計算索引'MCI'的函數。代碼是

mci<-function(Data){ 
    N<-dim(Data)[1] 
    L<-dim(Data)[2] 
    r <- rowSums(Data)  
    i.sc <- colSums(Data)  
    r.matrix <- matrix(r,N,1) %*% matrix(1,1,L) 
    p.cor.i <- (i.sc/N)     
    p.cor.i.matrix <- t(matrix(p.cor.i,L,N)) 
    gutt <- r.matrix - t(matrix(1:L,L,N)) 
    gutt<<- ifelse(gutt<0,0,1) 
    antigutt <- (L-r.matrix) - t(matrix(1:L,L,N)) 
    antigutt<<- ifelse(antigutt<0,1,0) 
    Covgp<-diag(cov(t(gutt),t(p.cor.i.matrix))) 
    Covdp<-diag(cov(t(Data),t(p.cor.i.matrix))) 
    Covagp<-diag(cov(t(antigutt),t(p.cor.i.matrix))) 
    MCI <- (Covgp-Covdp)/(Covgp-Covagp) 
    return(MCI) 
} 

上的數據集是這樣的：

V1 V2 V3 V4 V5 V6 
1 1 1 1 1 0 1 
2 0 0 0 1 0 0 
3 1 1 0 1 1 1 
4 1 1 1 1 0 1 
5 1 0 1 0 0 0 
6 1 1 1 1 1 1 
7 1 1 1 0 1 1 
8 0 1 1 1 1 1 
9 1 1 0 1 1 1 
10 0 1 1 0 0 0 

它返回：

[1] Inf Inf Inf Inf Inf Inf Inf Inf Inf Inf 

但如果我計算函數的每個部分獨立地結果將是：

[1] -Inf NaN NaN -Inf 0.3333333 NaN NaN NaN NaN 0.0000000 

爲什麼發生這種情況，我如何在將來的編程中阻止這種情況？

Answer 1

這些線的問題：

gutt<<- ifelse(gutt<0,0,1) 

antigutt<<- ifelse(antigutt<0,1,0) 

你分配給在封閉環境中的對象，而不是修改函數的環境中的對象（名稱相同）。當在這些行中的第一行中分配給gutt時，它不影響函數中的gutt，然後您稍後使用幾行。

變化這兩行改爲：

gutt<- ifelse(gutt<0,0,1) 

antigutt<- ifelse(antigutt<0,1,0) 

和功能與運行代碼行由行一致。