def all_scatters(X,names,casename):
'''
scatter plots of all variables against one another in the X matrix
'''
from scat_plotter import scat_plot
nv = len(names)
k=1
for i in np.arange(nv):
for j in np.arange(nv):
k+=1
scat_plot(k,X[:,i],X[:,j],names[i] + ' ' + names[j],names[i],names[j],
casename + "_" + names[i] + "_" + names[j] + '.png',False)
os.mkdir(newpath)
newpath += '/'
xname = ['SULFATE','TOC','pH','Wetland %']
for i in xrange(RHS.shape[1]):
xvals = RHS[:,i]
yvals = HG_reg
title = 'Scatter ' + xname[i]
if (cstudy < 9999):
outfile = xname[i] + '_' + currstud + '_scatter.png'
else:
outfile = xname[i] + '_scatter.png'
outfile = newpath + outfile
showme = False
sp.scat_plot(i,xvals,yvals,title,xname[i],outfile,HG_RAT_FLAG,showme)
if (cstudy < 9999):
outHGfile = 'methyl_to_total_HG'+ '_' + currstud + '.png'
else:
outHGfile = 'methyl_to_total_HG.png'
outHGfile = newpath + outHGfile
sp.scat_plot(10,np.arange(1,len(HG_reg)+1),HG_reg,hglabel,'meas. number',outHGfile,HG_RAT_FLAG,showme)
HGr = rob.conversion.py2ri(HG_reg)
RHSr = rob.conversion.py2ri(RHS)
#
CalData = rlc.TaggedList([HGr,RHSr],tags=('hg','rhs'))
CalData = rob.DataFrame(CalData)
