comment calculer intégrale avec R

Question

Ici, je voudrais avoir votre aide sur la mise en œuvre du calcul de l'intégrale sur deux vecteurs. J'ai vérifié les pages sur le calcul intégral par rapport à R. Mais, j'ai peu de formation sur les mathématiques, donc je ne peux toujours pas le faire moi-même.

Mes objectifs est de mettre en œuvre l'idée de cette phrase: « Si vous indiquerez les estimations du taux de la position, la carte génétique est juste l'intégrale de ce complot. » Cela signifie que j'ai des variables (taux, positions), chaque position a son propre taux. Je veux calculer l'intégrale des taux pour chaque position. Ici, la position augmente de façon monotone.

Cette tâche ne devrait pas être si complexe pour ceux qui ont de bons antécédents en calcul mathématique. Alors, pourriez-vous s'il vous plaît me donner des instructions/instructions à ce sujet?

Merci d'avance.

# here I make dummy data 

position <- c(2,34,58) 
rate <- c(14, 20, 5) 

Answer 1

En mathématiques, une intégrale est l'aire sous la courbe. Dans votre exemple, vous voulez que l'aire sous la courbe soit définie par la position et le taux.

position <- c(2,34,58) 
rate <- c(14, 20, 5) 

plot(position, rate, type="l", ylim=c(0, 25)) 

Vous pouvez calculer la surface sous la courbe à la main, en utilisant la règle trapézoïdale:

32*17 + 24*12.5 = 844 

Ou, pour le faire programme:

AUC <- function(x, y){ 
    sum(diff(x)*rollmean(y,2)) 
} 

AUC(position, rate) 
[1] 844