Utilisation de la fonction SI dans la R

Je suis en cours d'exécution d'un court if fonction dans R, mais je reçois le message d'avertissement suivant:

In if ((runif(50, 0, 1) < 0.69)) { :  the condition has length > 1 and only the first element will be used`.

Mon rudimentaire saisir de R m'amène à croire que runif génère un vecteur mais if donne une valeur unique, donc je pense que c'est le problème. Est-il une simple substitution de si ici?

Aussi, je veux que le produit d'une matrice de combinaison des deux arguments mais je n'étais pas sûr si c'était correct pour mettre 50 dans le rnorm fonction pour les deux scénarios.

Test <-
 if((runif(50, 0, 1)<0.69)) {
 rnorm(50, 25, 4)
 } else {
 rnorm(50, 28, 4.3)
}

Que voulez-vous obtenir à la fin ? Un vecteur de 50 articles ?

OriginalL'auteur YesSure | 2012-08-01

if-statement r

8

if() s'attend à ce que est dans les parenthèses pour s'évaluer à une logique de vecteur de longueur 1 (TRUE ou FALSE). Si le vecteur est supérieure à 1, alors que l'avertissement dit que seul le premier élément du vecteur est utilisé.

Une alternative si ifelse() où il est utilisé comme :
```
ifelse(test, TRUE_CASE, FALSE_CASE)
```
Pour votre exemple, nous avons:
```
set.seed(1)
ifelse(runif(50, 0, 1) < 0.69, rnorm(50, 25, 4), rnorm(50, 28, 4.3))
```
ou
```
set.seed(1)
test <- runif(50, 0, 1)
Tcase <- rnorm(50, 25, 4)
Fcase <- rnorm(50, 28, 4.3)
ifelse(test, Tcase, Fcase)
```
Est OK comme ifelse() génère les deux vecteurs de TRUE_CASE et FALSE_CASE avant de faire les comparaisons.

Exemple de sortie est:
```
R> ifelse(test < 0.69, Tcase, Fcase)
 [1] 24.77549 24.37682 19.11699 28.31967 26.67177 25.55472 27.41873 26.55069
 [9] 24.78478 19.49176 23.34002 23.42284 24.76275 29.40010 29.14852 24.34191
[17] 33.19383 32.98973 27.22665 34.82338 30.40149 26.45833 28.07413 24.55062
[25] 28.52443 26.59242 22.55189 26.36448 28.67952 30.73209 32.92160 23.53111
[33] 20.82346 27.27888 35.23336 34.60647 26.01493 27.75896 25.20201 22.02691
[41] 26.31093 17.78017 26.79981 25.61301 33.69045 25.82438 22.16021 27.44291
[49] 27.22791 27.56918
```
Une autre façon est de le faire est de générer seulement les valeurs requises
```
want <- test < 0.69
res <- test ## copy vector of correct length
res[want] <- rnorm(S <- sum(want), 25, 4)
res[!want] <- rnorm(length(res) - S, 28, 4.3)

R> res
 [1] 27.85067 24.70574 24.84946 22.04696 22.27336 36.03795 29.82793 23.70292
 [9] 25.24064 22.64442 27.12598 18.92642 26.22623 18.85420 26.97382 23.79610
[17] 32.55148 31.81162 22.88688 25.33725 37.48624 22.39162 24.77241 17.34256
[25] 29.70633 18.34011 23.14588 20.53632 26.90338 21.99672 33.34867 25.06958
[33] 19.85480 18.43758 21.87467 26.80075 27.37908 24.92576 28.89241 23.72773
[41] 37.92431 21.28255 28.45495 19.05016 20.69923 29.96509 29.00012 22.51493
[49] 27.66824 26.56380
```
Qui donne évidemment des numéros différents pour les précédents, mais seulement parce que nous avons généré moins de nombres aléatoires. Je doute que cela sera plus efficace que ifelse() pour la plupart des problèmes.

Excellent, Merci. J'avais juste trébuché à travers ifelse
Gavin Juste pour la clarification complète de ce ifelse est en train de faire: prendre la première des 50 numéros générés de manière aléatoire entre 0 et 1 à l'aide de runif, si cette valeur est inférieure à 0,69$, alors la valeur est générée à partir de la distribution normale (50, 25, 4) et si elle est plus grande, alors la valeur est générée à partir de (50, 28, 4.3). Ensuite, le résultat final est tout simplement un ammalgamation des valeurs dans une liste ou une matrice de 50 valeurs. Mon résumé correct pour l'intervention?
Oui, sauf pour le dernier bit; le résultat est un nombre vecteur de 50 valeurs. J'ai ajouté un exemple de la sortie.
Il est plus efficace de générer des 50 valeurs en une seule fois, pour chaque scénario, que d'avoir de nombreux appels à rnorm() déterminé sur la base de la clause. Il y a d'autres façons de le faire.
Veuillez ne pas faire des choses comme ça. DONC, ce n'est pas seulement un moyen pour vous d'obtenir des Réponses, mais est également une ressource pour les autres personnes qui pourraient rencontrer des problèmes similaires. Il n'y a rien de mal avec le démarrage d'une nouvelle Question, un lien vers celui-ci pour des détails pertinents. Je vais jeter un oeil à le Modifier par la suite.

OriginalL'auteur Gavin Simpson
1

Vous avez besoin any ou all, qui convertir un vecteur de partitions logiques contenues dans une seule logique (et, en fait un vecteur de longueur 1).

OriginalL'auteur Ari B. Friedman

Vous devez vous connecter pour publier un commentaire.