Skip to content

Instantly share code, notes, and snippets.

@AmritPatel

AmritPatel/code.R

Created May 20, 2014 12:47
Show Gist options
  • Star 0 You must be signed in to star a gist
  • Fork 0 You must be signed in to fork a gist
  • Save AmritPatel/90cad2e868937a744147 to your computer and use it in GitHub Desktop.
Save AmritPatel/90cad2e868937a744147 to your computer and use it in GitHub Desktop.
Download ZIP
NCS Database Plotting
Raw
code.R
library(ggplot2)
library(rCharts)
Fuel.Types <- read.csv("C:/Users/axp13/Desktop/Fuel Types.csv")
View(Fuel.Types)
Region.Specific.Uncertainties <- read.csv("C:/Users/axp13/Desktop/Region Specific Uncertainties.csv")
View(Region.Specific.Uncertainties)
merged <- merge(Fuel.Types, Region.Specific.Uncertainties, by="ID")
View(merged)
r1 <- rPlot(Total_rack_mfg_dk ~ Total_fuel_mfg_dk, data = merged, type = 'point', color = 'Manufacturer')
r1
Raw
index.html
<!doctype HTML>
<meta charset = 'utf-8'>
<html>
<head>
<link rel='stylesheet' href="http://netdna.bootstrapcdn.com/bootswatch/2.3.1/cosmo/bootstrap.min.css">
<link rel="stylesheet" href="http://netdna.bootstrapcdn.com/twitter-bootstrap/2.3.1/css/bootstrap-responsive.min.css" >
<link rel='stylesheet' href="http://getbootstrap.com/2.3.2/assets/js/google-code-prettify/prettify.css">
<link rel='stylesheet' href="http://aozora.github.io/bootplus/assets/css/docs.css">
<script src='http://ramnathv.github.io/rCharts/libraries/widgets/polycharts/js/polychart2.standalone.js' type='text/javascript'></script>
<style>
.rChart {
display: block
margin: auto auto;
width: 100%;
height: 400px;
}
/*
body {
margin-top: 60px;
}
*/
</style>
</head>
<body>
<div class='container'>
<div class='row'>
<div class='span8'>
<div class="bs-docs-example">
<div id='chart227c6c8e1432' class='rChart polycharts'>
</div>
<br/>
<pre><code class='r'>library(ggplot2)
library(rCharts)
Fuel.Types &lt;- read.csv(&quot;C:/Users/axp13/Desktop/Fuel Types.csv&quot;)
View(Fuel.Types)
Region.Specific.Uncertainties &lt;- read.csv(&quot;C:/Users/axp13/Desktop/Region Specific Uncertainties.csv&quot;)
View(Region.Specific.Uncertainties)
merged &lt;- merge(Fuel.Types, Region.Specific.Uncertainties, by=&quot;ID&quot;)
View(merged)
r1 &lt;- rPlot(Total_rack_mfg_dk ~ Total_fuel_mfg_dk, data = merged, type = 'point', color = 'Manufacturer')
r1
</code></pre>
</div>
</div>
</div>
</div>
<script type='text/javascript'>
var chartParams = {
"dom": "chart227c6c8e1432",
"width": 700,
"height": 400,
"layers": [
{
"x": "Total_fuel_mfg_dk",
"y": "Total_rack_mfg_dk",
"data": {
"ID": [ 1, 2, 3, 6, 7, 9, 10, 11, 12, 13, 14, 15, 22, 23, 24, 25, 26, 28, 29, 31, 32, 33, 34, 36, 37, 38, 39, 40, 42, 43, 44, 45 ],
"ML": [ "ML003691703", "ML003691703", "ML12188A098", "ML12188A098", "ML12188A098", "ML12188A099", "ML12188A099", "ML12188A099", "ML12188A099", "ML003691703", "ML003691703", "ML003691703", "ML112360231", "ML112360231", "ML112360231", "ML112360231", "ML112360231", "ML112360231", "ML12061A290", "ML12061A290", "ML103330094", "ML103330094", "ML103330094", "ML103330094", "ML103330094", "ML103330094", "ML103330094", "ML12132A416", "ML11291A036", "ML11291A036", "ML110560337", "ML110560337" ],
"Manufacturer": [ "WEC", "WEC", "GNF", "GNF", "GNF", "GNF", "GNF", "GNF", "GNF", "B&W", "ENC", "ENC", "WEC", "WEC", "WEC", "WEC", "ENC", "ENC", "CE", "AREVA", "GNF", "GNF", "GNF", "GNF", "AREVA", "AREVA", "AREVA", "CE", "ENC", "CE", "WEC", "WEC" ],
"Type": [ "OFA 14x14", "STANDARD 14x14", "GNF2 10x10", "9x9 (unspecified, but includes GE11)", "10x10 (unspecified)", "7x7 (unspecified)", "8x8 (unspecified)", "9x9 (unspecified, but includes GE11)", "10x10 (unspecified, includes GE14)", "14x14 (unspecified)", "SOLID 14x14", "ANNULAR 14x14", "OFA 14x14", "HIGH BURNUP-OFA 14x14", "400V+", "422V+", "STANDARD", "TOPROD", "15x15 (unspecified)", "15x15 (unspecified, Batch R+)", "GE6", "GE7", "GE11", "GNF2", "ANF 8x8", "ANF 9x9", "ATRIUM10", "16x16 (unspecified)", "14x14 (unspecified)", "14x14 (unspecified)", "STANDARD 15x15", "OFA 15x15" ],
"FreshLim": [ true, false, true, false, false, false, false, false, false, false, false, false, false, false, false, true, false, false, false, true, false, false, false, false, false, false, false, true, true, false, true, false ],
"BurnedLim": [ false, true, true, false, false, false, false, false, true, false, false, false, false, false, false, true, false, false, false, true, false, false, false, false, false, false, false, true, true, false, true, false ],
"Notes.x": [ "", "", "", "", "", "", "", "", "", "", "", "", "", "", "", "OFA is more reactive early in life, but the higher UO2 density assumed is considered to bound this reactivity difference.", "", "", "", "", "", "", "", "", "", "", "", "SFP only.", "", "In some cases with fresh or low BU, CE 14x14 is slightly more limiting. This was incorporated in the analysis via a bias.", "", "" ],
"U.235_dk": [ 0.00319, 0.00458, 0.00436, 0.00319, 0.00324, 0.00537, 0.00078, 0.00825, 0.0019, 0.00617, 0.00873, 0.00167, 0.00561, 0.00573, 0.00561, 0.00491, 0.0091, 0.00491, 0.00871, 0.00244, 0.00191, 0.00244, 0.00223, 0.00261, 0.00222, 0.00573, 0.00491, 0.0091, 0.00871, 0.00345, 0.00315, 0.00373 ],
"U.235_change_.": [ 1, 1, 1, 1.79, 1.79, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1 ],
"U_density_dk": [ 0.00221, 0.00326, 0.00394, 0.0005, 0.00041, 0.00134, 0.00256, 0.00402, 0.00317, 0.0037, 0.00192, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null ],
"U_density_change_.": [ 2.1, 2.1, 2.1, 1.36, 1.36, 2.05, 2.05, 2.05, 2.05, 2.05, 2.05, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null ],
"Pellet_diameter_dk": [ 0.00073, 0.00151, 0.00191, null, null, null, null, null, null, null, null, 0.00061, 0.00091, 0.00095, 0.00091, 0.00104, 0.00053, 0.00104, 0.00063, 0.00078, 0, 0.00078, 0.00073, 0.00068, 0, 0.00095, 0.00104, 0.00053, 0.00063, 0.00085, 0.0009, 0.00059 ],
"Pellet_diameter_change_.": [ 0.218, 0.218, 0.218, null, null, null, null, null, null, null, null, 0.137, 0.137, 0.137, 0.137, 0.137, 0.137, 0.137, 0.137, 0.137, 0.137, 0.137, 0.137, 0.137, 0.137, 0.137, 0.137, 0.137, 0.137, 0.137, 0.137, 0.137 ],
"Fuel_rod_pitch_dk": [ null, null, null, null, null, null, null, null, null, null, null, 0.01191, 0.01286, 0.01313, 0.01286, 0.013, 0.01242, 0.013, 0.0128, 0.01485, 0.01541, 0.01485, 0.01439, 0.01585, 0.01598, 0.01313, 0.013, 0.01242, 0.0128, 0.01459, 0.01324, 0.01529 ],
"Fuel_rod_pitch_change_.": [ null, null, null, null, null, null, null, null, null, null, null, -1.35, -1.35, -1.35, -1.35, -1.35, -1.35, -1.35, -1.35, -1.35, -1.35, -1.35, -1.35, -1.35, -1.35, -1.35, -1.35, -1.35, -1.35, -1.35, -1.35, -1.35 ],
"Clad_diameter_dk": [ 0.00163, 0.00319, 0.00479, 0.00152, 0.00165, null, null, null, null, null, null, 0.00184, 0.00132, 0.00221, 0.00132, 0.00139, 0.0013, 0.00139, 0.00194, 0.00166, 0.00125, 0.00166, 0.00157, 0.00151, 0.00122, 0.00221, 0.00139, 0.00129, 0.00194, 0.00168, 0.00246, 0.00078 ],
"Clad_diameter_change_.": [ -0.592, -0.592, -0.592, -13, -13, null, null, null, null, null, null, -0.711, -0.711, -0.711, -0.711, -0.711, -0.711, -0.711, -0.711, -0.711, -0.711, -0.711, -0.711, -0.711, -0.711, -0.711, -0.711, -0.711, -0.711, -0.711, -0.711, -0.711 ],
"Total_fuel_mfg_dk": [ 0.00427, 0.00664, 0.00782, 0.007, 0.00897, 0.00554, 0.00268, 0.00918, 0.0037, 0.00719, 0.00894, 0.01216, 0.01299, 0.01339, 0.01299, 0.01319, 0.01257, 0.01319, 0.01298, 0.01503, 0.01551, 0.01503, 0.01455, 0.016, 0.01607, 0.01339, 0.01319, 0.01257, 0.01298, 0.01478, 0.01355, 0.01536 ],
"Rack_cell_pitch_dk": [ null, null, 0.0139, 0.00281, 0.00284, 0.01748, 0.01181, 0.00496, 0.00023, 0.00408, 0.00497, 0.00228, 0.00608, 0.00614, 0.00608, 0.00498, 0.00757, 0.00498, 0.00753, 0.00554, 0.00565, 0.00554, 0.00452, 0.00761, 0.00739, 0.00614, 0.00498, 0.00757, 0.00753, 0.00601, 0.00457, 0.00733 ],
"Rack_cell_pitch_change_.": [ null, null, -0.853, -0.398, -0.398, -0.976, -0.976, -0.387, -0.387, -0.387, -0.388, -0.632, -0.632, -0.632, -0.632, -0.632, -0.632, -0.632, -0.632, -0.632, -0.632, -0.632, -0.632, -0.632, -0.632, -0.632, -0.632, -0.632, -0.632, -0.632, -0.632, -0.632 ],
"Rack_thickness_dk": [ 0.0019, 0.0009, 0.00265, 0.00099, 0.00066, 0.00744, 0.00197, 0.00361, 0.00345, 0.00305, 0.00353, 0.00327, 0.00533, 0.00539, 0.00533, 0.00481, 0.00526, 0.00481, 0.00542, 0.00471, 0.00463, 0.00471, 0.0046, 0.00544, 0.00502, 0.00539, 0.00481, 0.00526, 0.00542, 0.00513, 0.00434, 0.00509 ],
"Rack_thickness_change_.": [ -4.44, -9, -9, -33.3, -33.3, 50, 50, -6.7, -6.7, -6.7, -6.7, -11.1, -11.1, -11.1, -11.1, -11.1, -11.1, -11.1, -11.1, -11.1, -11.1, -11.1, -11.1, -11.1, -11.1, -11.1, -11.1, -11.1, -11.1, -11.1, -11.1, -11.1 ],
"Flux_trap_water_gap_dk": [ null, 0.00085, null, null, null, null, null, null, null, null, null, 0.00132, 0.00162, 0.00191, 0.00162, 0.00168, 0.00153, 0.00168, 0.00174, 0.00128, 0.00113, 0.00128, 0.00119, 0.00092, 0.00097, 0.00191, 0.00168, 0.00153, 0.00174, 0.0013, 0.00096, 0.00107 ],
"Flux_trap_water_gap_change_.": [ null, -35.3, null, null, null, null, null, null, null, null, null, -9.52, -9.52, -9.52, -9.52, -9.52, -9.52, -9.52, -9.52, -9.52, -9.52, -9.52, -9.52, -9.52, -9.52, -9.52, -9.52, -9.52, -9.52, -9.52, -9.52, -9.52 ],
"Poison_width_dk": [ null, null, null, 0.00017, 0.00199, 0.00337, 0.00235, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null ],
"Poison_width_change_.": [ null, null, null, -0.5, -1.53, -1.56, -1.56, null, null, null, null, null, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0 ],
"Total_rack_mfg_dk": [ 0.00259, 0.00599, 0.01415, 0.005275, 0.00353, 0.022676, 0.014302, 0.00617, 0.00347, 0.00526, 0.00618, 0.0042, 0.00822, 0.00839, 0.00822, 0.00712, 0.00934, 0.00712, 0.00944, 0.00738, 0.00739, 0.00738, 0.00656, 0.0094, 0.00899, 0.00839, 0.00712, 0.00934, 0.00944, 0.00801, 0.00638, 0.00899 ],
"IFBA_loading_dk": [ null, null, null, null, null, null, null, null, null, null, 0.00292, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null ],
"IFBA_loading_change_.": [ null, null, null, null, null, null, null, null, null, null, null, null, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0 ],
"Gd_concentration_dk": [ null, null, null, 0.00602, 0.00819, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null ],
"Gd_change_.": [ null, null, null, -7.5, -7.5, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null ],
"Boraflex_modeling_uncertainty": [ null, null, null, null, 0.0208, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null ],
"Burnup_measurement_uncertainty_dk": [ null, null, null, null, null, null, null, null, null, null, null, 0, 0, 0, 0, 0, 0, 0, 0, 0.00523, 0.00578, 0.00523, 0.00513, 0.0089, 0.00832, 0, 0, 0, 0, 0.00335, 0.00313, 0.00779 ],
"Depletion_uncertainty_dk": [ null, null, null, 0.00649, 0.00782, 0.00122, 0, 0.00942, 0, 0.00702, 0.0072, 0, 0, 0, 0, 0, 0, 0, 0, 0.00756, 0.00754, 0.00756, 0.00737, 0.0131, 0.01406, 0, 0, 0, 0, 0.00362, 0.00461, 0.00969 ],
"Fission_product_uncertainty_dk": [ null, null, null, 0.00176, 0.00204, null, null, null, null, null, null, 0, 0, 0, 0, 0, 0, 0, 0, 0.00358, 0.00321, 0.00358, 0.00244, 0.00447, 0.00494, 0, 0, 0, 0, 0.00213, 0.00189, 0.00414 ],
"Eccentric_position_dk": [ null, null, null, null, null, null, null, null, null, null, null, 0.00373, 0.00676, 0.00699, 0.00676, 0.0045, 0.00805, 0.0045, 0.01102, 0.00726, 0.00758, 0.00726, 0.00696, 0.0122, 0.00878, 0.00699, 0.00619, 0.00805, 0.01102, 0.00741, 0.00619, 0.0089 ],
"Uncertainty_in_bias_contributors": [ null, null, null, 0.00223, 0.00204, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null ],
"Uncertainty_in_mfg_tolerance_contributors": [ null, null, null, 0.00195, 0.00174, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null ],
"Notes.y": [ "", "", "", "", "", "Total rack mfg uncertainty includes several Boral related uncertainities.", "Total rack mfg uncertainty includes several Boral related uncertainities.", "", "", "", "", "Does not include fuel pellet density because design is already 1% above nominal. Fuel mfg uncertainty includes GT/IT related uncertainties.", "Does not include fuel pellet density because design is already 1% above nominal. Fuel mfg uncertainty includes GT/IT related uncertainties.", "Does not include fuel pellet density because design is already 1% above nominal. Fuel mfg uncertainty includes GT/IT related uncertainties.", "Does not include fuel pellet density because design is already 1% above nominal. Fuel mfg uncertainty includes GT/IT related uncertainties.", "Does not include fuel pellet density because design is already 1% above nominal. Fuel mfg uncertainty includes GT/IT related uncertainties.", "Does not include fuel pellet density because design is already 1% above nominal. Fuel mfg uncertainty includes GT/IT related uncertainties.", "Does not include fuel pellet density because design is already 1% above nominal. Fuel mfg uncertainty includes GT/IT related uncertainties.", "Does not include fuel pellet density because design is already 1% above nominal. Fuel mfg uncertainty includes GT/IT related uncertainties.", "Does not include fuel pellet density because design is already 1% above nominal. Fuel mfg uncertainty includes GT/IT related uncertainties.", "Does not include fuel pellet density because design is already 1% above nominal. Fuel mfg uncertainty includes GT/IT related uncertainties.", "Does not include fuel pellet density because design is already 1% above nominal. Fuel mfg uncertainty includes GT/IT related uncertainties.", "Does not include fuel pellet density because design is already 1% above nominal. Fuel mfg uncertainty includes GT/IT related uncertainties.", "Does not include fuel pellet density because design is already 1% above nominal. Fuel mfg uncertainty includes GT/IT related uncertainties.", "Does not include fuel pellet density because design is already 1% above nominal. Fuel mfg uncertainty includes GT/IT related uncertainties.", "Does not include fuel pellet density because design is already 1% above nominal. Fuel mfg uncertainty includes GT/IT related uncertainties.", "Does not include fuel pellet density because design is already 1% above nominal. Fuel mfg uncertainty includes GT/IT related uncertainties.", "Does not include fuel pellet density because design is already 1% above nominal. Fuel mfg uncertainty includes GT/IT related uncertainties.", "Does not include fuel pellet density because design is already 1% above nominal. Fuel mfg uncertainty includes GT/IT related uncertainties.", "Does not include fuel pellet density because design is already 1% above nominal. Fuel mfg uncertainty includes GT/IT related uncertainties.", "Does not include fuel pellet density because design is already 1% above nominal. Fuel mfg uncertainty includes GT/IT related uncertainties.", "Does not include fuel pellet density because design is already 1% above nominal. Fuel mfg uncertainty includes GT/IT related uncertainties." ]
},
"facet": null,
"type": "point",
"color": "Manufacturer"
}
],
"facet": [],
"guides": [],
"coord": [],
"id": "chart227c6c8e1432"
}
_.each(chartParams.layers, function(el){
el.data = polyjs.data(el.data)
})
var graph_chart227c6c8e1432 = polyjs.chart(chartParams);
</script>
</body>
<!-- Google Prettify -->
<script src="http://cdnjs.cloudflare.com/ajax/libs/prettify/188.0.0/prettify.js"></script>
<script
src='https://google-code-prettify.googlecode.com/svn-history/r232/trunk/src/lang-r.js'>
</script>
<script>
var pres = document.getElementsByTagName("pre");
for (var i=0; i < pres.length; ++i) {
pres[i].className = "prettyprint linenums";
}
prettyPrint();
</script>
</html>
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment