Since my amateur discursion into stochasticity appears to have flushed out all of the mathematically savvy, I’m going pose a real life statistics question for you. I have some data that are non-normally distributed (in fact, they really don’t seem to fit any distribution well–and, yes, I’ve tried various transformations and the data still don’t fit anything). If the data were normally distributed, I would like to perform ANOVA to partition the sources of variance (using percent sums of squares).
Since the data aren’t normally distributed, ANOVA is not the right test to use. Because I’m using a full-factorial model with four factors, using Friedman’s test or any of the one-way non-parametric tests doesn’t seem to apply either. I should add that, while my data are continuous, I can transform the data into ordinal and/or nominal categories if that would help. While logistic regression could circumvent the non-normality problem, I’m unaware of how one partitions variance with that test (if there’s a way to do so, please let me know).
Any ideas? If there are any statistics programs or R packages out there that can handle this, please let me know. Don’t be shy…
-
Recent Posts
Top Posts & Pages
- The Four (Anti-) Freedoms
- Links 5/26/23
- You're Not the Weirdo If You're Still Taking Precautions Against (Long) COVID
- Regarding the Debt Ceilling, This Time Biden Really Is The Green Lantern
- New Democrats Are So Exhausting
- In Case You Missed It...
- Links 5/27/23
- Some Thoughts on Biden and Aging
- Professional Democrats Are Bad at Their Jobs: The 14th Amendment Edition
- The CDC EIS Meeting Was Not a 'Superspreader Event', It Was the New Normal
Meta
Mike the Mad Biologist
Email Subscription
Join 4,242 other subscribersCategories
- 'Snuggles' McCain
- 3-D Printing
- A Fungus Among Us
- A Message to You
- A Message to You, Rudy
- A Moose Is Loose
- Acinetobacter
- Afghanistan
- Agriculture
- Alba Gu Brath
- Algae
- Aliens! I Seen 'Em!
- ALL UR SAMPLEZ R BELONGING TO ME!
- Ammosexuals
- Amoeba
- Animalcules
- Anthrax
- Antibiotics
- Archery
- Architecture
- Art Deco
- Art History
- Arthropods
- Ascidians
- Astronomy
- Authoritarianism
- Autism
- Automobiles
- Aves
- Axolotl
- Baby
- Balimore!
- BANG! BANG!
- Basic Human Decency
- Beer
- Behavior
- Bidness
- Big Shitpile
- Bioinformatics
- Biotech
- Blastocyst Liberation
- Bloggity Blog
- Books
- Booze
- Boston
- BPSDB
- BRAINS! I EAT THEM!!
- Bullshit As a Load Bearing Structure
- Bureaucracy
- Burma
- Campylobacter
- Camus
- Cancer
- Career Stuff
- Carnivora
- CDC
- Census
- CIA Spy Shit
- Civil Liberties
- Civility Trolls
- Cloning
- Coelenterata
- Coins
- Compulsive Centrist Disorder
- Compulsive Contrarian Disorder
- Computers
- Condolences
- Conservatives
- COVID-19
- COVID-19 Benchmark
- Creationism
- Crime
- Crime Pays?
- Danger Ain't Cool!
- Data Stuff
- Daylight Saving Time
- DC
- Death
- Democrats
- Demonic Possession
- Dinosaurs
- Dogs
- Drugs
- E. coli
- Earthquakes
- Easter Island
- Ecology
- Economics
- Edge Cases
- Education
- Elephants
- Energy
- Engineering
- Environment
- Ethics
- Evo Devo
- Evo Psych
- Evolution
- Experts
- Feminism
- Fish
- Food
- Football
- Framing
- Froggies
- Fucking Morons
- Fucking Royals
- Funding
- Gambling
- Gastropods
- Genetics
- Genomics
- Get a Pair, Will Ya?
- Global Warming
- Gravity, How Does It Work?
- Group Selection
- Gut Yontif
- Healthcare
- HIV
- Homelessness
- Housing
- Human Genetics
- Humor
- Hunger
- Hypotheses
- I Am Aware of All Internet Traditions
- I For One Welcome Our Arachnoid Overlords
- I For One Welcome Our Austerity Overlords
- I For One Welcome Our Cephalopod Overlords
- I For One Welcome Our Silicon Valley Overlords
- I For One Welcome Our Socialist Overlords
- Ick
- Immigration
- Inflation
- Influenza
- Insects
- Internet
- iPhone, uPhone, We All Like the iPhone
- IQ
- Iran
- Iraq
- Ireland
- Jewish Stuff
- Jobs
- John "Snuggles" McCain
- Jump Around Everybody Get Down
- Jury Duty
- Katrina
- Korea
- KPC
- Lady Bits
- Lakoff
- Legislation
- Legos
- Let's Talk About Sex
- Let's Talk About Sex, Baby
- LGBT
- Libraries
- Libya
- Likelihood
- Little Lord Pontchartrain
- Lobstah!
- Lotsa Links
- Lucy with the Football
- Machiavelli
- Makin' Babies
- Management
- Maps
- Massachusetts
- Mathematics
- Meetings
- Memes
- Microbiology
- Microbiome
- Middle East
- Military
- Model Systems
- Movies
- MRSA
- Museums etc.
- Music
- Mutation
- NASA
- NDM-1
- New
- New Orleans
- News Media
- NIH
- Ninnies
- NIST
- NSF
- Nukular Stuff
- Obesity
- Occupy Wall Street
- Oh Canada!
- OH NOES! MAH PENIS!
- Oil
- Olympics
- OOH! SHINY PEBBLE!
- Orangutan
- Paleontology
- Pee in a Cup!
- PhD School
- Philosophy Shit
- Pig Screwer
- Pizza!
- Placebo
- Plants
- Plasmids
- Please. Just Have One Chicken Finger
- Polling
- Poop
- Posters
- Progressives
- Propaganda
- Public Health
- Public Speaking
- published
- Publishing
- Racism
- Rape
- Reading
- Religion
- Reptiles
- Resistance
- Resistance Rebellion And Death
- Rhinos
- Robots
- Rodents
- Romney
- Ron Paul
- Rudy
- Run Ya Bastids!
- Salmonella
- Secession
- Secret Spy Shit
- Sept. 11
- Sexist Pigs
- Sharks
- Shut Up and Watch the Damn Movie!
- Sleep
- Social Security
- Speciation
- Sportsball!
- Statehood
- Statistics
- Stem Cells
- Stupid Internet Quiz
- Syria
- Tasmanian Devil
- Taxes
- Team Social Justice Warrior
- Teeth
- Terrorism
- The Freedom To Hunt Your Slaves
- The Identity That Shall Not Be Named
- The Rule of Law
- The Soviet
- The War of Southern Aggression
- The War on Art
- The War on Science
- There Be Dragons!
- Things That Make Ya Go Boom
- Tolerance
- Torture
- Transportation
- TWEET!!!
- Uncategorized
- Unions
- Urban Planning
- USPS
- Vaccination
- Vampires
- VIM
- Violence
- Viruses
- Vision
- Volcanoes How Do They Work?
- Vomit
- Voting
- We're Really Fucked
- Weekly Roundup
- WhatEVAH!
- Why Is Arianna Huffington a Murderous Idiot?
- Worms
- Ya Bastids!
- Year In Review
- Year of Jubilee
- Yersinia
- You're Having a Bad Hair Life
- You've Got Mail But Not For Long!
- Yurp
- Zombies
Archives
- May 2023
- April 2023
- March 2023
- February 2023
- January 2023
- December 2022
- November 2022
- October 2022
- September 2022
- August 2022
- July 2022
- June 2022
- May 2022
- April 2022
- March 2022
- February 2022
- January 2022
- December 2021
- November 2021
- October 2021
- September 2021
- August 2021
- July 2021
- June 2021
- May 2021
- April 2021
- March 2021
- February 2021
- January 2021
- December 2020
- November 2020
- October 2020
- September 2020
- August 2020
- July 2020
- June 2020
- May 2020
- April 2020
- March 2020
- February 2020
- January 2020
- December 2019
- November 2019
- October 2019
- September 2019
- August 2019
- July 2019
- June 2019
- May 2019
- April 2019
- March 2019
- February 2019
- January 2019
- December 2018
- November 2018
- October 2018
- September 2018
- August 2018
- July 2018
- June 2018
- May 2018
- April 2018
- March 2018
- February 2018
- January 2018
- December 2017
- November 2017
- October 2017
- September 2017
- August 2017
- July 2017
- June 2017
- May 2017
- April 2017
- March 2017
- February 2017
- January 2017
- December 2016
- November 2016
- October 2016
- September 2016
- August 2016
- July 2016
- June 2016
- May 2016
- April 2016
- March 2016
- February 2016
- January 2016
- December 2015
- November 2015
- October 2015
- September 2015
- August 2015
- July 2015
- June 2015
- May 2015
- April 2015
- March 2015
- February 2015
- January 2015
- December 2014
- November 2014
- October 2014
- September 2014
- August 2014
- July 2014
- June 2014
- May 2014
- April 2014
- March 2014
- February 2014
- January 2014
- December 2013
- November 2013
- October 2013
- September 2013
- August 2013
- July 2013
- June 2013
- May 2013
- April 2013
- March 2013
- February 2013
- January 2013
- December 2012
- November 2012
- October 2012
- September 2012
- August 2012
- July 2012
- June 2012
- May 2012
- April 2012
- March 2012
- February 2012
- January 2012
- December 2011
- November 2011
- October 2011
- September 2011
- August 2011
- July 2011
- June 2011
- May 2011
- April 2011
- March 2011
- February 2011
- January 2011
- December 2010
- November 2010
- October 2010
- September 2010
- August 2010
- July 2010
- June 2010
- May 2010
- April 2010
- March 2010
- February 2010
- January 2010
- December 2009
- November 2009
- October 2009
- September 2009
- August 2009
- July 2009
- June 2009
- May 2009
- April 2009
- March 2009
- February 2009
- January 2009
- December 2008
- November 2008
- October 2008
- September 2008
- August 2008
- July 2008
- June 2008
- May 2008
- April 2008
- March 2008
- February 2008
- January 2008
- December 2007
- November 2007
- October 2007
- September 2007
- August 2007
- July 2007
- June 2007
- May 2007
- April 2007
- March 2007
- February 2007
- January 2007
- December 2006
- November 2006
- October 2006
- September 2006
- August 2006
- July 2006
- June 2006
- February 2006
- September 2005
- August 2005
- March 2005
- January 2005
- December 2004
- November 2004
I remember vaguely that the output for logistic regression (in Stata, SAS, SPSS etc) has the amount of variation ‘explained’ by the model called ‘deviance’ The amount of variation due to a certain factor is hence the difference between deviance with the factor in the model and without. It is also used for the chi-square test for significance of the factor.
R-wise, maybe the examples in Frank Harrell’s Design::lrm package would help:
http://lib.stat.cmu.edu/S/Harrell/help/Design/html/lrm.html
With a logistic model the variance isn’t a good measure of fit because the domains of the errors is restricted. Maybe your data is nasty enough that the optimizing some function of the model errors doesn’t give you meaningful results, and any partitioning that error measure won’t tell you anything significant. You need to have some idea of the error distribution to make meaningful inferences about apportioning that error to factors.
I had a similar problem.
I’m using JMP (the full-on, expensive version).
JMP lets you decree that your continuous numbers are ordinal, without having to change the actual numbers. (I found it worked better when I truncated my numbers to 3 or 4 sig figs, though.) Then it lets you run a general linear model, giving you logistic regression output (instead of parametric). The GLM can handle ties in the ordinal data, multiple input variables, interactions, and even a little nesting.
It gives you the likelihood ratio tests for each effect — like effects tests in ANOVA — and provides parameter estimates (with P-values and 95%CIs) for each level in each effect. Ergo, the analysis does partition the error for you.
It also has the refs for the analyses it’s doing, so you can both figure out how the logistic GLM handles the error, and cite the authorities to back it up.
Finally, for those of us with advisors unfamiliar with logistic regression, you can run a (technically invalid) parametric ANOVA that is exactly congruent, just by switching your response variable back to “continuous.” When they give the same answer, this congruent-but-invalid analysis can be useful for educating advisors.
I can email you an example output if you’re interested.
PS
Mike, you’ve got two copies of the original post.
And they’ve got different comments.
Read both!