BULLETINS OF THE Zoological Society of San Diego No. 18 1. TAIL-LENGTH DIFFERENCES IN SNAKES WITH NOTES ON SEXUAL DIMORPHISM AND THE COEFFICIENT OF DIVERGENCE. 2. A GRAPHIC METHOD OF SHOWING RELATIONSHIPS. By L. M. KLAUBER Consulting Curator of Reptiles, Zoological Society of San Diego SAN DIEGO, CALIFORNIA September 1, 1943 Digitized by the Internet Archive in 2017 with funding from IMLS LG-70-15-0138-15 https://archive.org/details/bulletinsofzoolo1819unse ZOOLOGICAL SOCIETY OF SAN DIEGO BOARD OF DIRECTORS 1943-1944 T. M. Russell, President C. L. Cotant, First Vice-President Gordon Gray, Second Vice-President Fred Kunzel, Secretary Robert J. Sullivan, Treasurer F. F. Annable Mrs. Robert P. Scripts Dr. Thos. O. Burger L. M. Klauber W. D. Crandall F. T. Olmstead Milton G. Wegeforth EXECUTIVE COMMITTEE T. M. Russell, Chairman C. L. Cotant Gordon Gray Fred Kunzel Robert Sullivan F. T. Olmstead FINANCE COMMITTEE C. F. Cotant, Chairman F. F. Annable F. T. Olmstead Robert J. Sullivan hospital and research COMMITTEE Rawson J. Pickard, M.D., Chairman W. C. Crandall Joshua F. Baily, Sc. D. Howard A. Ball, M.D. Thos. O. Burger, M.D. Denis L. Fox, Ph.D. George F. Kilgore, M.D. Hall G. Holder, M.D. Eaton M. MacKay, M.D. Francis M. Smith, M.D. Q. M. Stephen-Hassard, D.D.S. Claude E. ZoBell, Ph.D. Theodore D. Beckwith, Ph.D. James D. Edgar, M.D. Meyer Wiener, M.D. A. M. Burlingame C. B. Perkins F. M. Klauber LEGAL COMMITTEE Fred Kunzel Gordon Gray building and maintenance COMMITTEE F. T. Olmstead, Chairman George Ray William P. Kemper Robert J. Sullivan EXHIBIT AND EQUIPMENT COMMITTEE F. F. Annable, Chairman Dr. T. O. Burger Gordon Gray Mrs. Robert P. Scripps F. M. Klauber Earl Warren GROUNDS COMMITTEE Mrs. Robt. P. Scripps, Chairman Mrs. Seth Anderton Fred Kunzel F. F. Annable Allen Perry EDUCATION COMMITTEE W. C. Crandall, Chairman James H. House Mrs. W. H. Porterfield W. E. Harper PUBLICATIONS COMMITTEE Robert J. Sullivan, Chairman Mrs. Robert P. Scripps C. F. Cotant MARINE MUSEUM (STAR OF INDIA) Capt. Robert Baker, Chairman Fred Kunzel Gerald MacMullf.n Bates Harper D. P. Marvin W. C. Crandall BULLETINS OF THE ZOOLOGICAL SOCIETY OF SAN DIEGO No. 18 I. TAIL-LENGTH DIFFERENCES IN SNAKES, WITH NOTES ON SEXUAL DIMORPHISM AND THE COEFFICIENT OF DIVERGENCE. II. A GRAPHIC METHOD OF SHOWING RELATIONSHIPS. BY L. M. Klauber Consulting Curator of Reptiles, Zoological Society of San Diego SAN DIEGO, CALIFORNIA SEPTEMBER 1, 1943 TABLE OF CONTENTS Page I. Tail-Length Differences in Snakes, with Notes on Sexual Dimorphism and the Coefficient of Divergence 5 Introduction 5 Statistical Formulas 5 Spurious Correlation 9 The Coefficient of Divergence as a Measure of Sexual Dimorphism in Scutellation 11 Methods and Precautions in Tail-length Studies 15 Variation of Tail-length within Fdomogeneous Populations 21 Correlation of Tail Length and Body Thickness 39 Difference Problems 40 Sexual Dimorphism 46 Rattlesnake Tail Proportionalities 51 Acknowledgments 5 6 Summary and Conclusions 57 Bibliography 5 9 II. A Graphic Method of Showing Relationships 61 Klauber: Tail-length Differences in Snakes 5 TAIL-LENGTH DIFFERENCES IN SNAKES, WITH NOTES ON SEXUAL DIMORPHISM AND THE COEFFICIENT OF DIVERGENCE Introduction In the course of certain studies of herpetological correlations the pro- portionate tail lengths of snakes were investigated. This led to a detour of such extent that it is deemed advisable, in the interest of clarity and balanced treatment, to offer this discussion of tail length separately. At the same time it appears opportune to give examples of the use of the coefficient of divergence as a measure of sexual dimorphism or ontogenetic differentiation. The tail-length ratios of snakes often prove useful in diagnosis and systematics. The proportion seems to be a rather stable character, so that when differences between related forms do appear, they are likely to be of importance. However, it is a character with respect to which the determination of the significance of differences is not simple. In contrast with characters of lepidosis, which are subject to individual and terri- torial variations, and sometimes sexual dimorphism, the tail ratio has all of these and usually an ontogenetic variation as well, that is, a change in the ratio of the tail length to the length of body as a snake grows. As is always the case when comparisons are to be made, we must make sure that our samples — that is, the groups of specimens to be compared — are homogeneous. Thus, if we are comparing two series with respect to a character in which sexual dimorphism is present, and we fail to treat the sexes separately, we may find an apparent difference which really results from an accidental sexual unbalance in the samples, rather than a true difference between the forms being compared. And so, as it will be shown later that tail length is ontogenetically variable in most species of snakes, complete homogeneity can only be secured if we limit our samples to specimens of uniform age. The virtual impossibility of obtaining ade- quate series under such a restriction renders it necessary to make special statistical provisions for combining specimens of diverse ages. Failure to take this ontogenetic variation into consideration in making species com- parisons may lead to false conclusions. It is the purpose of this paper to discuss methods of combining specimens, and to show the extent of ontogenetic variation in several example species. Methods of evaluating differences are developed. Sexual dimorphism is treated. Statistical Formulas Problems such as those of tail proportionality, its variation within a homogeneous series of specimens, and differences between series, are sub- 6 Bulletin 18: Zoological Society of San Diego ject to both analytic and graphic attack. If the results are to be trust- worthy both methods are usually advisable, since they supplement each other in affording an understanding of the nature of the variation involved. Where, as in the present case, a preliminary survey of the most super- ficial nature indicates that ontogenetic variation is probably present, the problem becomes one of correlation — the correlation of tail length propor- tionality with age. This is not to say that the coefficient of correlation necessarily affords the best measure of the concomitant relationship; on the contrary, in most morphological surveys of this kind, where the cor- relation between a body part and the whole body, or between two body parts is under investigation, the correlation is sure to be high. In such cases the direction of the regression line, and the extent and nature of the scatter of the individual specimens about that line, will be of greater in- terest and importance than the numerical value of the correlation co- efficient. The methods of calculating coefficients of correlation, of determining regression lines and errors of estimate, will be found fully detailed in every text book of statistical methods, and therefore will not be discussed here, ffowever, certain formulas involving the relationship between two parts and a whole — the body length, tail length, and length over-all of a snake, for example — are not so readily available, or, if given in a text, may be in a form not directly applicable to the present problem. These formulas are therefore given here, although it should be understood that they in- volve no originality whatever. The symbols I shall use are as follows: L, B, T, represent length over-all, body, and tail lengths, respectively. Thus, in each specimen, L = B + T. Ml, Mi-„ Mr are the sample means of the same quantities, and (Tl, , Mr — Ml — Mt and a' = Mx — b'M p or the form b' — (AL2 — AT2) / 2Ap2 - /2 may be used. These equations involve only factors derived directly from the original statistics of T and L, or which can be determined by the use of the equations already given. If the relationship between T and L is linear, that between T and B will likewise be linear. Further, if there is a constant term in the T on L equation, there must also be one in the T on B equation, provided the correlation is high. The latter statement is equivalent to saying that if the ratio of the tail length to length over-all changes as a snake grows, then the ratio of the tail length to body length does likewise. Klauber: Tail-length Differences in Snakes 9 If the correlation Ttl is quite high, approximations to the regression constants of the T on B equation may be had directly from the following: a' = a / ( 1 - b) and b' — b / ( 1 — b) With respect to the coefficient of divergence, its use will be enhanced if we know its standard error. In problems involving large samples, when determining the significance of the difference between two means, Mx - My, it is satisfactory to compute the ratio of this difference to its standard error (o 1% * V 1 ) » *\ \ © 0 • • \© 0 'i o \ n§\ » 4 V< ° £ ° O ju 1 ’V® \* Sf CX° \ • ; • ^ s o \ ois !o vO *!• *% • ‘ Vo c > _c k- o BOREGO SERIES TIONSHIP OF LENGTH OVER-ALL TO TAIL LENGTH MALES • FEMALES o ft, • • Vi • 1 0 \ V \ \ N \ V, V \\ \\ \ \ \ \ — \ \ \ \ \ >, < -I d o : \ \ \ V \ oooooooo cor-"> n OjC co -o c « O c u (— l P3 ri J ^ ° CO X ° w "O O M « O • -h rQ “ “■£ ” « ¥ a g -a £j < rS C* T3 V-i crj >% *-» - 3 -g ^ C -n O C 3 < « C/5 a> *_» J-. t/5 60 G a> 0 oP U O w Zr> C in CJ u! *5 *<3 a> -Q g 3 ^ £ ^ Mh O 1—H ON o OO O r-H r-H l\ o *n~s r-H t-H r-H i\ ON T-H Kn*\ NO N~ CM XN OO lx fA CM Kns N" S *S~\ C*^ ON o ‘✓“N ON N NO N" NO N" OO N" r\ *n~\ OO r\ rx Kns NO T-H r<^ CM oo T— < 1”H r<5 r<5 tN * r^> rO N" CN t— H r-H t-H CN CN r^i Wt-n r*~\ CN O o o o O o o o o O o o o o o o o o o o o o o o o o o E Kn~\ (N *V^N o t^N o o IjTN o o o o o o o o o o Kn~\ o CN o o Kn-\ Kns o E (N m CN ON OO [\ CM r-H o N" Kn~\ ON CM CN NO NO IX CN NO CN t-H Kns r-H r— H t-H T— H t-H ON oo NO CN O O ON r-H ON o CN N" CN NO m IX IX CN On t-H L«“n N" Po O o o (N OO tn 0“N o o NO hN OO N“ l\ N" NO ON t-H N" i\ t/-N U«-N oo fx IX Wr\ OO CN fx IX IX (N CN r^i rO <^5 r^5 CN CN • • • ' CN CN in~\ Po o ON OO ^t" t\ N" ON r\ OO i\ ON r**i ON ON ON K h^i <^N N 3 Vh CD <3 5h o c 03 $H o3 £ o3 C _T ri C r*. o3 ^ a, . o3 r-c S ^ ^ £ 3 g «UeqUh4 c ri -u> 03 CJ 3 ^ d -U "03 w N C/d 4-> o ^ *>T ?s ^ - ^ & <5 <3 V. <3 k^5 £ ^ 40 ^ O -C' ^ Kn} <3 <3 '■o • *s O S3 S3 — «^ ^ • **•* ■ **>■• UO^ o •>* S *i? rS O ^ u o •K* Co • *v» TV* o e» C- <5 v v cs -ie -~r .A ^ u u ^ u PQ CD WD C n u CD &0 ° > o in in % £> S: <=o ^ <3 O V- 'i 6 *K» ^~* o tyj r -5 C^5 C^; in} * *'■>« * **«• * ^ •Cl, -cl, -Cl o o o Co Co Co •Ni •% * *N ^ Hv* in} in} 0-3 >S CS « £ § £ *3 >* OD °? ^ Co O » 'tr O S-1 ?s Co hj rt X a» rt C o N ’ C < r3 5 CJ N «rC < -^ < ° O - U c3 - . o3 , l rt D • 40 p .t! 40 fl c S *j g w o 3 £ cj ^ 0 ^ W U bo <^5 ^ ‘ O ’S- Co * N, • >v. ^ iao ^ . in} in} ,_ CO d •cp, tu 'S r>« —o o N • O, V to C <5 -8, -Cp o o o v -Plo rt £ o oC ri O on ri ri X C o o c; in) in n £ ££ U £ n ** Co Lo ^ be •N *»v* <-o 3J t>4 •S «o Lo ^ bo S3 •* s H v -2 0; 8 ^ u -3S o .in} 3j «-o 3S >o r- &c ^ to to 3a 3j 8 t 3a >o f—< h~t Except Laguna Island. Klauber: Tail-length Differences in Snakes 3 5 data on three series, T. o. biscut atus, T. hammondii, and T. o. ordinoides (vicinity of Portland), are carried forward from Tables 6 and 9. A considerable consistency in adult proportionalities will be found evident in the last two columns of Table 11, both between territorial groups of the same race, and between races. It is of interest to note these adult tail-length proportionalities in the light of the relationships deduced by Dr. Fitch in his monograph.18 On page 8 I have given an equation showing how the mean of the tail ratios Mt/l (or Mt/b), which involves many individual computa- tions, may be derived from the mean of the tails divided by the mean of the lengths over-all, or body lengths (MT/ML or MT/MP>) from other statistics readily available, provided the coefficients of correlation have been calculated. In order to ascertain the importance of the differences between these statistics, I have made four example computations with the following results, which are given in terms of the multiplier to be applied to Mt Mr in order to derive Mt/l: Multiplier Species Male Female P. d. perkinsi (S. D. Co.) 0.9849 0.9926 S. o. annulata 0.9977 0.9969 It will be observed that in these four examples, where the samples are fairly large and the correlation high, the mean proportion differs from the proportion of the means by a maximum of 1.5 per cent. The latter pro- portion therefore may well serve as an approximation. It is to be remem- bered that neither of these statistics is particularly useful, since both depend too much on the ontogenetic distribution of the samples. Should a single statement of tail proportionality be desired to represent a species in mak- ing interspecific comparisons, and it is not possible to restrict the state- ment to some particular age (juvenile or adult, for example) then I should recommend the use of the proportion at the median length, as determined from the regression equation. But for most taxonomic studies adult pro- portionalities are to be preferred, as will be discussed hereafter. With respect to the figures given in Tables 5 to 9, it should be stated that all of the derived statistics were calculated to a further degree of accuracy than shown, usually to 6 figures after the decimal point. This was done so that the derived statistics would not lose in accuracy through rounding off; however, these additional figures are not entered in the tables, as they would lead only to confusion and difficulty in making comparisons, and would give an entirely unwarranted impression with re- spect to the accuracy of the results. This matter of dropping figures is 18 A Biogeographical Study of the Ordinoides Artenkreis of Garter Snakes (Genus Thamnophis) . Univ. Calif. Pubs, in Zook, Vol. 44, pp. 1-1 5 0, 1940. TABLE 11 ■M Can l-H NO NO IX lx 1—4 OO 1—4 r-H l-H X "O T1 ^ PI 0 r-H XI rn OO ON ON ON n O p ^ o -C cs O ri to 4-J 4—1 CN ON 0 CN Can NO ON 0 X 1—4 00 Can l-H pi n ~ H d "D W3 P C H3 O C q 2 ^ pq ^ c c n t-q o o o o »an NO 0 O O O O O 0 O 0 O O 0 O 0 0 O O O »AN Can 0 O Can (50 fi -n t\ X ON oo o o ON A C ^ O 4-T C/5 C on 0> O • — ^ CJ ^ CJ U* OO 0 0 0 ON 0 X ON Ca Cn 0 ON r^i ca> rA, rA> q js .s I « .3 — z ^ h'N 0-H I n ON NO IX NO O ON X ox x- 0N ON X X" oo vo ^ ^ K CA\ i— I CAN CN 1—1 \0 IA X ON T VAN 1-H NO 1— • I l-H OO Ca o oo N“ N~ ^ VO (N OO O (N r-H l-H OO oo xj- NO NO Can \0 *•— ! |\ NO t— 4 < c3 * q> rs X3 - -d rt £ o3 4— > CO 05 C/} 2 lT ^ t*> o rS ri *-. CO A ,nj <73 V3 H— i . ^ — c3 £ Z Z i— I C/5 05 Co Co Cb q C/5 O a CO © © g •© * ’■'* 22S ^ ^ ^ © 3 ° -£5 >-£5 •ci, -ci, o o X X O -ON X O'. rt X -O* p "X ?N Co Cb 5 Co * »S4 ■Cl, <3 N « o 2 VJ I o ^ o Su b 0 u CD c o | c^Z J ■ 4 * "V* • 1 ^ ^ 2: •S -p o 2 o o o U o £ ’u o "■a c QJ X 3 d o Sd CD V, TO c rt J-i O d mh o o rs u pt z cT o rt «> X O 05 PQ r3 "ID r3 C 4H ^5 '5 Jh 0 d 0 X a 00 • — 1 as CJ . N PO r3 U ^ > £ co ‘C < So C^ * *>k • 1X4 o o « <3 • ^ - 1= ^ ^ n. ^ o o 03 ^3 o ^o ^3 o r5 «o co ^ c^ rs ^2 ^ Ooo . rr *< Z- • ^s 1= ^3 V v- o O o co c^i "ts X o Co Co s o "ts X ■C5 ■~p "Cli -Cin o o « o »o o NO ^t" *^4 Ky-\ » O rO ON O ON o 04 04 a\ o ON ■3 04 04 cK 1—^ bC w o Q c o Is. OO 04 rt 04 ^o ON OO rO o rO CO OO t\ ls~s r\ OO rO fs rO ON o NO o t\ 04 ON o Ol OO 04 NO Ui NO 04 NO O O NO NO -\ Is rO ON O y—< OO 04 O OO o ro Ns OD rt C rt 04 no ^ Ns 04 ON ^ ^ ^ r^N On O ON > P< Ph > t. o ^t- O NO |\ O Ns t— < i o O 1-H '/'V U'-S OO o 04 04 ON V^-S V^N (N OO (N ^ O ^KKCNON^O^K '^-CNOrxOrsIr^OOO <^o Ns o NO NO O to c hJ rt H M-l o Cm JS uo C o 4— < rt "w a> on o rt C < 6 ^ U o to ■'1" Tt" i— i o vo o O K ^ 04 ^ ' — 1 O t-A rA ^ * fO NO 04 04 NO 04 04 l^N o NO OO NO O ro rO N 04 tin V NO ON O O ON 04 Os ro Ns rO ON O O OO NO O O O cO NO ON ^ ON 04 ^ K \o 04 Ns O © *^s ^ o rO ON NO OO 04 o N Ui ON W^s NO t}- »A\ o ON o Kr\ 04 04 04 O NO 3 rt § c • rO OO 1-H N S g N V. fsj ON OO 04 ^ Va Ns OO NO ON ON o A 6 < OO NO N Ns o r-H 1— H o K ON O 04 o 04 T-H o t~H •A OO £ .2 4-> c 22 2" s-i c/5 o 'S 4-4 E C o *-i c 4_r .2 c C_> a! £ - E QJ E “ c M « ° C_) u C C .2 .2 * OO ’ (/) OO OO CD oj > •c ° > -c K* i-> M— I 5? O G CO 4-> r3 CT3 p /-s 2 O ’S « ^.§ .95), or highly significantly (O i> .99) are clearly indicated. A third scheme is to go one step farther and utilize Fisher’s method of converting the probability values P into chi-square values at two degrees of freedom, these chi-square results being used as the co-ordinates. The additive quality of chi-square suggests the utility of this scheme; and it may be useful because chi-square tends to spread the higher values of Q (lower values of P). But this very effect is a disadvantage in some problems, since some of the co-ordinates will be so large that the others must be unduly compressed. In cases wherein large series of specimens are available, not only of obesus and tenuis, but of a as well, even small differences in characters will be found to be statistically significant, that is, both Oo and Qt will closely approach unity in almost every character. The result, if prob- abilities be the co-ordinates, is a chart in which all points cluster about the upper right hand corner and conclusions are drawn with difficulty. There is a rather important difference between problems of small samples and those of large. In the former, means, and therefore mean differences, are known only approximately — that is, the samples may or may not accurately define or represent the populations, and it is the populations and not the samples with which we are really concerned." The confidence 1 A further modification may be made by finding the value of P which corresponds to each value of t in a t- table, and retranslating this to the value of t on the normal curve, thus correcting for small samples. 2 The samples comprise the specimens available for study; the populations are the creatures remaining in the wild, from which the samples were collected. 64 Bulletin 18: Zoological Society of San Diego to be assigned to a mean may be deduced from its standard error. Inevitably the sampling variation of the difference between two means is of first importance in judging the reality of the difference. This reality or sig- nificance is dependent on the extent of the difference, and also on the in- ternal dispersions of the character under consideration, in each of the two populations between which the relationship is being judged. Thus Q be- comes a measure of both extent and validity, and is a necessary compro- mise when only small samples are available. But where the samples are relatively large, say about 100 specimens or more, and internal population dispersions are of the order of magnitude usually found in these biological problems, the means are known within a narrow confidence range, and we are then more interested in the direct measure of differences than in a statement of their significances. In this case a method of reducing diverse characters to some comparable standard, independent of the unit of meas- urement, is still necessary. Such a criterion is furnished by the coefficient of divergence, that is, the difference between two means, divided by half their sum. This comprises the fourth scale of co-ordinates which may be suggested. But it should be remembered that the coefficient of diverg- ence, in itself, contains no indication of the significance of a difference, that is, whether it may not be the result of the chance composition of the available samples; it serves merely as a yard-stick to place a variety of characters on a comparable basis of measurement. It may be thought that there is something fundamentally inaccurate in using for co-ordinates a statistic which is based entirely on means, thus ignoring the intrapopu- lation dispersions of the character under consideration. But investigation indicates that there is only a slight correlation between intrapopulation dispersion and interpopulation differences. The use of the coefficients of divergence as co-ordinates occasionally tends to crowd the points of some characters because a few others take high values — a situation which also occurs with / or chi-square. In such situa- tions the use of log-log paper may be recommended. In calculating probabilities (if x be a single specimen) one may determine the position of x directly from the place where it falls in the arrays of each of the major populations, using twice the percentage values of the tails beyond that point as the co-ordinates. In case there is more than one specimen in the class in which x falls, the tail should be considered as be- ginning at the midpoint of the class. For example, assume this distribution in the subcaudals of obesus: Number of subcaudals 20 21 22 23 24 2 5 26 Total Number of specimens 1 13 27 47 41 8 3 140 Klauber: Graphic Method of Showing Relationships 65 Suppose that x has 25 subcaudals. Then, splitting the 2 5 class in obesus, we take the tail of the distribution to be 4 + 3= 7, and find P0 = 2 X 7/140 = .10, orQo = .90 While this modification of the probability method is simple and avoids the necessity of computing the means and standard deviations of each character of obesus and tenuis, the result depends rather too much on the accidental distributions of the samples of obesus and tenuis available to us. Such a method, however, may be justified if the distributions of several of the characters are badly skewed. The use of co-ordinates based on probabilities permits the inclusion of any character the x- value of which may be expressed as a probability (as, for example, by a four-fold contingency table) ; it is not restricted to the difference between two means. If we have more than one specimen of x available, two courses are open. We can plot a O-value separately for each character of each specimen, thus finding as many points for each character as there are specimens. In plotting, suitable colors or symbols may be used to identify either speci- mens, characters, or both. Or, we may calculate the means and standard deviations of each character of the specimens comprising x; and then, by one of the formulas for the significance of the difference between two means, find the values of t, Q, or chi-square corresponding to the dif- ference between the mean of x and the mean of obesus, and between the means of * and tenuis. Thus, as was the case when v was only a single specimen, we shall have a pair of co-ordinates defining a single point for each character. This latter is the only practical alternative if x contains more than three or four specimens. Of course, the fourth scale suggested uses means in any case, but is not to be considered adequate if v is a small sample. As a first illustrative problem I shall take a single rattlesnake specimen from Shoup, Lemhi County, Idaho. This is on the upper Salmon River. It is known that rattlesnakes are found along this river down to its con- fluence with the Snake River, and that the lower reaches of the Salmon and beyond are inhabited by Crotalus viridis oreganus. In the opposite direction Crotalus viridis viridis has been collected across the Continental Divide (near the Bitter Root Range) in the Valley of the Beaverhead River, Montana. The Divide can be crossed at a number of points at about 8000 feet elevation; and, as rattlers have been observed in this vicinity at least up to 7000 feet, there is little doubt that they range across the Divide. The problem is whether the specimen from Shoup more nearly resembles C. v. oreganus down the river to the west, or C. v. viridis across the mountains to the east. We shall take, as representative populations of these two subspecies, the two nearest available series; these are a collection of oreganus from 66 Bulletin 18: Zoological Society of San Diego Pateros, in eastern Washington, and a series of viridis from Montana. The essential data on these series are set forth in Table 1, ten characters being given. Only males are included in the statistics of those characters wherein sexual dimorphism is present, because our specimen from Shoup is a male. Table 2 gives the derived differences, and the values of Q calculated there- from. Chi-square results are also included for comparative purposes. TABLE 1 Character Data The Relationship of a Specimen from Shoup, Idaho, to C. v. viridis, and C. v. oreganus Montana viridis Pateros oreganus Specimen — - — — - from Character N M (7 N M (7 Shoup Scale rows 174 25.443 0.913 615 25.511 0.873 27 Ventrals 70 174.943 4.363 326 173.313 3.125 175 Subcaudals 72 24.93 1 1.417 324 22.985 1.446 26 Supralabials 340 14.406 1.185 1229 15.218 0.832 14/2 Infralabials 345 15.149 1.093 1230 16.178 0.892 14/2 Scales on crown 151 18.728 4.516 267 24.985 4.839 18 Intersupraoculars 151 2.775 0.799 278 5.813 0.986 2 Body blotches 171 45.819 3.557 616 33.195 2.293 48 Tail rings 70 9.786 1.512 326 5.497 0.774 10 Width of button, mm. 25 5.248 0.2 56 252 5.611 0.365 4.8 * Mean of the two sides. TABLE 2 Difference Data The Relationship of a Specimen from Shoup, Idaho, to C. v. viridis and C. v. oreganus Differences from viridis Differences from oreganus Character D t Q Chi-square D t Q Chi-square Scale rows 1.557 1.705 .912 4.860 1.489 1.706 .912 4.861 Ventrals 0.057 0.013 .010 0.020 1.687 0.540 .41 1 1.059 Subcaudals 1.069 0.754 .546 1.579 3.01 5 2.085 .963 6.594 Supralabials 0.094 0.079 .063 0.130 0.718 0.863 .612 1.893 Infralabials 0.649 0.594 .447 1.185 1.678 1.881 .940 5.627 Scales on crown 0.728 0.161 .128 0.274 6.985 1.443 .851 3.808 Intersupraoculars 0.775 0.970 .668 2.205 3.813 3.867 .999+ 18.226 Body blotches 2.181 0.613 .460 1.232 14.805 6.457 .999+ 45.916 Tail rings 0.214 0.142 .112 0.238 4.503 5.818 .999+ 37.858 Width of button, mm. 0.448 1.750 .907 4.750 0.811 2.222 .974 7.299 Klauber: Graphic Method of Showing Relationships 67 1.0 .8 to D Z < O Id cr o to O tr UJ h- < CL 2 O o: z UJ a: UJ u — .4 O u O >■ f- CD < CD O T2 • TAIL RIN SS • BOD • INFR — Y BLOTCHE • SUBC ^LABIALS • — S IN AUDALS FERSUPRAO RATT CULARS LE BUTTOr' SCALE ROW c S • • SCALES ON CROWf' SUP • RALABIALS VENTRAL 9 S FIG. 1 RELATIONSHIP OF A SPECIMEN FROM SHOUP, IDAHO TO C.v.viridis AND C.v. oreganus (basis of comparison, q) .2 .3 .4 .5 .6 .7 .8 PROBABILITY OF DIFFERENCE FROM MONTANA VIRIDIS Using the paired values of O as co-ordinates, the graphic results are shown in Fig. 1. There can be no question that the specimen from Shoup more closely follows viridh than oreganus. With the single exception of the scale rows, all characters fall on the viridh side of a diagonal from 0-0 to 1-1. Our specimen v differs from oreganus significantly in no less than five characters; in three of these the differences are highly significant. It is of interest to note that this conclusion is substantiated by other non- measurable differences, such as the width of the postocular light stripe, the nature of the blotch borders, etc. Other specimens of rattlers from the Lemhi area are available; they have been omitted from considera- tion in this study in order to reduce our problem to an example in which x comprises a single specimen. Parenthetically, it may be stated that the additional specimens verify the conclusion that the Lemhi snakes are nearer viridh than oreganus, although with some suggestions of intergrada- tion with both. There seems little doubt that a full series collected down the Salmon would complete the link. 68 Bulletin 18: Zoological Society of San Diego As a second illustrative problem I shall consider the relationships of Crotalus viridis abyssus of the Grand Canyon in Arizona, with its neigh- bors to the north and south. This will serve as an example wherein a few specimens of each of three forms are available, as is typical of many taxonomic problems. To the south of the Grand Canyon, between El Tovar and Williams, lies the Coconino Plateau. This is inhabited by a form of Crotalus viridis nuntius, a somewhat larger and differently colored snake from the typical nuntius of the vicinity of Winslow. Intergradation of nuntius and abyssus seems evident along the south rim of the Canyon. Most of the available specimens of this Coconino Plateau form of nuntius were collected some miles south of the rim, at Anita and Valle. North of the Canyon occurs Crotalus viridis lutosus, the Great Basin Rattler. While considerable areas of the north rim may be too high for rattlers, there is a ready means of access between the territory of lutosus and the Canyon via the gorge cut by Kanab Creek. Fortunately, we have available a small series of specimens from the upper Kanab in western Kane TABLE 3 Character Data The Relationship of C. v. abyssus with C. v. nuntius and C. v. lutosus Coconino nuntius Grand Canyon abyssus Kanab lutosus Character N M jj. cr N M o' N M a* Scale rows 23 24.217 0.930 16 25.375 0.780 18 25.667 1.154 Vent rals males 13 167.33 8 1.678 11 177.818 2.249 11 180.3 5 5 2.307 females 10 172.900 1.921 4 180.750 1.480 7 184.571 4.370 Subcaudals males 13 24.308 1.726 11 25.545 1.827 11 24.091 1.564 females 10 18.700 1.418 5 20.200 0.748 7 20.143 1.884 Supralabials 46 14.935 0.919 32 15.781 0.892 35 15.486 0.873 Infralabials 46 15.478 1.246 32 16.281 1.123 35 15.629 0.928 Scales on crown 23 3.739 0.529 15 6.067 0.998 18 5.778 1.272 Intersupraoculars 23 23.435 4.726 15 36.067 5.423 18 31.278 9.056 Body blotches 23 41.000 2.904 13 42.385 2.923 18 39.722 3.228 Tail rings males 13 8.923 0.917 11 8.635 1.494 12 7.167 0.898 females 10 7.000 0.894 4 6.750 0.435 6 5.833 1.067 Width of button, mm... 4 4.600 0.123 5 5.280 0.232 5 5.680 0.117 * Standard deviation of sample, not estimated standard deviation of the population. Correction was made elsewhere in the computations. Klauber: Graphic Method of Showing Relationships 69 TABLE 4 Difference Data * The Relationship of C. v. abyssus with C. v. nun tius and C. v. lutosus Character Differences from nuntius Differences from lutosus D t Q CD% D t Q CD% Scale rows 1.158 3.97 .999 + 4.7 0.292 0.82 .58 1.1 Ventrals males 10.280 9.65 .999 + 6.0 2.537 2.49 .98 1.4 females . 7.850 6.80 .999 + 4.4 3.821 1.54 .84 2.1 Subcaudals males 1.237 1.63 .88 5.0 1.454 1.92 .93 5.9 females 1.500 2.06 .94 7.7 0.057 0.06 .05 0.3 Supra labials 0.846 3.99 .999 + 5.5 0.295 1.35 .82 1.9 Infralabials 0.803 2.88 .995 5.1 0.652 2.56 .99 4.1 Scales on crown 2.328 9.05 .999 + 47.4 0.289 0.70 .51 4.9 Intrasupraoculars 12.632 6.77 .999 + 42.4 4.789 1.63 .89 14.2 Body blotches 1.385 0.92 .63 3.3 2.663 2.28 .97 6.5 Tail rings males 0.288 0.55 .42 3.3 1.468 2.75 .99 18.5 females 0.250 0.49 .37 3.6 0.917 1.46 .82 14.6 Width of button, mm.. 0.680 4.69 .998 13.8 0.400 3.08 .98 7.3 County, Utah. Our problem then, is to determine the degree of relation- ship between abyssus and Coconino nuntius, on the one hand, and abyssus and the Kanab lutosus on the other. I have omitted from the abyssus sample, specimens from along the south rim that might be intergrades; thus, the sample comprises only true abyssus from within the Canyon. To make the problem more typical, not all available nuntius from the Coconino Plateau have been included, since such large series are rarely at hand in taxonomic studies of this kind. Tables 3 and 4 present the statistics of the problem, while Fig. 2 gives the graphic result. The significances of the differences between means have been computed by the null method. Because of dimorphism the sexes have been segregated in ventrals, subcaudals, and tail rings. It is not to be deemed inconsistent that the sexes do not agree in the significance of re- sults; it should be remembered that Q depends on N as well as on means and standard deviations. The results in pictorial form show conclusively that the Canyon speci- mens are more closely related to the population to the north than that to the south. In several characters (male ventrals, infralabials, and rattle buttons) abyssus differs significantly from both adjacent populations. (Note how the points representing these characters crowd into the upper right-hand corner.) The relationship with nuntius (as measured by prob- 70 Bulletin 18: Zoological Society of San Diego • SUBC AUDALS 9 SCALES — ON CROWr • Si sc 1 r*~* — * ALE ROWS | |/ i t SUPRALABIALS' / VENTRALS 9 INTERSUPRAOCULARS 1 _ L •• y '• RA VENTRAL. TTLE BUT ' INFRALAf SUBCA o' roN' x 3IALS " UDALS • 4 BODY BL DTCHES / TAIL UNGS J _ © • TAIL RIN ;s 9 FIG. 2 RELATIONSHIP OF C.v. abyssus TO C.v. nuntius AND C.v. lutosus (basis of comparison, q.) .8 D Z O cr ^ .6 Id O Z UJ o: w .5 u_ Ll. Q U- O CD < CD .3 O cr CL .3 .4 .5 .6 .7 PROBABILITY OF DIFFERENCE FROM LUTOSIJS .8 .9 1.0 ability) is closer than with lutosus in four characters (three of which are independent3) : Male subcaudals, body blotches, and the tail rings of both sexes. It will be observed that the link with nuntius is closer in all pattern characteristics; to this extent pattern gives a different answer than scale counts. I should consider the latter the more important. As a final example I take one in which large series of all three forms to be compared are available, so that even relatively small differences are likely to be statistically significant, yet may be unimportant in extent. In this example the populations compared are territorially more scattered than those previously discussed, as is generally the case in such taxonomic problems. 3 It is obvious that two independent characters should be given greater weight, in reaching a decision, than two which are correlated. This matter of character correlation will be the subject of a paper now in course of preparation. Klauber: Graphic Method of Showing Relationships 71 The diamondback rattlesnakes, C rot aim cinereous , C. lucasensis and C. ruber are obviously related. Superficially, lucasensis appears to be inter- mediate, as might be expected from the evidence of other elements of the Cape fauna, which were derived from mainland Mexico rather than from California. I shall endeavor to ascertain the comparative relationship of lucasensis with ruber and with cinereous. For the lucasensis material I shall use all the available non-island specimens from southern Lower Cali- fornia. The ruber material is from southern California, and northern and central Lower California. The cinereous series is from southern and western Arizona, southeastern California, and northeastern Lower California. Since all series are large, the significances of differences have been cal- culated by dividing the differences by their standard errors and taking the values of Q directly from a table of the normal curve. Several additional characters have been introduced, including head length, and tail propor- tionality. I have discussed elsewhere the method of determining the sig- nificances of differences in these ontogenetically variable characters at standard body lengths. The statistics of this problem are set forth in Tables 5 and 6; the graphic presentation in Fig. 3. We should first verify the validity of placing lucasensis between cinereous and ruber. We find, from a study of Table 5, that lucasensis falls between cinereous and ruber in 7 independent characters; it falls outside of cinereous in 2%, and outside of ruber in 5Yg. Most of the characters in which lucasensis is more highly differentiated from cinereous than is ruber are head scales — labials, scales on the crown, scales around the rostral, and loreals. I deem these less important than scale rows, ventrals, and body proportions, and therefore conclude that lucasensis is intermediate. That both were derived from cinereous seems clear; but whether separately, or lucasensis first, and then ruber as an offshoot of the latter, is somewhat uncertain. The fact that island forms are rather close to one or the other, suggests a separate development, followed by a gradual invasion of inter- mediate territories. Whether there is intergradation at present is not definitely known, but seems probable, from the few scattered specimens available between Concepcion Bay and La Paz. There is no indication of intergradation of either lucasensis or ruber with cinereous. Fig. 3 indicates without question that lucasensis is more closely allied to ruber than to cinereous, thus solving our problem. In only 8 out of 21 characters (not all independent) are the points closer to cinereous than ruber and all of these fall near the diagonal which represents the line of equal separation. On the other hand, 7 of the 13 characters which favor ruber as the closer relative, are much nearer to the ruber base than that of cinereous. This cinereous-lucasensis-rubcr example illustrates the inadequacy of using co-ordinates based on Q where large numbers of specimens are available, for many O-values closely approach unity. Values of t are better Character Data The Relationship of C. lucasensis with C. cinereous and C. ruber 72 Bulletin 18: Zoological Society of San Diego b 04 o r\ 04 ON oo ON ’F O ON ee, oo O oo NO r-H NO tX OO Tt~ o O NO ee> t-H o Lex ON oo O © NO tx ee> NO NO 04 r-H NO t-H ON r-H T-H K OO NO oo o ON r*N 04 KS"\ © OO NO T-H r-H i\ rr\ ee, ee> rO rX i— » — 1 04 r-H 1 — 1 r-H r-H ON On ON rO Tf- r\ eO Ce On o ^Ct" ON E*N 04 ON NO OO O r^i NO Cr-\ 04 o OO ON O K T-H o <^N T-H On o NO OO ON »— • «o ON o OO OO ON NO o NO KS~\ <^N NO oo 04 Cr\ r-H o to T-H O — 1 d d o 04 o o 04 C/-N N|- EE-> Ctn •r* 04 (N| oo ON 04 04 r^\ 04 een> o o t\ o ^1" e l< OO 04 Ctn d NO K Nt" Cr-N NO rO T-H ON EE% oo OO t-H NO 04 t\ 04 oo ON 04 04 r-H T-H (N 04 iy~\ ON (X (X ON o t-H 1-1 r-H o Ct-n rO ON I\ l\ <^N 04 NE EET o o E*N 04 ON tT ON t\ OO rO rO ^cj- ON EE% NO NO NO NO rO T-H r-H r-H T-H NO NO r t Jj rti 2 E £ 03 ’TD G 03 CJ _C c n CO CD 03 E CO o u G O 03 15 03 03 U co 03 CE ^ G O 60 O 2 o S n O CL h 3 w E § C3 Vh C "TD £ «N c3 £ 4-> o Ih a G "T" o O o o r— H o o r-H o NO o o o G— n r-H HO OO r— H r-H oo ’-4 <73 CD re CJ CJ G5 mal £ c <73 CD r-Q n u CX r3 U Mh oo CJ <73 CO u CJ 4-* U o CO CJ <73 E C/5 CD c3 £ r o _C Ijt ^ <73 £ C/5 — ' 4-> <73 *“3 £-5 £ r-H O jy „ SJ rt 4-1 2 3 ’T3 £ O U O •- C73 74 Bulletin 18: Zoological Society of San Diego criteria, but I think the extent of the differences, when reduced to a comparable basis, as afforded by the use of the coefficient of divergence,1 gives the better illustration. However, if one desires to incorporate the effect of internal dispersion, then t should be employed. In many genera of snakes there will be fewer countable characters than in the rattlers. Scale rows and labials, for example, are often invariant. On the other hand, use can be made of characters which I have not employed in these rattler examples, such as internasals, oculars, temporals, genials, the ventral count at locations where scale rows are dropped, or any other counts in which the forms being compared are not constant within a population. If there is enough variation in any character to produce an K x 2 table, a value of Q may be deduced. There are also the variables of head and tail ratios, although, in using them, care must be taken to avoid the fictitious effects of ontogenetic variation. 4 The percentages by which the character means of lucasensis are above or below th corresponding means of cinereous and ruber may also be used as co-ordinates. 100 7 / / gs VENTRA 4 VENT -s d RALS • s > • CAL AIL _ES AROUN LENGTH D ROSTRA ? Dl L IN VIDE FRA D F LA IRST ~ BIALS / / / / / / / 0 SI BCAU daC'i / d / / / FIG. 3 RELATIONSHIP OF TO C. cinereous AND (basis of comparison, C. lucasensis / / / / / / T AIL Rifs GS 9 ® INFRAL ABIALS A C. ruber COEFFICIENT s UPRALABI OF DIVERGENCE IN PER CENT) ALS# 80 60 40 20 10 z Ld u CC Ld CL (T bJ CO D CC o tr I.J z Ld o cr Ld L_ o z Id U Li_ !l_ LU O o .2 .4 .6 .8 I. 2 4 6 8 10 20 40 60 80 100 200 COEFFICIENT OF DIVERGENCE FROM CINEREOUS PER CENT HEAD LENGTH Klauber: Graphic Method of Showing Relationships 75 It will be observed that in only one of these examples have I taken account of the direction of differences. Even when our x lies territorially between K. obesus and K. tenuis, its characters may not always fall be- tween then means. Quantitatively this is provided for in the values of Q or the coefficient of divergence. If it be desired to distinguish between the characters in x which fall between obesus and tenuis, and those which fall outside of the one or the other, this can be done by the use of colors, symbols, 01 a diagram designed to show the directions of differences, as well as their amounts. Such a diagram can be prepared by employing all of the four quadrants of analytic geometry, instead of only the first quadrant as in the previous prob- lems, wherein the directions of differences were ignored. As before, we plot differences from obesus as abscissas and differences from tenuis as ordinates, but in addition we follow the standard practices of analytic geometry, plot- tmg positive values to the right and above the origin, and negative values to the left and below. Employing the usual designation of the quadrants, the first being that to the upper right of the origin, the second upper left, the third lower left, and the fourth lower right, we note that the first quadrant would contain points representing characters in which a is higher than either obesus or tenuis, the second quadrant, points in which a is higher than tenuis but lower than obesus, the third, points in which a is below either, and the fourth, characters in which a is below tenuis but higher than obesus. It is evident that when a lies between tenuis and obesus in any character, the co-ordinate representing that character will fall in the second or fourth quadrant, depending on whether tenuis or obesus is the higher. If a point falls in the first or third quadrant, .a lies outside the range between obesus and tenuis. While this type of diagram has the advantage of indicating directional trends, its evaluation, in judging the cumulative effect of all characters to determine the relative closeness of a to obesus and to tenuis, is not so easy as with the simpler single-quadrant method. In this more elaborate directional analysis we must draw two limiting diagonals. Characters falling within the range of 45 to 135 and 225° to 315° indicate a closer relationship of a to obesus, while points lying between 13 5° to 225°, and from 315 through zero to 45" show a to be closer to tenuis. This may be stated more simply, but less definitely, by saying that, if the character points of a tend to adhere closely to the vertical axis, the re- lationship of a is nearer to obesus; whereas, if the majority are nearer the horizontal axis, a is more like tenuis. The principal advantage of this type of analysis is that it will indicate whether, in the majority of characters, a is intermediate between obesus and tenuis or falls outside the means of both. Such a type of directional presentation is thought particularly suit- able for a study of hybrids. It also suggests the matter of diagrams to illustrate dines, a complicated problem beyond the scope of the present dis- cussion. 76 Bulletin 18: Zoological Society of San Diego Bibliography Fisher, R. A. 1932. Statistical Methods for Research Workers, Fourth Edition. Edin- burgh. Section 21.1. (Will be found in subsequent editions under the same section number. In the eighth edition, 1941, pp. 97-99). 1936. The Use of Multiple Measurements in Taxonomic Problems. Annals of Eugenics, Vol. 7, Part 2, pp. 179-18 8. Ginsburg, Isaac 193 8. Arithmetical Definition of the Species, Subspecias and Race Con- cept, with a Proposal for a Modified Nomenclature. Containing a Simple Method for the Comparison of Related Populations. Zoo- logica, Vol. 23, pp. 253-286. 1939. The Measure of Population Divergence and Multiplicity of Characters. Journal Washington Academy of Sciences. Vol. 29, pp. 317-330. 1940. Divergence and Probability in Taxonomy. Zoologica, Vol. 2 5, pp. 15-31. Pearson, Karl 1926. On the Coefficient of Racial Likeness. Biometrica, Vol. 18, pp. 105-117. Seltzer, Carl C. 1937. A Critique of the Coefficient of Racial Likeness. American Journal of Physical Anthropology, Vol. 20, pp. 101-109. Tippett, L. FT C. 1937. The Methods of Statistics, Second Edition. London, pp. 102-104. Wallis, W. Allen 1942. Compounding Probabilities from Independent Significance Tests. Econometrica, Vol. 10, pp. 229-248.