2017-12-03 104 views
0

我試圖使用外部函數獲取3D表面。並且我定義的函數在data.frame中引用列。以下是可重現的簡要代碼。使用函數包括數據幀引用的外部函數錯誤

的test.txt是

a b c 
    1 0 0 
    2 0 0 
    3 0 0 
    6 0 0 
    4 0 0 

和代碼顯示爲錯誤

test<- read.table(file="test.txt",header=TRUE) 

testf<-function(q,w){ 
test$b<-q 
test$c<-w 
sum(test$b)+sum(test$c) 
} 

x<-seq(1,10,by = 1) 
y<-seq(0,30,by = 1) 

z<-outer(x,y,testf) 

結果是

Error in `$<-.data.frame`(`*tmp*`, "b", value = c(1, 2, 3, 4, 5, 6, 7, : 
replacement has 310 rows, data has 5 

5.stop(sprintf(ngettext(N, "replacement has %d row, data has %d", 
"replacement has %d rows, data has %d"), N, nrows), domain = NA) 

4.`$<-.data.frame`(`*tmp*`, "b", value = c(1, 2, 3, 4, 5, 6, 7, 
8, 9, 10, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 1, 2, 3, 4, 5, 6, 7, 
8, 9, 10, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 1, 2, 3, 4, 5, 6, 7, 
8, 9, 10, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 1, 2, 3, 4, 5, 6, 7, ... 

3.`$<-`(`*tmp*`, "b", value = c(1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 
1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 
1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 
1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, ... 

2. FUN(X, Y, ...) 

1. outer(x, y, testf) 

功能本身似乎因爲testf(1,1)的作品沒有問題

> testf(1,1) 
[1] 10 

任何意見將不勝感激!

+0

test'你定義了'5'行,但你想用'10'行數字矢量替換它的data.frame'的方式。這正是你想要做的嗎?因爲行數應該相等。 – SirSaleh

+0

是的,我要做的是獲取每個用戶定義的函數結果,這是通過更改列數據基於數據幀計算出來的(在實際代碼中,我使用MAPLLY從數據框列中獲取輸入數據)。只有LOOP才能獲得z軸結果? –

+0

或者還有其他的功能與外部功能類似嗎? –

回答

0

我覺得for loop只產生Z軸值