Le traçage en échelle logarithmique en R

Je suis en train de tracer une échelle logarithmique, mais j'obtiens toujours cette erreur:

Error in plot.window(...) : invalid "log=1/h" specification.

Je ne suis pas sûr de ce que je fais de mal. Ci-dessous mon code:

   #function
function(stepsize, temp_val, counter) {
  while(counter < 0) {
    counter <- counter + stepsize
    px_norm <- dnorm(counter, mean = 0, sd = .04)
    temp_val <- temp_val + px_norm }
  temp_val <- 2*temp_val
  temp_val <- temp_val *(stepsize/2)
  print(temp_val, digits = 12)
}

#Initial step size
h <- .01 
while (h > .00001) {
  x <- calc_error(h, 0, -5) #Gives me a result around .5
  err <- x - (exp(-.02)*0.5)
  plot(1/h, err, log = "1/h")
  h <- h/10 }

Fondamentalement, dans ce court laps de fonction que j'essaie de montrer que l'étape de la taille augmente, l'erreur de la solution diminue. Cependant, je vais avoir de la difficulté à tracé de cette. Toute aide sera appréciée. Grâce

Ce n'est pas un exemple reproductible. # est utilisé pour les commentaires dans R, pas //. Qu'est-ce que calc_error? nous n'avons pas cette fonction.
Cette fonction fonctionne bien, il donne un résultat autour de .5. Je suis juste à avoir des problèmes de traçage
nous n'avons aucune foutue idée de ce que la fonction est. Nous ne pouvons pas lire dans votre esprit. sans calc_error nous ne pouvons pas générer les données pour tracer quoi que ce soit. Aussi, aucun sens de l'intrigue de chaque point à l'intérieur de la boucle.
la fonction peut bien fonctionner, mais nous n'avons pas la fonction calc_error. Vous devrez fournir le code pour cette fonction, si elle est écrite par vous ou le package dans lequel il se trouve. Sinon nous ne pouvons pas reproduire votre code.
Ok, je vais le poster

InformationsquelleAutor Josh | 2012-09-23

13

Cela devrait être un moyen facile de créer de la parcelle:
```
h <- 10^-seq(2, 4)
err <- lapply(h, function(x) calc_error(x, 0, -5) - (exp(-.02) * .5))
plot(1/h, err, log = "x")
```
- Hey merci pour le post. J'ai cherché en ligne pour la fonction lapply, mais je suis un peu confus comment il fonctionne/quand vous l'utiliser. Aussi, dans l'appel lapply(h, function(x), calc_error(x, 0, -5) - [...] pourquoi utilisez-vous des "x" et pas de "h" depuis que nous avons sont en train de changer les valeurs de x, ou est-il quelque chose comme "passage par valeur"
- La fonction est appliquée à toutes les valeurs du vecteur h. Une seule valeur est indiquée par x depuis que j'ai défini la fonction comme une fonction de x (function(x)). Fondamentalement, h est l'ensemble, mais la fonction est appliquée à chaque élément séparément.
- Aussi, lorsque j'utilise le lapply fonction, il me montre le calcul à chaque étape, est de toute façon il le cacher?
- Cela est dû à la print commande dans votre fonction. Il suffit de retirer la ligne print(temp_val, digits = 12) et vous ne verrez pas la sortie.
- Est-il une commande qui me permet de modifier une fonction? Comme éditer ou quelque chose? Le problème est que j'ai vraiment une drôle de version de R sur Ubuntu et j'utilise le terminal d'accès (en tapant R). Je ne sais pas comment je peux accéder aux fonctions que j'ai déjà créé, donc je dois en créer de nouveaux.
- Si vous n'êtes pas à l'aise avec la ligne de commande, je vous recommande de prendre un coup d'oeil à RStudio ( rstudio.org ), qui est une belle interface qui est disponible sur toutes les grandes plates-formes.
InformationsquelleAutor Sven Hohenstein
7
```
plot(1/h, err, log = "1/h")
```
journal doit être l'axe de votre parcelle de terrain que vous voulez dans l'espace du journal, mais pas les données.

c'est à dire
```
plot(1/h, err, log = "x")
```
permettra de tracer votre axe x dans l'espace du journal
- Salut, je n'arrive pas à voir un "journal" graphique, il y a juste une ligne droite sur l'axe des x, mais l'axe des x par incréments de 100 à 100000
- Ressemble plus à vous est voyant un "journal" graphique, mais que vous avez fait quelques mauvais choix pour votre y de l'échelle. Mais tout comme avec le calc_error question, je ne suis pas sûr de la façon dont vous vous attendez de nous aider sans vraiment avoir de vos données en face de nous.
InformationsquelleAutor AGS

Je suppose que vous voulez quelque chose comme ceci:

i <- 1
h <- 0.01
err <- vector(0,mode="numeric")

while (h[i] > .00001) {
  x <- rnorm(1,mean=0.5,sd=0.05) #use calc_error instead
  err <- c(err,x - (exp(-.02)*0.5))
  i <- i+1
  h <- c(h,h[i-1]/10) }

plot(1/h[-length(h)],err,log="x")

Sven code est mieux, bien sûr (ou au moins dans l'esprit de R), mais vous êtes probablement plus familier avec une boucle while.

InformationsquelleAutor Roland

Vous devez vous connecter pour publier un commentaire.