[PX] Scala 最強伝説

「～最強伝説」というと、Opera だったり、黒魔だったり赤魔だったりするわけですが、いずれにしてもその内実は、嫉妬・羨望・侮蔑・自虐が綯い交ぜになったものであるのが常です。

…ま、意味はよくわかりませんが、XDS の出力した XDS_ASCII.HKL を、Scala で処理してみましたよ、ということです。

波長 1.5 Å での Thaumatin 結晶の回折データを XDS でINTEGRATE・CORRECTし、その XDS_ASCII.HKL を、前回エントリで書いた「pointless」を使った方法でマルチレコードMTZにしました。

それを Scala の「on rotation axis with secondary beam absorption」でスケーリングした統計値が以下になります。

  N 1/d^2 Dmin(A) Rmrg  Rfull   Rcum  Ranom  Nanom    Av_I  SIGMA I/sigma   sd Mn(I/sd)  Nmeas   Nref  Ncent FRCBIAS  Nbias
  1 0.0399  5.01  0.035  0.035  0.035  0.015    802     1915    117  16.4     92  77.0    31741    1313   509    -         0
  2 0.0797  3.54  0.032  0.032  0.033  0.010   1665     3843    203  18.9    177  91.5    58940    2184   519    -         0
  3 0.1196  2.89  0.042  0.042  0.036  0.011   2221     2349    148  15.8    132  75.7    75281    2733   512    -         0
  4 0.1594  2.50  0.066  0.066  0.042  0.016   2672     1205    111  10.8    100  54.4    88125    3194   522    -         0
  5 0.1993  2.24  0.085  0.085  0.048  0.019   3083      966    113   8.5    102  44.7    97705    3597   514    -         0
  6 0.2391  2.05  0.108  0.108  0.055  0.023   3423      791    118   6.7    107  36.3   104587    3934   511    -         0
  7 0.2790  1.89  0.161  0.161  0.062  0.032   3741      538    118   4.6    111  25.9   109812    4259   519    -         0
  8 0.3188  1.77  0.268  0.268  0.071  0.051   4041      319    118   2.7    116  16.1   113646    4551   510    -         0
  9 0.3587  1.67  0.447  0.447  0.082  0.088   4316      209    130   1.6    133   9.9   116315    4824   513    -         0
 10 0.3985  1.58  0.620  0.620  0.089  0.154   4058      145    123   1.2    130   5.6    77988    4458   427    -         0

以前、「XDS + XSCALE」および「iMosflm + Scala」で処理したときの統計値をアップしていましたが、数値で見ても比較しにくいので、グラフにしてみましょう。

びッみょう～～に、「XDS + Scala」が一番低い R_merge を出してるような気がしますね…

…最大分解能変えて、shelxC/D/E をかけてみますか…

最大分解能	CCallの最頻値	BestCCall	Contrast(ori)	Contrast(inv)	Maxσ(mode)
2.0	10.500	12.92	0.402	0.376	2.671
2.1	11.522	13.82	0.391	0.388	1.976
2.2	13.675	20.13	0.735	0.440	5.068
2.4	15.541	27.25	0.737	0.435	8.450
2.6	17.500	30.92	0.736	0.444	10.380
2.8	20.625	35.13	0.737	0.444	10.442
3.0	23.615	40.44	0.738	0.443	10.511
3.2	26.559	31.96	0.423	0.413	2.922
3.4	29.463	40.39	0.618	0.418	6.601
3.6	31.625	47.38	0.743	0.436	8.420
3.8	34.524	38.87	0.408	0.357	1.929
4.0	36.457	49.19	0.518	0.417	6.022
4.2	40.476	51.53	0.504	0.414	3.786
4.4	43.304	49.28	0.452	0.415	1.913

と、こんな感じでした。

高分解能の限界は、2.2 Åと上がってはいるんですが、低分解能側の限界が 3.6 ～4.0 Åとちょっとあやふやです。
Try数の多寡の問題があるので、一概に「 XDS + Scala 最強！」とは言えませんが、そこがまた、いかにも「最強伝説」らしいっちゃ、らしいところですねｗ

---

永らく続いた [PX] も、とりあえずネタが切れたんで、これで一時休止ですかね。