[R] Problems with plot function

[Uwe Ligges]

KOITA Lassana - STAC/ACE wrote:

> 
> 
> 
> Hello  all  R   users,
> My simulation function works correctly, but I have problems with plot
> function. You will find the following code using it.
> Thank you for your help
> ##################################################"
> 
> simulation <- function(k, n){
> 
> conc <- seq(0,10,by=0.5)
> #choixg <- seq(1, length(conc))
> choixg <- rep(0,length(conc))
> for (i in 1:length(conc)){
>     choixg[i] <- (k + conc[i])^2/((k+conc[i])^n + (k+1)^n)
> 
>     }
>    return(choixg)
> 
> }
> simulation(5,1)


Please read the manuals!
The objects "conc" and "choixg" ar local to your function "simulation"...

If you return

   return(list(choixg=choixg, conc=conc))

from your function, then you can plot as follows:

   simResult <- simulation(5,1)

   with(simResult,
     plot(conc, choixg, main ="fonction de choix",
       col= "blue", pch=20, xlab = " concentration",
       ylab="proba de choisir la gauche"))


Uwe Ligges


> plot(conc, choixg, main ="fonction de choix", col= "blue", pch=20,
> xlab = " concentration", ylab="proba de choisir la gauche")
> ##########################################################
> 
> Lassana KOITA
> Service Technique de l'Aviation Civile (STAC)
> Direction Générale de l'Aviation Civile (DGAC)
> Tel: 01 49 56 80 60
> Fax: 01 49 56 82 14
> http://www.stac.aviation-civile.gouv.fr
> 
> ______________________________________________
> R-help at stat.math.ethz.ch mailing list
> https://stat.ethz.ch/mailman/listinfo/r-help
> PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html