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)