public void run(String arg) {
ImageWindow iw = WindowManager.getCurrentWindow();
pw = jutils.getPW4SelCopy(iw);
String title = pw.getTitle();
float[][] yvals = pw.getYValues();
float[][] xvals = pw.getXValues();
int length = yvals[0].length;
if (pw.getShowErrors()) errs = pw.getErrors(0, false);
int[] colors = pw.getColors();
colors[0] = 0;
ScriptEngineManager manager = new ScriptEngineManager();
engine = manager.getEngineByName("js");
ce = (Compilable) engine;
// hitcounter=0;
c2 = 0.0f;
iterations = 0;
checkc2 = false;
double[] stats = new double[3];
tempx = new float[length];
tempdata = new float[length];
System.arraycopy(xvals[0], 0, tempx, 0, length);
System.arraycopy(yvals[0], 0, tempdata, 0, length);
pw.addPoints(tempx, new float[tempx.length], false);
series = pw.getNpts().length - 1;
double[] params = new double[10];
int[] fixes = {0, 0, 0, 1, 1, 1, 1, 1, 1, 1};
init_options(params, fixes);
if (!init_functions()) {
return;
}
while (showoptions(params, fixes)) {
NLLSfit_v2 fitclass;
if (checkc2) {
fitclass = new NLLSfit_v2(this, 0);
} else {
fitclass = new NLLSfit_v2(this, 0.0001, 50, 0.1);
}
float[] fit = fitclass.fitdata(params, fixes, constraints, yvals[0], weights, stats, true);
pw.updateSeries(fit, series, false);
c2 = (float) stats[1];
iterations = (int) stats[0];
}
IJ.log("Chi Squared = " + (float) stats[1]);
IJ.log("Iterations = " + (int) stats[0]);
for (int i = 0; i < 10; i++) {
IJ.log("P" + (i + 1) + " = " + (float) params[i] + " fixed = " + fixes[i]);
}
IJ.log("AIC = " + (float) stats[2]);
// IJ.log("hits = "+hitcounter);
set_options(params, fixes);
}
void init(String[] names1, float[][][] pch1, int psfflag1) {
setLayout(null);
names = names1;
pch = pch1;
psfflag = psfflag1;
ncurves = pch.length;
nparams = 11;
xpts = pch[0].length;
ypts = pch[0][0].length;
checkarray = new Checkbox[ncurves];
include = new boolean[ncurves];
namearray = new TextField[ncurves + 1];
int1array = new TextField[ncurves + 1];
intensity1 = new double[ncurves + 1];
e1array = new TextField[ncurves + 1];
n1array = new TextField[ncurves + 1];
bright1 = new double[ncurves + 1];
number1 = new double[ncurves + 1];
int2array = new TextField[ncurves + 1];
intensity2 = new double[ncurves + 1];
e2array = new TextField[ncurves + 1];
n2array = new TextField[ncurves + 1];
bright2 = new double[ncurves + 1];
number2 = new double[ncurves + 1];
eccarray = new TextField[ncurves + 1];
brightcc = new double[ncurves + 1];
eminccarray = new TextField[ncurves + 1];
brightmincc = new double[ncurves + 1];
c2array = new TextField[ncurves + 1];
c2 = new double[ncurves + 1];
nmeas = new int[ncurves + 1];
avg = new float[xpts][ypts];
indices = new int[ncurves];
beta = 0.05;
getintbright();
for (int i = 0; i < ncurves; i++) {
include[i] = true;
indices[i] = i;
}
updateavg();
int starty = 60;
int startx = 10;
int yinc = 25;
for (int i = 0; i <= ncurves; i++) {
if (i != ncurves) {
checkarray[i] = new Checkbox("", include[i]);
checkarray[i].setBounds(startx, starty + i * yinc, 20, 20);
checkarray[i].addItemListener(this);
add(checkarray[i]);
}
namearray[i] = new TextField(names[i]);
namearray[i].setBounds(startx + 30, starty + i * yinc, 200, 20);
add(namearray[i]);
int1array[i] = new TextField("" + (float) intensity1[i]);
int1array[i].setBounds(startx + 30 + 210, starty + i * yinc, 40, 20);
add(int1array[i]);
e1array[i] = new TextField("" + (float) bright1[i]);
e1array[i].setBounds(startx + 30 + 210 + 50, starty + i * yinc, 40, 20);
add(e1array[i]);
n1array[i] = new TextField("" + (float) number1[i]);
n1array[i].setBounds(startx + 30 + 210 + 50 + 50, starty + i * yinc, 40, 20);
add(n1array[i]);
int2array[i] = new TextField("" + (float) intensity2[i]);
int2array[i].setBounds(startx + 30 + 210 + 50 + 50 + 50, starty + i * yinc, 40, 20);
add(int2array[i]);
e2array[i] = new TextField("" + (float) bright2[i]);
e2array[i].setBounds(startx + 30 + 210 + 50 + 50 + 50 + 50, starty + i * yinc, 40, 20);
add(e2array[i]);
n2array[i] = new TextField("" + (float) number2[i]);
n2array[i].setBounds(startx + 30 + 210 + 50 + 50 + 50 + 50 + 50, starty + i * yinc, 40, 20);
add(n2array[i]);
eccarray[i] = new TextField("" + (float) brightcc[i]);
eccarray[i].setBounds(
startx + 30 + 210 + 50 + 50 + 50 + 50 + 50 + 50, starty + i * yinc, 40, 20);
add(eccarray[i]);
eminccarray[i] = new TextField("" + (float) brightmincc[i]);
eminccarray[i].setBounds(
startx + 30 + 210 + 50 + 50 + 50 + 50 + 50 + 50 + 50, starty + i * yinc, 40, 20);
add(eminccarray[i]);
c2[i] = 0.0;
c2array[i] = new TextField("" + (float) c2[i]);
c2array[i].setBounds(
startx + 30 + 210 + 50 + 50 + 50 + 50 + 50 + 50 + 50 + 50, starty + i * yinc, 80, 20);
add(c2array[i]);
}
namelabel = new Label("Filename");
namelabel.setBounds(startx + 30, starty - 25, 100, 20);
add(namelabel);
intlabel = new Label("<Ig>");
intlabel.setBounds(startx + 30 + 210, starty - 25, 40, 20);
add(intlabel);
brightlabel = new Label("<eg>");
brightlabel.setBounds(startx + 30 + 210 + 50, starty - 25, 40, 20);
add(brightlabel);
nlabel = new Label("<Ng>");
nlabel.setBounds(startx + 30 + 210 + 50 + 50, starty - 25, 40, 20);
add(nlabel);
int2label = new Label("<Ir>");
int2label.setBounds(startx + 30 + 210 + 50 + 50 + 50, starty - 25, 40, 20);
add(int2label);
bright2label = new Label("<er>");
bright2label.setBounds(startx + 30 + 210 + 50 + 50 + 50 + 50, starty - 25, 40, 20);
add(bright2label);
n2label = new Label("<Nr>");
n2label.setBounds(startx + 30 + 210 + 50 + 50 + 50 + 50 + 50, starty - 25, 40, 20);
add(n2label);
brightcclabel = new Label("<ecc>");
brightcclabel.setBounds(startx + 30 + 210 + 50 + 50 + 50 + 50 + 50 + 50, starty - 25, 40, 20);
add(brightcclabel);
brightccminlabel = new Label("min");
brightccminlabel.setBounds(
startx + 30 + 210 + 50 + 50 + 50 + 50 + 50 + 50 + 50, starty - 25, 40, 20);
add(brightccminlabel);
c2label = new Label("chi^2");
c2label.setBounds(
startx + 30 + 210 + 50 + 50 + 50 + 50 + 50 + 50 + 50 + 50, starty - 25, 80, 20);
add(c2label);
int buttonsx = startx + 30 + 210 + 50 + 50 + 50 + 50 + 50 + 50 + 50 + 50 + 90;
fitavgbutton = new Button("Fit Avg");
fitavgbutton.setBounds(buttonsx, starty - 25, 100, 40);
fitavgbutton.addActionListener(this);
add(fitavgbutton);
fitglobalbutton = new Button("Fit Global");
fitglobalbutton.setBounds(buttonsx, starty - 25 + 50, 100, 40);
fitglobalbutton.addActionListener(this);
add(fitglobalbutton);
clearparamsbutton = new Button("Reset Fit Params");
clearparamsbutton.setBounds(buttonsx, starty - 25 + 50 + 50, 100, 40);
clearparamsbutton.addActionListener(this);
add(clearparamsbutton);
checkc2 = false;
fitclass = new NLLSfit(this, 0.0001, 50, 0.1);
globalfitclass = new NLLSglobalfit(this, 0.0001, 50, 0.1);
pchfunc = new pch2D((int) ((double) xpts * 1.5), (int) ((double) ypts * 1.5), psfflag);
avgfit = new float[xpts][ypts];
fit = new float[ncurves][xpts][ypts];
xvals = new float[ncurves][xpts][ypts];
for (int i = 0; i < ncurves; i++) {
for (int j = 0; j < xpts; j++) {
for (int k = 0; k < ypts; k++) {
xvals[i][j][k] = (float) k;
fit[i][j][k] = 1.0f;
}
}
}
globalc2label = new Label("Global chi^2 = " + (float) 0.0);
globalc2label.setBounds(buttonsx, starty - 25 + 50 + 50 + 50, 140, 20);
add(globalc2label);
dispcurvelabel = new Label("Display Fit #");
dispcurvelabel.setBounds(buttonsx, starty - 25 + 50 + 50 + 50 + 30, 70, 20);
add(dispcurvelabel);
dispcurvechoice = new Choice();
for (int i = 0; i < ncurves; i++) {
dispcurvechoice.add("" + (i + 1));
}
dispcurve = 0;
dispcurvechoice.select(0);
dispcurvechoice.setBounds(buttonsx + 80, starty - 25 + 50 + 50 + 50 + 30, 40, 20);
dispcurvechoice.addItemListener(this);
add(dispcurvechoice);
betalabel = new Label("Bleedthrough f");
betalabel.setBounds(buttonsx, starty - 25 + 50 + 50 + 50 + 30 + 30, 100, 20);
add(betalabel);
betaval = new TextField("" + (float) beta);
betaval.setBounds(buttonsx + 110, starty - 25 + 50 + 50 + 50 + 30 + 30, 40, 20);
betaval.addActionListener(this);
add(betaval);
beta = Double.parseDouble(betaval.getText());
updatebeta();
undobutton = new Button("Undo Global Fit");
undobutton.setBounds(buttonsx, starty - 25 + 50 + 50 + 50 + 30 + 30 + 50, 100, 40);
undobutton.addActionListener(this);
add(undobutton);
geterrorsbutton = new Button("Get Errors");
geterrorsbutton.setBounds(buttonsx, starty - 25 + 50 + 50 + 50 + 30 + 30 + 50 + 50, 100, 40);
geterrorsbutton.addActionListener(this);
add(geterrorsbutton);
copylabel = new Label("copyright 2009 Jay Unruh ([email protected]) non-profit use only");
copylabel.setBounds(10, 790, 400, 20);
add(copylabel);
n_b_label = new Label("N and B Analysis");
n_b_label.setBounds(250, 10, 100, 20);
add(n_b_label);
pwavg = new PlotWindow3D("Avg", "kg", "kr", "Frequency", avg, 0);
pwavg.setLogAxes(false, false, true);
pwavg.draw();
pwavg.addPoints(avgfit, true, 0);
float[] temp = pwavg.getLimits();
temp[4] = 1.0f;
pwavg.setLimits(temp);
float[][] temppch = new float[xpts][ypts];
for (int i = 0; i < xpts; i++) {
System.arraycopy(pch[dispcurve][i], 0, temppch[i], 0, ypts);
}
pwfit = new PlotWindow3D("Selected Curve", "kg", "kr", "Frequency", temppch, 0);
pwfit.setLogAxes(false, false, true);
pwfit.draw();
pwfit.addPoints(fit[dispcurve], true, 0);
float[] temp2 = pwfit.getLimits();
temp2[4] = 1.0f;
pwfit.setLimits(temp2);
resetparams();
repaint();
}
public void run(String arg) {
Frame[] niframes = WindowManager.getNonImageWindows();
String[] titles = new String[niframes.length + 1];
for (int i = 0; i < niframes.length; i++) {
titles[i] = niframes[i].getTitle();
}
titles[niframes.length] = "Clipboard";
GenericDialog gd = new GenericDialog("Windows");
boolean importfile = false;
gd.addCheckbox("Import from file?", importfile);
gd.addChoice("Windows", titles, titles[0]);
boolean hasxvals = false;
gd.addCheckbox("X Vals Column?", hasxvals);
boolean multix = false;
gd.addCheckbox("Multi_X_Columns?", multix);
boolean skipendzeros = false;
gd.addCheckbox("Skip_end_zeros?", skipendzeros);
String[] delimiters = {"Tab", "Comma", "Space"};
gd.addChoice("Delimiter", delimiters, delimiters[0]);
gd.showDialog();
if (gd.wasCanceled()) {
return;
}
importfile = gd.getNextBoolean();
int index = gd.getNextChoiceIndex();
hasxvals = gd.getNextBoolean();
multix = gd.getNextBoolean();
skipendzeros = gd.getNextBoolean();
int delimindex = gd.getNextChoiceIndex();
if (multix) hasxvals = true;
String textdata = "";
if (importfile) {
OpenDialog od = new OpenDialog("Open File", "", ".txt");
String directory = od.getDirectory();
String name = od.getFileName();
if (name == null) {
return;
}
try {
File infile = new File(directory + name);
BufferedReader b = new BufferedReader(new FileReader(infile));
textdata = (new jdataio()).readstringfile(b);
b.close();
} catch (IOException e) {
return;
}
} else {
if (index == niframes.length) {
// here we get the data from the clipboard
Transferable t = Toolkit.getDefaultToolkit().getSystemClipboard().getContents(null);
try {
if (t != null && t.isDataFlavorSupported(DataFlavor.stringFlavor)) {
textdata = (String) t.getTransferData(DataFlavor.stringFlavor);
}
} catch (UnsupportedFlavorException e) {
} catch (IOException e) {
}
if (textdata.equals("")) {
IJ.error("Error copying from clipboard.");
return;
}
} else {
if (niframes[index] instanceof Editor) {
Editor tw = (Editor) niframes[index];
textdata = tw.getText();
} else {
if (niframes[index] instanceof TextWindow) {
TextWindow tw = (TextWindow) niframes[index];
textdata = tw.getTextPanel().getText();
} else {
IJ.showMessage("Not a valid text window");
return;
}
}
}
}
if (textdata == null) {
IJ.showMessage("Error in Obtaining String");
return;
}
if (textdata.indexOf("\r") >= 0) {
textdata = textdata.replace('\r', '\n');
}
char[] delims = {'\t', ',', ' '};
delimit_string ds = new delimit_string(delims[delimindex]);
String[] rows = ds.getrows(textdata);
int lines = rows.length;
int columns = ds.getnumcolumns(rows[0]);
int ycolumns = columns;
if (hasxvals) {
if (multix) {
ycolumns /= 2;
} else {
ycolumns--;
}
}
if (multix) {
float[][] ydata = new float[ycolumns][lines];
float[][] xdata = new float[ycolumns][lines];
for (int i = 0; i < lines; i++) {
float[] temp = ds.delim2float(rows[i], columns);
for (int j = 0; j < ycolumns; j++) {
ydata[j][i] = temp[2 * j + 1];
xdata[j][i] = temp[2 * j];
}
}
int[] npts = new int[ycolumns];
for (int i = 0; i < ycolumns; i++) {
npts[i] = lines;
}
if (skipendzeros) {
for (int i = 0; i < ycolumns; i++) {
int counter = lines - 1;
while ((xdata[i][counter] == 0.0f || Float.isNaN(xdata[i][counter])) && counter > 0) {
xdata[i][counter] = 0.0f;
ydata[i][counter] = 0.0f;
npts[i]--;
counter--;
}
}
}
(new PlotWindow4("Text Plot", "x", "y", xdata, ydata, npts)).draw();
} else {
float[][] tempydata = new float[ycolumns][lines];
float[] tempxdata = new float[lines];
float[][] xdata = null;
float[][] ydata = null;
int startcolumn = 0;
if (hasxvals) startcolumn = 1;
for (int i = 0; i < lines; i++) {
float[] temp = ds.delim2float(rows[i], columns);
if (hasxvals) {
tempxdata[i] = temp[0];
} else {
tempxdata[i] = (float) (i + 1);
}
for (int j = 0; j < ycolumns; j++) {
tempydata[j][i] = temp[j + startcolumn];
}
}
int[] npts = new int[ycolumns];
npts[0] = lines;
if (skipendzeros) {
int maxpts = 0;
for (int i = 0; i < ycolumns; i++) {
int counter = lines - 1;
npts[i] = lines;
while ((tempydata[i][counter] == 0.0f || Float.isNaN(tempydata[i][counter]))
&& counter > 0) {
npts[i]--;
counter--;
}
if (npts[i] > maxpts) maxpts = npts[i];
IJ.log("" + npts[i]);
}
ydata = new float[ycolumns][maxpts];
xdata = new float[ycolumns][maxpts];
for (int i = 0; i < ycolumns; i++) {
// npts[i]=npts[0];
System.arraycopy(tempxdata, 0, xdata[i], 0, npts[i]);
System.arraycopy(tempydata[i], 0, ydata[i], 0, npts[i]);
}
} else {
ydata = tempydata;
xdata = new float[ycolumns][];
for (int i = 0; i < ycolumns; i++) {
npts[i] = npts[0];
xdata[i] = tempxdata.clone();
}
}
(new PlotWindow4("Text Plot", "x", "y", xdata, ydata, npts)).draw();
}
}