j'ai cherché sur le forum mais il semble que je n'ai pas trouvé ma réponse. Ma question est assez simple mais j'ai toujours du mal avec les boucles ...
J'étudie de la survie et j'ai une courbe pour une relation entre une survie et un indice de condition corporelle de mes individus (= BCI et j'ai 4405 individus donc 4405 valeurs).
Voici un exemple des données que j'ai :
Code : Tout sélectionner
> test <- cbind(BCI,S_logit,S,phi)
> test[1:10,]
BCI S_logit S phi
[1,] -1.1958264 -1.1707333 0.2367225 0.3101394
[2,] -1.1723595 -1.1317988 0.2438293 0.3224527
[3,] -1.1408950 -1.0795952 0.2535826 0.3397330
[4,] -1.1223771 -1.0488717 0.2594418 0.3503328
[5,] -1.0766264 -0.9729654 0.2742898 0.3779606
[6,] -1.0378945 -0.9087043 0.2872651 0.4030461
[7,] -1.0283422 -0.8928558 0.2905208 0.4094847
[8,] -0.9934615 -0.8349843 0.3025922 0.4338813
[9,] -0.9845606 -0.8202167 0.3057177 0.4403362
[10,] -0.9825113 -0.8168165 0.3064398 0.4418360
Je voudrais maintenant créer une matrice regroupant la valeur de 2 dérivés partielles pour chaque valeur de BCI
j'ai fait une formule pour calculer chaque dérivée
Code : Tout sélectionner
DerPart_1 <- phi/(1+phi) - (phi^2)/((1+phi)^2)
DerPart_2 <- (BCI*phi)/(1+phi) - (BCI*(phi^2))/((1+phi)^2)
Et je souhaiterai accoler ces 2 valeurs dans une matrice 1 ligne 2 col car j'dois ensuite faire du calcul matriciel.
J'ai essayé naivement ceci:
Code : Tout sélectionner
DerPart <- matrix(NA, nrow = 1, ncol = 2)
for (i in 1:4405){
DerPart[i] <- c(DerPart_1[i],DerPart_2[i], nrow=1,ncol=2)
}
Mais bien sur, ca ne fonctionne pas ^^
j'obtiens le message d'erreur suivant:
Code : Tout sélectionner
> for (i in 1:4405){
+ DerPart[i] <- c(DerPart_1[i],DerPart_2[i], nrow=1,ncol=2)
+ }
Il y a eu 50 avis ou plus (utilisez warnings() pour voir les 50 premiers)
> warnings()
Messages d'avis :
1: In DerPart[i] <- c(DerPart_1[i], DerPart_2[i], nrow = 1, ... :
le nombre d'objets à remplacer n'est pas multiple de la taille du remplacement
2: In DerPart[i] <- c(DerPart_1[i], DerPart_2[i], nrow = 1, ... :
le nombre d'objets à remplacer n'est pas multiple de la taille du remplacement
3: In DerPart[i] <- c(DerPart_1[i], DerPart_2[i], nrow = 1, ... :
Si vous avez des propositions et des corrections, j'en serai ravi !
Merci bien !
Guillaume