Re: [R] Problems with plot function

From: Petr Pikal <petr.pikal_at_precheza.cz>
Date: Wed 12 Oct 2005 - 16:48:09 EST


Hi

Just a small modification. You has to return also conc.

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(data.frame(choixg, conc))

}
mydf<-simulation(5,1)

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

HTH
Petr

On 11 Oct 2005 at 18:54, KOITA Lassana - STAC/ACE wrote:

To:             	r-help@stat.math.ethz.ch
From:           	"KOITA Lassana - STAC/ACE" <lassana.koita@aviation-civile.gouv.fr>
Date sent:      	Tue, 11 Oct 2005 18:54:18 +0200
Subject:        	[R] Problems with plot function

>
>
>
>
> 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)
>
> 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@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@precheza.cz



R-help@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 Received on Wed Oct 12 17:00:04 2005

This archive was generated by hypermail 2.1.8 : Sun 23 Oct 2005 - 18:47:05 EST