[R] normalised curve fitting with error bars

John Kane jrkrideau at yahoo.ca
Fri Jul 3 18:54:37 CEST 2009


Reading in your data
======================================================================
xx <- t(read.table(textConnection(x), header=FALSE, as.is=TRUE)); xx


mm <- apply(xx, 1, mean)
mhi <- apply(xx,1, max)
mlow <- apply(xx, 1, min)

plot(mm, ylim=c(min(xx)-1, max(xx)+1))
points(mhi, col="red")
points(mlow, col= "blue")

arrows(1:9, mlow, 1:9, mhi, length=0)

leg.text <- c("max", "mean", "min")
legend(7,25, leg.text, col=c("red", "black", "green"), pch= 1, cex=.75 )
==========================================================================
I think this gives you what you asked for but I don't understand why you are calling it "normalized curve fitting' since I don't see anything that you are doing to normalize anything.



--- On Fri, 7/3/09, Pooja Jain <pcxpj1 at nottingham.ac.uk> wrote:

> From: Pooja Jain <pcxpj1 at nottingham.ac.uk>
> Subject: [R] normalised curve fitting with error bars
> To: "r-help" <r-help at r-project.org>
> Received: Friday, July 3, 2009, 11:56 AM
> Dear List,
> 
> 
>  My data consist of nine columns and about 50,000 rows. It
> looks like this.
> 
> -9.0225    3.46464   
> 2.80926    -0.3847   
> 3.73735    1.1058   
> -2.98936    1.38901   
> -8.1846
> -2.4315    -5.1189   
> 1.8225    3.3798   
> 1.7874    4.693   
> -3.9286    1.4266    5.7849
> -3.4894    -4.0305   
> 3.7879    3.5195   
> 2.9186    2.8685   
> -6.126    4.978    4.9381
> 4.5282    3.62558   
> -3.0455    4.6518   
> 1.39746    0.68652   
> 3.5708    -3.6404    -4.2963
> -1.3183    0.6752   
> -4.0382    -2.5386   
> -0.6459    1.0689   
> -0.6392    -6.4141    -4.101
> -1.3735    3.1098   
> -2.8291    -5.2548   
> -3.3798    1.3959   
> -1.8605    0.1522    2.1818
> -1.0488    0.1071   
> -3.7154    -1.3748   
> -5.6218    -0.9989   
> -2.6763    -4.6548    0.5449
> 7.948    2.0673   
> 3.8729    8.0537   
> 2.79    -3.6963   
> 8.4584    -1.5122    -6.3354
> 4.5827    6.1787   
> -3.1787    1.4554   
> 5.6973    5.3386   
> 10.6077    -0.0424    3.3653
> 4.1653    8.0266   
> 1.9509    4.2077   
> 5.4182    4.1797   
> 7.9248    1.1502    0.753
> 3.6046    3.6743   
> 11.8299    9.5704   
> 10.7384    8.675   
> 4.9277    13.6898    13.3279
> 18.4431    12.3946   
> 22.1126    22.7109   
> 19.4623    17.3565   
> 15.27    17.5922    19.5873
> 
> 
> The values in each of the nine columns are almost normally
> distributed but have slightly different heights of nine
> normalised bell shaped plots.
> 
> I was trying to plot each row as a point with error bars.
> The central point will be the mean of the nine values in
> each of rows and the max and min of the row will be the
> upper and lower bound of the error bar.  I expect a
> normally distributed plot with spikes around. Can any one
> please help me plotting data in this way ?
> 
> Many thanks for any input.
> 
> -Pooja
> 
> This message has been checked for viruses but the contents
> of an attachment
> may still contain software viruses, which could damage your
> computer system:
> you are advised to perform your own checks. Email
> communications with the
> University of Nottingham may be monitored as permitted by
> UK legislation.
> 
> ______________________________________________
> R-help at r-project.org
> mailing list
> https://stat.ethz.ch/mailman/listinfo/r-help
> PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
> and provide commented, minimal, self-contained,
> reproducible code.
> 


      __________________________________________________________________
Make your browsing faster, safer, and easier with the new Internet Explorer® 8. Optimi
etexplorer/




More information about the R-help mailing list