<p><div><div><div><div data-width="100" data-show-faces="false" data-send="true" data-layout="button_count" data-height="21" data-href="http://www.r-bloggers.com/analyzing-public-health-data-in-r/"></div></div><div><a data-url="http://www.r-bloggers.com/analyzing-public-health-data-in-r/" data-counturl="http://www.r-bloggers.com/analyzing-public-health-data-in-r/" data-count="horizontal" data-text="Analyzing Public Health Data in R" class="twitter-share-button" data-via="rbloggers" href="https://twitter.com/share"></a></div><div></div></div><div></div></div>
<div>
(This article was first published on <strong><a href="http://www.arilamstein.com/blog/2016/06/20/public-health-r/"> R – AriLamstein.com</a></strong>, and kindly contributed to <a rel="nofollow" href="http://www.r-bloggers.com/">R-bloggers)</a> </div>
<p><em>Today’s post is by <a target="_blank" href="http://www.arilamstein.com/?feed-stats-url=aHR0cDovL3Rob21hc3lva290YS5jb20v&feed-stats-url-post-id=2729">Thomas Yokota</a>, an epidemiologist in Hawaii. I’ve been corresponding with Thomas via email and telephone for a while. I asked Thomas if he could write an introduction to how R, mapping and open data are used in the public health community. This is his reply.</em></p>
<div>
<div>
<p>Bonus: Download the code from this post!</p>
</div>
<div><a target="_blank" href="http://www.arilamstein.com/blog/2016/06/20/public-health-r/#content_upgrade_popup_2764">Click Here</a></div>
</div>
<h2>What is public health?</h2>
<p>Public health is a broad and sometimes overreaching field. Consequently, I am asked, “what <em>IS</em> public health?” more so than “what school did you go to (in Hawaii)?” I also work for a hospital system where I am often mistaken as someone knowledgeable of clinical care. With that said, I’d like to start off with a quote by C.E. Taylor:</p>
<blockquote><p>“The focus of public health is the community. The ‘patient’ is a whole population unit. The shift in professional orientation which occurs as the unit of attention moves from the individual to the group must be clearly recognized and explicitly stated because it has led to many misunderstandings in the past.”</p>
</blockquote>
<p>Once it’s clearly understood that public health is <em>not</em> brain surgery and – according to the CDC – more about <a target="_blank" href="http://www.arilamstein.com/?feed-stats-url=aHR0cDovL3d3dy5jZGMuZ292L3BocHIvem9tYmllcy5odG0%3D&feed-stats-url-post-id=2729">zombie invasions</a>, the next response I usually get is, “so what does public health do if it doesn’t cure diabetes?” Such a loaded question. In a nutshell, the <a target="_blank" href="http://www.arilamstein.com/?feed-stats-url=aHR0cDovL3d3dy5uYXRpb25hbGFjYWRlbWllcy5vcmcvaG1kL1JlcG9ydHMvMTk4OC9UaGUtRnV0dXJlLW9mLVB1YmxpYy1IZWFsdGguYXNweA%3D%3D&feed-stats-url-post-id=2729">four core functions</a> of public health are as follows:</p>
<ul>
<li>Assessment</li>
<li>Policy development</li>
<li>Assurance</li>
<li>Communication</li>
</ul>
<p>In other words, public health professionals rely on studies and surveys, such as the <a target="_blank" href="http://www.arilamstein.com/?feed-stats-url=aHR0cDovL3d3dy5jZGMuZ292L2JyZnNzLw%3D%3D&feed-stats-url-post-id=2729">Behavioral Risk Factors Surveillance System </a>(BRFSS), to understand health trends. Consequently, campaigns and/or policies are developed to break trends when needed.</p>
<h2>So why R?</h2>
<p>R is a great asset for making powerful <a target="_blank" href="http://www.arilamstein.com/?feed-stats-url=aHR0cDovL2dncGxvdDIub3JnLw%3D%3D&feed-stats-url-post-id=2729">graphs</a> and <a target="_blank" href="http://www.arilamstein.com/?feed-stats-url=aHR0cDovL2Jsb2cucmV2b2x1dGlvbmFuYWx5dGljcy5jb20vMjAxNC8wMS9lYXN5LWRhdGEtbWFwcy13aXRoLXItdGhlLWNob3JvcGxldGhyLXBhY2thZ2UtLmh0bWw%3D&feed-stats-url-post-id=2729">maps</a>, and can complement current SAS/STATA/SPSS/etc. workflows. Proprietary statistical programs will always be public health’s frenemies in the area of research and complex survey analysis – it’s inevitable that we must sacrifice funding health fair incentives to pay for license fees – but let’s not write off R too soon. R is also desirable for pubic health data analysts (i.e., epidemiologists) who work in community organizations.</p>
<h2>BRFSS</h2>
<p>In laymen terms, many public health organizations and professionals cite the BRFSS when interested in health risk behaviors, health access and chronic disease prevalence. Introductions and examples have already been created on BRFSS for R, and have been<a target="_blank" href="http://www.arilamstein.com/?feed-stats-url=aHR0cDovL3d3dy5yLWJsb2dnZXJzLmNvbS9hbmFseXplLXRoZS1iZWhhdmlvcmFsLXJpc2stZmFjdG9yLXN1cnZlaWxsYW5jZS1zeXN0ZW0tYnJmc3Mtd2l0aC1yLWFuZC1tb25ldGRiLw%3D%3D&feed-stats-url-post-id=2729"> featured on R-bloggers</a>. Anthony Damico maintains a <a target="_blank" href="http://www.arilamstein.com/?feed-stats-url=aHR0cHM6Ly9naXRodWIuY29tL2FqZGFtaWNvL2FzZGZyZWUvdHJlZS9tYXN0ZXIvQmVoYXZpb3JhbCUyMFJpc2slMjBGYWN0b3IlMjBTdXJ2ZWlsbGFuY2UlMjBTeXN0ZW0%3D&feed-stats-url-post-id=2729">GitHub repository</a> on downloading and storing complex surveys. I recommend his solution for those who are new to R and/or complex survey analysis and because there is a significant speed increase when working on analyzing the BRFSS alongside <a target="_blank" href="http://www.arilamstein.com/?feed-stats-url=aHR0cHM6Ly93d3cubW9uZXRkYi5vcmcvSG9tZQ%3D%3D&feed-stats-url-post-id=2729">MonetDB</a>. I also recommend Dr. Thomas Lumley’s book <em>Complex Survey Analysis</em> to understand complex surveys.</p>
<h3>Analyzing BRFSS in R</h3>
<p>In this example, I am taking a couple of shortcuts to import the data and to test out the <a target="_blank" href="http://www.arilamstein.com/?feed-stats-url=aHR0cHM6Ly9naXRodWIuY29tL2dlcmduZXNzL3Nydnly&feed-stats-url-post-id=2729">srvyr </a>package because… pipes. It also helps me to cut down some of the code and instructions entailed to setting up R and MonetDB – it’s also been done before. It is important to note, however, that  <a target="_blank" href="http://www.arilamstein.com/?feed-stats-url=aHR0cDovL3Itc3VydmV5LnItZm9yZ2Uuci1wcm9qZWN0Lm9yZy9zdXJ2ZXkv&feed-stats-url-post-id=2729">Survey</a> package is doing all the heavy lifting.</p>
<h4>load dependencies and data</h4>
<p>I use <a target="_blank" href="http://www.arilamstein.com/?feed-stats-url=aHR0cHM6Ly9jcmFuLnItcHJvamVjdC5vcmcvd2ViL3BhY2thZ2VzL3BhY21hbi9pbmRleC5odG1s&feed-stats-url-post-id=2729">pacman </a>package when handling all additional package loading and installs. I find it to be more convenient than <code>library()</code>. I am also borrowing Anthony’s download script to pull the BRFSS from CDC before loading it into R.</p>
<pre>
install.packages('pacman')
pacman::p_load(RCurl, foreign, downloader, survey, srvyr, ggplot2, dplyr)

source_url(&quot;https://raw.githubusercontent.com/ajdamico/asdfree/master/Download%20Cache/download%20cache.R&quot;, prompt=F, echo=F)
# download ez-pz brought to you by anthony joseph damico [ajdamico@gmail.com]

tf &lt;- tempfile(); td &lt;- tempdir()
xpt &lt;- &quot;http://www.cdc.gov/brfss/annual_data/2014/files/LLCP2014XPT.ZIP&quot;
download_cached(xpt, tf, mode='wb')
local.fn &lt;- unzip(tf, exdir=td)

brfss14 &lt;- read.xport(local.fn)
save(brfss14, file=&quot;brfss14.rda&quot;)
# load(&quot;brfss14.rda&quot;)

state.id &lt;- read.csv(&quot;https://raw.githubusercontent.com/tyokota/complexsurvey/master/brfss_state.csv&quot;, stringsAsFactors=F)
ageg65yr.id &lt;- read.csv(&quot;https://raw.githubusercontent.com/tyokota/complexsurvey/master/ageg65yr.csv&quot;, stringsAsFactors=F)
</pre>
<h4>create the survey object</h4>
<p>After importing the data, I use the <a target="_blank" href="http://www.arilamstein.com/?feed-stats-url=aHR0cDovL3Itc3VydmV5LnItZm9yZ2Uuci1wcm9qZWN0Lm9yZy9zdXJ2ZXkv&feed-stats-url-post-id=2729">Survey</a> package to utilize the BRFSS sampling design and weighting so that we can create estimates that represent the population. The <a target="_blank" href="http://www.arilamstein.com/?feed-stats-url=aHR0cDovL3d3dy5jZGMuZ292L2JyZnNzL2FubnVhbF9kYXRhLzIwMTQvcGRmL3dlaWdodGluZy1kYXRhLnBkZg%3D%3D&feed-stats-url-post-id=2729">raking weighting methodology</a> should be understood before proceeding. In this example, I will look at BMI categories as inferred by the following variable and question: BMI5CAT-<em>Computed body mass index categories</em>. I am re-coding BMI5CAT so that I can create estimates for overweight and/or obese prevalence.</p>
<pre>
brfss.design14 &lt;- brfss14 %>%
  as_survey_design(ids=X_PSU, weight=X_LLCPWT, nest=TRUE, strata=X_STSTR, variables= c(X_BMI5CAT, X_MRACE1, X_STATE))

brfss.design14 &lt;- brfss.design14 %>%
  mutate(X_BMI5CAT2 = car::recode(X_BMI5CAT, &quot;c(3,4)=1; NA=NA;  else=0&quot;),
         X_BMI5CAT=as.factor(X_BMI5CAT))
</pre>
<h4>state-by-state comparison</h4>
<p>Now that we have created the survey design object, we can use R to create estimates from the data. Before doing this, however, it is advised that the reader studies the BRFSS code book. Another good source to learning how to use the BRFSS is to study past reports and case studies. You can always verify your query against the <a target="_blank" href="http://www.arilamstein.com/?feed-stats-url=aHR0cDovL3d3dy5jZGMuZ292L2JyZnNzL2JyZnNzcHJldmFsZW5jZS9pbmRleC5odG1s&feed-stats-url-post-id=2729">CDC BRFSS Prevalence &amp; Trends Tool</a>.</p>
<pre>
options(survey.lonely.psu = &quot;certainty&quot;)

brfss.design14.0 %>%
  filter(X_STATE==15) %>%
  group_by(X_STATE, X_BMI5CAT) %>%
  summarize(prevalence = survey_mean(na.rm=T, vartype = c(&quot;ci&quot;)),
            N = survey_total(na.rm=T)) # matches CDC tool

brfss.design14 &lt;- brfss.design14 %>%
  mutate(X_BMI5CAT2=as.factor(X_BMI5CAT2))

BMI5CAT2.1 &lt;- brfss.design14 %>%
  group_by(X_STATE, X_BMI5CAT2) %>%
  summarize(prevalence = survey_mean(na.rm=T),
            N = survey_total(na.rm=T)) %>%
  mutate(X_STATE = as.character(X_STATE)) %>%
  left_join(state.id, by=c(&quot;X_STATE&quot;=&quot;VALUE&quot;))
</pre>
<p>Most commonly, BRFSS is used to make state-by-state comparisons. In the past (and sometimes currently), state-by-state comparisons were often ranked. This made for great click-bait news headlines. Most recently, choropleth maps have been used to visualize prevalence rates across the nation.</p>
<pre>
BMI5CAT2.1 %>%
  filter(X_BMI5CAT2==1) %>%
  ggplot(aes(x=reorder(STATE, prevalence), y=prevalence)) +
  scale_y_continuous(labels=scales::percent) +
  geom_bar(stat=&quot;identity&quot;) +
  ylab(&quot;YES (%)&quot;) +
  xlab(&quot;STATE&quot;) +
  ggtitle(&quot;Computed body mass index categories (overweight or obese) by state&quot;) +
  theme(axis.text.x = element_text(angle = 90, hjust = 1),
        panel.background = element_rect(fill = 'white' )) +
  guides(fill=FALSE)
</pre>
<p><img alt="state overweight and obese" src="http://i1.wp.com/raw.githubusercontent.com/tyokota/complexsurvey/master/state_overweightobese.png?w=676&ssl=1" data-recalc-dims="1"></p>
<pre>
BMI5CAT2.2 &lt;- BMI5CAT2.1 %>%
  filter(X_BMI5CAT2==1) %>%
  select(region=STATE, value=prevalence) %>%
  mutate(region = tolower(region))

choroplethr::state_choropleth(BMI5CAT2.2, title=&quot;Computed body mass index categories (overweight or obese) choropleth map&quot;,
                              num_colors=9)
</pre>
<p><img alt="state overweight and obese map" src="http://i1.wp.com/raw.githubusercontent.com/tyokota/complexsurvey/master/chorplethr.png?w=676&ssl=1" data-recalc-dims="1"></p>
<blockquote><p>I’m not a big fan of this map projection, and hope there will be an option in the future to change the projection.</p>
</blockquote>
<h4>deep dive Hawaii – stratify by race</h4>
<p>Before we start dropping leaflets on healthy eating across the Bible Belt, why don’t we dig into the data a bit more. Let’s take my home state for example. Hawaii appears to have fewer obesity and overweight prevalence than most of the nation – California, drool at our sculpted natural beach bodies. Seems like there’s no issues here…</p>
<pre>
brfss14a &lt;- brfss14 %>%
  filter(X_STATE==15)

brfss.design14a &lt;- brfss14a %>%
  as_survey_design(ids=X_PSU, weight=X_LLCPWT, nest=TRUE, strata=X_STSTR, variables= c(X_BMI5CAT, X_MRACE1, X_STATE, X_AGEG5YR))

brfss.design14a &lt;- brfss.design14a %>%
  mutate(X_BMI5CAT2 = car::recode(X_BMI5CAT, &quot;c(3,4)=1; NA=NA;  else=0&quot;),
         X_BMI5CAT2a=car::recode(X_BMI5CAT2, &quot;1='Obese/overweight'; NA=NA; else='Not Obese/overweight'&quot;),
         X_BMI5CAT2a=factor(X_BMI5CAT2a, levels=c('Obese/overweight', 'Not Obese/overweight'), ordered=TRUE),
         X_MRACE1a=car::recode(X_MRACE1, &quot;1='White'; 2='Black'; 3='AIAN'; 4='Asian'; 5='NHOPI';  NA=NA; c(6,7)='other/multiracial'; else=NA&quot;),
         X_AGEG5YR=as.factor(X_AGEG5YR))

BMI5CAT2.3 &lt;- brfss.design14a %>%
  group_by(X_STATE, X_MRACE1a, X_BMI5CAT2a) %>%
  summarize(prevalence = survey_mean(na.rm=T),
            N = survey_total(na.rm=T),
            prevalence_ci = survey_mean(na.rm=T, vartype = c(&quot;ci&quot;))) %>%
  mutate(X_STATE = as.character(X_STATE)) %>%
  left_join(state.id, by=c(&quot;X_STATE&quot;=&quot;VALUE&quot;))

BMI5CAT2.3 %>%
  ggplot(aes(x = X_MRACE1a, y = prevalence, group = X_BMI5CAT2a)) +
    geom_bar(stat = &quot;identity&quot;, aes(fill = X_BMI5CAT2a), alpha = 0.3) +
  scale_y_continuous(labels=scales::percent) +
  ylab(&quot;Prevalence (%)&quot;) +
  xlab(&quot;Race&quot;) +
  ggtitle(&quot;Computed body mass index categories (overweight or obese) in Hawaii by race&quot;) +
  theme(axis.text.x = element_text(angle = 90, hjust = 1),
        panel.background = element_rect(fill = 'white')) +
  guides(fill = guide_legend(title = &quot;Categories&quot;))
</pre>
<p><img alt="state of HI overweight and obese by race" src="http://i2.wp.com/raw.githubusercontent.com/tyokota/complexsurvey/master/hawaiirace.png?w=676&ssl=1" data-recalc-dims="1"></p>
<p>Stratifying by race shows us a different picture: NHOPI [ed note: NHOPI stands for (N)ative (H)awaiian &amp; (O)ther (P)acific (I)slander] overweight/obese prevalence is at 80% (CI95%:  73-87%), easily topping all point estimates by state.</p>
<h4>deep dive Hawaii – stratify by race and age group</h4>
<p>Going one level deeper, we see that onset of overweight/obesity starts earlier for NHOPI and that prevalence across age groups is higher in comparison. The graph suggests that efforts to address overweight/obesity among NHOPI groups is needed in Hawaii; in addition, it also suggests that intervention at earlier stages would be beneficial for this group as well.</p>
<pre>
BMI5CAT2.4 &lt;- brfss.design14a %>%
  group_by(X_STATE, X_MRACE1a, X_AGEG5YR, X_BMI5CAT2a) %>%
  summarize(prevalence = survey_mean(na.rm=T),
            N = survey_total(na.rm=T),
            prevalence_ci = survey_mean(na.rm=T, vartype = c(&quot;ci&quot;))) %>%
  mutate(X_STATE = as.character(X_STATE)) %>%
  left_join(state.id, by=c(&quot;X_STATE&quot;=&quot;VALUE&quot;)) %>%
  mutate(X_STATE = as.character(X_STATE)) %>%
  left_join(ageg65yr.id %>% rename(X_AGEG5YR2=X_AGEG5YR), by=c(&quot;X_AGEG5YR&quot;=&quot;VALUE&quot;))

BMI5CAT2.4 %>%
  filter(X_MRACE1a==&quot;NHOPI&quot; | X_MRACE1a==&quot;White&quot; | X_MRACE1a==&quot;Asian&quot;, X_AGEG5YR2!=&quot;Dont know/Refused/Missing&quot;) %>%
  ggplot(aes(X_AGEG5YR2, prevalence)) +
    geom_bar(stat = &quot;identity&quot;, aes(fill = X_BMI5CAT2a), col = &quot;gray&quot;, alpha = 0.7) +
    facet_grid(~X_MRACE1a) +
    scale_y_continuous(labels=scales::percent) +
    ylab(&quot;Prevalence (%)&quot;) +
    xlab(&quot;Age Group&quot;) +
    ggtitle(&quot;Computed body mass index categories (overweight or obese) in Hawaii by race and age&quot;) +
    theme(axis.text.x = element_text(angle = 90, hjust = 1),
          panel.background = element_rect(fill = 'white' )) +
    guides(fill = guide_legend(title = &quot;Categories&quot;))
</pre>
<p><img alt="state of HI overweight and obese by age and race" src="http://i2.wp.com/raw.githubusercontent.com/tyokota/complexsurvey/master/hawaiiraceage.png?w=676&ssl=1" data-recalc-dims="1"></p>
<h4>Conclusion</h4>
<p>Relying on high-level estimates alone can sometimes mislead your agenda. For example, after stratifying by race, we saw that the shared benefits of living in Hawaii were not equally shared by all groups. Further stratifying by age shows that onset by race occurred earlier for NHOPIs. These cuts of the data can help to narrow down public health efforts and hopefully prevent the usual spray-and-pray approach to solving problems.</p>
<p>This example only scratches the surface; there are <a target="_blank" href="http://www.arilamstein.com/?feed-stats-url=aHR0cDovL3d3dy5yLWJsb2dnZXJzLmNvbS9zdGF0aXN0aWNhbGx5LXNpZ25pZmljYW50LXRyZW5kcy13aXRoLW11bHRpcGxlLXllYXJzLW9mLWNvbXBsZXgtc3VydmV5LWRhdGEv&feed-stats-url-post-id=2729">other methods that can be explored</a>. Analyzing data sets for public health such as BRFSS is more than just getting stuck on statistical models and debate on R vs. SAS vs. STATA. While it’s important that we continue to leverage statistical methods, analyzing data for public health is more than being a data monkey and computer babysitter. In recent years, with the onslaught of articles promoting “big data” and  “data is power”, I’ve experienced this elitism both in the workplace and in various communities. In my wiser years, I’ve realized that power doesn’t come from hogging the ability to analyze the data, but from contributing to the discussion as a team-player.</p>
<div>
<div>
<p>Bonus: Download the code from this post!</p>
</div>
<div><a target="_blank" href="http://www.arilamstein.com/blog/2016/06/20/public-health-r/#content_upgrade_popup_2764">Click Here</a></div>
</div>
<p> <img height="1" width="1" src="http://www.arilamstein.com/?feed-stats-post-id=2729">
<p>The post <a target="_blank" rel="nofollow" href="http://www.arilamstein.com/blog/2016/06/20/public-health-r/">Analyzing Public Health Data in R</a> appeared first on <a target="_blank" rel="nofollow" href="http://www.arilamstein.com/">AriLamstein.com</a>.</p>
<p><div><div><div><div data-width="55" data-show-faces="false" data-send="true" data-layout="box_count" data-height="62" data-href="http://www.r-bloggers.com/analyzing-public-health-data-in-r/"></div></div><div><a data-url="http://www.r-bloggers.com/analyzing-public-health-data-in-r/" data-counturl="http://www.r-bloggers.com/analyzing-public-health-data-in-r/" data-count="vertical" data-text="Analyzing Public Health Data in R" class="twitter-share-button" data-via="rbloggers" href="https://twitter.com/share"></a></div><div></div></div><div></div></div>
<div>
<div>To <strong>leave a comment</strong> for the author, please follow the link and comment on their blog: <strong><a href="http://www.arilamstein.com/blog/2016/06/20/public-health-r/"> R – AriLamstein.com</a></strong>.</div>
<hr>
<a rel="nofollow" href="http://www.r-bloggers.com/">R-bloggers.com</a> offers <strong><a rel="nofollow" href="http://feedburner.google.com/fb/a/mailverify?uri=RBloggers">daily e-mail updates</a></strong> about <a title="The R Project for Statistical Computing" rel="nofollow" href="http://www.r-project.org/">R</a> news and <a title="R tutorials" rel="nofollow" href="http://www.r-bloggers.com/?s=tutorial">tutorials</a> on topics such as: <a title="Data science" rel="nofollow" href="http://www.r-bloggers.com/?s=data%20science">Data science</a>, <a title="Big Data" rel="nofollow" href="http://www.r-bloggers.com/?s=Big%20Data">Big Data, <a title="R jobs" rel="nofollow" href="http://www.r-users.com/">R jobs</a>, visualization (<a title="ggplot and ggplot2 tutorials" rel="nofollow" href="http://www.r-bloggers.com/?s=ggplot2">ggplot2</a>, <a title="Boxplots using lattice and ggplot2 tutorials" rel="nofollow" href="http://www.r-bloggers.com/?s=boxplot">Boxplots</a>, <a title="Maps and gis" rel="nofollow" href="http://www.r-bloggers.com/?s=map">maps</a>, <a title="Animation in R" rel="nofollow" href="http://www.r-bloggers.com/?s=animation">animation</a>), programming (<a title="RStudio IDE for R" rel="nofollow" href="http://www.r-bloggers.com/?s=RStudio">RStudio</a>, <a title="Sweave and literate programming" rel="nofollow" href="http://www.r-bloggers.com/?s=sweave">Sweave</a>, <a title="LaTeX in R" rel="nofollow" href="http://www.r-bloggers.com/?s=LaTeX">LaTeX</a>, <a title="SQL and databases" rel="nofollow" href="http://www.r-bloggers.com/?s=SQL">SQL</a>, <a title="Eclipse IDE for R" rel="nofollow" href="http://www.r-bloggers.com/?s=eclipse">Eclipse</a>, <a title="git and github, Version Control System" rel="nofollow" href="http://www.r-bloggers.com/?s=git">git</a>, <a title="Large data in R using Hadoop" rel="nofollow" href="http://www.r-bloggers.com/?s=hadoop">hadoop</a>, <a title="Web Scraping of google, facebook, yahoo, twitter and more using R" rel="nofollow" href="http://www.r-bloggers.com/?s=Web+Scraping">Web Scraping</a>) statistics (<a title="Regressions and ANOVA analysis tutorials" rel="nofollow" href="http://www.r-bloggers.com/?s=regression">regression</a>, <a title="principal component analysis tutorial" rel="nofollow" href="http://www.r-bloggers.com/?s=PCA">PCA</a>, <a title="Time series" rel="nofollow" href="http://www.r-bloggers.com/?s=time+series">time series</a>, <a title="finance trading" rel="nofollow" href="http://www.r-bloggers.com/?s=trading">trading</a>) and more...
<img height="1" alt="" width="1" src="http://feeds.feedburner.com/~r/RBloggers/~4/3VbtN2isiX8"></a></div>