Discussion :

[lesuricate]

Bonjour à toutes et à tous,

Je suis bloqué sur une ligne de code en R.
J'aimerais pouvoir soustraire les valeurs contenues dans les vecteurs X et A qui sont tous les 2 issus de la fonction sapply. Ce sont donc 2 vecteurs composés de vecteurs de tailles variables.
Le vecteur X contient, dans chaque cellule, N variables suivant une loi gamma. N allant généralement de 0 jusqu'à 4.
Le vecteur A contient, dans chaque cellule, N constantes Alpha = 100.
Pourriez-vous m'aider à soustraire, dans chaque cellule de X, pour chaque variable de X, la valeur de la variable de X - 100.

Je vous mets un petit bout du code que j'ai fait pour vous aider à mieux comprendre ce qui pourrait me bloquer

Code :

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

1
2
3
4
5
6
7
8
9
n <-100
m <-10000
Lambda <- 1000
Alpha <- 100
 
N <- rpois(m,0.01*n)
X <- sapply(N , function(x)  rgamma(x,1000000/Lambda,scale = Lambda/1000))
A <- sapply(N , function(x)  rep(Alpha,x))
Y <- sapply(N , function(x) X[[x]] - A[[x]])

Je vous remercie d'avance pour votre aide

[lesuricate]

J'ai réussi à résoudre le problème,
il suffisait de poser :

Code :

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

X <- sapply(N , function(x)  if(x > 0) {rgamma(x,1000000/Lambda,scale = Lambda/1000)-Alpha} else {0})

J'ai néanmoins un nouveau problème :
j'aimerais calculer dans chaque cellule de X, pour chaque variable de la cellule X,
le max(X,Beta)où Beta <- 3000Pour se faire je simule W de cette manière

Code :

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

W <- sapply((1:m) , function(x) {for (i in 1:length(X[[x]])) {min(X[[x]][i],Beta)}})

mais cette commande ne marche pas :/
pourriez-vous m'aider svp
merci d'avance.