Here’s a slightly simpler version of delMissRows() that works just as well:

proc delMissRows(x);
 local mask;

 mask = sumr(x .== miss(0,0)); // number of missing entries per row
 mask = mask .< cols(x); // 1 if a row has some non-missing entries
 retp(selif(x,mask)); // returns only those rows where mask .== 1
endp;

As a workaround here is a procedure to remove all rows that only contain missing values.

new;
struct plotControl myPlot;
myplot = plotGetDefaults("xy");
load rawdataF;
PlotSetLegend(&myplot,"off");
x = seqa(2009-0.25*cols(rawdataF),.25,cols(rawdataF));
plotXY(myPlot, x, delMissRows(rawdataF)');


proc delMissRows(x);
	local mask, idx, xout, rtrim, xctr;
	
	//Which rows have as many missings as columns
	mask = sumr(x .== miss(0,0));
	
	//Preallocate output matrix
	rtrim = sumc(mask .== cols(x));
	xout = zeros(rows(x) - rtrim, cols(x));
	
	//Fill in xout
	xctr = 1;
	for i(1, rows(x), 1);
		if mask[i] != cols(x);
			xout[xctr, .] = x[i,.];
			xctr = xctr + 1;
		endif;
	endfor;
	
	retp(xout);	
endp;

Thanks for your answer. Your example displays perfectly for me.

I’ve tried to create a single example that reproduces the problem, using my original data saved in .fmt format as you suggested. (I can’t figure out how to attach the .fmt file to this forum post, so I’ll send that separately by e-mail.) Here’s a program that should reproduce the problem.

new;
struct plotControl myPlot;
myplot = plotGetDefaults("xy");
load rawdataF;
PlotSetLegend(&myplot,"off");
plotxy(myplot,seqa(2009-0.25*cols(rawdataF),.25,cols(rawdataF)),rawdataF');

This gives me a plot whose horizontal axis starts at 0, despite the fact that rawdataF is a 157×170 matrix.

Thanks in advance for your help.

Simon

struct plotControl myPlot;
myPlot = plotGetDefaults("xy");

x = seqa(1960, 0.25, 53*4);
y = cos(seqa(0.1, 0.1, rows(x)));
plotxy(myPlot, x,y);

Does this display correctly for you? Could you save your x and y variables as .fmt files with the save command and send them to support.