def getQuantAlphaData(request): samples = Sample.objects.all() samples.query = pickle.loads(request.session['selected_samples']) selected = samples.values_list('sampleid') qs1 = Sample.objects.all().filter(sampleid__in=selected) if request.is_ajax(): allJson = request.GET["all"] all = simplejson.loads(allJson) button = int(all["button"]) sig_only = int(all["sig_only"]) norm = int(all["normalize"]) taxaString = all["taxa"] taxaDict = simplejson.JSONDecoder(object_pairs_hook=multidict).decode(taxaString) metaString = all["meta"] metaDict = simplejson.JSONDecoder(object_pairs_hook=multidict).decode(metaString) metaDF = quantAlphaMetaDF(qs1, metaDict) myList = metaDF['sampleid'].tolist() mySet = list(set(myList)) taxaDF = taxaProfileDF(mySet) factor = 'none' if norm == 1: factor = 'none' elif norm == 2: factor = 'min' elif norm == 3: factor = '10th percentile' elif norm == 4: factor = '25th percentile' elif norm == 5: factor = 'median' final_fieldList = [] for key in metaDict: final_fieldList.append(metaDict[key]) normDF = normalizeAlpha(taxaDF, taxaDict, mySet, factor) finalDF = metaDF.merge(normDF, on='sampleid', how='outer') finalDF[[final_fieldList[0], 'count', 'rel_abund', 'rich', 'diversity']] = finalDF[[final_fieldList[0], 'count', 'rel_abund', 'rich', 'diversity']].astype(float) pd.set_option('display.max_rows', finalDF.shape[0], 'display.max_columns', finalDF.shape[1], 'display.width', 1000) finalDict = {} seriesList = [] xAxisDict = {} yAxisDict = {} grouped1 = finalDF.groupby(['rank', 'taxa_name', 'taxa_id']) for name1, group1 in grouped1: dataList = [] x = [] y = [] if button == 1: dataList = group1[[final_fieldList[0], 'count']].values.tolist() x = group1[final_fieldList[0]].values.tolist() y = group1['count'].values.tolist() elif button == 2: dataList = group1[[final_fieldList[0], 'rel_abund']].values.tolist() x = group1[final_fieldList[0]].values.tolist() y = group1['rel_abund'].values.tolist() elif button == 3: dataList = group1[[final_fieldList[0], 'rich']].values.tolist() x = group1[final_fieldList[0]].values.tolist() y = group1['rich'].values.tolist() elif button == 4: dataList = group1[[final_fieldList[0], 'diversity']].values.tolist() x = group1[final_fieldList[0]].values.tolist() y = group1['diversity'].values.tolist() if max(x) == min(x): stop = 0 else: stop = 1 slope, intercept, r_value, p_value, std_err = stats.linregress(x, y) p_value = "%0.3f" % p_value r_square = r_value * r_value r_square = "%0.4f" % r_square min_y = slope*min(x) + intercept max_y = slope*max(x) + intercept slope = "%.3E" % slope intercept = "%.3E" % intercept regrList = [] regrList.append([min(x), min_y]) regrList.append([max(x), max_y]) if sig_only == 0: seriesDict = {} seriesDict['type'] = 'scatter' seriesDict['name'] = name1 seriesDict['data'] = dataList seriesList.append(seriesDict) if stop == 0: regDict = {} elif stop == 1: regrDict = {} regrDict['type'] = 'line' name2 = list(name1) temp = 'R2: ' + str(r_square) + '; p-value: ' + str(p_value) + '<br>' + '(y = ' + str(slope) + 'x' + ' + ' + str(intercept) print temp name2.append(temp) print name2 regrDict['name'] = name2 regrDict['data'] = regrList seriesList.append(regrDict) if sig_only == 1: if p_value <= 0.05: seriesDict = {} seriesDict['type'] = 'scatter' name2 = list(name1) temp = 'R2: ' + str(r_square) + '; p-value: ' + str(p_value) + '<br>' + '(y = ' + str(slope) + 'x' + ' + ' + str(intercept) name2.append(temp) seriesDict['name'] = name2 seriesDict['data'] = dataList seriesList.append(seriesDict) regrDict = {} regrDict['type'] = 'line' regrDict['name'] = name1 regrDict['data'] = regrList seriesList.append(regrDict) xTitle = {} xTitle['text'] = final_fieldList[0] xAxisDict['title'] = xTitle yTitle = {} if button == 1: yTitle['text'] = 'Sequence Reads' elif button == 2: yTitle['text'] = 'Relative Abundance' elif button == 3: yTitle['text'] = 'Species Richness' elif button == 4: yTitle['text'] = 'Shannon Diversity' yAxisDict['title'] = yTitle finalDict['series'] = seriesList finalDict['xAxis'] = xAxisDict finalDict['yAxis'] = yAxisDict if not seriesList: finalDict['empty'] = 0 else: finalDict['empty'] = 1 finalDF.reset_index(drop=True, inplace=True) res_table = finalDF.to_html(classes="table display") res_table = res_table.replace('border="1"', 'border="0"') finalDict['res_table'] = str(res_table) res = simplejson.dumps(finalDict) return HttpResponse(res, content_type='application/json')
def getQuantAlphaData(request): samples = Sample.objects.all() samples.query = pickle.loads(request.session['selected_samples']) selected = samples.values_list('sampleid') qs1 = Sample.objects.all().filter(sampleid__in=selected) if request.is_ajax(): allJson = request.GET["all"] all = simplejson.loads(allJson) button = int(all["button"]) sig_only = int(all["sig_only"]) norm = int(all["normalize"]) taxaString = all["taxa"] taxaDict = simplejson.JSONDecoder( object_pairs_hook=multidict).decode(taxaString) metaString = all["meta"] metaDict = simplejson.JSONDecoder( object_pairs_hook=multidict).decode(metaString) metaDF = quantAlphaMetaDF(qs1, metaDict) myList = metaDF['sampleid'].tolist() mySet = list(set(myList)) taxaDF = taxaProfileDF(mySet) factor = 'none' if norm == 1: factor = 'none' elif norm == 2: factor = 'min' elif norm == 3: factor = '10th percentile' elif norm == 4: factor = '25th percentile' elif norm == 5: factor = 'median' final_fieldList = [] for key in metaDict: final_fieldList.append(metaDict[key]) normDF = normalizeAlpha(taxaDF, taxaDict, mySet, factor) finalDF = metaDF.merge(normDF, on='sampleid', how='outer') finalDF[[ final_fieldList[0], 'count', 'rel_abund', 'rich', 'diversity' ]] = finalDF[[ final_fieldList[0], 'count', 'rel_abund', 'rich', 'diversity' ]].astype(float) pd.set_option('display.max_rows', finalDF.shape[0], 'display.max_columns', finalDF.shape[1], 'display.width', 1000) finalDict = {} seriesList = [] xAxisDict = {} yAxisDict = {} grouped1 = finalDF.groupby(['rank', 'taxa_name', 'taxa_id']) for name1, group1 in grouped1: dataList = [] x = [] y = [] if button == 1: dataList = group1[[final_fieldList[0], 'count']].values.tolist() x = group1[final_fieldList[0]].values.tolist() y = group1['count'].values.tolist() elif button == 2: dataList = group1[[final_fieldList[0], 'rel_abund']].values.tolist() x = group1[final_fieldList[0]].values.tolist() y = group1['rel_abund'].values.tolist() elif button == 3: dataList = group1[[final_fieldList[0], 'rich']].values.tolist() x = group1[final_fieldList[0]].values.tolist() y = group1['rich'].values.tolist() elif button == 4: dataList = group1[[final_fieldList[0], 'diversity']].values.tolist() x = group1[final_fieldList[0]].values.tolist() y = group1['diversity'].values.tolist() if max(x) == min(x): stop = 0 else: stop = 1 slope, intercept, r_value, p_value, std_err = stats.linregress( x, y) p_value = "%0.3f" % p_value r_square = r_value * r_value r_square = "%0.4f" % r_square min_y = slope * min(x) + intercept max_y = slope * max(x) + intercept slope = "%.3E" % slope intercept = "%.3E" % intercept regrList = [] regrList.append([min(x), min_y]) regrList.append([max(x), max_y]) if sig_only == 0: seriesDict = {} seriesDict['type'] = 'scatter' seriesDict['name'] = name1 seriesDict['data'] = dataList seriesList.append(seriesDict) if stop == 0: regDict = {} elif stop == 1: regrDict = {} regrDict['type'] = 'line' name2 = list(name1) temp = 'R2: ' + str(r_square) + '; p-value: ' + str( p_value) + '<br>' + '(y = ' + str( slope) + 'x' + ' + ' + str(intercept) print temp name2.append(temp) print name2 regrDict['name'] = name2 regrDict['data'] = regrList seriesList.append(regrDict) if sig_only == 1: if p_value <= 0.05: seriesDict = {} seriesDict['type'] = 'scatter' name2 = list(name1) temp = 'R2: ' + str(r_square) + '; p-value: ' + str( p_value) + '<br>' + '(y = ' + str( slope) + 'x' + ' + ' + str(intercept) name2.append(temp) seriesDict['name'] = name2 seriesDict['data'] = dataList seriesList.append(seriesDict) regrDict = {} regrDict['type'] = 'line' regrDict['name'] = name1 regrDict['data'] = regrList seriesList.append(regrDict) xTitle = {} xTitle['text'] = final_fieldList[0] xAxisDict['title'] = xTitle yTitle = {} if button == 1: yTitle['text'] = 'Sequence Reads' elif button == 2: yTitle['text'] = 'Relative Abundance' elif button == 3: yTitle['text'] = 'Species Richness' elif button == 4: yTitle['text'] = 'Shannon Diversity' yAxisDict['title'] = yTitle finalDict['series'] = seriesList finalDict['xAxis'] = xAxisDict finalDict['yAxis'] = yAxisDict if not seriesList: finalDict['empty'] = 0 else: finalDict['empty'] = 1 finalDF.reset_index(drop=True, inplace=True) res_table = finalDF.to_html(classes="table display") res_table = res_table.replace('border="1"', 'border="0"') finalDict['res_table'] = str(res_table) res = simplejson.dumps(finalDict) return HttpResponse(res, content_type='application/json')
def getCatAlphaData(request): ### get sample list from cookie samples = Sample.objects.all() samples.query = pickle.loads(request.session['selected_samples']) selected = samples.values_list('sampleid') qs1 = Sample.objects.all().filter(sampleid__in=selected) if request.is_ajax(): allJson = request.GET["all"] all = simplejson.loads(allJson) button = int(all["button"]) sig_only = int(all["sig_only"]) norm = int(all["normalize"]) selectAll = int(all["selectAll"]) metaString = all["meta"] ### function to merge values on common keys metaDict = simplejson.JSONDecoder(object_pairs_hook=multidict).decode(metaString) ### function to create a meta variable DataFrame metaDF = catAlphaMetaDF(qs1, metaDict) myList = metaDF['sampleid'].tolist() mySet = list(set(myList)) ### function to create a taxa DataFrame taxaDF = taxaProfileDF(mySet) ### function to merge values on common keys taxaString = all["taxa"] ### this taxaDict is from the dynatree (ajax call) taxaDict = simplejson.JSONDecoder(object_pairs_hook=multidict).decode(taxaString) # change dict if selectAll levels is on (avoids loading entire tree first) if selectAll == 1: taxaDict = {} qs1 = Profile.objects.all().filter(sampleid__in=mySet).values_list('kingdomid', flat='True').distinct() taxaDict['Kingdom'] = qs1 elif selectAll == 2: taxaDict = {} qs1 = Profile.objects.all().filter(sampleid__in=mySet).values_list('phylaid', flat='True').distinct() taxaDict['Phyla'] = qs1 elif selectAll == 3: taxaDict = {} qs1 = Profile.objects.all().filter(sampleid__in=mySet).values_list('classid', flat='True').distinct() taxaDict['Class'] = qs1 elif selectAll == 4: taxaDict = {} qs1 = Profile.objects.all().filter(sampleid__in=mySet).values_list('orderid', flat='True').distinct() taxaDict['Order'] = qs1 elif selectAll == 5: taxaDict = {} qs1 = Profile.objects.all().filter(sampleid__in=mySet).values_list('familyid', flat='True').distinct() taxaDict['Family'] = qs1 elif selectAll == 6: taxaDict = {} qs1 = Profile.objects.all().filter(sampleid__in=mySet).values_list('genusid', flat='True').distinct() taxaDict['Genus'] = qs1 elif selectAll == 7: taxaDict = {} qs1 = Profile.objects.all().filter(sampleid__in=mySet).values_list('speciesid', flat='True').distinct() taxaDict['Species'] = qs1 factor = 'none' if norm == 1: factor = 'none' elif norm == 2: factor = 'min' elif norm == 3: factor = '10th percentile' elif norm == 4: factor = '25th percentile' elif norm == 5: factor = 'median' ### function to normalize the number of sequence reads per sample normDF = normalizeAlpha(taxaDF, taxaDict, mySet, factor) finalDF = metaDF.merge(normDF, on='sampleid', how='outer') finalDF[['count', 'rel_abund', 'rich', 'diversity']] = finalDF[['count', 'rel_abund', 'rich', 'diversity']].astype(float) pd.set_option('display.max_rows', finalDF.shape[0], 'display.max_columns', finalDF.shape[1], 'display.width', 1000) final_fieldList = [] for key in metaDict: final_fieldList.append(key) finalDict = {} result = "" seriesList = [] xAxisDict = {} yAxisDict = {} ### group DataFrame by each taxa level selected grouped1 = finalDF.groupby(['rank', 'taxa_name', 'taxa_id']) equal_error = 'no' ### group DataFrame by each meta variable selected for name1, group1 in grouped1: trtList = [] valList = [] grouped2 = pd.DataFrame() if button == 1: grouped2 = group1.groupby(final_fieldList)['count'] elif button == 2: grouped2 = group1.groupby(final_fieldList)['rel_abund'] elif button == 3: grouped2 = group1.groupby(final_fieldList)['rich'] ### for taxa with only 1 species all values will be '1' and cause an anova error if group1['rich'].sum() == group1['rich'].count(): equal_error = 'yes' elif button == 4: grouped2 = group1.groupby(final_fieldList)['diversity'] for name2, group2 in grouped2: if isinstance(name2, unicode): trt = name2 else: trt = ' & '.join(list(name2)) trtList.append(trt) valList.append(list(group2.T)) ### One-way ANOVA with some error checking D = Anova1way() if equal_error == 'no': try: D.run(valList, conditions_list=trtList) anova_error = 'no' except: D['p'] = 1 anova_error = 'yes' else: D['p'] = 1 anova_error = 'yes' ### select only significant ANOVAs for output (graph & text area) if sig_only == 1: if D['p'] <= 0.05: result = result + '===============================================\n' result = result + 'Taxa level: ' + str(name1[0]) + '\n' result = result + 'Taxa name: ' + str(name1[1]) + '\n' if button == 1: result = result + 'Dependent Variable: Sequence Reads' + '\n' elif button == 2: result = result + 'Dependent Variable: Relative Abundance' + '\n' elif button == 3: result = result + 'Dependent Variable: Species Richness' + '\n' elif button == 4: result = result + 'Dependent Variable: Shannon Diversity' + '\n' indVar = ' x '.join(final_fieldList) result = result + 'Independent Variable: ' + str(indVar) + '\n' if equal_error == 'yes' or anova_error == 'yes': result = result + 'Analysis cannot be performed...' + '\n' else: result = result + str(D) + '\n' result = result + '===============================================\n' result = result + '\n\n\n\n' dataList = [] grouped2 = group1.groupby(final_fieldList).mean() if button == 1: dataList.extend(list(grouped2['count'].T)) elif button == 2: dataList.extend(list(grouped2['rel_abund'].T)) elif button == 3: dataList.extend(list(grouped2['rich'].T)) elif button == 4: dataList.extend(list(grouped2['diversity'].T)) seriesDict = {} seriesDict['name'] = name1 seriesDict['data'] = dataList seriesList.append(seriesDict) xTitle = {} xTitle['text'] = indVar xAxisDict['title'] = xTitle xAxisDict['categories'] = trtList yTitle = {} if button == 1: yTitle['text'] = 'Sequence Reads' elif button == 2: yTitle['text'] = 'Relative Abundance' elif button == 3: yTitle['text'] = 'Species Richness' elif button == 4: yTitle['text'] = 'Shannon Diversity' yAxisDict['title'] = yTitle ### select all ANOVAs for output (graph & text area) if sig_only == 0: result = result + '===============================================\n' result = result + 'Taxa level: ' + str(name1[0]) + '\n' result = result + 'Taxa name: ' + str(name1[1]) + '\n' if button == 1: result = result + 'Dependent Variable: Sequence Reads' + '\n' elif button == 2: result = result + 'Dependent Variable: Relative Abundance' + '\n' elif button == 3: result = result + 'Dependent Variable: Species Richness' + '\n' elif button == 4: result = result + 'Dependent Variable: Shannon Diversity' + '\n' indVar = ' x '.join(final_fieldList) result = result + 'Independent Variable: ' + str(indVar) + '\n' if equal_error == 'yes' or anova_error == 'yes': result = result + 'Analysis cannot be performed...' + '\n' else: result = result + str(D) + '\n' result = result + '===============================================\n' result = result + '\n\n\n\n' dataList = [] grouped2 = group1.groupby(final_fieldList).mean() if button == 1: dataList.extend(list(grouped2['count'].T)) elif button == 2: dataList.extend(list(grouped2['rel_abund'].T)) elif button == 3: dataList.extend(list(grouped2['rich'].T)) elif button == 4: dataList.extend(list(grouped2['diversity'].T)) seriesDict = {} seriesDict['name'] = name1 seriesDict['data'] = dataList seriesList.append(seriesDict) xTitle = {} xTitle['text'] = indVar xAxisDict['title'] = xTitle xAxisDict['categories'] = trtList yTitle = {} if button == 1: yTitle['text'] = 'Sequence Reads' elif button == 2: yTitle['text'] = 'Relative Abundance' elif button == 3: yTitle['text'] = 'Species Richness' elif button == 4: yTitle['text'] = 'Shannon Diversity' yAxisDict['title'] = yTitle finalDict['series'] = seriesList finalDict['xAxis'] = xAxisDict finalDict['yAxis'] = yAxisDict finalDict['text'] = result if not seriesList: finalDict['empty'] = 0 else: finalDict['empty'] = 1 finalDF.reset_index(drop=True, inplace=True) res_table = finalDF.to_html(classes="table display") res_table = res_table.replace('border="1"', 'border="0"') finalDict['res_table'] = str(res_table) res = simplejson.dumps(finalDict) return HttpResponse(res, content_type='application/json')
def getCatAlphaData(request): ### get sample list from cookie samples = Sample.objects.all() samples.query = pickle.loads(request.session['selected_samples']) selected = samples.values_list('sampleid') qs1 = Sample.objects.all().filter(sampleid__in=selected) if request.is_ajax(): allJson = request.GET["all"] all = simplejson.loads(allJson) button = int(all["button"]) sig_only = int(all["sig_only"]) norm = int(all["normalize"]) selectAll = int(all["selectAll"]) metaString = all["meta"] ### function to merge values on common keys metaDict = simplejson.JSONDecoder( object_pairs_hook=multidict).decode(metaString) ### function to create a meta variable DataFrame metaDF = catAlphaMetaDF(qs1, metaDict) myList = metaDF['sampleid'].tolist() mySet = list(set(myList)) ### function to create a taxa DataFrame taxaDF = taxaProfileDF(mySet) ### function to merge values on common keys taxaString = all["taxa"] ### this taxaDict is from the dynatree (ajax call) taxaDict = simplejson.JSONDecoder( object_pairs_hook=multidict).decode(taxaString) # change dict if selectAll levels is on (avoids loading entire tree first) if selectAll == 1: taxaDict = {} qs1 = Profile.objects.all().filter(sampleid__in=mySet).values_list( 'kingdomid', flat='True').distinct() taxaDict['Kingdom'] = qs1 elif selectAll == 2: taxaDict = {} qs1 = Profile.objects.all().filter(sampleid__in=mySet).values_list( 'phylaid', flat='True').distinct() taxaDict['Phyla'] = qs1 elif selectAll == 3: taxaDict = {} qs1 = Profile.objects.all().filter(sampleid__in=mySet).values_list( 'classid', flat='True').distinct() taxaDict['Class'] = qs1 elif selectAll == 4: taxaDict = {} qs1 = Profile.objects.all().filter(sampleid__in=mySet).values_list( 'orderid', flat='True').distinct() taxaDict['Order'] = qs1 elif selectAll == 5: taxaDict = {} qs1 = Profile.objects.all().filter(sampleid__in=mySet).values_list( 'familyid', flat='True').distinct() taxaDict['Family'] = qs1 elif selectAll == 6: taxaDict = {} qs1 = Profile.objects.all().filter(sampleid__in=mySet).values_list( 'genusid', flat='True').distinct() taxaDict['Genus'] = qs1 elif selectAll == 7: taxaDict = {} qs1 = Profile.objects.all().filter(sampleid__in=mySet).values_list( 'speciesid', flat='True').distinct() taxaDict['Species'] = qs1 factor = 'none' if norm == 1: factor = 'none' elif norm == 2: factor = 'min' elif norm == 3: factor = '10th percentile' elif norm == 4: factor = '25th percentile' elif norm == 5: factor = 'median' ### function to normalize the number of sequence reads per sample normDF = normalizeAlpha(taxaDF, taxaDict, mySet, factor) finalDF = metaDF.merge(normDF, on='sampleid', how='outer') finalDF[['count', 'rel_abund', 'rich', 'diversity' ]] = finalDF[['count', 'rel_abund', 'rich', 'diversity']].astype(float) pd.set_option('display.max_rows', finalDF.shape[0], 'display.max_columns', finalDF.shape[1], 'display.width', 1000) final_fieldList = [] for key in metaDict: final_fieldList.append(key) finalDict = {} result = "" seriesList = [] xAxisDict = {} yAxisDict = {} ### group DataFrame by each taxa level selected grouped1 = finalDF.groupby(['rank', 'taxa_name', 'taxa_id']) equal_error = 'no' ### group DataFrame by each meta variable selected for name1, group1 in grouped1: trtList = [] valList = [] grouped2 = pd.DataFrame() if button == 1: grouped2 = group1.groupby(final_fieldList)['count'] elif button == 2: grouped2 = group1.groupby(final_fieldList)['rel_abund'] elif button == 3: grouped2 = group1.groupby(final_fieldList)['rich'] ### for taxa with only 1 species all values will be '1' and cause an anova error if group1['rich'].sum() == group1['rich'].count(): equal_error = 'yes' elif button == 4: grouped2 = group1.groupby(final_fieldList)['diversity'] for name2, group2 in grouped2: if isinstance(name2, unicode): trt = name2 else: trt = ' & '.join(list(name2)) trtList.append(trt) valList.append(list(group2.T)) ### One-way ANOVA with some error checking D = Anova1way() if equal_error == 'no': try: D.run(valList, conditions_list=trtList) anova_error = 'no' except: D['p'] = 1 anova_error = 'yes' else: D['p'] = 1 anova_error = 'yes' ### select only significant ANOVAs for output (graph & text area) if sig_only == 1: if D['p'] <= 0.05: result = result + '===============================================\n' result = result + 'Taxa level: ' + str(name1[0]) + '\n' result = result + 'Taxa name: ' + str(name1[1]) + '\n' if button == 1: result = result + 'Dependent Variable: Sequence Reads' + '\n' elif button == 2: result = result + 'Dependent Variable: Relative Abundance' + '\n' elif button == 3: result = result + 'Dependent Variable: Species Richness' + '\n' elif button == 4: result = result + 'Dependent Variable: Shannon Diversity' + '\n' indVar = ' x '.join(final_fieldList) result = result + 'Independent Variable: ' + str( indVar) + '\n' if equal_error == 'yes' or anova_error == 'yes': result = result + 'Analysis cannot be performed...' + '\n' else: result = result + str(D) + '\n' result = result + '===============================================\n' result = result + '\n\n\n\n' dataList = [] grouped2 = group1.groupby(final_fieldList).mean() if button == 1: dataList.extend(list(grouped2['count'].T)) elif button == 2: dataList.extend(list(grouped2['rel_abund'].T)) elif button == 3: dataList.extend(list(grouped2['rich'].T)) elif button == 4: dataList.extend(list(grouped2['diversity'].T)) seriesDict = {} seriesDict['name'] = name1 seriesDict['data'] = dataList seriesList.append(seriesDict) xTitle = {} xTitle['text'] = indVar xAxisDict['title'] = xTitle xAxisDict['categories'] = trtList yTitle = {} if button == 1: yTitle['text'] = 'Sequence Reads' elif button == 2: yTitle['text'] = 'Relative Abundance' elif button == 3: yTitle['text'] = 'Species Richness' elif button == 4: yTitle['text'] = 'Shannon Diversity' yAxisDict['title'] = yTitle ### select all ANOVAs for output (graph & text area) if sig_only == 0: result = result + '===============================================\n' result = result + 'Taxa level: ' + str(name1[0]) + '\n' result = result + 'Taxa name: ' + str(name1[1]) + '\n' if button == 1: result = result + 'Dependent Variable: Sequence Reads' + '\n' elif button == 2: result = result + 'Dependent Variable: Relative Abundance' + '\n' elif button == 3: result = result + 'Dependent Variable: Species Richness' + '\n' elif button == 4: result = result + 'Dependent Variable: Shannon Diversity' + '\n' indVar = ' x '.join(final_fieldList) result = result + 'Independent Variable: ' + str(indVar) + '\n' if equal_error == 'yes' or anova_error == 'yes': result = result + 'Analysis cannot be performed...' + '\n' else: result = result + str(D) + '\n' result = result + '===============================================\n' result = result + '\n\n\n\n' dataList = [] grouped2 = group1.groupby(final_fieldList).mean() if button == 1: dataList.extend(list(grouped2['count'].T)) elif button == 2: dataList.extend(list(grouped2['rel_abund'].T)) elif button == 3: dataList.extend(list(grouped2['rich'].T)) elif button == 4: dataList.extend(list(grouped2['diversity'].T)) seriesDict = {} seriesDict['name'] = name1 seriesDict['data'] = dataList seriesList.append(seriesDict) xTitle = {} xTitle['text'] = indVar xAxisDict['title'] = xTitle xAxisDict['categories'] = trtList yTitle = {} if button == 1: yTitle['text'] = 'Sequence Reads' elif button == 2: yTitle['text'] = 'Relative Abundance' elif button == 3: yTitle['text'] = 'Species Richness' elif button == 4: yTitle['text'] = 'Shannon Diversity' yAxisDict['title'] = yTitle finalDict['series'] = seriesList finalDict['xAxis'] = xAxisDict finalDict['yAxis'] = yAxisDict finalDict['text'] = result if not seriesList: finalDict['empty'] = 0 else: finalDict['empty'] = 1 finalDF.reset_index(drop=True, inplace=True) res_table = finalDF.to_html(classes="table display") res_table = res_table.replace('border="1"', 'border="0"') finalDict['res_table'] = str(res_table) res = simplejson.dumps(finalDict) return HttpResponse(res, content_type='application/json')