Discussion :

[vince3001]

Bonjour,

J'ai créé deux tableaux de données par des procédés différents.
Ces tableaux ont en apparence les mêmes caractéristiques.
J'ai aussi créé une fonction Logit qui agit sur une ligne d'un tableau.
Je fais appel à la fonction apply pour affecter ma fonction Logit à toutes les lignes de mon tableau.

Problème : pour un tableau cela marche, mais pas pour le second.
Cela m'est incompréhensible.


Voici le code en question :

Code :

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

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
rm(list=ls())#vider la memoire
 
#Premier tableau
Y<-as.matrix(sample(1:0,248,replace=T))
N<<-NROW(Y)#variable globale
X0<-c(rep(1,N))
X1<-sample(1:10,248,replace=T)
X2<-sample(1:10,248,replace=T)
X<-cbind(X0,X1,X2)
P<<-NCOL(X)-1
 
#second tableau
Y<-infert$case
N<<-NROW(Y)#variable globale
X0<-c(rep(1,N))
X1<-infert$induced
X2<-infert$spontaneous
X<-cbind(X0,X1,X2)
P<<-NCOL(X)-1
 
#fonction Logit
Logit<-function(i,beta,X){
temp<-exp(X[i,]%*%beta)
return(temp/(1+temp))
}
 
#apply problematique a essayer avec les deux différents tableaux X
beta<-t(t(rnorm(P+1,mean=0,sd=0.7))) #vecteur colonne de taille P+1
apply(as.matrix(c(1:N)),MARGIN=1,FUN=Logit,beta=beta,X=X)

Je vous remercie infiniment de votre collaboration.

[A. D.]

Bonjour,

Dans la mesure où nous ne disposons pas de votre jeu de données appelé "infert", il ne nous est pas possible de voir où cela coince avec le second tableau.
Vous pouvez mettre un extrait du jeu de données utilisé dans votre message (si le problème apparaît sur cet extrait) ou mettre votre jeu de données en pièce jointe du message.

De plus, quel est le problème rencontré avec le second tableau?

Bonne continuation


Cordialement,

A.D.

[vince3001]

merci de cette première réponse.

En piece jointe le jeu de données (disponible sous R).
Le message d'erreur obtenu pour le second tableau est le suivant :

Code :

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

1
2
Erreur dans X[i, ] %*% beta : arguments inadéquats
Calls: apply -> FUN