[R] plotting gaussian data
Jacques VESLOT
jacques.veslot at good.ibl.fr
Tue Jun 13 10:41:44 CEST 2006
curve(coef$A * dnorm(x, coef$mu, coef$sig), -4, 4)
-------------------------------------------------------------------
Jacques VESLOT
CNRS UMR 8090
I.B.L (2ème étage)
1 rue du Professeur Calmette
B.P. 245
59019 Lille Cedex
Tel : 33 (0)3.20.87.10.44
Fax : 33 (0)3.20.87.10.31
http://www-good.ibl.fr
-------------------------------------------------------------------
H. Paul Benton a écrit :
> Ok I guess it's time to ask.
> So I want to plot my data. It's my data from a frequency table, "temp". My
> formula is just a Gaussian eq. I have done the nls function to get my
> parameters and now I want to do the whole plot (...) and then lines(..)
> This is what I have done.
>
>
>>temp
>
> bin x
> 1 -4.0 0
> 2 -3.9 0
> 3 -3.8 0
> 4 -3.7 0
> 5 -3.6 0
> 6 -3.5 0 .... and so on
>
>>fo
>
> x ~ (A/(sig * sqrt(2 * pi))) * exp(-1 * ((bin - mu)^2/(2 * sig^2)))
>
>>fo.v
>
> x ~ (335.48/(0.174 * sqrt(2 * pi))) * exp(-1 * ((bin - (-0.0786))^2/(2 *
> 0.174^2)))
>
>>coef
>
> $A
> [1] 335.4812
>
> $mu
> [1] -0.07863746
>
> $sig
> [1] 0.1746473
>
> plot(fo, temp, coef)
> Error in eval(expr, envir, enclos) : object "sig" not found
>
>>plot(fo, coef, temp)
>
> Error in eval(expr, envir, enclos) : object "x" not found
>
>>plot(fo, temp, list=coef)
>>
>
>
> If someone could point me in the right direction I would be very grateful.
>
> Cheers,
>
> ______________________________________________
> 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
>
More information about the R-help
mailing list