Corrected p-value during correlation analysis
2 次查看(过去 30 天)
显示 更早的评论
Hi everyone,
I am interested in doing multiple corrections while doing some correlation analysis.
I have brain thickness data of 50 subject's 14 brain networks: 7 networks on left hemisphere and 7 networks on right hemisphere i.e. X = 50x7x2. I would like to see whether thickness of each network of these 50 subjects is correlated with age i.e. my age variable is a vector of size 50 values i.e. Y = 50x1.
Could you please tell me how and which MATLAB program can I use to estimate 14 correlation coefficient values and corrected p-value of each of these correlations?
Thanks a lot.
0 个评论
回答(2 个)
Walter Roberson
2017-2-20
T1 = reshape(thicknesses, [], 14);
T2 = reshape(ages, [], 1);
T3 = [T1, T2];
c = cov(T3);
correlations_you_care_about = reshape( c(end,1:end-1), 7, 2);
Star Strider
2017-2-20
I would do a regression. See the ‘Statistical Analysis’ section in ‘Associations between age and gray matter volume in anatomical brain networks in middle-aged to older adults.’ Aging Cell. 2014 Dec;13(6):1068-74.
If you have the Statistics and Machine Learning Toolbox, you can do the sort of categorical regression they discuss, with multiple comparisons. The most obvious reason to adapt their statistical analysis procedure is that these authors have published in a refereed journal, and you would be using an accepted technique.
There may be other papers with statistical analysis techniques better suited to your experimental design. Search PubMed with the search string brain & volume & age for all 12276 similar articles.
I wish you success in your research.
另请参阅
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!