[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