rickets

Well, damn.

I've been very inactive as of late. Mostly I've been working in a location in lab that is far away from my "SG computer," which makes it less and less convenient to pop over for a few posts. This week, I'm back at my normal computer, but I'm busy with things like,

"corrdata(:,n:n+1)=inf;

newfile=load(['c:\ascii\' input('Input IV data file (in single quotes)>')]); %load in IV data

[r,c]=size(newfile);
corrdata(1:r,n:n+1)=newfile;
leaksub=polyval(fitleak,corrdata(1:r,n));
corrdata(1:r,n+1)=corrdata(1:r,n+1)-leaksub;

volt=[-100:10:10]';
[row,col]=size(volt);
num=r/row %figures out N value for each potential, assuming the N is the same for each.

a=0;
while corrdata((1+num*a),(n+1))~=inf
transpose((b+1),1:num)=(corrdata((1+num*a)num*(a+1)),n+1))';%transpose data

a=a+1;
b=b+1;%row location in transpose
end

transpose2=transpose';
means(:,m)=mean(((transpose(((num*c)+1)num*(c+1)),1:num))),2);
serr=(std(transpose2(1:num,((num*c)+1)(num*(c+1))./(sqrt(num))))))';
sem(:,m)=serr;

if m==1
symbol='o';
elseif m==2
symbol='d';
elseif m==3
symbol='^';
elseif m==4
symbol='s';
elseif m==5
symbol='*';
elseif m==6
symbol='p'
else
disp('Error: cannot exceed 6 concentrations')
end
figure
plot(volt,means,symbol)
hold
errorbar(volt,means(:,m),sem(:,m))."

This puts me in a state of mind, where I can think of little or nothing to say that would be of any interest to anybody but me and (maybe) my boss.

Hope you're all well. I'll do my best to stop by everybody's journals soon.

Mike

PS I just posted this and realized that some of my code was replaced with frowning faces. Very appropriate.