SAXS

Small Angle X-ray Scattering Initiative for Europe :: Forum
It is currently 2018.05.26 23:28

All times are UTC + 2 hours




Post new topic Reply to topic  [ 9 posts ] 
Author Message
PostPosted: 2016.11.18 18:39 
Offline
Member

Joined: 2016.11.18 18:26
Posts: 6
In order to compare theoretical scattering from xtal structure with experimental data from solution I used CRYSOL. Although I received .fit output file with two curves, the chi^2 value is not reported. What is the way to have it printed? The log file in attachment.

thank you in advance,


Attachments:
dimer_model00.log [3 KiB]
Downloaded 52 times
Top
 Profile  
 
 Post subject: chi-square is too large
PostPosted: 2016.11.18 19:20 
Offline
Administrator
User avatar

Joined: 2007.08.03 18:55
Posts: 629
Location: EMBL Hamburg, Germany
You may check 'crysol_summary.txt' and the header of the .fit file but basically your chi-square is way too large to make sense, this is why you see '*******' instead of a number.
You can recalculate chi-square of the fit using datcmp:
Code:
$ datcmp --test=CHI-SQUARE dimer_model00.fit


Top
 Profile  
 
PostPosted: 2016.11.22 13:49 
Offline
Member

Joined: 2016.11.18 18:26
Posts: 6
Thank you for the answer. Do you know what is the range of Chi^2 "too high to be displayed"?. When I use datcmp the calculated Chi^2 is -1, what is a ridiculous value.

I am aware that the fit is far from perfect, because xtal structure is lacking 52 residues at N-terminus (out of 133 contributing to scattering in solution); still I would like to have the numerical value to show the improvement after including the missing part in the model.


Attachments:
dimer_model00.fit [36.77 KiB]
Downloaded 44 times
Top
 Profile  
 
PostPosted: 2016.11.22 22:15 
Offline
Administrator

Joined: 2007.08.10 11:09
Posts: 406
aleful wrote:
Thank you for the answer. Do you know what is the range of Chi^2 "too high to be displayed"?. When I use datcmp the calculated Chi^2 is -1, what is a ridiculous value.

Not ridiculous. Lemme guess, no error estimates? Then that's what you get. Get accurate errors, and datcmp gives you an accurate chi^2.

With respect to the range: depending on the number of points of experimental data, the chi^2 value should be in the rule-of-thumb range of 0.9 to 1.1 (assuming ~1500 experimental points, a slightly wider range if less, more narrow if more). Of course, this assumes correct errors. Without correct errors, all other bets are off, too. And with accurate errors "*****" would indicate, roughly translated, "too embarrassingly high to even show the exact numerical value as it has no further meaning here", without errors it may mean anything. Probably division-by-zero. Or Not-A-Number. Or maybe the winning numbers of next weeks lottery.


Top
 Profile  
 
PostPosted: 2016.11.23 19:10 
Offline
Member

Joined: 2016.11.18 18:26
Posts: 6
Thank you for your answer. By saying ridiculous I meant I expect any value > 0. Especially that I have error estimates in my .dat file (column entitled sig(q), if I understand correctly). However, I cannot evaluate how accurate they are, but they are always at least one order of magnitude smaller than the corresponding intensity. The .dat file in the attachmnent.


Attachments:
dimer_edited.dat [194.97 KiB]
Downloaded 39 times
Top
 Profile  
 
PostPosted: 2016.11.24 10:01 
Offline
Active member

Joined: 2015.03.25 09:03
Posts: 86
Location: EMBL Hamburg
What version of ATSAS are you using? A recent version should copy the sig(q) column into the .fit file so that it can be used to calculate the chi^2 - but your .fit file has no sigma column.


Last edited by ckerr on 2016.11.24 10:02, edited 1 time in total.
ATSAS version not CRYSOL version


Top
 Profile  
 
 Post subject: ATSAS version?
PostPosted: 2016.11.24 10:19 
Offline
Administrator
User avatar

Joined: 2007.08.03 18:55
Posts: 629
Location: EMBL Hamburg, Germany
Which ATSAS version do you have? Please check what you get if you execute
Code:
crysol --version
Starting from ATSAS version 2.7.0 CRYSOL does not skip the experimental errors from the fit file by default.


Top
 Profile  
 
PostPosted: 2016.11.24 16:37 
Offline
Member

Joined: 2016.11.18 18:26
Posts: 6
crysol, ATSAS 2.7.1 (r6669)

it was downloaded and installed not more than a month ago


Top
 Profile  
 
PostPosted: 2016.11.24 17:52 
Offline
Active member

Joined: 2015.03.25 09:03
Posts: 86
Location: EMBL Hamburg
The latest version is 2.7.2 so I suggest that you upgrade. Then try running crysol and datcmp again and tell us if the problem still appears.

UPD: ATSAS 2.8.0 is out now.


Top
 Profile  
 
Display posts from previous:  Sort by  
Post new topic Reply to topic  [ 9 posts ] 

All times are UTC + 2 hours


Who is online

Users browsing this forum: No registered users and 1 guest


You cannot post new topics in this forum
You cannot reply to topics in this forum
You cannot edit your posts in this forum
You cannot delete your posts in this forum
You cannot post attachments in this forum

Search for:
Jump to:  
Powered by phpBB® Forum Software © phpBB Group