RNS_Figure_CumulativeTherapySeizureCorr.m

conn = database('RNS','USERNAME','PASSWORD', 'Vendor', 'Microsoft SQL Server', 'Server', 'IP ADDRESS', 'AuthType', 'Server', 'PortNumber', 1433);
sqlQuery = 'select rt.*, es.electrographic_seizures_calc_total, es.electrographic_seizures_calc_per_wk from rns_dm.sm_pe_nshd_weighted_responsive_therapies rt join rns_dm.sm_pe_nshd_weighted_electrographic_seizures es on rt.rns_deid_id = es.rns_deid_id and rt.programming_dt = es.programming_dt';
curs = exec(conn,sqlQuery);
curs = fetch(curs);

sData = cell2struct(curs.Data, columnnames(curs, true)', 2);

pt_id = extractfield(sData, 'rns_deid_id');
x = extractfield(sData, 'electrographic_seizures_calc_per_wk');
r = extractfield(sData, 'rt_total');

plot(x,r,'.', 'MarkerSize', 20);
lsline;
[rho,p]=corr(reshape(x,[],1),reshape(r,[],1));

ylabel('Cumulative Charge Density per Programming Epoch (mC/cm^2)');
xlabel('Electrographic Seizures per Week');