Discussion :

[CupOfAppleTea]

Bonjour,

Je dispose d'un jeu de données de 9731 colonnes (individus) et 451 lignes (observations). Chaque colonne représente une même mesure effectuée pour un individu différent. Chaque ligne correspond à une longueur d'onde à laquelle la mesure a été effectuée. J'ai donc pour chaque individu (colonne) un spectre sur le domaine de longueur d'onde étudié (domaine de longueur d'onde identique pour tous les individus).

Je souhaite calculer l'intégrale de la courbe décrite par chaque colonne et obtenir en sortie un tableau d'une ligne et 9731 colonnes avec pour chaque colonne l'intégrale de la courbe correspondante.

J'ai donc tenté cela mais ça ne fonctionne pas, je n'ai en sortie qu'une ligne remplie de NA, excepté pour la dernière colonne où j'ai bien la dernière valeur calculée.

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
c=read.table("DiffCarre.txt", header=TRUE)
 
for ( i in 2:9731) {
	AUC=trapz(c[,1],c[,i])
        c2[1,i]=AUC
}

Merci pour votre aide !

[CupOfAppleTea]

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
c2=matrix(ncol=9730,nrow=1)
 
for ( i in 2:9731) {
	AUC=trapz(c[,1],c[,i])
	c2[1,i-1]=AUC
}
c2

Voila!