<p><div><div><div><div data-show-faces="false" data-height="21" data-width="100" data-layout="button_count" data-send="true" data-href="https://www.r-bloggers.com/how-happy-is-your-country%e2%80%8a-%e2%80%8ahappy-planet-index-visualized/"></div></div><div></div></div><div></div></div>
<div>
(This article was first published on <strong><a href="https://datascienceplus.com/how-happy-is-your-country-visualized/"> R Programming – DataScience+</a></strong>, and kindly contributed to <a rel="nofollow" href="https://www.r-bloggers.com/">R-bloggers)</a> </div>
<p><a rel="nofollow" href="http://happyplanetindex.org/" target="_blank">The Happy Planet Index (HPI)</a> is an index of human well-being and environmental impact that was introduced by <a rel="nofollow" href="http://neweconomics.org/" target="_blank">NEF</a>, a UK-based economic think tank promoting social, economic and environmental justice. It ranks 140 countries according to “what matters most — sustainable wellbeing for all”.</p>
<p>This is how HPI is calculated:<br>
<a rel="nofollow" href="https://i0.wp.com/datascienceplus.com/wp-content/uploads/2017/11/1_qiUXIMC7p_iuH5_YRq9W9Q.png?ssl=1" target="_blank"><img sizes="(max-width: 490px) 100vw, 490px" src="https://i1.wp.com/datascienceplus.com/wp-content/uploads/2017/11/1_qiUXIMC7p_iuH5_YRq9W9Q-490x60.png?w=450&ssl=1" alt="" class="wp-image-9397" data-recalc-dims="1"></a><br>
It’s tells us “how well nations are doing at achieving long, happy, sustainable lives”. The index is weighted to give progressively higher scores to nations with lower ecological footprints.</p>
<p>I downloaded the 2016 dataset from the <a rel="nofollow" href="http://happyplanetindex.org/countries" target="_blank">HPI website</a>. Inspired by “<a rel="nofollow" href="https://datascienceplus.com/web-scraping-and-applied-clustering-global-happiness-and-social-progress-index/" target="_blank">Web Scraping and Applied Clustering Global Happiness and Social Progress Index</a>” written by Dr. Mesfin Gebeyaw, I am interested to find correlations among happiness, wealth, life expectancy, footprint and so on, and then put these 140 countries into different clusters, according to the above measures. I wonder whether the findings will surprise me.</p>
<p>Note: for those who want to see the results right now, I have created a Tableau story, that can be accessed from <a rel="nofollow" href="https://public.tableau.com/profile/susan.li#!/vizhome/HappyPlanetIndex_2/2016WorldHappinessStory" target="_blank">here</a>.</p>
<p><b>Load the packages</b></p>
<pre>
library(dplyr)
library(plotly)
library(stringr)
library(cluster)
library(FactoMineR)
library(factoextra)
library(ggplot2)
library(reshape2)
library(ggthemes)
library(NbClust)
</pre>
<p><b>Data Preprocessing</b></p>
<pre>
library(xlsx)
hpi &lt;- read.xlsx('hpi-data-2016.xlsx',sheetIndex = 5, header = TRUE)

# Remove the unnecessary columns
hpi &lt;- hpi[c(3:14)]
# remove footer
hpi &lt;- hpi[-c(141:158), ]
# rename columns
hpi &lt;- hpi[,c(grep('Country', colnames(hpi)), grep('Region', colnames(hpi)), grep('Happy.Planet.Index', colnames(hpi)), grep('Average.Life..Expectancy', colnames(hpi)), grep('Happy.Life.Years', colnames(hpi)), grep('Footprint..gha.capita.', colnames(hpi)), grep('GDP.capita...PPP.', colnames(hpi)), grep('Inequality.of.Outcomes', colnames(hpi)), grep('Average.Wellbeing..0.10.', colnames(hpi)), grep('Inequality.adjusted.Life.Expectancy', colnames(hpi)), grep('Inequality.adjusted.Wellbeing', colnames(hpi)), grep('Population', colnames(hpi)))]
names(hpi) &lt;- c('country', 'region','hpi_index', 'life_expectancy', 'happy_years', 'footprint', 'gdp', 'inequality_outcomes', 'wellbeing', 'adj_life_expectancy', 'adj_wellbeing', 'population')
# change data type
hpi$country &lt;- as.character(hpi$country)
hpi$region &lt;- as.character(hpi$region)
</pre>
<p><b>The structure of the data</b></p>
<pre>
str(hpi)
<em>'data.frame':	140 obs. of 12 variables:
$ country : chr &quot;Afghanistan&quot; &quot;Albania&quot; &quot;Algeria&quot; &quot;Argentina&quot; ##...
$ region : chr &quot;Middle East and North Africa&quot; &quot;Post-communist&quot; ##&quot;Middle East and North Africa&quot; &quot;Americas&quot; ...
$ hpi_index : num 20.2 36.8 33.3 35.2 25.7 ...
$ life_expectancy : num 59.7 77.3 74.3 75.9 74.4 ...
$ happy_years : num 12.4 34.4 30.5 40.2 24 ...
$ footprint : num 0.79 2.21 2.12 3.14 2.23 9.31 6.06 0.72 5.09 7.44 ...
$ gdp : num 691 4247 5584 14357 3566 ...
$ inequality_outcomes: num 0.427 0.165 0.245 0.164 0.217 ...
$ wellbeing : num 3.8 5.5 5.6 6.5 4.3 7.2 7.4 4.7 5.7 6.9 ...
$ adj_life_expectancy: num 38.3 69.7 60.5 68.3 66.9 ...
$ adj_wellbeing : num 3.39 5.1 5.2 6.03 3.75 ...
$ population : num 29726803 2900489 37439427 42095224 2978339 ...
</em>
</pre>
<p><b>The summary</b></p>
<pre>
summary(hpi[, 3:12])
<em> hpi_index life_expectancy happy_years footprint Min. :12.78 Min. :48.91 Min. : 8.97 Min. : 0.610 1st Qu.:21.21 1st Qu.:65.04 1st Qu.:18.69 1st Qu.: 1.425 Median :26.29 Median :73.50 Median :29.40 Median : 2.680 Mean :26.41 Mean :70.93 Mean :30.25 Mean : 3.258 3rd Qu.:31.54 3rd Qu.:77.02 3rd Qu.:39.71 3rd Qu.: 4.482 Max. :44.71 Max. :83.57 Max. :59.32 Max. :15.820 gdp inequality_outcomes wellbeing adj_life_expectancy
Min. : 244.2 Min. :0.04322 Min. :2.867 Min. :27.32 1st Qu.: 1628.1 1st Qu.:0.13353 1st Qu.:4.575 1st Qu.:48.21 Median : 5691.1 Median :0.21174 Median :5.250 Median :63.41 Mean : 13911.1 Mean :0.23291 Mean :5.408 Mean :60.34 3rd Qu.: 15159.1 3rd Qu.:0.32932 3rd Qu.:6.225 3rd Qu.:72.57 Max. :105447.1 Max. :0.50734 Max. :7.800 Max. :81.26 adj_wellbeing population Min. :2.421 Min. :2.475e+05 1st Qu.:4.047 1st Qu.:4.248e+06 Median :4.816 Median :1.065e+07 Mean :4.973 Mean :4.801e+07 3rd Qu.:5.704 3rd Qu.:3.343e+07 Max. :7.625 Max. :1.351e+09
</em>
</pre>
<pre>
ggplot(hpi, aes(x=gdp, y=life_expectancy)) + 
  geom_point(aes(size=population, color=region)) + coord_trans(x = 'log10') +
  geom_smooth(method = 'loess') + ggtitle('Life Expectancy and GDP per Capita in USD log10') + theme_classic()
</pre>
<p>Gives this plot:<br>
<a rel="nofollow" href="https://i1.wp.com/datascienceplus.com/wp-content/uploads/2017/11/life-1.png?ssl=1" target="_blank"><img sizes="(max-width: 490px) 100vw, 490px" src="https://i1.wp.com/datascienceplus.com/wp-content/uploads/2017/11/life-1-490x228.png?w=450&ssl=1" alt="" class="wp-image-9400" data-recalc-dims="1"></a><br>
After log transformation, the relationship between GDP per capita and life expectancy is more clear and looks relatively strong. These two variables are correlated. The Pearson correlation between this two variable is reasonably high, at approximate 0.62.</p>
<pre>
cor.test(hpi$gdp, hpi$life_expectancy)
<em>Pearson's product-moment correlation
data: hpi$gdp and hpi$life_expectancy
t = 9.3042, df = 138, p-value = 2.766e-16
alternative hypothesis: true correlation is not equal to 0
95 percent confidence interval:
0.5072215 0.7133067
sample estimates:
cor 0.6208781
</em>
</pre>
<pre>
ggplot(hpi, aes(x=life_expectancy, y=hpi_index)) + 
  geom_point(aes(size=population, color=region)) + geom_smooth(method = 'loess') + ggtitle('Life Expectancy and Happy Planet Index Score') + theme_classic()
</pre>
<p>Gives this plot:<br>
<a rel="nofollow" href="https://i0.wp.com/datascienceplus.com/wp-content/uploads/2017/11/life-2.png?ssl=1" target="_blank"><img sizes="(max-width: 490px) 100vw, 490px" src="https://i0.wp.com/datascienceplus.com/wp-content/uploads/2017/11/life-2-490x228.png?w=450&ssl=1" alt="" class="wp-image-9401" data-recalc-dims="1"></a><br>
Many countries in Europe and Americas end up with middle-to-low HPI index probably because of their big carbon footprints, despite the long life expectancy.</p>
<pre>
ggplot(hpi, aes(x=gdp, y=hpi_index)) + geom_point(aes(size=population, color=region)) + geom_smooth(method = 'loess') + ggtitle('GDP per Capita(log10) and Happy Planet Index Score') + coord_trans(x = 'log10')
</pre>
<p>Gives this plot:<br>
<a rel="nofollow" href="https://i0.wp.com/datascienceplus.com/wp-content/uploads/2017/11/life-3.png?ssl=1" target="_blank"><img sizes="(max-width: 490px) 100vw, 490px" src="https://i1.wp.com/datascienceplus.com/wp-content/uploads/2017/11/life-3-490x228.png?w=450&ssl=1" alt="" class="wp-image-9402" data-recalc-dims="1"></a></p>
<p>Money can’t buy happiness. The correlation between GDP and Happy Planet Index score is indeed very low, at about 0.11.</p>
<pre>
cor.test(hpi$gdp, hpi$hpi_index)
<em>Pearson's product-moment correlation
data: hpi$gdp and hpi$hpi_index
t = 1.3507, df = 138, p-value = 0.179
alternative hypothesis: true correlation is not equal to 0
95 percent confidence interval:
-0.05267424 0.27492060
sample estimates:
cor 0.1142272 </em>
</pre>
<h3>Scale the data</h3>
<p>An important step of meaningful clustering consists of transforming the variables such that they have mean zero and standard deviation one.</p>
<pre>
hpi[, 3:12] &lt;- scale(hpi[, 3:12])
summary(hpi[, 3:12])
<em>hpi_index life_expectancy happy_years footprint Min. :-1.86308 Min. :-2.5153 Min. :-1.60493 Min. :-1.1493 1st Qu.:-0.71120 1st Qu.:-0.6729 1st Qu.:-0.87191 1st Qu.:-0.7955 Median :-0.01653 Median : 0.2939 Median :-0.06378 Median :-0.2507 Mean : 0.00000 Mean : 0.0000 Mean : 0.00000 Mean : 0.0000 3rd Qu.: 0.70106 3rd Qu.: 0.6968 3rd Qu.: 0.71388 3rd Qu.: 0.5317 Max. : 2.50110 Max. : 1.4449 Max. : 2.19247 Max. : 5.4532 gdp inequality_outcomes wellbeing adj_life_expectancy
Min. :-0.6921 Min. :-1.5692 Min. :-2.2128 Min. :-2.2192 1st Qu.:-0.6220 1st Qu.:-0.8222 1st Qu.:-0.7252 1st Qu.:-0.8152 Median :-0.4163 Median :-0.1751 Median :-0.1374 Median : 0.2060 Mean : 0.0000 Mean : 0.0000 Mean : 0.0000 Mean : 0.0000 3rd Qu.: 0.0632 3rd Qu.: 0.7976 3rd Qu.: 0.7116 3rd Qu.: 0.8221 Max. : 4.6356 Max. : 2.2702 Max. : 2.0831 Max. : 1.4059 adj_wellbeing population Min. :-2.1491 Min. :-0.2990 1st Qu.:-0.7795 1st Qu.:-0.2740 Median :-0.1317 Median :-0.2339 Mean : 0.0000 Mean : 0.0000 3rd Qu.: 0.6162 3rd Qu.:-0.0913 Max. : 2.2339 Max. : 8.1562
</em>
</pre>
<p>A simple correlation heatmap</p>
<pre>
qplot(x=Var1, y=Var2, data=melt(cor(hpi[, 3:12], use=&quot;p&quot;)), fill=value, geom=&quot;tile&quot;) +
  scale_fill_gradient2(limits=c(-1, 1)) + 
  theme(axis.text.x = element_text(angle = 90, hjust = 1)) + 
             labs(title=&quot;Heatmap of Correlation Matrix&quot;, 
                  x=NULL, y=NULL)
</pre>
<p>Gives this plot:<br>
<a rel="nofollow" href="https://i0.wp.com/datascienceplus.com/wp-content/uploads/2017/11/life-4.png?ssl=1" target="_blank"><img sizes="(max-width: 490px) 100vw, 490px" src="https://i1.wp.com/datascienceplus.com/wp-content/uploads/2017/11/life-4-490x228.png?w=450&ssl=1" alt="" class="wp-image-9403" data-recalc-dims="1"></a></p>
<h2>Principal Component Analysis (PCA)</h2>
<p><a rel="nofollow" href="https://datascienceplus.com/principal-component-analysis-unsupervised-learning/" target="_blank">PCA</a> is a procedure for identifying a smaller number of uncorrelated variables, called “principal components”, from a large set of data. The goal of principal components analysis is to explain the maximum amount of variance with the minimum number of principal components.</p>
<pre>
hpi.pca &lt;- PCA(hpi[, 3:12], graph=FALSE)
print(hpi.pca)
<em>**Results for the Principal Component Analysis (PCA)**
The analysis was performed on 140 individuals, described by 10 variables
*The results are available in the following objects:
name description 1 &quot;$eig&quot; &quot;eigenvalues&quot; 2 &quot;$var&quot; &quot;results for the variables&quot; 3 &quot;$var$coord&quot; &quot;coord. for the variables&quot; 4 &quot;$var$cor&quot; &quot;correlations variables - dimensions&quot;
5 &quot;$var$cos2&quot; &quot;cos2 for the variables&quot; 6 &quot;$var$contrib&quot; &quot;contributions of the variables&quot; 7 &quot;$ind&quot; &quot;results for the individuals&quot; 8 &quot;$ind$coord&quot; &quot;coord. for the individuals&quot; 9 &quot;$ind$cos2&quot; &quot;cos2 for the individuals&quot; 10 &quot;$ind$contrib&quot; &quot;contributions of the individuals&quot; 11 &quot;$call&quot; &quot;summary statistics&quot; 12 &quot;$call$centre&quot; &quot;mean of the variables&quot; 13 &quot;$call$ecart.type&quot; &quot;standard error of the variables&quot; 14 &quot;$call$row.w&quot; &quot;weights for the individuals&quot; 15 &quot;$call$col.w&quot; &quot;weights for the variables&quot; </em>
</pre>
<pre>
eigenvalues &lt;- hpi.pca$eig
head(eigenvalues)
<em>eigenvalue percentage of variance cumulative percentage of variance
comp 1 6.66741533 66.6741533 66.67415
comp 2 1.31161290 13.1161290 79.79028
comp 3 0.97036077 9.7036077 89.49389
comp 4 0.70128270 7.0128270 96.50672
comp 5 0.24150648 2.4150648 98.92178
comp 6 0.05229306 0.5229306 99.44471
</em>
</pre>
<p><b>Interpretation:</b><br>
* The proportion of variation retained by the principal components was extracted above.<br>
* Eigenvalues is the amount of variation retained by each PC. The first PC corresponds to the maximum amount of variation in the data set. In this case, the first two principal components are worthy of consideration because <a rel="nofollow" href="http://www.rc.usf.edu/~jdorio/FA/Cliff%20(1988)%20The%20Eigenvalues-Greater-Than-One%20Rule%20and%20the%20Reliability%20of%20Components.pdf" target="_blank">A commonly used criterion for the number of factors to rotate is the eigenvalues-greater-than-one rule proposed by Kaiser (1960)</a>.</p>
<pre>
fviz_screeplot(hpi.pca, addlabels = TRUE, ylim = c(0, 65))
</pre>
<p>Gives this plot:<br>
<a rel="nofollow" href="https://i1.wp.com/datascienceplus.com/wp-content/uploads/2017/11/life-5.png?ssl=1" target="_blank"><img sizes="(max-width: 490px) 100vw, 490px" src="https://i2.wp.com/datascienceplus.com/wp-content/uploads/2017/11/life-5-490x228.png?w=450&ssl=1" alt="" class="wp-image-9404" data-recalc-dims="1"></a><br>
The scree plot shows us which components explain most of the variability in the data. In this case, almost 80% of the variances contained in the data are retained by the first two principal components.</p>
<pre>
head(hpi.pca$var$contrib)
<em>Dim.1 Dim.2 Dim.3 Dim.4 Dim.5
hpi_index 3.571216 50.96354921 5.368971166 2.1864830 5.28431372
life_expectancy 12.275001 2.29815687 0.002516184 18.4965447 0.31797242
happy_years 14.793710 0.01288175 0.027105103 0.7180341 0.03254368
footprint 9.021277 24.71161977 2.982449522 0.4891428 7.62967135
gdp 9.688265 11.57381062 1.003632002 2.3980025 72.49799232
inequality_outcomes 13.363651 0.30494623 0.010038818 9.7957329 2.97699333
</em>
</pre>
<p>* Variables that are correlated with PC1 and PC2 are the most important in explaining the variability in the data set.<br>
* The contribution of variables was extracted above: The larger the value of the contribution, the more the variable contributes to the component.</p>
<pre>
fviz_pca_var(hpi.pca, col.var=&quot;contrib&quot;,
             gradient.cols = c(&quot;#00AFBB&quot;, &quot;#E7B800&quot;, &quot;#FC4E07&quot;),
             repel = TRUE 
             )
</pre>
<p>Gives this plot:<br>
<a rel="nofollow" href="https://i0.wp.com/datascienceplus.com/wp-content/uploads/2017/11/life-6.png?ssl=1" target="_blank"><img sizes="(max-width: 490px) 100vw, 490px" src="https://i0.wp.com/datascienceplus.com/wp-content/uploads/2017/11/life-6-490x228.png?w=450&ssl=1" alt="" class="wp-image-9405" data-recalc-dims="1"></a><br>
This highlights the most important variables in explaining the variations retained by the principal components.</p>
<h3>Group countries by wealth, development, carbon emissions, and happiness</h3>
<p>When using clustering algorithms, k must be specified. I use the following method to help to find the best k.</p>
<pre>
number &lt;- NbClust(hpi[, 3:12], distance=&quot;euclidean&quot;,
               min.nc=2, max.nc=15, method='ward.D', index='all', alphaBeale = 0.1)
<em>*** : The Hubert index is a graphical method of determining the number of clusters.
In the plot of Hubert index, we seek a significant knee that corresponds to a significant increase of the value of the measure i.e the significant peak in Hubert
index second differences plot. *** : The D index is a graphical method of determining the number of clusters. In the plot of D index, we seek a significant knee (the significant peak in Dindex
second differences plot) that corresponds to a significant increase of the value of
the measure. ******************************************************************* * Among all indices: * 4 proposed 2 as the best number of clusters * 7 proposed 3 as the best number of clusters * 1 proposed 5 as the best number of clusters * 5 proposed 6 as the best number of clusters * 3 proposed 10 as the best number of clusters * 3 proposed 15 as the best number of clusters ***** Conclusion ***** * According to the majority rule, the best number of clusters is 3 </em>
</pre>
<p>I will apply K=3 in the following steps:</p>
<pre>
set.seed(2017)
pam &lt;- pam(hpi[, 3:12], diss=FALSE, 3, keep.data=TRUE)
fviz_silhouette(pam)
<em>cluster size ave.sil.width
1 1 43 0.46
2 2 66 0.32
3 3 31 0.37
</em>
</pre>
<p>Number of countries assigned in each cluster</p>
<pre>
hpi$country[pam$id.med]
<em>[1] &quot;Liberia&quot; &quot;Romania&quot; &quot;Ireland&quot;
</em>
</pre>
<p>This prints out one typical country represents each cluster.</p>
<pre>
fviz_cluster(pam, stand = FALSE, geom = &quot;point&quot;,
             ellipse.type = &quot;norm&quot;)
</pre>
<p>Gives this plot:<br>
<a rel="nofollow" href="https://i2.wp.com/datascienceplus.com/wp-content/uploads/2017/11/life-8.png?ssl=1" target="_blank"><img sizes="(max-width: 490px) 100vw, 490px" src="https://i1.wp.com/datascienceplus.com/wp-content/uploads/2017/11/life-8-490x228.png?w=450&ssl=1" alt="" class="wp-image-9406" data-recalc-dims="1"></a><br>
It is always a good idea to look at the cluster results, see how these three clusters were assigned.</p>
<p><b>A World map of three clusters</b></p>
<pre>
hpi['cluster'] &lt;- as.factor(pam$clustering)
map &lt;- map_data(&quot;world&quot;)
map &lt;- left_join(map, hpi, by = c('region' = 'country'))
ggplot() + geom_polygon(data = map, aes(x = long, y = lat, group = group, fill=cluster, color=cluster)) +
  labs(title = &quot;Clustering Happy Planet Index&quot;, subtitle = &quot;Based on data from:http://happyplanetindex.org/&quot;, x=NULL, y=NULL) + theme_minimal()
</pre>
<p>Gives this plot:<br>
<a rel="nofollow" href="https://i1.wp.com/datascienceplus.com/wp-content/uploads/2017/11/life-9.png?ssl=1" target="_blank"><img sizes="(max-width: 490px) 100vw, 490px" src="https://i0.wp.com/datascienceplus.com/wp-content/uploads/2017/11/life-9-490x228.png?w=450&ssl=1" alt="" class="wp-image-9407" data-recalc-dims="1"></a></p>
<h2>Summary</h2>
<p>The Happy Planet index has been criticized for weighting the ecological footprint too heavily; and the ecological footprint is a controversial concept. In addition, the Happy Planet Index has been misunderstood as a measure of personal “happiness”, when in fact, it is a measure of the “happiness” of the planet.</p>
<p>Nevertheless, the Happy Planet Index has been a consideration in the political area. For us, it is useful because it combines well being and environmental aspects, and it is simple and understandable. Also, it is available online, so we can create a story out of it.</p>
<p>Source code that created this post can be found <a rel="nofollow" href="https://github.com/susanli2016/Data-Analysis-with-R/blob/master/happy_planet_index.Rmd" target="_blank">here</a>. I am happy to hear any feedback or questions.</p>
<ul></ul>
<p><strong>
<p>Related Post</p>
<p></p></strong>
<ol>
<li><a rel="nofollow" href="https://datascienceplus.com/exploring-clustering-and-mapping-torontos-crimes/" title="Permanent Link toExploring, Clustering, and Mapping Toronto’s Crimes" target="_blank">Exploring, Clustering, and Mapping Toronto’s Crimes</a></li>
<li><a rel="nofollow" href="https://datascienceplus.com/spring-budget-2017-circle-visualisation/" title="Permanent Link toSpring Budget 2017: Circle Visualisation" target="_blank">Spring Budget 2017: Circle Visualisation</a></li>
<li><a rel="nofollow" href="https://datascienceplus.com/qualitative-research-in-r/" title="Permanent Link toQualitative Research in R" target="_blank">Qualitative Research in R</a></li>
<li><a rel="nofollow" href="https://datascienceplus.com/multi-dimensional-reduction-and-visualisation-with-t-sne/" title="Permanent Link toMulti-Dimensional Reduction and Visualisation with t-SNE" target="_blank">Multi-Dimensional Reduction and Visualisation with t-SNE</a></li>
<li><a rel="nofollow" href="https://datascienceplus.com/comparing-trump-and-clintons-facebook-pages-during-the-us-presidential-election-2016/" title="Permanent Link toComparing Trump and Clinton’s Facebook pages during the US presidential election, 2016" target="_blank">Comparing Trump and Clinton’s Facebook pages during the US presidential election, 2016</a></li>
</ol>
	<p>
<div>
<div>To <strong>leave a comment</strong> for the author, please follow the link and comment on their blog: <strong><a href="https://datascienceplus.com/how-happy-is-your-country-visualized/"> R Programming – DataScience+</a></strong>.</div>
<hr>
<a rel="nofollow" href="https://www.r-bloggers.com/">R-bloggers.com</a> offers <strong><a rel="nofollow" href="https://feedburner.google.com/fb/a/mailverify?uri=RBloggers">daily e-mail updates</a></strong> about <a rel="nofollow" href="https://www.r-project.org/" title="The R Project for Statistical Computing">R</a> news and <a rel="nofollow" href="https://www.r-bloggers.com/search/tutorial" title="R tutorials">tutorials</a> on topics such as: <a rel="nofollow" href="https://www.r-bloggers.com/search/data%20science" title="Data science">Data science</a>, <a rel="nofollow" href="https://www.r-bloggers.com/search/Big%20Data" title="Big Data">Big Data, <a rel="nofollow" href="https://www.r-users.com/" title="R jobs">R jobs</a>, visualization (<a rel="nofollow" href="https://www.r-bloggers.com/search/ggplot2" title="ggplot and ggplot2 tutorials">ggplot2</a>, <a rel="nofollow" href="https://www.r-bloggers.com/search/boxplot" title="Boxplots using lattice and ggplot2 tutorials">Boxplots</a>, <a rel="nofollow" href="https://www.r-bloggers.com/search/map" title="Maps and gis">maps</a>, <a rel="nofollow" href="https://www.r-bloggers.com/search/animation" title="Animation in R">animation</a>), programming (<a rel="nofollow" href="https://www.r-bloggers.com/search/RStudio" title="RStudio IDE for R">RStudio</a>, <a rel="nofollow" href="https://www.r-bloggers.com/search/sweave" title="Sweave and literate programming">Sweave</a>, <a rel="nofollow" href="https://www.r-bloggers.com/search/LaTeX" title="LaTeX in R">LaTeX</a>, <a rel="nofollow" href="https://www.r-bloggers.com/search/SQL" title="SQL and databases">SQL</a>, <a rel="nofollow" href="https://www.r-bloggers.com/search/eclipse" title="Eclipse IDE for R">Eclipse</a>, <a rel="nofollow" href="https://www.r-bloggers.com/search/git" title="git and github, Version Control System">git</a>, <a rel="nofollow" href="https://www.r-bloggers.com/search/hadoop" title="Large data in R using Hadoop">hadoop</a>, <a rel="nofollow" href="https://www.r-bloggers.com/search/Web+Scraping" title="Web Scraping of google, facebook, yahoo, twitter and more using R">Web Scraping</a>) statistics (<a rel="nofollow" href="https://www.r-bloggers.com/search/regression" title="Regressions and ANOVA analysis tutorials">regression</a>, <a rel="nofollow" href="https://www.r-bloggers.com/search/PCA" title="principal component analysis tutorial">PCA</a>, <a rel="nofollow" href="https://www.r-bloggers.com/search/time+series" title="Time series">time series</a>, <a rel="nofollow" href="https://www.r-bloggers.com/search/trading" title="finance trading">trading</a>) and more...
<img width="1" alt="" src="http://feeds.feedburner.com/~r/RBloggers/~4/IrS2KXlGM3s" height="1"></a></div>