Differences

This shows you the differences between two versions of the page.

--- kvalobs:kvoss:system:qc2:8hdk376snf09zj37dk82s92:ncdf [2009-06-20 23:36:06]
paule
+++ kvalobs:kvoss:system:qc2:8hdk376snf09zj37dk82s92:ncdf [2022-05-31 09:29:32] (current)
@@ Line 15: / Line 15: @@
   $ R CMD INSTALL ncdf
   $ history
+  ...
+  ...
+  # there are perhaps even easier ways to do this:
+  $ install.packages("gplots", dependencies = TRUE)
+  # no source was linked ... but the mirror choose automatically opens
+  # and one just selects a site and the install happens!
+  ...
+  $
 </code>
@@ Line 22: / Line 32: @@
   $ R
    Method='unspecified'
+   Shortname='outfile'
    library(ncdf)
    nd = open.ncdf("interpolations.nc")
@@ Line 33: / Line 44: @@
    yyy=yy[!is.na(yy) & yy>=0 & yy < 100]
    xxx=xx[!is.na(yy) & yy>=0 & yy < 100]
-   plot(xxx,yyy,pch=46,xlabel='Original Value',ylabel='Interpoltion',main=Method)
+   jpeg(filename=paste(Shortname,'jpg',sep='.'))
-   cor(xxx,yyy)
+   plot(xxx,yyy,pch=46,xlab='Original Value',ylab='Interpoltion',main=Method)
+   XC=cor(xxx,yyy)
    ff <- lm(xxx ~ yyy)
    ff$coefficients[1]
@@ Line 40: / Line 52: @@
    z=ff$coefficients[1]*xxx + ff$coefficients[2]
    #lines(xxx,z)
-   text(10,95,"Linear Regression",pos=2)
+   text(0,95,"Linear Regression",pos=4)
-   text(10,90,"Gradient",pos=2)
+   text(0,90,paste("Gradient =",format(ff$coefficients[1],digits=5)),pos=4)
-   text(10,85,"Intercept",pos=2)
+   text(0,85,paste("Intercept =",format(ff$coefficients[2],digits=5)),pos=4)
-   text(10,75,"Correlation coefficient",Pos=4)
+   text(0,75,paste("Correlation coefficient =",format(XC,digits=5)),pos=4)
+   dev.off()
 </code>
+Alternatively, place the code in a file and run from the command line:
+<code>
+ $ Rscript plot_netcdf.R "Test Plot" OUTPUT
+</code>
+that will produce the file OUTPUT.jpg with the plot title "Test Plot".
+  * {{kvalobs:kvoss:system:qc2:8hdk376snf09zj37dk82s92:plot_netcdf.pdf|plot_netcdf.R}}