##-------------------------------------------## ## WARNING: ## ## Number of residues unspecified ## ##-------------------------------------------## Effective parameters: #phil __ON__ scaling { input { asu_contents { sequence_file = None n_residues = None n_bases = None n_copies_per_asu = 2 } xray_data { file_name = "/net/chevy/raid1/afonine/work/crom/final/fig_22_bad_lig_1SE6/autobuild/AutoBuild_run_1_/refinement_PHX.mtz" obs_labels = FOBS calc_labels = None unit_cell = 59.554 79.206 87.43 90 92.26 90 space_group = "P 1 21 1" high_resolution = 1.6096 low_resolution = None reference { data { file_name = None labels = None unit_cell = None space_group = None } structure { file_name = None } } } parameters { reporting { verbose = 1 log = "/net/chevy/raid1/afonine/work/crom/final/fig_22_bad_lig_1SE6/autobuild/AutoBuild_run_1_/refinement_PHX.mtz_xtriage.log" loggraphs = False } merging { n_bins = 10 skip_merging = False } misc_twin_parameters { missing_symmetry { sigma_inflation = 1.25 } twinning_with_ncs { perform_analyses = False n_bins = 7 } twin_test_cuts { low_resolution = 10 high_resolution = None isigi_cut = 3 completeness_cut = 0.85 } apply_basic_filters_prior_to_twin_analysis = True } } optional { hklout = None hklout_type = mtz sca *Auto label_extension = "massaged" aniso { action = *remove_aniso None final_b = *eigen_min eigen_mean user_b_iso b_iso = None } outlier { action = *extreme basic beamstop None parameters { basic_wilson { level = 1e-06 } extreme_wilson { level = 0.01 } beamstop { level = 0.001 d_min = 10 } } } symmetry { action = detwin twin *None twinning_parameters { twin_law = None fraction = None } } } } gui { result_file = None output_dir = None job_title = None } } #phil __END__ Symmetry, cell and reflection file content summary Miller array info: /net/chevy/raid1/afonine/work/crom/final/fig_22_bad_lig_1SE6/autobuild/AutoBuild_run_1_/refinement_PHX.mtz:FOBS,SIGFOBS Observation type: xray.amplitude Type of data: double, size=84260 Type of sigmas: double, size=84260 Number of Miller indices: 84260 Anomalous flag: False Unit cell: (59.554, 79.206, 87.43, 90, 92.26, 90) Space group: P 1 21 1 (No. 4) Systematic absences: 0 Centric reflections: 2332 Resolution range: 35.8944 1.6096 Completeness in resolution range: 0.802293 Completeness with d_max=infinity: 0.802201 Wavelength: 1.0000 ################################################################################ # Basic statistics # ################################################################################ =================== Solvent content and Matthews coefficient ================== Number of residues unknown, assuming 50% solvent content Best guess : 754 residues in the ASU Caution: this estimate is based on the distribution of solvent content across structures in the PDB, but it does not take into account the resolution of the data (which is strongly correlated with solvent content) or the physical properties of the model (such as oligomerization state, et cetera). If you encounter problems with molecular replacement and/or refinement, you may need to consider the possibility that the ASU contents are different than expected. Number of copies per asymmetric unit provided Will use user specified value of 2.0 ======================== Data strength and completeness ======================= ----------Completeness at I/sigma cutoffs---------- The following table lists the completeness in various resolution ranges, after applying a I/sigI cut. Miller indices for which individual I/sigI values are larger than the value specified in the top row of the table, are retained, while other intensities are discarded. The resulting completeness profiles are an indication of the strength of the data. **********WARNING********** Please be aware that the input data were given as amplitudes and squared for the purposes of this analysis, therefore the numbers displayed here are less reliable than the values calculated from the raw intensities. ---------------------------------------------------------------------------------------- | Completeness and data strength | |--------------------------------------------------------------------------------------| | Res. range | I/sigI>1 | I/sigI>2 | I/sigI>3 | I/sigI>5 | I/sigI>10 | I/sigI>15 | |--------------------------------------------------------------------------------------| | 35.90 - 3.97 | 96.8 | 96.7 | 96.7 | 96.4 | 95.6 | 94.5 | | 3.97 - 3.15 | 77.8 | 77.7 | 77.5 | 77.1 | 75.8 | 74.1 | | 3.15 - 2.75 | 99.4 | 99.1 | 98.6 | 97.1 | 92.9 | 89.1 | | 2.75 - 2.50 | 99.0 | 98.2 | 97.1 | 95.1 | 89.2 | 83.1 | | 2.50 - 2.32 | 98.5 | 97.1 | 95.7 | 92.4 | 84.0 | 76.2 | | 2.32 - 2.18 | 54.8 | 53.8 | 52.7 | 50.5 | 45.1 | 38.7 | | 2.18 - 2.08 | 96.8 | 94.5 | 91.4 | 86.0 | 72.5 | 60.3 | | 2.08 - 1.98 | 94.8 | 91.3 | 87.5 | 79.8 | 62.8 | 48.2 | | 1.98 - 1.91 | 73.6 | 70.5 | 67.1 | 60.7 | 44.8 | 30.1 | | 1.91 - 1.84 | 67.4 | 60.9 | 55.4 | 45.1 | 26.9 | 15.7 | | 1.84 - 1.78 | 77.5 | 68.7 | 59.9 | 44.5 | 22.2 | 11.2 | | 1.78 - 1.73 | 66.4 | 57.4 | 48.7 | 34.8 | 14.4 | 6.3 | | 1.73 - 1.69 | 54.1 | 45.6 | 37.6 | 25.3 | 9.0 | 3.5 | | 1.69 - 1.65 | 43.5 | 36.0 | 28.6 | 18.2 | 5.7 | 1.9 | ---------------------------------------------------------------------------------------- The completeness of data for which I/sig(I)>3.00, exceeds 85 % for resolution ranges lower than 1.98A. The data are cut at this resolution for the potential twin tests and intensity statistics. ----------Low resolution completeness analyses---------- The following table shows the completeness of the data to 5.0 A. Poor low-resolution completeness often leads to map distortions and other difficulties, and is typically caused by problems with the crystal orientation during data collection, overexposure of frames, interference with the beamstop, or omission of reflections by data-processing software. --------------------------------------------------------- | Resolution range | N(obs)/N(possible) | Completeness | --------------------------------------------------------- | 35.8953 - 10.6883 | [305/378] | 0.807 | | 10.6883 - 8.5210 | [318/361] | 0.881 | | 8.5210 - 7.4548 | [357/361] | 0.989 | | 7.4548 - 6.7782 | [362/362] | 1.000 | | 6.7782 - 6.2952 | [361/361] | 1.000 | | 6.2952 - 5.9258 | [355/355] | 1.000 | | 5.9258 - 5.6302 | [361/361] | 1.000 | | 5.6302 - 5.3859 | [350/351] | 0.997 | | 5.3859 - 5.1792 | [351/351] | 1.000 | | 5.1792 - 5.0010 | [355/356] | 0.997 | --------------------------------------------------------- ----------Completeness (log-binning)---------- The table below presents an alternative overview of data completeness, using the entire resolution range but on a logarithmic scale. This is more sensitive to missing low-resolution data (and is complementary to the separate table showing low-resolution completeness only). -------------------------------------------------- | Resolution | Reflections | Completeness | -------------------------------------------------- | 35.8944 - 14.8137 | 101/135 | 74.8% | | 14.7624 - 12.2142 | 99/116 | 85.3% | | 12.1664 - 9.9933 | 168/203 | 82.8% | | 9.9852 - 8.2092 | 348/376 | 92.6% | | 8.2041 - 6.7403 | 649/653 | 99.4% | | 6.7356 - 5.5371 | 1187/1187 | 100.0% | | 5.5336 - 4.5462 | 2096/2110 | 99.3% | | 4.5457 - 3.7343 | 3043/3790 | 80.3% | | 3.7333 - 3.0664 | 5920/6820 | 86.8% | | 3.0664 - 2.5185 | 12226/12260 | 99.7% | | 2.5184 - 2.0683 | 18869/22049 | 85.6% | | 2.0683 - 1.6096 | 39554/55317 | 71.5% | -------------------------------------------------- ----------Analysis of resolution limits---------- Your data have been examined to determine the resolution limits of the data along the reciprocal space axes (a*, b*, and c*). These are expected to vary slightly depending on unit cell parameters and overall resolution, but should never be significantly different for complete data. (This is distinct from the amount of anisotropy present in the data, which changes the effective resolution but does not actually exclude reflections.) overall d_min = 1.610 d_min along a* = 1.653 d_min along b* = 1.650 d_min along c* = 1.747 max. difference between axes = 0.097 Resolution limits are within expected tolerances. ===================== Absolute scaling and Wilson analysis ==================== ----------Maximum likelihood isotropic Wilson scaling---------- ML estimate of overall B value of /net/chevy/raid1/afonine/work/crom/final/fig_22_bad_lig_1SE6/autobuild/AutoBuild_run_1_/refinement_PHX.mtz:FOBS,SIGFOBS: 11.06 A**(-2) Estimated -log of scale factor of /net/chevy/raid1/afonine/work/crom/final/fig_22_bad_lig_1SE6/autobuild/AutoBuild_run_1_/refinement_PHX.mtz:FOBS,SIGFOBS: -0.70 The overall B value ("Wilson B-factor", derived from the Wilson plot) gives an isotropic approximation for the falloff of intensity as a function of resolution. Note that this approximation may be misleading for anisotropic data (where the crystal is poorly ordered along an axis). The Wilson B is strongly correlated with refined atomic B-factors but these may differ by a significant amount, especially if anisotropy is present. ----------Maximum likelihood anisotropic Wilson scaling---------- ML estimate of overall B_cart value: 7.42, 0.00, -0.70 14.12, 0.00 13.98 Equivalent representation as U_cif: 0.09, -0.00, -0.00 0.18, 0.00 0.18 Eigen analyses of B-cart: ------------------------------------------------ | Eigenvector | Value | Vector | ------------------------------------------------ | 1 | 14.125 | (0.00, 1.00, 0.00) | | 2 | 14.053 | (-0.11, 0.00, 0.99) | | 3 | 7.346 | ( 0.99, 0.00, 0.11) | ------------------------------------------------ ML estimate of -log of scale factor: -0.68 ----------Anisotropy analyses---------- For the resolution shell spanning between 1.64 - 1.61 Angstrom, the mean I/sigI is equal to 4.07. 49.2 % of these intensities have an I/sigI > 3. When sorting these intensities by their anisotropic correction factor and analysing the I/sigI behavior for this ordered list, we can gauge the presence of 'anisotropy induced noise amplification' in reciprocal space. The quarter of Intensities *least* affected by the anisotropy correction show : 4.58e+00 Fraction of I/sigI > 3 : 5.53e-01 ( Z = 2.74 ) The quarter of Intensities *most* affected by the anisotropy correction show : 2.55e+00 Fraction of I/sigI > 3 : 3.08e-01 ( Z = 8.23 ) Z-scores are computed on the basis of a Bernoulli model assuming independence of weak reflections with respect to anisotropy. ----------Wilson plot---------- The Wilson plot shows the falloff in intensity as a function in resolution; this is used to calculate the overall B-factor ("Wilson B-factor") for the data shown above. The expected plot is calculated based on analysis of macromolecule structures in the PDB, and the distinctive appearance is due to the non-random arrangement of atoms in the crystal. Some variation is natural, but major deviations from the expected plot may indicate pathological data (including ice rings, detector problems, or processing errors). ----------Mean intensity analyses---------- Inspired by: Morris et al. (2004). J. Synch. Rad.11, 56-59. The following resolution shells are worrisome: ----------------------------------------------------------------- | Mean intensity by shell (outliers) | |---------------------------------------------------------------| | d_spacing | z_score | completeness | / | |---------------------------------------------------------------| | 2.247 | 14.32 | 0.34 | 3.531 | | 2.033 | 6.51 | 0.99 | 0.823 | | 2.016 | 5.46 | 0.99 | 0.848 | | 2.000 | 6.28 | 0.98 | 0.825 | | 1.984 | 7.04 | 0.97 | 0.808 | | 1.969 | 6.77 | 0.96 | 0.810 | | 1.954 | 8.70 | 0.94 | 0.756 | | 1.924 | 12.13 | 0.42 | 2.082 | | 1.910 | 15.91 | 0.45 | 2.812 | | 1.897 | 4.59 | 0.40 | 0.817 | | 1.883 | 7.46 | 0.92 | 0.782 | | 1.870 | 8.30 | 0.94 | 0.767 | | 1.857 | 7.15 | 0.93 | 0.803 | | 1.844 | 7.85 | 0.89 | 0.786 | | 1.832 | 8.54 | 0.90 | 0.775 | | 1.820 | 4.99 | 0.87 | 0.856 | | 1.808 | 5.15 | 0.86 | 0.850 | ----------------------------------------------------------------- Possible reasons for the presence of the reported unexpected low or elevated mean intensity in a given resolution bin are : - missing overloaded or weak reflections - suboptimal data processing - satellite (ice) crystals - NCS - translational pseudo symmetry (detected elsewhere) - outliers (detected elsewhere) - ice rings (detected elsewhere) - other problems Note that the presence of abnormalities in a certain region of reciprocal space might confuse the data validation algorithm throughout a large region of reciprocal space, even though the data are acceptable in those areas. ----------Possible outliers---------- Inspired by: Read, Acta Cryst. (1999). D55, 1759-1764 Acentric reflections: ----------------------------------------------------------------------------------------------------- | Acentric reflections | |---------------------------------------------------------------------------------------------------| | d_spacing | H K L | |E| | p(wilson) | p(extreme) | |---------------------------------------------------------------------------------------------------| | 3.057 | 19, 5, 2 | 3.75 | 7.82e-07 | 6.19e-02 | ----------------------------------------------------------------------------------------------------- p(wilson) : 1-(1-exp[-|E|^2]) p(extreme) : 1-(1-exp[-|E|^2])^(n_acentrics) p(wilson) is the probability that an E-value of the specified value would be observed if it were selected at random the given data set. p(extreme) is the probability that the largest |E| value is larger or equal than the observed largest |E| value. Both measures can be used for outlier detection. p(extreme) takes into account the size of the dataset. Centric reflections: ----------------------------------------------------------------------------------------------------- | Centric reflections | |---------------------------------------------------------------------------------------------------| | d_spacing | H K L | |E| | p(wilson) | p(extreme) | |---------------------------------------------------------------------------------------------------| | 2.053 | -29, 0, 3 | 4.16 | 3.12e-05 | 6.91e-02 | ----------------------------------------------------------------------------------------------------- p(wilson) : 1-(erf[|E|/sqrt(2)]) p(extreme) : 1-(erf[|E|/sqrt(2)])^(n_acentrics) p(wilson) is the probability that an E-value of the specified value would be observed when it would selected at random from the given data set. p(extreme) is the probability that the largest |E| value is larger or equal than the observed largest |E| value. Both measures can be used for outlier detection. p(extreme) takes into account the size of the dataset. ----------Ice ring related problems---------- The following statistics were obtained from ice-ring insensitive resolution ranges: mean bin z_score : 2.33 ( rms deviation : 2.02 ) mean bin completeness : 0.88 ( rms deviation : 0.19 ) The following table shows the Wilson plot Z-scores and completeness for observed data in ice-ring sensitive areas. The expected relative intensity is the theoretical intensity of crystalline ice at the given resolution. Large z-scores and high completeness in these resolution ranges might be a reason to re-assess your data processsing if ice rings were present. ------------------------------------------------------------- | d_spacing | Expected rel. I | Data Z-score | Completeness | ------------------------------------------------------------- | 3.897 | 1.000 | 0.85 | 0.29 | | 3.669 | 0.750 | 1.91 | 0.46 | | 3.441 | 0.530 | 2.47 | 1.00 | | 2.671 | 0.170 | 1.95 | 0.99 | | 2.249 | 0.390 | 14.32 | 0.34 | | 2.072 | 0.300 | 3.22 | 0.99 | | 1.948 | 0.040 | 4.27 | 0.74 | | 1.918 | 0.180 | 15.91 | 0.45 | | 1.883 | 0.030 | 8.30 | 0.94 | | 1.721 | 0.020 | 1.21 | 0.67 | ------------------------------------------------------------- Abnormalities in mean intensity or completeness at resolution ranges with a relative ice ring intensity lower than 0.10 will be ignored. At 2.25 A the z-score is more than 4.00 times the standard deviation of all z-scores, while at the same time, the occupancy does not go down. At 1.92 A the z-score is more than 4.00 times the standard deviation of all z-scores, while at the same time, the occupancy does not go down. There were 2 ice ring related warnings. This could indicate the presence of ice rings. ################################################################################ # Twinning and symmetry analyses # ################################################################################ ============================= Systematic absences ============================= ----------Table of systematic absence rules---------- The following table gives information about systematic absences allowed for the specified intensity point group. For each operator, the reflections are split in three classes: Absent : Reflections that are absent for this operator. Non Absent: Reflection of the same type (i.e. (0,0,l)) as above, but they should be present. Complement: All other reflections. For each class, the is reported, as well as the number of 'violations'. A 'violation' is designated as a reflection for which a I/sigI criterion is not met. The criteria are: Absent violation : I/sigI > 3.0 Non Absent violation : I/sigI < 3.0 Complement violation : I/sigI < 3.0 Operators with low associated violations for *both* absent and non absent reflections, are likely to be true screw axis or glide planes. Both the number of violations and their percentages are given. The number of violations within the 'complement' class, can be used as a comparison for the number of violations in the non-absent class. ---------------------------------------------------------------------------------------------------------------------------------------- | Operator | # absent | (violations) | # not absent | (violations) | # complete | (violations) | Score | ---------------------------------------------------------------------------------------------------------------------------------------- | 2_0 (b) | 0 | 0.00 (0, 0.0%) | 13 | 28.08 (1, 7.7%) | 50724 | 31.48 (1930, 3.8%) | 1.49e+00 | | 2_1 (b) | 0 | 0.00 (0, 0.0%) | 13 | 28.08 (1, 7.7%) | 50724 | 31.48 (1930, 3.8%) | 1.49e+00 | ---------------------------------------------------------------------------------------------------------------------------------------- ----------Space group identification---------- Analyses of the absences table indicates a number of likely space group candidates, which are listed below. For each space group, the number of absent violations are listed under the '+++' column. The number of present violations (weak reflections) are listed under '---'. The last column is a likelihood based score for the particular space group. Note that enantiomorphic spacegroups will have equal scores. Also, if absences were removed while processing the data, they will be regarded as missing information, rather then as enforcing that absence in the space group choices. ----------------------------------------------------------------------------------- | space group | n absent | _absent | _absent | +++ | --- | score | ----------------------------------------------------------------------------------- | P 1 2 1 | 0 | 0.00 | 0.00 | 0 | 1 | 0.000e+00 | | P 1 21 1 | 0 | 0.00 | 0.00 | 0 | 1 | 0.000e+00 | ----------------------------------------------------------------------------------- ----------List of individual systematic absences---------- Note: this analysis uses the original input data rather than the filtered data used for twinning detection; therefore, the results shown here may include more reflections than shown above. Also note that the input data were amplitudes, which means that weaker reflections may have been modified by French-Wilson treatment or discarded altogether, and the original intensities will not be recovered. P 1 2 1: no systematic absences possible P 1 21 1 (input space group): no absences found =============== Diagnostic tests for twinning and pseudosymmetry ============== Using data between 10.00 to 1.98 Angstrom. ----------Patterson analyses---------- Largest Patterson peak with length larger than 15 Angstrom: Frac. coord. : 0.212 0.500 0.001 Distance to origin : 41.574 Height relative to origin : 5.195 % p_value(height) : 8.830e-01 Explanation The p-value, the probability that a peak of the specified height or larger is found in a Patterson function of a macromolecule that does not have any translational pseudo-symmetry, is equal to 8.830e-01. p_values smaller than 0.05 might indicate weak translational pseudo symmetry, or the self vector of a large anomalous scatterer such as Hg, whereas values smaller than 1e-3 are a very strong indication for the presence of translational pseudo symmetry. ----------Wilson ratio and moments---------- Acentric reflections: /^2 :2.117 (untwinned: 2.000; perfect twin 1.500) ^2/ :0.775 (untwinned: 0.785; perfect twin 0.885) <|E^2 - 1|> :0.749 (untwinned: 0.736; perfect twin 0.541) Centric reflections: /^2 :2.936 (untwinned: 3.000; perfect twin 2.000) ^2/ :0.657 (untwinned: 0.637; perfect twin 0.785) <|E^2 - 1|> :1.014 (untwinned: 0.968; perfect twin 0.736) ----------NZ test for twinning and TNCS---------- The NZ test is diagnostic for both twinning and translational NCS. Note however that if both are present, the effects may cancel each other out, therefore the results of the Patterson analysis and L-test also need to be considered. Maximum deviation acentric : 0.030 Maximum deviation centric : 0.033 _acentric : +0.022 _centric : -0.017 ----------L test for acentric data---------- Using difference vectors (dh,dk,dl) of the form: (2hp, 2kp, 2lp) where hp, kp, and lp are random signed integers such that 2 <= |dh| + |dk| + |dl| <= 8 Mean |L| :0.502 (untwinned: 0.500; perfect twin: 0.375) Mean L^2 :0.334 (untwinned: 0.333; perfect twin: 0.200) The distribution of |L| values indicates a twin fraction of 0.00. Note that this estimate is not as reliable as obtained via a Britton plot or H-test if twin laws are available. Reference: J. Padilla & T. O. Yeates. A statistic for local intensity differences: robustness to anisotropy and pseudo-centering and utility for detecting twinning. Acta Crystallogr. D59, 1124-30, 2003. ================================== Twin laws ================================== ----------Twin law identification---------- Possible twin laws: ------------------------------------------------------------------------------- | Type | Axis | R metric (%) | delta (le Page) | delta (Lebedev) | Twin law | ------------------------------------------------------------------------------- | PM | 2-fold | 2.324 | 2.260 | 0.032 | h,-k,-l | ------------------------------------------------------------------------------- 0 merohedral twin operators found 1 pseudo-merohedral twin operators found In total, 1 twin operators were found Please note that the possibility of twin laws only means that the lattice symmetry permits twinning; it does not mean that the data are actually twinned. You should only treat the data as twinned if the intensity statistics are abnormal. ----------Twin law-specific tests---------- The following tests analyze the input data with each of the possible twin laws applied. If twinning is present, the most appropriate twin law will usually have a low R_abs_twin value and a consistent estimate of the twin fraction (significantly above 0) from each test. The results are also compiled in the summary section. WARNING: please remember that the possibility of twin laws, and the results of the specific tests, does not guarantee that twinning is actually present in the data. Only the presence of abnormal intensity statistics (as judged by the Wilson moments, NZ-test, and L-test) is diagnostic for twinning. ----------Analysis of twin law h,-k,-l---------- H-test on acentric data Only 50.0 % of the strongest twin pairs were used. mean |H| : 0.433 (0.50: untwinned; 0.0: 50% twinned) mean H^2 : 0.273 (0.33: untwinned; 0.0: 50% twinned) Estimation of twin fraction via mean |H|: 0.067 Estimation of twin fraction via cum. dist. of H: 0.045 Britton analyses Extrapolation performed on 0.00 < alpha < 0.495 Estimated twin fraction: 0.024 Correlation: 0.9979 R vs R statistics R_abs_twin = <|I1-I2|>/<|I1+I2|> (Lebedev, Vagin, Murshudov. Acta Cryst. (2006). D62, 83-95) R_abs_twin observed data : 0.443 R_sq_twin = <(I1-I2)^2>/<(I1+I2)^2> R_sq_twin observed data : 0.248 No calculated data available. R_twin for calculated data not determined. ======================= Exploring higher metric symmetry ====================== The point group of data as dictated by the space group is P 1 2 1 The point group in the niggli setting is P 1 2 1 The point group of the lattice is P 2 2 2 A summary of R values for various possible point groups follow. ---------------------------------------------------------------------------------------------- | Point group | mean R_used | max R_used | mean R_unused | min R_unused | BIC | choice | ---------------------------------------------------------------------------------------------- | P 1 2 1 | None | None | 0.443 | 0.443 | 2.749e+05 | | | P 2 2 2 | 0.443 | 0.443 | None | None | 4.453e+06 | | ---------------------------------------------------------------------------------------------- R_used: mean and maximum R value for symmetry operators *used* in this point group R_unused: mean and minimum R value for symmetry operators *not used* in this point group An automated point group suggestion is made on the basis of the BIC (Bayesian information criterion). The likely point group of the data is: P 1 2 1 Possible space groups in this point group are: Unit cell: (59.554, 79.206, 87.43, 90, 92.26, 90) Space group: P 1 21 1 (No. 4) Note that this analysis does not take into account the effects of twinning. If the data are (almost) perfectly twinned, the symmetry will appear to be higher than it actually is. ================== Twinning and intensity statistics summary ================== ----------Final verdict---------- The largest off-origin peak in the Patterson function is 5.19% of the height of the origin peak. No significant pseudotranslation is detected. The results of the L-test indicate that the intensity statistics behave as expected. No twinning is suspected. ----------Statistics independent of twin laws---------- /^2 : 2.117 (untwinned: 2.0, perfect twin: 1.5) ^2/ : 0.775 (untwinned: 0.785, perfect twin: 0.885) <|E^2-1|> : 0.749 (untwinned: 0.736, perfect twin: 0.541) <|L|>, : 0.502, 0.334 Multivariate Z score L-test: 1.973 The multivariate Z score is a quality measure of the given spread in intensities. Good to reasonable data are expected to have a Z score lower than 3.5. Large values can indicate twinning, but small values do not necessarily exclude it. ----------Statistics depending on twin laws---------- ----------------------------------------------------------------- | Operator | type | R obs. | Britton alpha | H alpha | ML alpha | ----------------------------------------------------------------- | h,-k,-l | PM | 0.443 | 0.024 | 0.045 | 0.022 | -----------------------------------------------------------------