Re: Add 1 more variable to print
* ela:
I tried to add one more variable to print by using our previous programmer's
program but failed. I tried using pointer variable but also failed. the
compilable codes are as follows and
lines following the comment
"/////////////// I want to print MCC so I add this variable here"
were added by me.
Please check the FAQ on how to post a question about Code That Does Not Work.
Sorry but unable to provide the ana files to test but I wish this problem
can still be solved by somebody.
What *is* "this problem"?
#include <stdio.h>
#include <math.h>
#include <stdlib.h>
#include <assert.h>
#include <map>
#include <vector>
#include <set>
using namespace std;
const char* anaFilename = "%d.txt";
Bad idea.
char buffer[65536];
char buffer2[65536];
Bad idea.
void wrongArgv(char **argv){
fprintf(stderr,
"Usage: %s [select how many] [total # model] \n",\
*argv);
Bad idea.
exit(1);
}
int main(int argc, char **argv){
map<char*, vector<double> > tptable;
Bad idea.
if(argc != 3)
wrongArgv(argv);
int select = atoi(argv[1]);
int all = atoi(argv[2]);
if(select<=0 || all<=0 || select > all)
wrongArgv(argv);
// use map for priority queue
map<double, int> mcc; //matthews correlation coefficient
Bad idea.
for(int i=0; i<all; ++i){
sprintf(buffer,anaFilename, i);
Unnecessary.
FILE *ana = fopen(buffer, "r");
for(int j=1; j<=17; ++j)
fgets(buffer, sizeof(buffer), ana);
double TP, TN, FP, FN, MCC;
Bad idea.
{
int zero, one;
fgets(buffer, sizeof(buffer), ana);
assert(3==sscanf(buffer," %d %lf %lf",&zero, &TN, &FP));
assert(0==zero);
fgets(buffer, sizeof(buffer), ana);
assert(3==sscanf(buffer," %d %lf %lf",&one, &FN, &TP));
assert(1==one);
/////////////// I want to print MCC so I add this variable here
MCC = (TP*TN-FP*FN)/ sqrt((TP+FP)*(TP+FN)*(TN+FP)*(TN+FN));
tptable["TP"].push_back(TP);
tptable["TN"].push_back(TN);
tptable["FP"].push_back(FP);
tptable["FN"].push_back(FN);
Bad idea.
/////////////// I want to print MCC so I add this variable here
tptable["MCC"].push_back(MCC);
}
assert(mcc.find(MCC)==mcc.end());
Unnecessary.
mcc[MCC] = i;
Bad idea.
fclose(ana);
}
set<int> top;
map<double,int>::reverse_iterator riter=mcc.rbegin();
Bad idea.
for(int i=0;
i<select;
++i, ++riter){
top.insert((*riter).second);
}
{
FILE *file = fopen("Eval.txt","w");
fprintf(file,"PPV\tNPV\tSP\tSE\tMCC\n");
for(int i=0; i<all; ++i){
double tp = tptable["TP"][i];
double fp = tptable["FP"][i];
double tn = tptable["TN"][i];
double fn = tptable["FN"][i];
Bad idea.
/////////////// I want to print MCC so I add this variable here
double tpmcc = tptable["MCC"][i];
fprintf(file,"%1.3f \t%1.3f\t%1.3f\t%1.3f\n",ppv,npv,sp,se,tpmcc);
I count four % and five args.
}
freopen("SEL.txt","w",file);
fprintf(file,"%d",*top.begin());
for(set<int>::iterator it = ++top.begin();
it!=top.end();
++it){
fprintf(file," %d",*it);
}
fclose(file);
}
return 0;
}
This code is horrible.
But presumably your problem is identified by the last of my in-code comments.
Cheers & hth.,
- Alf
--
Due to hosting requirements I need visits to <url: http://alfps.izfree.com/>.
No ads, and there is some C++ stuff! :-) Just going there is good. Linking
to it is even better! Thanks in advance!