Hi new here. Got this question for a brief to complete:
Run a regression of datsetX on b , c , d , e , and f . Report the estimated regression with robust standard errors.

I run the regression as normal assigning lm variable and running it, i get the call, the residuals, the coefficients but its specifically asking to report a robust standard error.

Do i then just run sandwich(regression) for the robust standard error? This returns me with the variables being the same in both rows and columns and showing me the intercept.

What you're getting is the matrix of the robust variances and covariances of the coefficient estimates. If you want the standard errors, take the square root of the diagonal elements. By default you're getting heteroskedasticity and autocorrelation robust standard errors. You may want to be robust only to the former. In that case use vcovHC(). Or try

Ok done this. Got the t-test with the standard errors.

I assume now (report the robust standard errors) that the answer would be just giving the standard errors of that for each variable. If so, done and thank you.

I've been looking everywhere and just couldn't understand this heteroskedasticity concept.

In the standard regression model, we assume the error terms all have the same variance. Heteroskedasticity just means the errors have differing variances.