首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 15 毫秒
1.
广西英罗港红树植物群落的非线性排序   总被引:1,自引:1,他引:1  
梁士楚  张炜银 《广西植物》2001,21(3):228-232
采用主分量分析 (PCA)、无偏主分量 (DPC)和非度量多维调节 (NMDS)等方法对广西英罗港 2 2个红树植物群落样地进行了排序。PCA和 DPC分析结果表明 ,取样数据具有明显的非线性结构。通过 NMDS分析 ,得到二维 NMDS排序格局 ,它能较好地反映了红树植物群落与环境因子之间的相互关系。  相似文献   

2.
Accuracy of estimated phylogenetic trees from molecular data   总被引:27,自引:0,他引:27  
The accuracies and efficiencies of three different methods of making phylogenetic trees from gene frequency data were examined by using computer simulation. The methods examined are UPGMA, Farris' (1972) method, and Tateno et al.'s (1982) modified Farris method. In the computer simulation eight species (or populations) were assumed to evolve according to a given model tree, and the evolutionary changes of allele frequencies were followed by using the infinite-allele model. At the end of the simulated evolution five genetic distance measures (Nei's standard and minimum distances, Rogers' distance, Cavalli-Sforza's f theta, and the modified Cavalli-Sforza distance) were computed for all pairs of species, and the distance matrix obtained for each distance measure was used for reconstructing a phylogenetic tree. The phylogenetic tree obtained was then compared with the model tree. The results obtained indicate that in all tree-making methods examined the accuracies of both the topology and branch lengths of a reconstructed tree (rooted tree) are very low when the number of loci used is less than 20 but gradually increase with increasing number of loci. When the expected number of gene substitutions (M) for the shortest branch is 0.1 or more per locus and 30 or more loci are used, the topological error as measured by the distortion index (dT) is not great, but the probability of obtaining the correct topology (P) is less than 0.5 even with 60 loci. When M is as small as 0.004, P is substantially lower. In obtaining a good topology (small dT and high P) UPGMA and the modified Farris method generally show a better performance than the Farris method. The poor performance of the Farris method is observed even when Rogers' distance which obeys the triangle inequality is used. The main reason for this seems to be that the Farris method often gives overestimates of branch lengths. For estimating the expected branch lengths of the true tree UPGMA shows the best performance. For this purpose Nei's standard distance gives a better result than the others because of its linear relationship with the number of gene substitutions. Rogers' or Cavalli-Sforza's distance gives a phylogenetic tree in which the parts near the root are condensed and the other parts are elongated. It is recommended that more than 30 loci, including both polymorphic and monomorphic loci, be used for making phylogenetic trees. The conclusions from this study seem to apply also to data on nucleotide differences obtained by the restriction enzyme techniques.  相似文献   

3.
Effectively summarizing complex community relationships is an important feature in studies such as biodiversity, global change, and invasion ecology. The reliability of such community summaries depends on the degree of sampling variability that is present in the data, the structure of the data, and the choice of ordination method, but the relative importance of these factors is not understood. We compared the validity of results from different ordination methods by applying five levels of sampling error to a simulated coenoplane model at two gradient lengths using two types of data (abundance and presence–absence). The multivariate methods we compared were correspondence analysis (CA), detrended correspondence analysis (DCA), non-metric multidimensional scaling (NMDS), principal component analysis (PCA) and principal coordinates analysis (PCoA). Our results showed CA and PCA using presence–absence data were the most successful methods regardless of sampling error and gradient length, closely followed by the other methods using presence–absence data. With abundance data, PCA and CA were the most successful approaches with the short and long gradients, respectively. Approaches based on PCoA and NMDS using abundance data did not perform well regardless of the choice of distance measure used in the analysis. Both of these methods, along with the PCA using abundance data, were strongly affected by the longer gradient, leading to more distorted results.  相似文献   

4.
《Ecological Informatics》2007,2(2):138-149
Ecological patterns are difficult to extract directly from vegetation data. The respective surveys provide a high number of interrelated species occurrence variables. Since often only a limited number of ecological gradients determine species distributions, the data might be represented by much fewer but effectively independent variables. This can be achieved by reducing the dimensionality of the data. Conventional methods are either limited to linear feature extraction (e.g., principal component analysis, and Classical Multidimensional Scaling, CMDS) or require a priori assumptions on the intrinsic data dimensionality (e.g., Nonmetric Multidimensional Scaling, NMDS, and self organized maps, SOM).In this study we explored the potential of Isometric Feature Mapping (Isomap). This new method of dimensionality reduction is a nonlinear generalization of CMDS. Isomap is based on a nonlinear geodesic inter-point distance matrix. Estimating geodesic distances requires one free threshold parameter, which defines linear geometry to be preserved in the global nonlinear distance structure. We compared Isomap to its linear (CMDS) and nonmetric (NMDS) equivalents. Furthermore, the use of geodesic distances allowed also extending NMDS to a version that we called NMDS-G. In addition we investigated a supervised Isomap variant (S-Isomap) and showed that all these techniques are interpretable within a single methodical framework.As an example we investigated seven plots (subdivided in 456 subplots) in different secondary tropical montane forests with 773 species of vascular plants. A key problem for the study of tropical vegetation data is the heterogeneous small scale variability implying large ranges of β-diversity. The CMDS and NMDS methods did not reduce the data dimensionality reasonably. On the contrary, Isomap explained 95% of the data variance in the first five dimensions and provided ecologically interpretable visualizations; NMDS-G yielded similar results. The main shortcoming of the latter was the high computational cost and the requirement to predefine the dimension of the embedding space. The S-Isomap learning scheme did not improve the Isomap variant for an optimal threshold parameter but substantially improved the nonoptimal solutions.We conclude that Isomap as a new ordination method allows effective representations of high dimensional vegetation data sets. The method is promising since it does not require a priori assumptions, and is computationally highly effective.  相似文献   

5.
Abstract: Previous ordination studies of land snail community composition have been limited to four or fewer habitat types from sites separated by no more than 300 km. To investigate the nature of large-scale patterns, North American land snail assemblages at 421 sites, representing 26 habitat types and covering a 1400 × 800 km area, were ordinated using global, nonmetric multi-dimensional scaling (NMDS). These data were then subjected to model-based cluster analysis and kmeans clustering to identify the main compositional groups and most important environmental covariables. Six primary compositional groups were identified. Three of these largely represent upland forest and rock outcrop sites, while the remaining largely represent either lowland forest, lowland grassland or upland grassland habitats. The geographical location and moisture level of sites also influences community composition. A strong compositional difference exists between sites having duff vs. turf soil surface layers. Only 8% of sites were improperly classified when soil surface architecture was used as the sole predictor variable. Fully 43% of taxa exhibited significant preferences towards one of these surface types, while only 15% of relatively common (10 + occurrence) taxa showed no preferences. Twelve groups of closely related taxa within the same genus had members that favoured different surface types, indicating that differential selection pressures have existed over evolutionary time scales. While turf faunas appeared unaffected by anthropogenic disturbance, duff faunas were strongly impacted, suggesting that their conservation will require protection of soil surface architecture.  相似文献   

6.
In 1978 and 1999, seven and eight species of Mytilidae (Mollusca: Bivalvia) were analyzed using gel electrophoresis. Mean heterozygosity per individual (Hobs and Hexp) and genetic distances (Rogers' DR, Nei's DN, and others) were estimated for 21 and 24 allozyme loci. Mytilus modiolus had the highest variation among the species examined. Genetic distances were lowest for the M. trossulus-M. galloprovincialis species pair: DR = 0.147, DN = 0.078. Overall, five species of the genera Mytilus and Crenomytilus were genetically closer to each other (DR = 0.147, DN = 0.078) than to the remaining three species of this group (DR = 0.807, DN = 2.243). The relationships among the species were examined using cluster analysis and parsimony methods. The densest clusters in the dendrograms consisted of (1) M. trossulus and M. galloprovincialis and (2) M. coruscus, M. californianus, and M. grayanus. These two clusters form a larger cluster (3), which comprises all representatives of the nominal genus Mutilus and C. grayanus. The Mytilus-Crenomytilus cluster is consecutively joined by Adula falcatoides, Mytilus modiolus, and Septifer keenae. According to Nei's genetic distances DN, the time of divergence between M. trossulus and M. galloprovincialis is 0.8-1.6 Myr; between M. californianus and C. grayanus, it is approximately 9 Myr; and between M. coruscus and the latter pair, it is 13 Myr before present. Two representatives of the Mytilus ex gr. edulis complex diverged from the Mytilus-Crenomytilus group of large-size Pacific species about 20 Myr ago. These results are in good agreement with paleontological data and indicate a relatively recent origin of the Mytilus ex gr. edulis complex. The results obtained can be used in systematics and phylogeny of modern Mytilidae.  相似文献   

7.
We conducted a 55-day long camera-trap survey in the Yungas subtropical forest in NW Argentina, to assess the effect of human accessibility, conservation status of the area, domestic animals and elevation on the diversity and composition of the large and medium-sized native terrestrial mammal assemblage. We deployed 24 camera-trap stations at distances of ~2 km from each other. The study area is covered by continuous forest and has its center in the small community of Acambuco, in the Acambuco Provincial Reserve. The main economic activity in the area is oil/gas exploitation. Local residents raise cattle, hunt and use timber and non-timber forest products. The human impact was indirectly measured with an accessibility cost model. We used a multiple regression ANCOVA to assess the effect of elevation (range: 628–1170 masl), accessibility, protection status (reserve vs not) and frequency of records of domestic animals on the native mammal species richness and on a nonmetric multidimensional scaling (NMDS) ordination based on the frequency of records of the native mammals recorded at >3 camera-trap stations. We recorded 15 species of native mammals. Native mammal species richness decreased with elevation. Elevation was correlated with NMDS axes. Other predictive variables had no effect on species richness or the NMDS ordination, probably as a result of the relatively narrow range of conditions assessed in this study. The effect of elevation on mammal assemblages should be considered in landscape planning processes aimed at promoting biodiversity conservation.  相似文献   

8.
利用AFLP技术对26个竹子种类进行了多样性分析, 以探索引物组合数量对准确研究竹子类群系统关系的影响。实验共随机选取10对AFLP引物, 并对所得10组AFLP标记数据随机组合后进行Nei氏遗传距离/UPGMA聚类分析。每对AFLP引物 扩增数据为一组, 随着用于聚类统计的AFLP标记数据随机组合数量的增加, 26个竹子种类的聚类关系趋向一致。这提示我们,在系统学研究中, 足够数量的引物组合是获得供试材料间准确聚类关系的基础, 应采用对各AFLP引物组合数据随机累加后进行聚类分析的方法, 以聚类关系为标准来确定用于分析供试品种的最少引物组合数量。  相似文献   

9.
The data obtained using electrophoresis strongly support the specific status of Atherina boyeri Risso, 1810 and A. presbyfer Cuvier, 1829, and thuscontradict the recently proposedsynonymyof the two species. Four populations of A. boyeri and six populations of A. presbyter were assayed for 11 enzymes and general protein using muscle and liver extracts. Eight of the 11 enzymes were shown to be polymorphic at the 95% level. Sixteen loci, encoding 40 putative alleles were consistently resolved in all 10 populations.
The two species were fixed for different alleles at the EST-3 locus. At the G3PDH locus. with the exception of two heterozygotes, all individuals of each species were also homozygous for different alleles. At the PGM locus the common allele was unique to each species.
The mean Nei's genetic distance ( ), over all loci, calculated between populations of A. boyeri (= 0.10 ± 0.06); between populations of A. presbyter ( = 0.02 0.02) and between populations of A. boyeri and A. presbyter ( = 0.42 0.09) indicated the separateness of the two species.
UPGMA cluster analysis based on genetic distances produced a dendrogram whose principal dichotomy resulted in the formation of two clusters. The ordination of populations in the UPGMA cluster analysis strongly reflected the geographic distribution of populations in both species.  相似文献   

10.
Abstract. This article investigates whether the Braun‐Blanquet abundance/dominance (AD) scores that commonly appear in phytosociological tables can properly be analysed by conventional multivariate analysis methods such as Principal Components Analysis and Correspondence Analysis. The answer is a definite NO. The source of problems is that the AD values express species performance on a scale, namely the ordinal scale, on which differences are not interpretable. There are several arguments suggesting that no matter which methods have been preferred in contemporary numerical syntaxonomy and why, ordinal data should be treated in an ordinal way. In addition to the inadmissibility of arithmetic operations with the AD scores, these arguments include interpretability of dissimilarities derived from ordinal data, consistency of all steps throughout the analysis and universality of the method which enables simultaneous treatment of various measurement scales. All the ordination methods that are commonly used, for example, Principal Components Analysis and all variants of Correspondence Analysis as well as standard cluster analyses such as Ward's method and group average clustering, are inappropriate when using AD data. Therefore, the application of ordinal clustering and scaling methods to traditional phytosociological data is advocated. Dissimilarities between relevés should be calculated using ordinal measures of resemblance, and ordination and clustering algorithms should also be ordinal in nature. A good ordination example is Non‐metric Multidimensional Scaling (NMDS) as long as it is calculated from an ordinal dissimilarity measure such as the Goodman & Kruskal γ coefficient, and for clustering the new OrdClAn‐H and OrdClAn‐N methods.  相似文献   

11.
The Agave deserti complex, comprising A. deserti, A. cerulata and A. subsimplex, represents a group of species and subspecies with a near allopatric distribution and clear differences in morphology. Genetic differentiation and taxonomic status with respect to spatial distribution of 14 populations of the complex were analyzed in an effort to understand the evolution and speciation process within the genus. Allelic frequencies, levels of genetic variation, expected heterozygosity (H(S)), proportion of polymorphic loci (P), and genetic differentiation (theta and Nei's genetic distance) were estimated using 41 putative RAPD loci. All three species show high levels of genetic variation (H(S)=0.12-0.29, P=63.4-95.1), and low genetic differentiation between populations and species (theta populations=0.14+/-0.02 (SE); G(st)=0.11+/-0.02). Accordingly, gene flow among populations was estimated as high by three different methods (N(m)=2.91-6.14). Nei's genetic distances between the three species were low compared to the values obtained from other Agavaceae, and there was no clear correlation with taxonomic divisions. In a UPGMA analysis, A. subsimplex and A. cerulata formed exclusive monospecific clusters, whereas the A. deserti populations appear in more than one cluster together with other species. The results were consistent with a pattern of genetic isolation by distance.  相似文献   

12.
Variation at four highly polymorphic allozyme loci (inorganic pyrophosphatase, peptidase, and two esterase loci) was examined in 25 settlements of the marine snail Littorina sitkana (Mollusca, Gastropoda). The sampling localities covered a wide part of the species range: from the Peter the Great Bay (the Sea of Japan) at the southwest to the Mednyi Island (Commander Islands) at the northeast. Like other littorines lacking the pelagic stage, L. sitkana was characterized by significant genetic differentiation (G(ST) for the pooled sample was 0.310). Cluster analysis and nonmetric multidimensional scaling conducted on a matrix of pairwise genetic distances between all of the settlements studied revealed four genetically different groups: southern Primorye, northern Prymorye, Sakhalin, and Kuril-Commanders. The population-genetic structure of the L. sitkana settlements is close to that described by the isolation-by-distance and stepping-stone models: the geographic and the genetic distances between the most settlements examined are distinctly correlated.  相似文献   

13.
We describe the mesoscale floristic patterns in the central Western Ghats of Karnataka, India, through combined analysis of woody species abundance and stand structure data from a network of ninety-six 1-ha sampling plots spread across 22,000 km2. A total of 61,906 individuals (≥10 cm gbh) comprising 400 plant species from 254 genera and 75 families were recorded. Euphorbiaceae, Rubiaceae, Lauraceae and Moraceae families constituted 23.5 percent of the total number of species encountered. The relative dominance of species was skewed with Poecilonueron indicum, Xylia xylocarpa, Terminalia tomentosa and Anogeissus latifolia being dominant in some plots. Correspondence analysis (CA) and a nonmetric multidimensional scaling (NMDS) of plots by species abundances data showed similar arching patterns, with significant correlation between the first axis of CA and NMDS (r=0.77). Hierarchical clustering of plot scores along the three first CA axes resulted in splitting the plots into five different categories that broadly reflect the major bioclimatic features of the region. A multiscale bootstrapping test indicated that categorization of the wettest (wet evergreen group 1 and 2) and driest (dry deciduous) groups were robust (P<0.05 with 1000 bootstraps), while the remaining two transitional groups were uncertain (P=0.12 and 0.26 for moist deciduous and semi-evergreen group, respectively). Principal component analysis revealed that plots with similar floristic composition can encompass contrastingly different physiognomic structures (canopy cover, canopy height and mean tree diameter) probably in relation to their levels of disturbance. Observed patterns in the floristic composition have been discussed in the light of the complex interaction between the bioclimatic and disturbance regimes that characterize the region.  相似文献   

14.
Aim We used modern pollen assemblages to develop a method for climate reconstruction that reduces the spatial autocorrelation of residuals and accounts for the strong topographic and climatic variation that occurs in British Columbia, Canada. Location British Columbia, Canada, including sites both on the mainland and on adjacent islands (Queen Charlotte Islands and Vancouver Island). Methods New pollen assemblages from surface‐sediment samples collected in British Columbia were combined with other published and unpublished samples (n = 284). Multivariate rank‐distances between sample sites and a randomized set of sites within the province were calculated for climate parameters to determine whether gaps in the current network of present‐day pollen sample sites exist. Lacustrine surface‐sediment pollen assemblages (n = 145) were ordinated using non‐metric multidimensional scaling (NMDS), and a generalized additive model (GAM) was used to reconstruct modern mean warmest month temperature (MWMT) and mean annual precipitation (MAP) from the NMDS ordination. The results were compared with standard climate reconstruction techniques, including the modern analogue technique, partial least squares, weighted averaging, weighted averaging–partial least squares and factor analysis. Results Reconstructions of MWMT and MAP using NMDS and GAM were comparable to those of existing models. When reconstructing MWMT, the NMDS/GAM method had a lower root‐mean‐squared error of prediction (RMSEP), lower spatial autocorrelation and higher correlation with observed temperature values than the other methods tested. When reconstructing MAP, the partial least squares method performed better than the NMDS/GAM method for RMSEP and correlation with observed values; however, the NMDS/GAM method had a lower spatial autocorrelation of residuals. Main conclusions NMDS reveals strong relationships among modern pollen assemblages, vegetation and climate parameters. Climate models using NMDS and GAM are comparable to other palaeoecological reconstruction models, but provide lower spatial autocorrelation of residuals for both parameters tested. An inverse distance‐weighted surface of multivariate rank‐climate distances generated from the network of pollen sample sites indicates that greater sampling intensity in north‐western and central‐interior British Columbia is required in order to obtain an accurate representation of climatic and vegetation diversity in the province.  相似文献   

15.
Five samples (30 individuals each) of the Mediterranean cyprinodontid fish, Aphanius fasciatus, were analysed by allozyme electrophoresis for a study concerning genetic structure and gene flow among populations. Forty-three loci, twelve of which were polymorphic at P0.99 level, were analysed. A. fasciatus showed low levels of genetic polymorphism, with expected heterozygosity values ranging from 0.027 (SE = 0.013) to 0.064 (SE = 0.023). The modified Rogers' genetic distances ranged from 0.053 to 0.202. The observed pattern of genetic differentiation among populations, obtained by multidimensional scaling of the modified Rogers' genetic distances, was consistent with geographic distribution. Coancestry coefficient ( = 0.302, SE = 0.045) indicated a very high degree of genetic subdivision among populations and the estimate of gene flow (the effective number of migrants per generation, Nm = 0.445), computed by Wright's method, showed very restricted gene flow among populations. The regression analysis carried out using log-transformed values of geographic distance against M^ (pairwise values of Nm) gave a regression coefficient, b = –1.03, indicating that a restricted migration rate among populations can occur according to the one-dimensional stepping-stone model. The amount and modality of gene flow in A. fasciatus inferred from the results of the present study are consistent with the absence of dispersal stages and with the fragmented nature of brackish-water habitats along the coasts.  相似文献   

16.
Randomly amplified polymorphic DNA (RAPD) markers have been used to study the genetic variation among androgenetic monoploids of diploid Solanum species. Cluster analysis of pairwise genetic distances was used to construct a genetic relationship among anther donor and anther-derived potato plants. The clustering based on Rogers' distances resembled classifications based on parental origins and hybrid combinations. Six of the 32 RAPD primers used resulted in the selective amplification of DNA fragments which were polymorphic between the two S. phureja parental clones, 1.22 and A95. It should be possible to construct a genetic linkage map, without making crosses, using monoploids derived from a single heterozygous diploid clone and RAPD markers.  相似文献   

17.
Genetic variation in four indigenous chicken breeds from the Veneto region of Italy was assessed using amplified fragment length polymorphism (AFLP) markers. A total of 99 individuals were analysed using three AFLP primer combinations that produced 70 polymorphisms. Four indigenous Veneto chicken breeds (Ermellinata, Padovana, Pépoi and Robusta) and a reference broiler line were included in the analysis. Breed-specific markers were identified in each breed. The expected heterozygosity did not differ significantly among the indigenous Veneto chicken breeds and the broiler line. The coefficient of gene variation (Gst) value across loci indicated that almost half of the total variability was observed among breeds. Nei's standard genetic distance between pairs of breeds showed that the distance between the broiler line and the Pépoi breed was greater than the distances between the broiler line and the other three chicken breeds. Cluster analysis based on standard genetic distances between breeds indicated that the Padovana and Pépoi breeds were closely related. Factorial analysis based on a binary matrix of the AFLP data showed a clear distinction of all breeds.  相似文献   

18.
Twenty-seven microsatellite loci were used to define genetic variation and relationships among eight Indian riverine buffalo breeds. The total number of alleles ranged from 166 in the Toda breed to 194 each in the Mehsana and the Murrah. Significant departures from the Hardy-Weinberg equilibrium were observed for 26 locus-breed combinations due to heterozygote deficiency. Breed differentiation was analysed by estimation of F(ST) index (values ranging from 0.75% to 6.00%) for various breed combinations. The neighbour-joining tree constructed from chord distances, multidimensional scaling (MDS) display of F(ST) values and Bayesian clustering approach consistently identified the Toda, Jaffarabadi, and Pandharpuri breeds as one lineage each, and the Bhadawari, Nagpuri, Surati, Mehsana and Murrah breeds as admixture. Analysis of molecular variance refuted the earlier classification of these breeds proposed on the basis of morphological and geographical parameters. The Toda buffaloes, reared by a tribe of the same name, represent an endangered breed from the Nilgiri hills in South India. Divergence time of the Toda buffaloes from the other main breeds, calculated from Nei's standard genetic distances based on genotyping data on seven breeds and 20 microsatellite loci, suggested separation of this breed approximately 1800-2700 years ago. The results of the present study will be useful for development of rational breeding and conservation strategies for Indian buffaloes.  相似文献   

19.
微卫星DNA标记探讨镜鲤的种群结构与遗传变异   总被引:17,自引:2,他引:15  
全迎春  李大宇  曹鼎辰  孙效文  梁利群 《遗传》2006,28(12):1541-1548
采用30个微卫星分子标记, 对5个镜鲤群体的观测杂合度(Ho)、期望杂合度(He)、多态信息含量(PIC)和有效等位基因数(Ae)等进行了遗传检测, 根据基因频率计算遗传相似系数和Nei氏标准遗传距离, 以c2检验估计Hardy-Weinberg平衡, 以近交系数(FST)和基因流(Nm)分析群体的遗传分化。同时, 使用PHYLIP3.63软件绘制基于Nei氏标准遗传距离的UPGMA聚类图, 并进行bootstrap自举检验验证进化树的可靠性。在德国镜鲤选育系(Scattered Cyprinus carpio L.)和来自4个不同养殖场(松浦、东岗、奉城和辽中)的德国镜鲤群体中共检测到7 083个扩增片段, 长度在102 ~ 446 bp之间, 在群体内扩增出等位基因1~16个不等, 共计356个等位基因。结果表明: (1)5个群体检测的有效等位基因数在1.07~12.30个不等, 平均多态信息含量为0.74、0.74、0.69、0.75和0.75, 无偏期望杂合度的平均值为0.74、0.78、0.70、0.76和0.78, 说明这几个群体属于高度多态, 遗传多样性水平较高。(2)群体间相似系数在0.52以上, 相似性较高。聚类分析显示, 东岗、奉城和辽中3个养殖场的德国镜鲤群体聚类成一个分支, 而德国镜鲤选育系与松浦群体聚类成另一分支。聚类的先后与它们在地理分布上距离远近有一定的相关性。(3)在与功能基因相关的多个微卫星基因座位上, 扩增产物呈现不同程度的缺失现象, 这些无效等位基因的产生可能与结构基因在育种中受到人工选择的影响较大有关。  相似文献   

20.
设为首页 | 免责声明 | 关于勤云 | 加入收藏

Copyright©北京勤云科技发展有限公司  京ICP备09084417号