Non-iterative, regression-based estimation of haplotype associations
|
|
- Kimberly Casey
- 5 years ago
- Views:
Transcription
1 Non-iterative, regression-based estimation of haplotype associations Benjamin French, PhD Department of Biostatistics and Epidemiology University of Pennsylvania National Cancer Center Republic of Korea 12 September 2012
2 nationalatlas.gov TM Where We Are STATES WASHINGTON CA NADA ALASKA ATLANTIC OCEAN HAWAII MONTANA NORTH DAKOTA L a k e S u p e r i o r PACIFIC OCEAN MAINE P A C I F I C O C E A N OREGON CALIFORNIA NEVADA IDAHO ARIZONA UTAH WYOMING COLORADO NEW MEXICO SOUTH DAKOTA NEBRASKA KANSAS OKLAHOMA MINNESOTA IOWA MISSOURI ARKANSAS WISCONSIN ILLINOIS MISSISSIPPI MICHIGAN a n L a k e M i c h i g INDIANA L a k e H u r o n KENTUCKY TENNESSEE ALABAMA OHIO L a k e E r i e GEORGIA L a k e O n t a r i o NEW YORK PENNSYLVANIA WEST VIRGINIA VIRGINIA NORTH CAROLINA SOUTH CAROLINA DC DELAWARE MARYLAND A T L A VERMONT NEW HAMPSHIRE CONNECTICUT NEW JERSEY MASSACHUSETTS RHODE ISLAND E A N N T I C O C TEXAS HAWAII A R C T I C O C E A N LOUISIANA R U S S I A P A C I F I C O C E A N mi km ALASKA CANADA ME G U L F O F M E X I C O FLORIDA BERING SEA XIC THE BAHAMAS GULF OF ALASKA O km 300 mi P A C I F I C O C E A N mi 200 km CU BA U.S. Department of the Interior U.S. Geological Survey O The National Atlas of the United States of America R states2.pdf INTERIOR-GEOLOGICAL SURVEY, RESTON, VIRGINIA-2003
3 Collaborators Nandita Mitra, PhD Department of Biostatistics and Epidemiology University of Pennsylvania Thomas P Cappola, MD Penn Cardiovascular Institute University of Pennsylvania Thomas Lumley, PhD Department of Statistics University of Auckland B French (bcfrench@upenn.edu) Estimating haplotype associations September / 31
4 Goals To integrate haplotypes into large association studies such that haplotype imputation is done once as a data-processing step Case-control studies (binary outcome) [French et al., 2006] Prospective studies (censored survival outcome) [French et al., 2012] To allow haplotype associations to be estimated in general-purpose statistical software (eg R) by researchers expert in the subject matter In world historical terms there is a lot to be said for keeping data analysis out of the hands of statisticians Thomas Lumley B French (bcfrench@upenn.edu) Estimating haplotype associations September / 31
5 Goals To integrate haplotypes into large association studies such that haplotype imputation is done once as a data-processing step Case-control studies (binary outcome) [French et al., 2006] Prospective studies (censored survival outcome) [French et al., 2012] To allow haplotype associations to be estimated in general-purpose statistical software (eg R) by researchers expert in the subject matter In world historical terms there is a lot to be said for keeping data analysis out of the hands of statisticians Thomas Lumley B French (bcfrench@upenn.edu) Estimating haplotype associations September / 31
6 Phase ambiguity Observed data is composed of a set of unphased genotypes Diplotype (pair of haplotypes) may be ambiguous; may not know which allele was transmitted from maternal or paternal chromosome Missing data problem; impute the unobserved diplotype Father A T Mother C G Child A C T G AT/CG or AG/CT B French (bcfrench@upenn.edu) Estimating haplotype associations September / 31
7 Haplotype imputation Expectation-maximization (EM) algorithm E: calculate expected phase given haplotype frequencies M: calculate MLEs for haplotype frequencies given phase Software: haplo.stats [Sinnwell and Schaid, 2012] Bayesian algorithm Observed genotype data combined with expected haplotype patterns Haplotypes estimated from posterior distribution Software: PHASE [Stephens and Donnelly, 2003] B French (bcfrench@upenn.edu) Estimating haplotype associations September / 31
8 Diplotype uncertainty Angiotensin II receptor type 1 (AGTR1) Haplotype Diplotype Label Haplotype frequency probability D TCCACGCATCTT F TCTGTGCATCTC C TCCACGCATCTC G TCTGTGCATCTT Rare TCCGCGCATCTC < < 0.01 Rare TCTATGCATCTT < B French (bcfrench@upenn.edu) Estimating haplotype associations September / 31
9 Target of inference If diplotypes were known, could fit a standard regression model including diplotypes D i and environmental exposures Z i (possibly with interaction) Logistic regression: Y i = {1 = case, 0 = control} logit P[Y i = 1] = α + D i β D + Z i β Z + Cox regression: Y i = min(t i, C i ); event indicator δ i = 1[Y i = T i ] log λ i (t) = log λ 0 (t) + D i β D + Z i β Z + in which β = {β D, β Z } represent association parameters of interest We integrate the set of imputed haplotypes into these regression models B French (bcfrench@upenn.edu) Estimating haplotype associations September / 31
10 Estimating associations Non-iterative weighted estimation [French et al., 2006] 1. Impute haplotypes and estimate population haplotype frequencies 2. Create multi-record data for each individual Design matrix: set of diplotypes consistent with observed genotype, possibly including environmental exposures Weights equal to conditional probability of each diplotype Weight A B C D E F G H I Rare Estimate associations using a weighted regression model Logistic regression for binary outcomes Cox regression for censored survival outcomes Robust or sandwich standard error estimator Account for uncertainty in phase B French (bcfrench@upenn.edu) Estimating haplotype associations September / 31
11 Estimating associations Non-iterative weighted estimation [French et al., 2006] 1. Impute haplotypes and estimate population haplotype frequencies 2. Create multi-record data for each individual Design matrix: set of diplotypes consistent with observed genotype, possibly including environmental exposures Weights equal to conditional probability of each diplotype Weight A B C D E F G H I Rare Estimate associations using a weighted regression model Logistic regression for binary outcomes Cox regression for censored survival outcomes Robust or sandwich standard error estimator Account for uncertainty in phase B French (bcfrench@upenn.edu) Estimating haplotype associations September / 31
12 Estimating associations Non-iterative weighted estimation Estimates are obtained from estimating functions Logistic regression for binary outcomes n U(β) = π id (p id, β)x id [Y id expit X id β] i=1 d d(g i ) Cox regression for censored survival outcomes n U(β) = i=1 d d(g i ) l id (β) β l id (β) = δ id π id (p id, β) X id β log exp X i dβ i R(t i ) Inference based on univariable or multivariable Wald tests Estimator is valid and efficient under the null hypothesis B French (bcfrench@upenn.edu) Estimating haplotype associations September / 31
13 Estimating associations Weighted haplotype combination [Rebbeck et al., 2009] 1. Impute haplotypes and estimate population haplotype frequencies 2. Create single-record data for each individual Design matrix: weighted combination of diplotypes consistent with observed genotype, with weights equal to conditional probability of each diplotype, possibly including environmental exposures A B C D E F G H I Rare Estimate associations using an unweighted regression model Logistic regression for binary outcomes Cox regression for censored survival outcomes Account for uncertainty in phase B French (bcfrench@upenn.edu) Estimating haplotype associations September / 31
14 Estimating associations Weighted haplotype combination [Rebbeck et al., 2009] 1. Impute haplotypes and estimate population haplotype frequencies 2. Create single-record data for each individual Design matrix: weighted combination of diplotypes consistent with observed genotype, with weights equal to conditional probability of each diplotype, possibly including environmental exposures A B C D E F G H I Rare Estimate associations using an unweighted regression model Logistic regression for binary outcomes Cox regression for censored survival outcomes Account for uncertainty in phase B French (bcfrench@upenn.edu) Estimating haplotype associations September / 31
15 Estimating associations Weighted haplotype combination Estimates are obtained from estimating functions Logistic regression for binary outcomes U(β) = n X i [Y i expit X i β] i=1 Cox regression for censored survival outcomes U(β) = n i=1 l i (β) β l i (β) = δ i X i β log i R(t i ) exp X i β Inference based on univariable or multivariable Wald tests Estimator is exactly valid under the null hypothesis B French (bcfrench@upenn.edu) Estimating haplotype associations September / 31
16 Simulation study In which situations do non-iterative methods provide reliable results? Angiotensin II receptor type 1 (AGTR1) Label Haplotype Frequency A ATTATGCATCTC B ATTATGTGATCC C TCCACGCATCTC D TCCACGCATCTT E TCTGTGCAACTT F* TCTGTGCATCTC G TCTGTGCATCTT H TTTACACATCTC I TTTACACATCTT * Referent n = 500, 25% censoring (independent) B French (bcfrench@upenn.edu) Estimating haplotype associations September / 31
17 Simulation results: Moderate effects Approximate bias Weighted estimation Weighted combination Known phase A 1/2.0 B 1/1.75 C 1.25 D 1.5 E 1/1.25 G 1/1.5 H 1.75 I 2.0 Hazard ratio Haplotype B French (bcfrench@upenn.edu) Estimating haplotype associations September / 31
18 Simulation results: Moderate effects % coverage of 95% confidence intervals Hazard Weighted Weighted Known Haplotype ratio estimation combination phase A 1/ B 1/ C D E 1/ G 1/ H I B French (bcfrench@upenn.edu) Estimating haplotype associations September / 31
19 Simulation results: Moderate effects Rejection rate (%) of two-sided hypothesis tests Hazard Weighted Weighted Known Haplotype ratio estimation combination phase A 1/ B 1/ C D E 1/ G 1/ H I B French (bcfrench@upenn.edu) Estimating haplotype associations September / 31
20 Simulation results: Strong SNP effects Approximate bias Weighted estimation Weighted combination Known phase A 1.0 B 1.0 C 1.0 D 4.0 E 4.0 G 4.0 H 1.0 I 4.0 Hazard ratio Haplotype B French (bcfrench@upenn.edu) Estimating haplotype associations September / 31
21 Simulation results: Strong SNP effects % coverage of 95% confidence intervals Hazard Weighted Weighted Known Haplotype ratio estimation combination phase A B C D E G H I B French (bcfrench@upenn.edu) Estimating haplotype associations September / 31
22 Simulation results: Strong SNP effects Rejection rate (%) of two-sided hypothesis tests Hazard Weighted Weighted Known Haplotype ratio estimation combination phase A B C D E G H I B French (bcfrench@upenn.edu) Estimating haplotype associations September / 31
23 Simulation results: Strong non-snp effects Approximate bias Weighted estimation Weighted combination Known phase A 1.0 B 1.0 C 4.0 D 1.0 E 4.0 G 4.0 H 4.0 I 1.0 Hazard ratio Haplotype B French (bcfrench@upenn.edu) Estimating haplotype associations September / 31
24 Simulation results: Strong non-snp effects % coverage of 95% confidence intervals Hazard Weighted Weighted Known Haplotype ratio estimation combination phase A B C D E G H I B French (bcfrench@upenn.edu) Estimating haplotype associations September / 31
25 Simulation results: Strong non-snp effects Rejection rate (%) of two-sided hypothesis tests Hazard Weighted Weighted Known Haplotype ratio estimation combination phase A B C D E G H I B French (bcfrench@upenn.edu) Estimating haplotype associations September / 31
26 Application HSPB7-CLCNKA haplotypes and adverse events in chronic heart failure Regulate renal potassium channels to control blood pressure SNP associated with heart failure in a large case-control study [Cappola et al., 2011] Genotypes available for 1150 genetically inferred Caucasians with heart failure enrolled in a prospective study 70% male; median age at study entry, 58 years 14 pre-selected SNPs Inferred 10 common haplotypes (frequency > 0.02) 65% had an unambiguous diplotype 90% had a highest posterior probability > B French (bcfrench@upenn.edu) Estimating haplotype associations September / 31
27 Application HSPB7-CLCNKA haplotypes and adverse events in chronic heart failure Outcome: time to all-cause mortality or cardiac transplantation Median follow-up, 3 years; maximum, 5 years 22% experienced an adverse event Cox regression framework Stratified by 4-level classification for disease severity Adjusted for gender, age, heart failure etiology, clinical site Time-varying covariate for age (exhibited non-proportional hazards) B French (bcfrench@upenn.edu) Estimating haplotype associations September / 31
28 Application results: Weighted estimation Label Haplotype Frequency HR (95% CI) P Q AGAGCGAGACGAGG (0.80, 1.77) 0.39 R AGAGCGAGGGAAGG (0.80, 1.35) 0.79 S AGAGCGGAGCAAGA (0.80, 1.77) 0.38 T AGCGAGAGGCAAGA (0.34, 0.88) 0.01 U GACGCGGAGCGCGG (0.53, 1.20) 0.28 V GGAACAAGGGAAGG (0.26, 0.92) 0.03 W GGAACAGAGCAAGA Referent X GGAACAGAGCAAGG (0.95, 1.95) 0.09 Y GGAGCAAGGCAAGG (0.78, 1.69) 0.49 Z GGCGCGGAGCAAGG (0.62, 1.92) 0.76 Overall 0.02 B French (bcfrench@upenn.edu) Estimating haplotype associations September / 31
29 Application results: Weighted estimation Label Haplotype Frequency HR (95% CI) P Q AGAGCGAGACGAGG (0.80, 1.77) 0.39 R AGAGCGAGGGAAGG (0.80, 1.35) 0.79 S AGAGCGGAGCAAGA (0.80, 1.77) 0.38 T AGCGAGAGGCAAGA (0.34, 0.88) 0.01 U GACGCGGAGCGCGG (0.53, 1.20) 0.28 V GGAACAAGGGAAGG (0.26, 0.92) 0.03 W GGAACAGAGCAAGA Referent X GGAACAGAGCAAGG (0.95, 1.95) 0.09 Y GGAGCAAGGCAAGG (0.78, 1.69) 0.49 Z GGCGCGGAGCAAGG (0.62, 1.92) 0.76 Overall 0.02 B French (bcfrench@upenn.edu) Estimating haplotype associations September / 31
30 Application results: Weighted combination Label Haplotype Frequency HR (95% CI) P Q AGAGCGAGACGAGG (0.76, 1.87) 0.44 R AGAGCGAGGGAAGG (0.81, 1.37) 0.73 S AGAGCGGAGCAAGA (0.74, 1.94) 0.47 T AGCGAGAGGCAAGA (0.32, 0.90) 0.02 U GACGCGGAGCGCGG (0.52, 1.21) 0.29 V GGAACAAGGGAAGG (0.21, 0.88) 0.02 W GGAACAGAGCAAGA Referent X GGAACAGAGCAAGG (0.95, 2.19) 0.09 Y GGAGCAAGGCAAGG (0.77, 1.72) 0.49 Z GGCGCGGAGCAAGG (0.64, 1.87) 0.73 Overall 0.03 B French (bcfrench@upenn.edu) Estimating haplotype associations September / 31
31 Application results: Weighted combination Label Haplotype Frequency HR (95% CI) P Q AGAGCGAGACGAGG (0.76, 1.87) 0.44 R AGAGCGAGGGAAGG (0.81, 1.37) 0.73 S AGAGCGGAGCAAGA (0.74, 1.94) 0.47 T AGCGAGAGGCAAGA (0.32, 0.90) 0.02 U GACGCGGAGCGCGG (0.52, 1.21) 0.29 V GGAACAAGGGAAGG (0.21, 0.88) 0.02 W GGAACAGAGCAAGA Referent X GGAACAGAGCAAGG (0.95, 2.19) 0.09 Y GGAGCAAGGCAAGG (0.77, 1.72) 0.49 Z GGCGCGGAGCAAGG (0.64, 1.87) 0.73 Overall 0.03 B French (bcfrench@upenn.edu) Estimating haplotype associations September / 31
32 R packages: Non-iterative weighted estimation haplo.ccs [French and Lumley, 2011] Weighted logistic regression for binary outcomes Depends on haplo.stats package to impute haplotypes Calls glm(..., family=quasibinomial(link=logit)) Includes GEE-type sandwich standard error estimator haplo.cph (in process) Weighted Cox regression for censored survival outcomes Will depend on haplo.stats package to impute haplotypes Will call cph(..., robust=true) from Design package Allow stratification and time-varying exposures B French (bcfrench@upenn.edu) Estimating haplotype associations September / 31
33 Conclusions Non-iterative estimation methods Incorporate all diplotypes consistent with the observed genotype, either weighted (weighted estimation) or averaged (weighted combination) by the conditional probability of each diplotype Provide valid tests for genetic associations and reliable estimates of modest genetic effects of common haplotypes Convenient regression framework Adjustment for or interaction with environmental exposures Stratification and time-varying exposures in Cox regression Straightforward to implement in R haplo.ccs for binary outcomes haplo.cph for censored survival outcomes See [French et al., 2006; 2012] for implementation in Stata B French (bcfrench@upenn.edu) Estimating haplotype associations September / 31
34 Target High Mostly known May not exist Relative risk Low Probably can t be found Our focus Rare Common Haplotype frequency B French (bcfrench@upenn.edu) Estimating haplotype associations September / 31
35 Limitations Our methods and/or software may not be applicable to Related individuals Rare haplotypes Small studies Longitudinal outcomes B French (bcfrench@upenn.edu) Estimating haplotype associations September / 31
36 References dbe.med.upenn.edu/biostat-research/bcfrench 1. Cappola TP, Matkovich SJ, et al Loss-of-function DNA sequence variant in the CLCNKA chloride channel implicates the cardio-renal axis in interindividual heart failure risk variation. Proc Natl Acad Sci U S A 108: French B, Lumley T haplo.ccs: Estimate haplotype relative risks in case-control data. R package French B, Lumley T, et al Simple estimates of haplotype relative risks in case-control data. Genet Epidemiol 30: French B, Lumley T, et al Non-iterative, regression-based estimation of haplotype associations with censored survival outcomes. Stat App Genet Mol Biol 11:4. 5. Rebbeck TR, Mitra N, et al Modification of ovarian cancer risk by BRAC1/2-interacting genes in a multi center cohort of BRAC1/2 mutation carriers. Cancer Res 71: Sinnwell JP, Schaid DJ haplo.stats: Statistical analysis of haplotypes with traits and covariates when linkage phase is ambiguous. R package version Stephens M, Donnelly P A comparison of Bayesian methods for haplotype reconstruction from population genotype data. Am J Hum Genet 73: B French (bcfrench@upenn.edu) Estimating haplotype associations September / 31
New Educators Campaign Weekly Report
Campaign Weekly Report Conversations and 9/24/2017 Leader Forms Emails Collected Text Opt-ins Digital Journey 14,661 5,289 4,458 7,124 317 13,699 1,871 2,124 Pro 13,924 5,175 4,345 6,726 294 13,086 1,767
More informationA. Geography Students know the location of places, geographic features, and patterns of the environment.
Learning Targets Elementary Social Studies Grade 5 2014-2015 A. Geography Students know the location of places, geographic features, and patterns of the environment. A.5.1. A.5.2. A.5.3. A.5.4. Label North
More informationStandard Indicator That s the Latitude! Students will use latitude and longitude to locate places in Indiana and other parts of the world.
Standard Indicator 4.3.1 That s the Latitude! Purpose Students will use latitude and longitude to locate places in Indiana and other parts of the world. Materials For the teacher: graph paper, globe showing
More informationJakarta International School 6 th Grade Formative Assessment Graphing and Statistics -Black
Jakarta International School 6 th Grade Formative Assessment Graphing and Statistics -Black Name: Date: Score : 42 Data collection, presentation and application Frequency tables. (Answer question 1 on
More informationChallenge 1: Learning About the Physical Geography of Canada and the United States
60ºN S T U D E N T H A N D O U T Challenge 1: Learning About the Physical Geography of Canada and the United States 170ºE 10ºW 180º 20ºW 60ºN 30ºW 1 40ºW 160ºW 50ºW 150ºW 60ºW 140ºW N W S E 0 500 1,000
More informationMultivariate Statistics
Multivariate Statistics Chapter 6: Cluster Analysis Pedro Galeano Departamento de Estadística Universidad Carlos III de Madrid pedro.galeano@uc3m.es Course 2017/2018 Master in Mathematical Engineering
More informationMultivariate Statistics
Multivariate Statistics Chapter 4: Factor analysis Pedro Galeano Departamento de Estadística Universidad Carlos III de Madrid pedro.galeano@uc3m.es Course 2017/2018 Master in Mathematical Engineering Pedro
More informationAdditional VEX Worlds 2019 Spot Allocations
Overview VEX Worlds 2019 Spot s Qualifying spots for the VEX Robotics World Championship are calculated twice per year. On the following table, the number in the column is based on the number of teams
More informationPreview: Making a Mental Map of the Region
Preview: Making a Mental Map of the Region Draw an outline map of Canada and the United States on the next page or on a separate sheet of paper. Add a compass rose to your map, showing where north, south,
More informationCorrection to Spatial and temporal distributions of U.S. winds and wind power at 80 m derived from measurements
JOURNAL OF GEOPHYSICAL RESEARCH, VOL. 109,, doi:10.1029/2004jd005099, 2004 Correction to Spatial and temporal distributions of U.S. winds and wind power at 80 m derived from measurements Cristina L. Archer
More informationCooperative Program Allocation Budget Receipts Southern Baptist Convention Executive Committee May 2018
Cooperative Program Allocation Budget Receipts May 2018 Cooperative Program Allocation Budget Current Current $ Change % Change Month Month from from Contribution Sources 2017-2018 2016-2017 Prior Year
More informationCooperative Program Allocation Budget Receipts Southern Baptist Convention Executive Committee October 2017
Cooperative Program Allocation Budget Receipts October 2017 Cooperative Program Allocation Budget Current Current $ Change % Change Month Month from from Contribution Sources 2017-2018 2016-2017 Prior
More informationCooperative Program Allocation Budget Receipts Southern Baptist Convention Executive Committee October 2018
Cooperative Program Allocation Budget Receipts October 2018 Cooperative Program Allocation Budget Current Current $ Change % Change Month Month from from Contribution Sources 2018-2019 2017-2018 Prior
More information, District of Columbia
State Capitals These are the State Seals of each state. Fill in the blank with the name of each states capital city. (Hint: You may find it helpful to do the word search first to refresh your memory.),
More informationMultivariate Analysis
Multivariate Analysis Chapter 5: Cluster analysis Pedro Galeano Departamento de Estadística Universidad Carlos III de Madrid pedro.galeano@uc3m.es Course 2015/2016 Master in Business Administration and
More informationMultivariate Statistics
Multivariate Statistics Chapter 3: Principal Component Analysis Pedro Galeano Departamento de Estadística Universidad Carlos III de Madrid pedro.galeano@uc3m.es Course 2017/2018 Master in Mathematical
More informationHourly Precipitation Data Documentation (text and csv version) February 2016
I. Description Hourly Precipitation Data Documentation (text and csv version) February 2016 Hourly Precipitation Data (labeled Precipitation Hourly in Climate Data Online system) is a database that gives
More informationQF (Build 1010) Widget Publishing, Inc Page: 1 Batch: 98 Test Mode VAC Publisher's Statement 03/15/16, 10:20:02 Circulation by Issue
QF 1.100 (Build 1010) Widget Publishing, Inc Page: 1 Circulation by Issue Qualified Non-Paid Circulation Qualified Paid Circulation Individual Assoc. Total Assoc. Total Total Requester Group Qualified
More informationSummary of Natural Hazard Statistics for 2008 in the United States
Summary of Natural Hazard Statistics for 2008 in the United States This National Weather Service (NWS) report summarizes fatalities, injuries and damages caused by severe weather in 2008. The NWS Office
More informationAbortion Facilities Target College Students
Target College Students By Kristan Hawkins Executive Director, Students for Life America Ashleigh Weaver Researcher Abstract In the Fall 2011, Life Dynamics released a study entitled, Racial Targeting
More informationOnline Appendix: Can Easing Concealed Carry Deter Crime?
Online Appendix: Can Easing Concealed Carry Deter Crime? David Fortunato University of California, Merced dfortunato@ucmerced.edu Regulations included in institutional context measure As noted in the main
More informationClub Convergence and Clustering of U.S. State-Level CO 2 Emissions
Methodological Club Convergence and Clustering of U.S. State-Level CO 2 Emissions J. Wesley Burnett Division of Resource Management West Virginia University Wednesday, August 31, 2013 Outline Motivation
More informationNorth American Geography. Lesson 2: My Country tis of Thee
North American Geography Lesson 2: My Country tis of Thee Unit Overview: As students work through the activities in this unit they will be introduced to the United States in general, different regions
More informationWhat Lies Beneath: A Sub- National Look at Okun s Law for the United States.
What Lies Beneath: A Sub- National Look at Okun s Law for the United States. Nathalie Gonzalez Prieto International Monetary Fund Global Labor Markets Workshop Paris, September 1-2, 2016 What the paper
More informationParametric Test. Multiple Linear Regression Spatial Application I: State Homicide Rates Equations taken from Zar, 1984.
Multiple Linear Regression Spatial Application I: State Homicide Rates Equations taken from Zar, 984. y ˆ = a + b x + b 2 x 2K + b n x n where n is the number of variables Example: In an earlier bivariate
More informationPrintable Activity book
Printable Activity book 16 Pages of Activities Printable Activity Book Print it Take it Keep them busy Print them out Laminate them or Put them in page protectors Put them in a binder Bring along a dry
More informationJAN/FEB MAR/APR MAY/JUN
QF 1.100 (Build 1010) Widget Publishing, Inc Page: 1 Circulation Breakdown by Issue Qualified Non-Paid Qualified Paid Previous This Previous This Total Total issue Removals Additions issue issue Removals
More informationMeteorology 110. Lab 1. Geography and Map Skills
Meteorology 110 Name Lab 1 Geography and Map Skills 1. Geography Weather involves maps. There s no getting around it. You must know where places are so when they are mentioned in the course it won t be
More information2005 Mortgage Broker Regulation Matrix
2005 Mortgage Broker Regulation Matrix Notes on individual states follow the table REG EXEMPTIONS LIC-EDU LIC-EXP LIC-EXAM LIC-CONT-EDU NET WORTH BOND MAN-LIC MAN-EDU MAN-EXP MAN-EXAM Alabama 1 0 2 0 0
More informationOffice of Special Education Projects State Contacts List - Part B and Part C
Office of Special Education Projects State Contacts List - Part B and Part C Source: http://www.ed.gov/policy/speced/guid/idea/monitor/state-contactlist.html Alabama Customer Specialist: Jill Harris 202-245-7372
More informationRELATIONSHIPS BETWEEN THE AMERICAN BROWN BEAR POPULATION AND THE BIGFOOT PHENOMENON
RELATIONSHIPS BETWEEN THE AMERICAN BROWN BEAR POPULATION AND THE BIGFOOT PHENOMENON ETHAN A. BLIGHT Blight Investigations, Gainesville, FL ABSTRACT Misidentification of the American brown bear (Ursus arctos,
More informationOsteopathic Medical Colleges
Osteopathic Medical Colleges Matriculants by U.S. States and Territories Entering Class 0 Prepared by the Research Department American Association of Colleges of Osteopathic Medicine Copyright 0, AAM All
More informationLecture 5: Ecological distance metrics; Principal Coordinates Analysis. Univariate testing vs. community analysis
Lecture 5: Ecological distance metrics; Principal Coordinates Analysis Univariate testing vs. community analysis Univariate testing deals with hypotheses concerning individual taxa Is this taxon differentially
More informationDOWNLOAD OR READ : USA PLANNING MAP PDF EBOOK EPUB MOBI
DOWNLOAD OR READ : USA PLANNING MAP PDF EBOOK EPUB MOBI Page 1 Page 2 usa planning map usa planning map pdf usa planning map Printable USA Blank Map, USA Blank Map PDF, Blank US State Map. Thursday, 19
More informationCrop Progress. Corn Mature Selected States [These 18 States planted 92% of the 2017 corn acreage]
Crop Progress ISSN: 00 Released October, 0, by the National Agricultural Statistics Service (NASS), Agricultural Statistics Board, United s Department of Agriculture (USDA). Corn Mature Selected s [These
More informationIntercity Bus Stop Analysis
by Karalyn Clouser, Research Associate and David Kack, Director of the Small Urban and Rural Livability Center Western Transportation Institute College of Engineering Montana State University Report prepared
More informationChapter. Organizing and Summarizing Data. Copyright 2013, 2010 and 2007 Pearson Education, Inc.
Chapter 2 Organizing and Summarizing Data Section 2.1 Organizing Qualitative Data Objectives 1. Organize Qualitative Data in Tables 2. Construct Bar Graphs 3. Construct Pie Charts When data is collected
More informationLecture 5: Ecological distance metrics; Principal Coordinates Analysis. Univariate testing vs. community analysis
Lecture 5: Ecological distance metrics; Principal Coordinates Analysis Univariate testing vs. community analysis Univariate testing deals with hypotheses concerning individual taxa Is this taxon differentially
More informationGrand Total Baccalaureate Post-Baccalaureate Masters Doctorate Professional Post-Professional
s by Location of Permanent Home Address and Degree Level Louisiana Acadia 19 13 0 3 0 3 0 0 0 Allen 5 5 0 0 0 0 0 0 0 Ascension 307 269 2 28 1 6 0 1 0 Assumption 14 12 0 1 0 1 0 0 0 Avoyelles 6 4 0 1 0
More informationOUT-OF-STATE 965 SUBTOTAL OUT-OF-STATE U.S. TERRITORIES FOREIGN COUNTRIES UNKNOWN GRAND TOTAL
Report ID: USSR8072-V3 Page No. 1 Jurisdiction: ON-CAMPUS IL Southern Illinois University - Carb 1 0 0 0 Black Hawk College Quad-Cities 0 0 1 0 John A Logan College 1 0 0 0 Rend Lake College 1 0 0 0 Aurora
More informationSUPPLEMENTAL NUTRITION ASSISTANCE PROGRAM QUALITY CONTROL ANNUAL REPORT FISCAL YEAR 2008
SUPPLEMENTAL NUTRITION ASSISTANCE PROGRAM QUALITY CONTROL ANNUAL REPORT FISCAL YEAR 2008 U.S. DEPARTMENT OF AGRICULTURE FOOD AND NUTRITION SERVICE PROGRAM ACCOUNTABILITY AND ADMINISTRATION DIVISION QUALITY
More informationMultivariate Classification Methods: The Prevalence of Sexually Transmitted Diseases
Multivariate Classification Methods: The Prevalence of Sexually Transmitted Diseases Summer Undergraduate Mathematical Sciences Research Institute (SUMSRI) Lindsay Kellam, Queens College kellaml@queens.edu
More informationGrand Total Baccalaureate Post-Baccalaureate Masters Doctorate Professional Post-Professional
s by Location of Permanent Home Address and Degree Level Louisiana Acadia 26 19 0 6 1 0 0 0 0 Allen 7 7 0 0 0 0 0 0 0 Ascension 275 241 3 23 1 6 0 1 0 Assumption 13 12 0 1 0 0 0 0 0 Avoyelles 15 11 0 3
More informationLABORATORY REPORT. If you have any questions concerning this report, please do not hesitate to call us at (800) or (574)
LABORATORY REPORT If you have any questions concerning this report, please do not hesitate to call us at (800) 332-4345 or (574) 233-4777. This report may not be reproduced, except in full, without written
More informationLABORATORY REPORT. If you have any questions concerning this report, please do not hesitate to call us at (800) or (574)
LABORATORY REPORT If you have any questions concerning this report, please do not hesitate to call us at (800) 332-4345 or (574) 233-4777. This report may not be reproduced, except in full, without written
More informationHigh School World History Cycle 2 Week 2 Lifework
Name: Advisory: Period: High School World History Cycle 2 Week 2 Lifework This packet is due Monday, November 7 Complete and turn in on Friday for 10 points of EXTRA CREDIT! Lifework Assignment Complete
More informationOffice of Budget & Planning 311 Thomas Boyd Hall Baton Rouge, LA Telephone 225/ Fax 225/
Louisiana Acadia 20 17 3 0 0 0 Allen 2 2 0 0 0 0 Ascension 226 185 37 2 1 1 Assumption 16 15 1 0 0 0 Avoyelles 20 19 1 0 0 0 Beauregard 16 11 4 0 0 1 Bienville 2 2 0 0 0 0 Bossier 22 18 4 0 0 0 Caddo 91
More informationAlpine Funds 2017 Tax Guide
Alpine s 2017 Guide Alpine Dynamic Dividend ADVDX 1/30/17 1/31/17 1/31/17 0.020000000 0.019248130 0.000000000 0.00000000 0.019248130 0.013842273 0.000000000 0.000000000 0.000751870 0.000000000 0.00 0.00
More informationGIS use in Public Health 1
Geographic Information Systems (GIS) use in Public Health Douglas Morales, MPH Epidemiologist/GIS Coordinator Office of Health Assessment and Epidemiology Epidemiology Unit Objectives Define GIS and justify
More informationAlpine Funds 2016 Tax Guide
Alpine s 2016 Guide Alpine Dynamic Dividend ADVDX 01/28/2016 01/29/2016 01/29/2016 0.020000000 0.017621842 0.000000000 0.00000000 0.017621842 0.013359130 0.000000000 0.000000000 0.002378158 0.000000000
More informationJAN/FEB MAR/APR MAY/JUN
QF 1.100 (Build 1010) Widget Publishing, Inc Page: 1 Circulation Breakdown by Issue Analyzed Nonpaid and Verified Paid Previous This Previous This Total Total issue Removals Additions issue issue Removals
More informationextreme weather, climate & preparedness in the american mind
extreme weather, climate & preparedness in the american mind Extreme Weather, Climate & Preparedness In the American Mind Interview dates: March 12, 2012 March 30, 2012. Interviews: 1,008 Adults (18+)
More informationKS PUBL 4YR Kansas State University Pittsburg State University SUBTOTAL-KS
Report ID: USSR8072-V3 Page No. 1 Jurisdiction: ON-CAMPUS IL PUBL TCH DeVry University Addison 1 0 0 0 Eastern Illinois University 1 0 0 0 Illinois State University 0 0 2 0 Northern Illinois University
More informationLast time: PCA. Statistical Data Mining and Machine Learning Hilary Term Singular Value Decomposition (SVD) Eigendecomposition and PCA
Last time: PCA Statistical Data Mining and Machine Learning Hilary Term 2016 Dino Sejdinovic Department of Statistics Oxford Slides and other materials available at: http://www.stats.ox.ac.uk/~sejdinov/sdmml
More informationPackage ZIM. R topics documented: August 29, Type Package. Title Statistical Models for Count Time Series with Excess Zeros. Version 1.
Package ZIM August 29, 2013 Type Package Title Statistical Models for Count Time Series with Excess Zeros Version 1.0 Date 2013-06-15 Author Ming Yang, Gideon K. D. Zamba, and Joseph E. Cavanaugh Maintainer
More informationNational Organization of Life and Health Insurance Guaranty Associations
National Organization of and Health Insurance Guaranty Associations November 21, 2005 Dear Chief Executive Officer: Consistent with prior years, NOLHGA is providing the enclosed data regarding insolvency
More informationMINERALS THROUGH GEOGRAPHY
MINERALS THROUGH GEOGRAPHY INTRODUCTION Minerals are related to rock type, not political definition of place. So, the minerals are to be found in a variety of locations that doesn t depend on population
More informationA GUIDE TO THE CARTOGRAPHIC PRODUCTS OF
A GUIDE TO THE CARTOGRAPHIC PRODUCTS OF THE FEDERAL DEPOSITORY LIBRARY PROGRAM (FDLP) This guide was designed for use as a collection development tool by map selectors of depository libraries that participate
More informationBlackRock Core Bond Trust (BHK) BlackRock Enhanced International Dividend Trust (BGY) 2 BlackRock Defined Opportunity Credit Trust (BHL) 3
MUNICIPAL FUNDS Arizona (MZA) California Municipal Income Trust (BFZ) California Municipal 08 Term Trust (BJZ) California Quality (MCA) California Quality (MUC) California (MYC) Florida Municipal 00 Term
More informationIntroducing North America
Introducing North America I. Quick Stats Includes U.S. & Canada U.S consists of 50 States Federal Government Democracy 4 th in world w/ land area 3 rd in population Economic leader of free world II. Major
More informationMO PUBL 4YR 2090 Missouri State University SUBTOTAL-MO
Report ID: USSR8072-V3 Page No. 1 Jurisdiction: ON-CAMPUS IL American Intercontinental Universit 0 0 1 0 Northern Illinois University 0 0 4 0 Southern Illinois Univ - Edwardsvil 2 0 2 0 Southern Illinois
More informationInsurance Department Resources Report Volume 1
2014 Insurance Department Resources Report Volume 1 201 Insurance Department Resources Report Volume One 201 The NAIC is the authoritative source for insurance industry information. Our expert solutions
More informationInfant Mortality: Cross Section study of the United State, with Emphasis on Education
Illinois State University ISU ReD: Research and edata Stevenson Center for Community and Economic Development Arts and Sciences Fall 12-15-2014 Infant Mortality: Cross Section study of the United State,
More informationAn Analysis of Regional Income Variation in the United States:
Modern Economy, 2017, 8, 232-248 http://www.scirp.org/journal/me ISSN Online: 2152-7261 ISSN Print: 2152-7245 An Analysis of Regional Income Variation in the United States: 1969-2013 Orley M. Amos Jr.
More informationFLOOD/FLASH FLOOD. Lightning. Tornado
2004 Annual Summaries National Oceanic and Atmospheric Administration National Environmental Satellite Data Information Service National Climatic Data Center FLOOD/FLASH FLOOD Lightning Tornado Hurricane
More informationRank University AMJ AMR ASQ JAP OBHDP OS PPSYCH SMJ SUM 1 University of Pennsylvania (T) Michigan State University
Rank University AMJ AMR ASQ JAP OBHDP OS PPSYCH SMJ SUM 1 University of Pennsylvania 4 1 2 0 2 4 0 9 22 2(T) Michigan State University 2 0 0 9 1 0 0 4 16 University of Michigan 3 0 2 5 2 0 0 4 16 4 Harvard
More informationPima Community College Students who Enrolled at Top 200 Ranked Universities
Pima Community College Students who Enrolled at Top 200 Ranked Universities Institutional Research, Planning and Effectiveness Project #20170814-MH-60-CIR August 2017 Students who Attended Pima Community
More informationUnited States Geography Unit 1
United States Geography Unit 1 I WANT YOU TO STUDY YOUR GEORGAPHY Name: Period: Due Date: Geography Key Terms Absolute Location: Relative Location: Demographic Map: Population Density: Sun-Belt: Archipelago:
More information112th U.S. Senate ACEC Scorecard
112th U.S. Senate ACEC Scorecard HR 658 FAA Funding Alaska Alabama Arkansas Arizona California Colorado Connecticut S1 Lisa Murkowski R 100% Y Y Y Y Y Y Y Y Y S2 Mark Begich D 89% Y Y Y Y N Y Y Y Y S1
More informationA Summary of State DOT GIS Activities
A Summary of State DOT GIS Activities Prepared for the 2006 AASHTO GIS-T Symposium Columbus, OH Introduction This is the 11 th year that the GIS-T Symposium has conducted a survey of GIS activities at
More informationMINERALS THROUGH GEOGRAPHY. General Standard. Grade level K , resources, and environmen t
Minerals through Geography 1 STANDARDS MINERALS THROUGH GEOGRAPHY See summary of National Science Education s. Original: http://books.nap.edu/readingroom/books/nses/ Concept General Specific General Specific
More informationOutline. Administrivia and Introduction Course Structure Syllabus Introduction to Data Mining
Outline Administrivia and Introduction Course Structure Syllabus Introduction to Data Mining Dimensionality Reduction Introduction Principal Components Analysis Singular Value Decomposition Multidimensional
More informationAll-Time Conference Standings
All-Time Conference Standings Pac 12 Conference Conf. Matches Sets Overall Matches Team W L Pct W L Pct. Score Opp Last 10 Streak Home Away Neutral W L Pct. Arizona 6 5.545 22 19.537 886 889 6-4 W5 4-2
More informationCCC-A Survey Summary Report: Number and Type of Responses
CCC-A Survey Summary Report: Number and Type of Responses Suggested Citation: American Speech-Language-Hearing Association. (2011). 2011 Membership survey. CCC-A survey summary report: Number and type
More informationOffice of Budget & Planning 311 Thomas Boyd Hall Baton Rouge, LA Telephone 225/ Fax 225/
Louisiana Acadia 25 19 4 2 0 0 Allen 8 7 1 0 0 0 Ascension 173 143 26 1 0 3 Assumption 14 12 2 0 0 0 Avoyelles 51 41 9 0 0 1 Beauregard 18 14 3 0 0 1 Bienville 5 0 4 0 1 0 Bossier 28 27 0 1 0 0 Caddo 95
More informationOctober 2016 v1 12/10/2015 Page 1 of 10
State Section S s Effective October 1, 2016 Overview The tables list the Section S items that will be active on records with a target date on or after October 1, 2016. The active on each item subset code
More informationThis project is supported by a National Crime Victims' Right Week Community Awareness Project subgrant awarded by the National Association of VOCA
TABLE OF CONTENTS ACTIVITY PAGES Coloring Pages Younger Children (6) Coloring Pages Older (2) Connect the Dots (1 English, 1 Spanish) Word Search 1 (1 English, 1 Spanish) Younger Children Word Search 2
More information2011 NATIONAL SURVEY ON DRUG USE AND HEALTH
2011 NATIONAL SURVEY ON DRUG USE AND HEALTH PERSON-LEVEL SAMPLING WEIGHT CALIBRATION Prepared for the 2011 Methodological Resource Book Contract No. HHSS283200800004C RTI Project No. 0211838.207.004 Phase
More informationA Standardized Digital Well-Record Database for the Glaciated U.S.
A Standardized Digital Well-Record Database for the Glaciated U.S. Randy Bayless and Les Arihood U.S. Geological Survey, Indiana Water Science Center, Indianapolis, IN Howard Reeves U.S. Geological Survey,
More informationPhysical Features of Canada and the United States
Physical Features of Canada and the United tates 170 ARCTIC OCA Aleutian s 1 1 Bering ea ALAKA Yukon R. Mt. McKinley (20,320 ft. 6,194 m) Gulf of Alaska BROOK RAG RAG Queen Charlotte s R Vancouver O C
More informationSt. Luke s College Institutional Snapshot
1. Student Headcount St. Luke s College Institutional Snapshot Enrollment # % # % # % Full Time 139 76% 165 65% 132 54% Part Time 45 24% 89 35% 112 46% Total 184 254 244 2. Average Age Average Age Average
More informationPhysical Features of Canada and the United States
I VIUAL Physical Features of Canada and the United tates 170 ARCTIC OCA Aleutian s 1 Bering ea ALAKA Yukon R. Mt. McKinley (20,320 ft. 6,194 m) Gulf of Alaska BROOK RAG RAG Queen Charlotte s R Vancouver
More informationDirect Selling Association 1
Total Number of Firms 1,100 104 58 46 14 41 22 27 44 60 1 U.S. Sales Volume & Growth 1 Average Retail Sales per Firm ($thousands) 2015 $32,836 $210,572 $320,820 $69,556 $1,041 $11,788 $100,216 $721,897
More informationEarl E. ~rabhl/ This report is preliminary and has not been reviewed for conformity with U.S. Geological Survey editoral standards GEOLOGICAL SURVEY
UNITED STATES DEPARTMENT OF THE INTERIOR GEOLOGICAL SURVEY Minimum landslide damage in the United States, 1973-1983 by Earl E. ~rabhl/ Open-File Report 84-486 This report is preliminary and has not been
More informationChapter 11 : State SAT scores for 1982 Data Listing
EXST3201 Chapter 12a Geaghan Fall 2005: Page 1 Chapter 12 : Variable selection An example: State SAT scores In 1982 there was concern for scores of the Scholastic Aptitude Test (SAT) scores that varied
More informationStem-and-Leaf Displays
3.2 Displaying Numerical Data: Stem-and-Leaf Displays 107 casts in your area? (San Luis Obispo Tribune, June 15, 2005). The responses are summarized in the table below. Extremely 4% Very 27% Somewhat 53%
More informationChapter 5 Linear least squares regression
Chapter 5 Linear least squares regression Consider first simple linear regression. For now, the model of interest is only the model of the conditional expectations; the distribution of the residuals is
More informationEvidence for increasingly extreme and variable drought conditions in the contiguous United States between 1895 and 2012
Evidence for increasingly extreme and variable drought conditions in the contiguous United States between 1895 and 2012 Sierra Rayne a,, Kaya Forest b a Chemologica Research, 318 Rose Street, PO Box 74,
More informationDivision I Sears Directors' Cup Final Standings As of 6/20/2001
1 Stanford (Cal.) 1359 2 90 9 75 20 62.5 0 0 0 0 0 0 3 75 1 100 5 60 8 76 4 80 0 0 2 70 2 California-Los Angeles 1138 0 0 5 78.5 17 67 0 0 0 0 0 0 2 90 9 50 5 60 2 0 33 48.5 2 70 1 100 3 Georgia 890.5
More informationState Section S Items Effective October 1, 2017
State Section S Items Effective October 1, 2017 Overview The tables list the Section S items that will be active on records with a target date on or after October 1, 2017. The active item on each item
More informationThe Heterogeneous Effects of the Minimum Wage on Employment Across States
The Heterogeneous Effects of the Minimum Wage on Employment Across States Wuyi Wang a, Peter C.B. Phillips b, Liangjun Su c a Institute for Economic and Social Research, Jinan University b Yale University,
More informationSNP Association Studies with Case-Parent Trios
SNP Association Studies with Case-Parent Trios Department of Biostatistics Johns Hopkins Bloomberg School of Public Health September 3, 2009 Population-based Association Studies Balding (2006). Nature
More informationFungal conservation in the USA
The following supplements accompany the article Fungal conservation in the USA Jessica L. Allen*, James C. Lendemer *Corresponding author: jlendemer@nybg.org Endangered Species Research 28: 33 42 (2015)
More informationConduent EDI Solutions, Inc. EDI Direct Pass Thru Fees
Conduent EDI Solutions, Inc. EDI Direct Pass Thru Fees Eligibility/Claim Status Payer $ Eligibility/Claim Status Payer $ Aetna Better Health (LA) 0.04 257 BCBS of Western New York (Healthnow) 0.07 255
More information(Specification B) (JUN H01) (JaN11GEOG101) General Certificate of Education Secondary Education Advanced Higher TierSubsidiary Examination
Centre Number Surname Candidate Number For Examiner s Use Other Names Candidate Signature Examiner s Initials Geography (Specification B) Unit 1 General Certificate of Education Secondary Education Advanced
More informationPrerequisite: STATS 7 or STATS 8 or AP90 or (STATS 120A and STATS 120B and STATS 120C). AP90 with a minimum score of 3
University of California, Irvine 2017-2018 1 Statistics (STATS) Courses STATS 5. Seminar in Data Science. 1 Unit. An introduction to the field of Data Science; intended for entering freshman and transfers.
More informationCrop / Weather Update
Crop / Weather Update Corn Crop Condition Percent of Acreage Rated Good or Excellent 85 80 75 70 65 60 55 50 45 As of September 9, USDA rates the crop at 68% good to excellent. The rating is up one point
More informationHP-35s Calculator Program Lambert 1
Convert Latitude and Longitude to Lambert Conformal Conic Projection Co-ordinates (SPCS) Programmer: Dr. Bill Hazelton Date: July, 2010. Version: 1.2 Line Instruction Display User Instructions N001 LBL
More informationGreen Building Criteria in Low-Income Housing Tax Credit Programs Analysis
Green Building Criteria in Low-Income Housing Tax Credit Programs 2010 Analysis www.globalgreen.org September 2010 Green Building Criteria in State Low-Income Housing Tax Credit Programs 2 Introduction
More informationPhysical Therapist Centralized Application Service. COURSE PREREQUISITES SUMMARY Admissions Cycle
Physical Therapist Centralized Application Service COURSE PREREQUISITES SUMMARY Report of the American Physical Therapy Association August 2011 2011-12 COURSE PREREQUISITES SUMMARY EXECUTIVE SUMMARY The
More information