Как получить обратную линейную интерполяцию в R

Question

Я прошел множество решений для этой проблемы, но не смог найти надежную автоматическую. Ниже вы найдете подробное самодостаточное описание.

Эти данные: data.txt

x y 
1 1 
2 2 
3 3 
4 2 
5 1 

Заговор как график рассеяния:

t=read.table("data.txt",header=T) 
plot(t$x,t$y,"l") 

Вы увидите пик, мой вопрос теперь: если предположить, я доволен линейная интерполяция, какова ширина в половине максимума этой «кривой»? Итак, для каких значений x0 из x I f (x0) = max (y)/2, где f - линейная интерполяция. Я пробовал с approxfun и некоторой плотностью ядра, но я не хочу сглаживать свои данные.

Любой вход очень приветствуется.

Answer 1

Есть, вероятно, много лучших способов сделать это, но вот одно решение. Было бы легче, если бы мы знали, как вы сделали интерполяцию в первую очередь.

# Extend your line for testing 
y<-c(1,2,3,2,1,2,3,4,5,4,3) 

# Split into linear segments. 
segments<-c(1,diff(diff(y)))!=0 
seg.points<-which(c(segments,TRUE)) 

# Calculate desired value 
val<-max(y)/2 

# Loop through and find value 
res<-c() 
for(i in 1:(length(seg.points)-1)) { 
    start<-y[seg.points[i]] 
    end<-y[seg.points[i+1]] 
    if ((val>=start & val<=end) | (val<=start & val >=end)) { 
    slope=(end-start)/(seg.points[i+1] - seg.points[i]) 
    res<-c(res,seg.points[i] + ((val - start)/slope)) 
    } 
} 
res 
# [1] 2.5 3.5 6.5