ggplot2 2.0.0

ggplot2 2.0.0


I have used the ggplot2 package for a while and really like this package. It's happy to see that Hadley Wickham has officially updated the ggplot2 to version 2.0.0. On the RStudio Blog, Hadley highlights several important changes:

• ggplot2 now has an official extension mechanism.
• There are a handful of new geoms, and updates to existing geoms.
• The default appearance has been thoroughly tweaked so most plots should look better.
• Facets have a much richer set of labelling options.
• The documentation has been overhauled to be more helpful, and require less integration across multiple pages.
• A number of older and less used features have been deprecated.

You can find the document/manua on the project website. Many times, I go to the dev website to find the latest document/vignettes (extension, aesthetic specifications, themes).

The R Graphics Cookbook by Winston Chang is a must-have book to learn and become an expert of ggplot2 user. You can find the codes here from the Cookbook-R, and the Google book here.

general linear models vs. generalized linear models

General linear models vs. generalized linear models

• Wikipedia: Comparison of general and generalized linear models (modified)

	General linear model	Generalized linear model
Typical estimation method	Least squares, best linear unbiased prediction	Maximum likelihood or Bayesian
Special cases	ANOVA, ANCOVA, MANOVA, MANCOVA, ordinary linear regression, mixed model, t-test, F-test	linear regression, logistic regression, Poisson regression, gamma regression
Function in R	lm()	glm()
Function in Matlab	mvregress()	glmfit()
Procedure in SAS	PROC GLM, PROC MIXED	PROC GENMOD (PROC LOGISTIC for logistic regression only), PROC GLIMMIX Comparing the MIXED and GLIMMIX
Command in Stata	regress	Stata Generalized linear models
Function in Mathematica	LinearModelFit	GeneralizedLinearModelFit
Command in EViews	ls

• Generalized linear models have the flexiblility for response variables that have other than a normal distribution. If a generalized linear model uses an identity link function and a normal family distribution, then this model is equivalent to a general linear model.
• Generalized linear mixed models have the flexibility to model random effects and correlated errors for nonmormal data.

non-probability sample

Non-Probability Sample


Definition

• Sage Publication (pdf)
• Wikipedia
• Laerd

Reflection

• Wretman (2010) Reflections on Probability vs Nonprobability Sampling (pdf)
• Doherty   (1994) Probability versus Non-Probability Sampling in Sample Surveys (pdf)

Video

• Washington Statistical Society Mini-Conference Presentations (9/9/2015, Vimeo)

All-cause mortality was increasing among US middle age whites

All-cause mortality was increasing among US middle age Whites


Title: Rising morbidity and mortality in midlife among white non-Hispanic Americans in the 21st century
Authors: Anne Case and Angus Deaton
Abstract: This paper documents a marked increase in the all-cause mortality of middle-aged white non-Hispanic men and women in the United States between 1999 and 2013. This change reversed decades of progress in mortality and was unique to the United States; no other rich country saw a similar turnaround. The midlife mortality reversal was confined to white non-Hispanics; black non-Hispanics and Hispanics at midlife, and those aged 65 and above in every racial and ethnic group, continued to see mortality rates fall. This increase for whites was largely accounted for by increasing death rates from drug and alcohol poisonings, suicide, and chronic liver diseases and cirrhosis. Although all education groups saw increases in mortality from suicide and poisonings, and an overall increase in external cause mortality, thosewith less education saw the most marked increases. Rising midlife mortality rates of white non-Hispanics were paralleled by increases in midlife morbidity. Self-reported declines in health, mental health, and ability to conduct activities of daily living, and increases in chronic pain and inability to work, as well as clinically measured deteriorations in liver function, all point to growing distress in this population. We comment on potential economic causes and consequences of this deterioration. Full text: PNAS

Related articles:

• Epimonitor: All Cause mortality rate takes surprising upward turn for middle age Whites in the US
• The Washington Post: Prestigious medical journals rejected stunning study on deaths among middle-aged whites
• The New York Times: Death rates rising for middle-aged White Americans, study finds
• Vox: A big study found spiking death rates for middle-aged white Americans. Critics say it's more complicated
• New York: Anne Case, Co-author of the White Mortality Paper, Responds to a New Critique of Its Approach to Gender (Updated)

How to recover a lost partition of a hard drive

How to recover a lost partition of a hard drive


There are two major reasons you might not see the drive letter of your computer: the logic drive letter lose  or partition table corrupted.

Try these steps first:

• Go to the 'cmd' window by holding the "Windows" key and press the "R" key
• Type and run 'diskmgmt.msc'
• "Disk Management" will be shown.
• If see a partition without a drive letter then right-click on it
• Select "Change Drive Letter and Paths..." 
• Click on "Add" button
• Select the drive letter and Click on OK.

If you cannot see the partition without a drive letter, the partition table may be corrupt, try the TestDisk after the original instruction here or abstracted steps below:

• Download the TestDisk
• Unzip and save it on the USB drive
• Run "testdisk_win"
• At the first window, select “No Log” and press the key
• Select which drive to analyse, choose “Proceed” and press key
• Select partition type (select "Intel" if it’s a PC) then press key
• Select “Analyse” then press key
• Select “Quick Search” at the next screen, then press  key
• Press key, if the partitions were created under Vista – press key if not.
• TestDisk should say “Structure OK”. If so, press key 
• Select “Write” and press key and press   key to confirm.
• "ok" to reboot the compute, press key
• Now, close TestDisk and RESTART the computer.

R documentation and Learning Resources

R! documentation and Learning Resources

• Blog: R commands and keyboard shortcuts
• Blog: R! Books
• Rdocumentation of Datacamp is a tool that helps you easily find and browse the documentation of all current and some past packages on CRAN.
• CRAN Task View organizes the packages into different groups such as Graphics, Survival Analysis, etc.
• R_note -- The Exploration of Statistical Software R
• R! Tips
• Resources of learning R, UCLA
• Quick-R is to help you quickly access this language in your work
• William N Venables, David M Smith, and the R Core Team. An Introduction to R (on-going updated)
• Daniel Navarro (2015). Learning Statistics with R
• Sharon Machlis (2013).Beginner's guide to R: Introduction (pdf), and (2016) Advanced Beginner's Guide to R (pdf). 60+ R resources to improve your data skills
• Emmanuel Paradis. R for Beginners
• Paul Johnson (2014). R Tips
• ZevRoss. Beautiful plotting in R: A ggplot2 Cheatsheet (pdf)
• Alboukadel Kassambara. ggplot2: Guide to Create Beautiful Graphics in R
• Jay Kerns (2010). Introduction to Probability and Statistics Using R
• Julian Faraway (2002) Practical Regression and Anova using R
• Benjamin Yakir (2011). Introduction to Statistical Thinking (With R, Without Calculus)
• Paul Helwson (2009). Multivariate Statistics with R
• Hadley Wickham. ggplot2 documents (Plotly ggplot2 docs with Plotly's open-source ggplotly converter)
• RStudio. Cheat Sheets, Webinars and Video
• Norman Matloff (2009). The Art of R Programming
• Gillespie (2016). Efficient R Programming
• Scherer: A ggplot2 tutorial for beautiful plotting in R
• Thomas Lumley (2011). Course Notes - Complex sampling and R use survey package. withReplicates {survey} computes variances by replicate weighting.
• Anthony Damico:  analyze survey data for free, unlocked public-use data sets, twotorials
• Aedin Culhane: Introduction to R, CDC, June 25-28th 2012 (it's one of my favorite training courses attended). Here is a Amazon list of R TextBooks by Aedin.
• Yiling Cheng (2016). Doing Basic Calculus using R! etc.
• Verzani (2002). simpleR - Using R for Introductory Statistics (pdf)
• Patrick Burns (1998). S Poetry
• Zhou (2010). Fun with the R Grid Package
• Torre-Reyna: R: A language and Environment for statistical computing (Full manual)
• Getting started in R~Stata notes on exploring data (pdf)
• Muenchen: R for SAS and SPSS Users (pdf). the list of SAS equivalent packages.
• Maindonald: Using R for Data Analysis and Graphics (pdf)
• Data Analysts Captivated by R’s Power from the New York Times.
• Kadane. Principles of Uncertainty (pdf)
• Kerns. Introduction to Probability and Statistics Using R (pdf)
• Gaston. Handling and Processing Strings in R
• Frank Harrell. Regression Modeling Strategies and the rms Package
• Effective Graphs with Microsoft R Open
• Roger Peng and Jeff Leek. Git and GitHub videos for beginners
• Roger Peng etc. Course materials of Data Science Specialization
• Kevin Markham (2014) Hands-on dplyr tutorial for faster data manipulation in R
• Hopper (2016).The Excel User's introduction to R!
• Rickert. Getting Started with Markov Chains
• Wikibook. R programming
• Walker (2014). International population pyramids with ggplot2
• Demos: A Language, not a Letter: Learning Statistics in R

Running and longevity

Running and Longevity

News

• DialyMail: Stop that binge jogging
• Long term study found slow joggers had the lowest rates of death
• Strenuous joggers were as likely to die as sedentary non-joggers
• Going jogging three times a week for no more than 2.4 hours is optimal
• Pace of the slow joggers corresponds to vigorous exercise and strenuous jogging corresponds to very vigorous exercise, researchers qualified
• ScienceDaily: Light jogging may be most optimal for longevity

Guidelines and Researches

• CDC: How much physical activity do you need?
• Schnohr (2015). Dose of Jogging and Long-Term Mortality: the Copenhagen City Heart Study.
• Mittlema (1993). Triggering of Acute Myocardial Infarction by Heavy Physical Exertion
• Cheng (2000): Physical activity and self-reported, physician-diagnosed osteoarthritis: is physical activity a risk factor?
• Cheng (2012): Prevalence of diagnosed arthritis and arthritis-attributable activity limitation among adults with and without diagnosed diabetes

The Roseto Mystery

The Roseto Mystery

A few days ago, I read a book, Outliers (2008), by Malcolm Gladwell. The book is quite insightful. All children, parents, grandparents (ha-ha), and educators can be benefited from reading it. The introduction chapter is about the Rosetans who living in the Pennsylvannia, which is more related to my work. The Rosetans are outliers who are eating a big portion of saturated fat including lard with some other bad habits such as smoking. However, “THESE PEOPLE WERE DYING OF OLD AGE. THAT’S IT.” The Gladwell showed some evidences of importance of a low stressful lifestyle, good supports from family, friends, and communities for the health and longevity.

Here are a few more stories related to the Rosetans who living in Roseto, Pennsylvania:

• The first Chapter of Outliers
• Joe (2011). The Roseto Effect – The Valley of the Roses
• Steward Wolf (1993). The Power of Clan
• Cal Gutkin (2009). Outliers: extended families, better health outcomes
• Rock Positano (2011). The Mystery of the Rosetan People

Cancer isn’t just bad luck

Cancer isn’t just bad luck
By Thomas Lumley

"
From Stuff, "Bad luck is responsible for two-thirds of adult cancer while the remaining cases are due to environmental risk factors and inherited genes, researchers from the Johns Hopkins Kimmel Cancer Center found."

...

So, in summary: the “two-thirds of cancers explained” is Just Wrong. Doing a mathematically correct calculation gives about one third. Doing a calculation that’s actually relevant to cancer in the population gives even smaller values. (update) That’s not to say that DNA replication errors are unimportant — the paper makes it clear that they are important.
"

The fulltext: Cancer isn't just bad luck