R help - Oct 2005 - problems with loop and plot function

Hi all R users,

I have problems with my second loop for drawing the three curves in the
same graphic. I need help please

Thank you in advance

#########################################################################

simulation <- function(k, n){

conc <- seq(0,100,by=0.5)
#choixg <- seq(1, length(conc))
choixg <- rep(0,length(conc))
for (i in 1:length(conc)){
    choixg[i] <- (k + conc[i])^n/((k+conc[i])^n + (k+1)^n)

    }
   #return(data.frame(choixg, conc))
   return(list(choixg=choixg, conc = conc))
}
#cbind(simResult$conc, simResult$choixg,
format(Sys.time(),"%H:%M:%S"))

for( n in c(1,3,5)) {
 x <- NULL; y <- NULL
 simResult <- simulation (5,n)
 conc <- simResult$conc
 choixg <- simResult$choixg
 prin("n");  print(c(k=5, n))
 points(conc, log10(1-choixg), main ="fonction de choix", col= n,
pch=20,
lwd = 3,
xlab = " concentration", ylab="proba de choisir la gauche",
type="l")
}


#########################################################################


Lassana KOITA
Etude S??curit?? et Exploitation a??roportuaires / Aerodrome Safety &
Statistical analysis
Service Technique de l'Aviation Civile (STAC) / Civil Aviation Technical
Department
Direction G??n??rale de l'Aviation Civile (DGAC) / French Civil Aviation
Tel: 01 49 56 80 60
Fax: 01 49 56 82 14
E-mail: Lassana.Koita at aviation-civile.gouv.fr

Hi

you forgot to plot something before adding points (and actually 
you wanted lines).

plot(conc, log10(1-choixg), main ="fonction de choix", col= n,
pch=20,
lwd = 3,
xlab = " concentration", ylab="proba de choisir la gauche",
type="n")

for( n in c(1,3,5)) {
 x <- NULL; y <- NULL
 simResult <- simulation (5,n)
 conc <- simResult$conc
 choixg <- simResult$choixg
 print("n");  print(c(k=5, n)) ### what is this???
 lines(conc, log10(1-choixg))
}

Isn't it time to read some introductory texts?
:-)


HTH
Petr

On 13 Oct 2005 at 12:47, KOITA Lassana - STAC/ACE wrote:

To:             	R-help at stat.math.ethz.ch
From:           	"KOITA Lassana - STAC/ACE" <lassana.koita at
aviation-civile.gouv.fr>
Date sent:      	Thu, 13 Oct 2005 12:47:27 +0200
Subject:        	[R] problems with loop and plot function
> 
> 
> 
> 
> Hi all R users,
> 
> I have problems with my second loop for drawing the three curves in
> the same graphic. I need help please
> 
> Thank you in advance
> 
> ######################################################################
> ###
> 
> simulation <- function(k, n){
> 
> conc <- seq(0,100,by=0.5)
> #choixg <- seq(1, length(conc))
> choixg <- rep(0,length(conc))
> for (i in 1:length(conc)){
>     choixg[i] <- (k + conc[i])^n/((k+conc[i])^n + (k+1)^n)
> 
>     }
>    #return(data.frame(choixg, conc))
>    return(list(choixg=choixg, conc = conc))
> }
> #cbind(simResult$conc, simResult$choixg,
> #format(Sys.time(),"%H:%M:%S"))
> 
> for( n in c(1,3,5)) {
>  x <- NULL; y <- NULL
>  simResult <- simulation (5,n)
>  conc <- simResult$conc
>  choixg <- simResult$choixg
>  prin("n");  print(c(k=5, n))
>  points(conc, log10(1-choixg), main ="fonction de choix", col= n,
>  pch=20,
> lwd = 3,
> xlab = " concentration", ylab="proba de choisir la
gauche", type="l")
> }
> 
> 
> ######################################################################
> ###
> 
> 
> Lassana KOITA
> Etude S??curit?? et Exploitation a??roportuaires / Aerodrome Safety &
> Statistical analysis Service Technique de l'Aviation Civile (STAC) /
> Civil Aviation Technical Department Direction G??n??rale de l'Aviation
> Civile (DGAC) / French Civil Aviation Tel: 01 49 56 80 60 Fax: 01 49
> 56 82 14 E-mail: Lassana.Koita at 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
Petr Pikal
petr.pikal at precheza.cz