Bonjour les ami(e)s,

J'espère que vous allez bien.

J'ai un problème avec une programma et votre aide sera grandement apprécié.

1) Premier point :

J'ai deux matrices dont j'aimerais faire le produit mais dans l'une d'elle il y a un inconnu x.

Exemple :

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
 
 
# Matrice A avec un inconu x
A <- cbind(rep(1, 1), 1, 1, "x", 0.83, 0.26, 0.19, 0.03, 0.39, 1) ; A
 
 
# Matrice B
B <- c(0,39, -0.02, 0.22, 0.96, 0.02, 0.06, 0.02, 0.01, 0.32, -0.96) ; B
 
matrix(B) # transpose the vector
 
B <- as.matrix(B) ; B
 
# Produit de A*B
 
AB <- A %*% B ; AB  
 
# J'ai cette erreur : Error in A %*% B : 
# nécessite des arguments numériques/complexes matrice/vecteur
 
# ce qui est normal
 
# J'essai donc ceci mais sans rien obtenir de concret
 
AB <- function(x) A(x) %*% B ; AB
 
# Je l'ai fait manuellement à la main et je dois obtenir AB = 0.96 * x - 0.18, une fonction de x

1) Deuxième point :

Je dois ensuite faire ce calcul :

Code : Sélectionner tout - Visualiser dans une fenêtre à part
1
2
3
4
5
6
7
8
9
 
 
 
Sigma = 4
# Voici ce que j'essai de calculer au final : 
 
x <- AB * pnorm(AB/sigma) + sigma * dnorm(AB/sigma)
 
#  pnorm(.) et dnorm(.) sont respectivement les fonctions de distribution cumulative et de densité de la loi normale

Donc au final je devrais avoir une équation non linéaire de la forme x = f(x) que je vais ensuite résoudre.

Pouvez vous m'aider ?

Je ne suis pas un expert en R et je n'arrive pas a faire le produit de matrice avec inconnu a plus forte raison calculer les fonctions de distribution cumulative et de densité de la loi normale avec l'inconnu.



GRAND merci d'avance.