Doing regression only with correlation matrix, means, and SDs [duplicate]












1












$begingroup$



This question already has an answer here:




  • Is there a way to use the covariance matrix to find coefficients for multiple regression?

    1 answer




I was wondering how mathematically is it possible to run a full regression analysis between 3 predictors (x1 x2 x3) and a dependent variable (y) by only knowing the: Means, Ns, SDs, and the Correlations between all these 4 variables (without the original data)?



I highly appreciate an R demonstration.



ns    <- c(273, 273, 273, 273)

means <- c(15.4, 7.1, 3.5, 6.2)

sds   <- c(3.4, 0.9, 1.5, 1.4)



r <- matrix( c(

1.0,  .57,  -.4,  .48,  

.57,  1.0, -.61,  .66,   

-.4, -.61,  1.0, -.68,  

.48,  .66, -.68,   1.0), 4)



rownames(r) <- colnames(r) <- c('y', paste0('x', 1:3))





r regression multiple-regression regression-coefficients







share|cite|improve this question











$endgroup$



marked as duplicate by whuber♦ r
Users with the  r badge can single-handedly close r questions as duplicates and reopen them as needed.

StackExchange.ready(function() {
if (StackExchange.options.isMobile) return;

$('.dupe-hammer-message-hover:not(.hover-bound)').each(function() {
var $hover = $(this).addClass('hover-bound'),
$msg = $hover.siblings('.dupe-hammer-message');

$hover.hover(
function() {
$hover.showInfoMessage('', {
messageElement: $msg.clone().show(),
transient: false,
position: { my: 'bottom left', at: 'top center', offsetTop: -7 },
dismissable: false,
relativeToBody: true
});
},
function() {
StackExchange.helpers.removeMessages();
}
);
});
});
 Dec 5 '18 at 14:34


This question has been asked before and already has an answer. If those answers do not fully address your question, please ask a new question.














  • 1




    $begingroup$
    You cannot run a truly "full" regression analysis with just these statistics, because you will not be able to construct residuals and perform regression diagnostics that depend on them. You are limited to making and testing parameter estimates.
    $endgroup$
    – whuber♦
    Dec 5 '18 at 14:35










  • $begingroup$
    Structural equation modeling is the most obvious answer in my mind, which is more or less the math provided in the answer below. Most SEM programs will be able to accept your inputs (all need to be able to estimate a covariance matrix in the end).
    $endgroup$
    – Matt Barstead
    Dec 10 '18 at 2:29










  • $begingroup$
    @whuber, do I need to first convert my correlation matrix into a var-covariance matrix using r_x_iy_i * sd_x_i*sd_y_i and then go from there?
    $endgroup$
    – rnorouzian
    Dec 10 '18 at 4:52












  • $begingroup$
    That's one approach, because it reduces your problem to one with a known, explicit solution.
    $endgroup$
    – whuber♦
    Dec 10 '18 at 15:13
















1












$begingroup$



This question already has an answer here:




  • Is there a way to use the covariance matrix to find coefficients for multiple regression?

    1 answer




I was wondering how mathematically is it possible to run a full regression analysis between 3 predictors (x1 x2 x3) and a dependent variable (y) by only knowing the: Means, Ns, SDs, and the Correlations between all these 4 variables (without the original data)?



I highly appreciate an R demonstration.



ns    <- c(273, 273, 273, 273)

means <- c(15.4, 7.1, 3.5, 6.2)

sds   <- c(3.4, 0.9, 1.5, 1.4)



r <- matrix( c(

1.0,  .57,  -.4,  .48,  

.57,  1.0, -.61,  .66,   

-.4, -.61,  1.0, -.68,  

.48,  .66, -.68,   1.0), 4)



rownames(r) <- colnames(r) <- c('y', paste0('x', 1:3))





r regression multiple-regression regression-coefficients







share|cite|improve this question











$endgroup$



marked as duplicate by whuber♦ r
Users with the  r badge can single-handedly close r questions as duplicates and reopen them as needed.

StackExchange.ready(function() {
if (StackExchange.options.isMobile) return;

$('.dupe-hammer-message-hover:not(.hover-bound)').each(function() {
var $hover = $(this).addClass('hover-bound'),
$msg = $hover.siblings('.dupe-hammer-message');

$hover.hover(
function() {
$hover.showInfoMessage('', {
messageElement: $msg.clone().show(),
transient: false,
position: { my: 'bottom left', at: 'top center', offsetTop: -7 },
dismissable: false,
relativeToBody: true
});
},
function() {
StackExchange.helpers.removeMessages();
}
);
});
});
 Dec 5 '18 at 14:34


This question has been asked before and already has an answer. If those answers do not fully address your question, please ask a new question.














  • 1




    $begingroup$
    You cannot run a truly "full" regression analysis with just these statistics, because you will not be able to construct residuals and perform regression diagnostics that depend on them. You are limited to making and testing parameter estimates.
    $endgroup$
    – whuber♦
    Dec 5 '18 at 14:35










  • $begingroup$
    Structural equation modeling is the most obvious answer in my mind, which is more or less the math provided in the answer below. Most SEM programs will be able to accept your inputs (all need to be able to estimate a covariance matrix in the end).
    $endgroup$
    – Matt Barstead
    Dec 10 '18 at 2:29










  • $begingroup$
    @whuber, do I need to first convert my correlation matrix into a var-covariance matrix using r_x_iy_i * sd_x_i*sd_y_i and then go from there?
    $endgroup$
    – rnorouzian
    Dec 10 '18 at 4:52












  • $begingroup$
    That's one approach, because it reduces your problem to one with a known, explicit solution.
    $endgroup$
    – whuber♦
    Dec 10 '18 at 15:13














1












1








1


1



$begingroup$



This question already has an answer here:




  • Is there a way to use the covariance matrix to find coefficients for multiple regression?

    1 answer




I was wondering how mathematically is it possible to run a full regression analysis between 3 predictors (x1 x2 x3) and a dependent variable (y) by only knowing the: Means, Ns, SDs, and the Correlations between all these 4 variables (without the original data)?



I highly appreciate an R demonstration.



ns    <- c(273, 273, 273, 273)

means <- c(15.4, 7.1, 3.5, 6.2)

sds   <- c(3.4, 0.9, 1.5, 1.4)



r <- matrix( c(

1.0,  .57,  -.4,  .48,  

.57,  1.0, -.61,  .66,   

-.4, -.61,  1.0, -.68,  

.48,  .66, -.68,   1.0), 4)



rownames(r) <- colnames(r) <- c('y', paste0('x', 1:3))





r regression multiple-regression regression-coefficients







share|cite|improve this question











$endgroup$





This question already has an answer here:




  • Is there a way to use the covariance matrix to find coefficients for multiple regression?

    1 answer




I was wondering how mathematically is it possible to run a full regression analysis between 3 predictors (x1 x2 x3) and a dependent variable (y) by only knowing the: Means, Ns, SDs, and the Correlations between all these 4 variables (without the original data)?



I highly appreciate an R demonstration.



ns    <- c(273, 273, 273, 273)

means <- c(15.4, 7.1, 3.5, 6.2)

sds   <- c(3.4, 0.9, 1.5, 1.4)



r <- matrix( c(

1.0,  .57,  -.4,  .48,  

.57,  1.0, -.61,  .66,   

-.4, -.61,  1.0, -.68,  

.48,  .66, -.68,   1.0), 4)



rownames(r) <- colnames(r) <- c('y', paste0('x', 1:3))




This question already has an answer here:




  • Is there a way to use the covariance matrix to find coefficients for multiple regression?

    1 answer







r regression multiple-regression regression-coefficients




r regression multiple-regression regression-coefficients






share|cite|improve this question















share|cite|improve this question













share|cite|improve this question




share|cite|improve this question








edited Dec 10 '18 at 4:15







rnorouzian

















asked Dec 5 '18 at 3:49









rnorouzianrnorouzian

7081821




7081821




marked as duplicate by whuber♦ r
Users with the  r badge can single-handedly close r questions as duplicates and reopen them as needed.

StackExchange.ready(function() {
if (StackExchange.options.isMobile) return;

$('.dupe-hammer-message-hover:not(.hover-bound)').each(function() {
var $hover = $(this).addClass('hover-bound'),
$msg = $hover.siblings('.dupe-hammer-message');

$hover.hover(
function() {
$hover.showInfoMessage('', {
messageElement: $msg.clone().show(),
transient: false,
position: { my: 'bottom left', at: 'top center', offsetTop: -7 },
dismissable: false,
relativeToBody: true
});
},
function() {
StackExchange.helpers.removeMessages();
}
);
});
});
 Dec 5 '18 at 14:34


This question has been asked before and already has an answer. If those answers do not fully address your question, please ask a new question.









marked as duplicate by whuber♦ r
Users with the  r badge can single-handedly close r questions as duplicates and reopen them as needed.

StackExchange.ready(function() {
if (StackExchange.options.isMobile) return;

$('.dupe-hammer-message-hover:not(.hover-bound)').each(function() {
var $hover = $(this).addClass('hover-bound'),
$msg = $hover.siblings('.dupe-hammer-message');

$hover.hover(
function() {
$hover.showInfoMessage('', {
messageElement: $msg.clone().show(),
transient: false,
position: { my: 'bottom left', at: 'top center', offsetTop: -7 },
dismissable: false,
relativeToBody: true
});
},
function() {
StackExchange.helpers.removeMessages();
}
);
});
});
 Dec 5 '18 at 14:34


This question has been asked before and already has an answer. If those answers do not fully address your question, please ask a new question.










  • 1




    $begingroup$
    You cannot run a truly "full" regression analysis with just these statistics, because you will not be able to construct residuals and perform regression diagnostics that depend on them. You are limited to making and testing parameter estimates.
    $endgroup$
    – whuber♦
    Dec 5 '18 at 14:35










  • $begingroup$
    Structural equation modeling is the most obvious answer in my mind, which is more or less the math provided in the answer below. Most SEM programs will be able to accept your inputs (all need to be able to estimate a covariance matrix in the end).
    $endgroup$
    – Matt Barstead
    Dec 10 '18 at 2:29










  • $begingroup$
    @whuber, do I need to first convert my correlation matrix into a var-covariance matrix using r_x_iy_i * sd_x_i*sd_y_i and then go from there?
    $endgroup$
    – rnorouzian
    Dec 10 '18 at 4:52












  • $begingroup$
    That's one approach, because it reduces your problem to one with a known, explicit solution.
    $endgroup$
    – whuber♦
    Dec 10 '18 at 15:13














  • 1




    $begingroup$
    You cannot run a truly "full" regression analysis with just these statistics, because you will not be able to construct residuals and perform regression diagnostics that depend on them. You are limited to making and testing parameter estimates.
    $endgroup$
    – whuber♦
    Dec 5 '18 at 14:35










  • $begingroup$
    Structural equation modeling is the most obvious answer in my mind, which is more or less the math provided in the answer below. Most SEM programs will be able to accept your inputs (all need to be able to estimate a covariance matrix in the end).
    $endgroup$
    – Matt Barstead
    Dec 10 '18 at 2:29










  • $begingroup$
    @whuber, do I need to first convert my correlation matrix into a var-covariance matrix using r_x_iy_i * sd_x_i*sd_y_i and then go from there?
    $endgroup$
    – rnorouzian
    Dec 10 '18 at 4:52












  • $begingroup$
    That's one approach, because it reduces your problem to one with a known, explicit solution.
    $endgroup$
    – whuber♦
    Dec 10 '18 at 15:13








1




1




$begingroup$
You cannot run a truly "full" regression analysis with just these statistics, because you will not be able to construct residuals and perform regression diagnostics that depend on them. You are limited to making and testing parameter estimates.
$endgroup$
– whuber♦
Dec 5 '18 at 14:35




$begingroup$
You cannot run a truly "full" regression analysis with just these statistics, because you will not be able to construct residuals and perform regression diagnostics that depend on them. You are limited to making and testing parameter estimates.
$endgroup$
– whuber♦
Dec 5 '18 at 14:35












$begingroup$
Structural equation modeling is the most obvious answer in my mind, which is more or less the math provided in the answer below. Most SEM programs will be able to accept your inputs (all need to be able to estimate a covariance matrix in the end).
$endgroup$
– Matt Barstead
Dec 10 '18 at 2:29




$begingroup$
Structural equation modeling is the most obvious answer in my mind, which is more or less the math provided in the answer below. Most SEM programs will be able to accept your inputs (all need to be able to estimate a covariance matrix in the end).
$endgroup$
– Matt Barstead
Dec 10 '18 at 2:29












$begingroup$
@whuber, do I need to first convert my correlation matrix into a var-covariance matrix using r_x_iy_i * sd_x_i*sd_y_i and then go from there?
$endgroup$
– rnorouzian
Dec 10 '18 at 4:52






$begingroup$
@whuber, do I need to first convert my correlation matrix into a var-covariance matrix using r_x_iy_i * sd_x_i*sd_y_i and then go from there?
$endgroup$
– rnorouzian
Dec 10 '18 at 4:52














$begingroup$
That's one approach, because it reduces your problem to one with a known, explicit solution.
$endgroup$
– whuber♦
Dec 10 '18 at 15:13




$begingroup$
That's one approach, because it reduces your problem to one with a known, explicit solution.
$endgroup$
– whuber♦
Dec 10 '18 at 15:13










1 Answer
1






active

oldest

votes


















4












$begingroup$


  1. Based on that correlation matrix, you can estimate the standardized regression coefficients (excerpt intercept) by following. (suppose the first column is for y)


$$left(begin{matrix} 1.00 & -0.61 & 0.66\
-.61 & 1.00 & -0.68\
.66 & -.68 & 1.00end{matrix}right)^{-1}left(begin{matrix} 0.57\
-.40\
.48end{matrix}right)$$




  1. Combining standard deviations you can convert standardized regression coefficients into general regression coefficients.


  2. Using the information about means, you can get the estimate of intercept.







share|cite|improve this answer











$endgroup$









  • 1




    $begingroup$
    I do not know how to use R.
    $endgroup$
    – user158565
    Dec 5 '18 at 5:00






  • 1




    $begingroup$
    How does collinearity influence this?
    $endgroup$
    – PascalVKooten
    Dec 5 '18 at 9:25










  • $begingroup$
    Full collinearity ==> not exist of the inverse of that matrix in the answer, and the generalized inverse should be used and there are infinite number of estimates. Partial collinearity ==> the inverse of that matrix is unstable, i.e., the little change in X can result in tremendous change in estimate.
    $endgroup$
    – user158565
    Dec 5 '18 at 14:20


















1 Answer
1






active

oldest

votes








1 Answer
1






active

oldest

votes









active

oldest

votes






active

oldest

votes









4












$begingroup$


  1. Based on that correlation matrix, you can estimate the standardized regression coefficients (excerpt intercept) by following. (suppose the first column is for y)


$$left(begin{matrix} 1.00 & -0.61 & 0.66\
-.61 & 1.00 & -0.68\
.66 & -.68 & 1.00end{matrix}right)^{-1}left(begin{matrix} 0.57\
-.40\
.48end{matrix}right)$$




  1. Combining standard deviations you can convert standardized regression coefficients into general regression coefficients.


  2. Using the information about means, you can get the estimate of intercept.







share|cite|improve this answer











$endgroup$









  • 1




    $begingroup$
    I do not know how to use R.
    $endgroup$
    – user158565
    Dec 5 '18 at 5:00






  • 1




    $begingroup$
    How does collinearity influence this?
    $endgroup$
    – PascalVKooten
    Dec 5 '18 at 9:25










  • $begingroup$
    Full collinearity ==> not exist of the inverse of that matrix in the answer, and the generalized inverse should be used and there are infinite number of estimates. Partial collinearity ==> the inverse of that matrix is unstable, i.e., the little change in X can result in tremendous change in estimate.
    $endgroup$
    – user158565
    Dec 5 '18 at 14:20
















4












$begingroup$


  1. Based on that correlation matrix, you can estimate the standardized regression coefficients (excerpt intercept) by following. (suppose the first column is for y)


$$left(begin{matrix} 1.00 & -0.61 & 0.66\
-.61 & 1.00 & -0.68\
.66 & -.68 & 1.00end{matrix}right)^{-1}left(begin{matrix} 0.57\
-.40\
.48end{matrix}right)$$




  1. Combining standard deviations you can convert standardized regression coefficients into general regression coefficients.


  2. Using the information about means, you can get the estimate of intercept.







share|cite|improve this answer











$endgroup$









  • 1




    $begingroup$
    I do not know how to use R.
    $endgroup$
    – user158565
    Dec 5 '18 at 5:00






  • 1




    $begingroup$
    How does collinearity influence this?
    $endgroup$
    – PascalVKooten
    Dec 5 '18 at 9:25










  • $begingroup$
    Full collinearity ==> not exist of the inverse of that matrix in the answer, and the generalized inverse should be used and there are infinite number of estimates. Partial collinearity ==> the inverse of that matrix is unstable, i.e., the little change in X can result in tremendous change in estimate.
    $endgroup$
    – user158565
    Dec 5 '18 at 14:20














4












4








4





$begingroup$


  1. Based on that correlation matrix, you can estimate the standardized regression coefficients (excerpt intercept) by following. (suppose the first column is for y)


$$left(begin{matrix} 1.00 & -0.61 & 0.66\
-.61 & 1.00 & -0.68\
.66 & -.68 & 1.00end{matrix}right)^{-1}left(begin{matrix} 0.57\
-.40\
.48end{matrix}right)$$




  1. Combining standard deviations you can convert standardized regression coefficients into general regression coefficients.


  2. Using the information about means, you can get the estimate of intercept.







share|cite|improve this answer











$endgroup$




  1. Based on that correlation matrix, you can estimate the standardized regression coefficients (excerpt intercept) by following. (suppose the first column is for y)


$$left(begin{matrix} 1.00 & -0.61 & 0.66\
-.61 & 1.00 & -0.68\
.66 & -.68 & 1.00end{matrix}right)^{-1}left(begin{matrix} 0.57\
-.40\
.48end{matrix}right)$$




  1. Combining standard deviations you can convert standardized regression coefficients into general regression coefficients.


  2. Using the information about means, you can get the estimate of intercept.








share|cite|improve this answer














share|cite|improve this answer



share|cite|improve this answer








edited Dec 5 '18 at 5:58

























answered Dec 5 '18 at 4:29









user158565user158565

5,3921518




5,3921518








  • 1




    $begingroup$
    I do not know how to use R.
    $endgroup$
    – user158565
    Dec 5 '18 at 5:00






  • 1




    $begingroup$
    How does collinearity influence this?
    $endgroup$
    – PascalVKooten
    Dec 5 '18 at 9:25










  • $begingroup$
    Full collinearity ==> not exist of the inverse of that matrix in the answer, and the generalized inverse should be used and there are infinite number of estimates. Partial collinearity ==> the inverse of that matrix is unstable, i.e., the little change in X can result in tremendous change in estimate.
    $endgroup$
    – user158565
    Dec 5 '18 at 14:20














  • 1




    $begingroup$
    I do not know how to use R.
    $endgroup$
    – user158565
    Dec 5 '18 at 5:00






  • 1




    $begingroup$
    How does collinearity influence this?
    $endgroup$
    – PascalVKooten
    Dec 5 '18 at 9:25










  • $begingroup$
    Full collinearity ==> not exist of the inverse of that matrix in the answer, and the generalized inverse should be used and there are infinite number of estimates. Partial collinearity ==> the inverse of that matrix is unstable, i.e., the little change in X can result in tremendous change in estimate.
    $endgroup$
    – user158565
    Dec 5 '18 at 14:20








1




1




$begingroup$
I do not know how to use R.
$endgroup$
– user158565
Dec 5 '18 at 5:00




$begingroup$
I do not know how to use R.
$endgroup$
– user158565
Dec 5 '18 at 5:00




1




1




$begingroup$
How does collinearity influence this?
$endgroup$
– PascalVKooten
Dec 5 '18 at 9:25




$begingroup$
How does collinearity influence this?
$endgroup$
– PascalVKooten
Dec 5 '18 at 9:25












$begingroup$
Full collinearity ==> not exist of the inverse of that matrix in the answer, and the generalized inverse should be used and there are infinite number of estimates. Partial collinearity ==> the inverse of that matrix is unstable, i.e., the little change in X can result in tremendous change in estimate.
$endgroup$
– user158565
Dec 5 '18 at 14:20




$begingroup$
Full collinearity ==> not exist of the inverse of that matrix in the answer, and the generalized inverse should be used and there are infinite number of estimates. Partial collinearity ==> the inverse of that matrix is unstable, i.e., the little change in X can result in tremendous change in estimate.
$endgroup$
– user158565
Dec 5 '18 at 14:20



-

Popular posts from this blog

Plaza Victoria

Image
Plaza Victoria Plaza Victoria en 2017 Otros nombres Plaza de La Victoria Localización El Almendral, Valparaíso, Chile Vías adyacentes Chacabuco, Molina, Plaza Victoria y Edwards Creación 1841 Metro Bellavista Ubicación 33°02′46″S 71°37′11″O  /  -33.04624167, -71.61980833 Coordenadas: 33°02′46″S 71°37′11″O  /  -33.04624167, -71.61980833 [editar datos en Wikidata] La Plaza Victoria , también llamada Plaza de la Victoria , [ 1 ] ​ es una plaza ubicada en el sector El Almendral del plan de la ciudad de Valparaíso, Chile, de importancia patrimonial, que incluye valiosas estatuas y especies vegetales. [ 2 ] ​ La plaza se creó como tal en la primera mitad del siglo XIX, tomando previamente los nombres de Plaza Nueva y Plaza de Orrego . Si bien desde sus inicios tuvo relevancia histórica, no fue sino a partir de fines de los años 1860 cuando comenzó a cobrar mayor importancia y se convirtió en uno de los principales centros sociales de l...
Read more

Brian Clough

Image
Brian Clough Datos personales Nombre completo Brian Howard Clough Apodo(s) Old Big 'Ead, Cloughie [ 1 ] ​ Nacimiento Middlesbrough, Inglaterra, Reino Unido 21 de marzo de 1935 Nacionalidad(es) Fallecimiento Derby, Reino Unido 20 de septiembre de 2004 (69 años) Altura 1,78 [ 1 ] ​ metros Carrera Deporte Fútbol Debut deportivo 1955 (Como jugador) 1965 (Como entrenador) (Middlesbrough F. C. (Como jugador) Hartlepool United F. C. (Como entrenador) ) Posición Delantero Goles en clubes 269 [ 2 ] ​ Retirada deportiva 1964 (Como jugador) 1993 (Como entrenador) (Sunderland A. F. C. (Como jugador) Nottingham Forest (Como entrenador) ) Carrera internacional Part. 2 [editar datos en Wikidata] Brian Howard Clough , OBE [ 3 ] ​ (Middlesbrough, Inglaterra, Reino Unido, 21 de marzo de 1935 – Derby, Inglaterra, Reino Unido, 20 de septiembre de 2004) fue un jugador y entrenador británico de fútbol. Es conocido por...
Read more

Cáceres

Image
Para otros usos de este término, véase Cáceres (desambiguación). Cáceres Municipio y ciudad Bandera Escudo Collage de Cáceres Cáceres Ubicación de Cáceres en España. Cáceres Ubicación de Cáceres en la provincia de Cáceres. País  España • Com. autónoma  Extremadura • Provincia  Cáceres • Partido judicial Cáceres Ubicación 39°28′23″N 6°22′16″O  /  39.473055555556, -6.3711111111111 Coordenadas: 39°28′23″N 6°22′16″O  /  39.473055555556, -6.3711111111111 • Altitud 457 [ 1 ] ​ msnm (mín.:211 [ 2 ] ​, máx.:708 [ 2 ] ​) Superficie 1750,33 km² Núcleos de población Cáceres Estación Arroyo-Malpartida Rincón de Ballesteros Valdesalor Fundación 34 a. C. Población 96 068 hab. (2018) • Densidad 54,8 hab./km² Gentilicio Cacereño/a [ 3 ] ​ Código postal 10001-10005, 10195 Alcaldesa (2015) María Elena Nevado del Campo [ 4 ] ​ Presupuesto 69.045.004€ [ 5 ] ​  (año ...
Read more