Применение функции в ddply

Question

Я пытаюсь применить групповую функцию, используя ddply. Мои данные выглядит следующим образом:

https://echange-fichiers.inra.fr/get?k=b1jD63CWkT93hDsbZ0g

Следующая функция вычисляет значение х при у = 0

intercept=function (x,y){ 

    if (length(y[y==0])==0){ 
     x1=min(x[y >0]) 
     y1=min(y[y >0]) 

     x2=max(x[y <0]) 
     y2=max(y[y <0]) 

     p =(y2-y1)/(x2-x1) 

     x0= -y1/p + x1 

    }else{ 
     x0=x[y==0] 

    } 
return(x0) 
} 

тогда я пытаюсь применить это группой с ddply:

zzz=ddply(Data,.(Genotype,T), summarise 
     InterDays=(intercept(Data$Days,Data$Diff)) 
) 

Моя грамматика должна быть неправильной, поскольку она не возвращает результат с аналогичным форматом, чем:

zzz=ddply(DataAll,.(Genotype,T), summarise 
     InterDays=mean(DataAll$Days) 
) 

Я думаю, что что-то не так с тем, как я называю свои переменные функции, но я не вижу, где ... Любая помощь приветствуется. Все самое лучшее, Винсент

Answer 1

Я думаю, что вы должны назвать ddply таким образом:

ddply(Data,.(Genotype,T), summarise, 
     InterDays=(intercept(Days,Diff)) 
) 

IIUC, если вы обратитесь к функции перехвата с Data$Days и Data$Diff в качестве аргументов, целые векторы будут переданы в функция каждый раз. С Days и Diff используются только значения, соответствующие текущему подмножеству Data.