# using aggregate( )

painting=read.csv("paint.csv",header=T)
painting

aggregate(reflect~brand,data=painting,mean)
aggregate(reflect~brand,data=painting,var)

# use stargazer( )
# make sure the package is intalled and loaded

stargazer(painting$reflect,title="Summary Statistics for Paint Brand",type="text",flip=T)