Summary Inspired by the MaCzek Visual Basic program we pprovide an R package, RMaCzek, that produces Czekanowski’s diagrams. Our package permits any seriation and distance method the user provides. In this paper we focus on the “OLO” and “QAP_2SUM” methods from the seriation package. We illustrate the possibilities of our package with three anthropological studies, one socioeconomic study and a phylogenetically motivated simulation study.
{"title":"“Old Techniques for New Times”: the RMaCzek package for producing Czekanowski’s diagrams","authors":"K. Bartoszek, Albin Västerlund","doi":"10.2478/bile-2020-0008","DOIUrl":"https://doi.org/10.2478/bile-2020-0008","url":null,"abstract":"Summary Inspired by the MaCzek Visual Basic program we pprovide an R package, RMaCzek, that produces Czekanowski’s diagrams. Our package permits any seriation and distance method the user provides. In this paper we focus on the “OLO” and “QAP_2SUM” methods from the seriation package. We illustrate the possibilities of our package with three anthropological studies, one socioeconomic study and a phylogenetically motivated simulation study.","PeriodicalId":8933,"journal":{"name":"Biometrical Letters","volume":"6 1","pages":"89 - 118"},"PeriodicalIF":0.0,"publicationDate":"2020-12-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"87652234","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Achilleas Anastasiou, A. Karagrigoriou, A. Katsileros
Summary The normal distribution is considered to be one of the most important distributions, with numerous applications in various fields, including the field of agricultural sciences. The purpose of this study is to evaluate the most popular normality tests, comparing the performance in terms of the size (type I error) and the power against a large spectrum of distributions with simulations for various sample sizes and significance levels, as well as through empirical data from agricultural experiments. The simulation results show that the power of all normality tests is low for small sample size, but as the sample size increases, the power increases as well. Also, the results show that the Shapiro–Wilk test is powerful over a wide range of alternative distributions and sample sizes and especially in asymmetric distributions. Moreover the D’Agostino–Pearson Omnibus test is powerful for small sample sizes against symmetric alternative distributions, while the same is true for the Kurtosis test for moderate and large sample sizes.
{"title":"Comparative evaluation of goodness of fit tests for normal distribution using simulation and empirical data","authors":"Achilleas Anastasiou, A. Karagrigoriou, A. Katsileros","doi":"10.2478/bile-2020-0015","DOIUrl":"https://doi.org/10.2478/bile-2020-0015","url":null,"abstract":"Summary The normal distribution is considered to be one of the most important distributions, with numerous applications in various fields, including the field of agricultural sciences. The purpose of this study is to evaluate the most popular normality tests, comparing the performance in terms of the size (type I error) and the power against a large spectrum of distributions with simulations for various sample sizes and significance levels, as well as through empirical data from agricultural experiments. The simulation results show that the power of all normality tests is low for small sample size, but as the sample size increases, the power increases as well. Also, the results show that the Shapiro–Wilk test is powerful over a wide range of alternative distributions and sample sizes and especially in asymmetric distributions. Moreover the D’Agostino–Pearson Omnibus test is powerful for small sample sizes against symmetric alternative distributions, while the same is true for the Kurtosis test for moderate and large sample sizes.","PeriodicalId":8933,"journal":{"name":"Biometrical Letters","volume":"47 1","pages":"237 - 251"},"PeriodicalIF":0.0,"publicationDate":"2020-12-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"90498456","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Summary A field study was conducted at the Department of Agronomy of Poznań University of Life Sciences to determine the effect of the depth of NP fertilizer application in maize cultivation on the dynamics of initial maize growth, expressed in dry matter of a single plant at two juvenile maize stages. The adopted assumptions were verified on the basis of a four-year field experiment using four depths of NP fertilizer application, two nitrogen fertilizers and two nitrogen dose application dates. Thermal conditions in the early maize growing season had a significant effect on maize response to the depth of application of a phosphorus starting dose. Row fertilization (regardless of the depth of fertilizer application) was more effective than broadcast fertilization at both studied developmental stages.
{"title":"Assessment of the impact of NP fertilizer application depth on the rate of initial dry matter accumulation of maize (Zea mays L.)","authors":"K. Ambroży-Deręgowska, P. Szulc, I. Mejza","doi":"10.2478/bile-2020-0016","DOIUrl":"https://doi.org/10.2478/bile-2020-0016","url":null,"abstract":"Summary A field study was conducted at the Department of Agronomy of Poznań University of Life Sciences to determine the effect of the depth of NP fertilizer application in maize cultivation on the dynamics of initial maize growth, expressed in dry matter of a single plant at two juvenile maize stages. The adopted assumptions were verified on the basis of a four-year field experiment using four depths of NP fertilizer application, two nitrogen fertilizers and two nitrogen dose application dates. Thermal conditions in the early maize growing season had a significant effect on maize response to the depth of application of a phosphorus starting dose. Row fertilization (regardless of the depth of fertilizer application) was more effective than broadcast fertilization at both studied developmental stages.","PeriodicalId":8933,"journal":{"name":"Biometrical Letters","volume":"1 1","pages":"253 - 260"},"PeriodicalIF":0.0,"publicationDate":"2020-12-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"76035213","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Malwina Janiszewska, A. Szczepańska-Álvarez, Emilia E. Zawieja
Summary The aim of this paper is to study the association between two random vectors related to two groups of characteristics. To analyze the multivariate association, the ρV coefficient and distance correlation are used. Two methods (classical and recent) are compared and illustrated with real data.
{"title":"Application of the ρV coefficient and distance correlation to the analysis of multivariate association","authors":"Malwina Janiszewska, A. Szczepańska-Álvarez, Emilia E. Zawieja","doi":"10.2478/bile-2020-0009","DOIUrl":"https://doi.org/10.2478/bile-2020-0009","url":null,"abstract":"Summary The aim of this paper is to study the association between two random vectors related to two groups of characteristics. To analyze the multivariate association, the ρV coefficient and distance correlation are used. Two methods (classical and recent) are compared and illustrated with real data.","PeriodicalId":8933,"journal":{"name":"Biometrical Letters","volume":"1 1","pages":"119 - 129"},"PeriodicalIF":0.0,"publicationDate":"2020-12-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"90161141","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
P. Fasahat, M. Aghaeezadeh, M. Kakueinezhad, Leila Jabbari
Summary The evaluation of sugar beet genotypes under different climate conditions is a principal goal of breeding programs. In most studies, environment has a high influence on the qualitative and quantitative traits of sugar beet. Therefore, data collected from different environments may contribute to more accurate genotype selection. In this study, the effect of different environments on sugar beet genotypes’ performance was evaluated using a meta-analysis method based on Hedges’ technique. Data were collected from 149 trials conducted in twelve regions in Iran over 15 years (2003–18). For all trials, the value of the traits was weighted, and subsequently the effect size, reaction ratio and confidence interval were estimated. Among the studied environments, Khoy had a positive effect on root yield, sugar content, sugar yield and white sugar yield. As could be expected, the effect of environment on final yield formation was high, so that the Shiraz environment had a negative effect on root yield and sugar yield. Overall, the ranking of environments based on the meta-analysis results was quite different from that obtained by comparison of mean results.
{"title":"A meta-analysis of genotype × environment interaction on sugar beet performance","authors":"P. Fasahat, M. Aghaeezadeh, M. Kakueinezhad, Leila Jabbari","doi":"10.2478/bile-2020-0014","DOIUrl":"https://doi.org/10.2478/bile-2020-0014","url":null,"abstract":"Summary The evaluation of sugar beet genotypes under different climate conditions is a principal goal of breeding programs. In most studies, environment has a high influence on the qualitative and quantitative traits of sugar beet. Therefore, data collected from different environments may contribute to more accurate genotype selection. In this study, the effect of different environments on sugar beet genotypes’ performance was evaluated using a meta-analysis method based on Hedges’ technique. Data were collected from 149 trials conducted in twelve regions in Iran over 15 years (2003–18). For all trials, the value of the traits was weighted, and subsequently the effect size, reaction ratio and confidence interval were estimated. Among the studied environments, Khoy had a positive effect on root yield, sugar content, sugar yield and white sugar yield. As could be expected, the effect of environment on final yield formation was high, so that the Shiraz environment had a negative effect on root yield and sugar yield. Overall, the ranking of environments based on the meta-analysis results was quite different from that obtained by comparison of mean results.","PeriodicalId":8933,"journal":{"name":"Biometrical Letters","volume":"44 1","pages":"221 - 236"},"PeriodicalIF":0.0,"publicationDate":"2020-12-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"89071547","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
I. Mejza, K. Ambroży-Deręgowska, K. Ozawa, S. Mejza, Shinji Kuriki
Summary We consider a new method of constructing non-orthogonal (incomplete) split-split-plot designs (SSPDs) for three (A, B, C) factor experiments. The final design is generated by some resolvable incomplete block design (for the factor A) and by square lattice designs for factors B and C using a modified Kronecker product of those designs (incidence matrices). Statistical properties of the constructed designs are investigated under a randomized-derived linear model. This model is strictly connected with a four-step randomization of units (blocks, whole plots, subplots, sub-subplots inside each block). The final SSPD has orthogonal block structure (OBS) and satisfies the general balance (GB) property. The statistical analysis of experiments performed in the SSPD is based on the analysis of variance often used for multistratum experiments. We characterize the SSPD with respect to the stratum efficiency factors for the basic estimable treatment contrasts. The structures of the vectors defining treatment contrasts are also given.
{"title":"Constructing non-orthogonal split-split-plot designs using some resolvable block designs","authors":"I. Mejza, K. Ambroży-Deręgowska, K. Ozawa, S. Mejza, Shinji Kuriki","doi":"10.2478/bile-2020-0012","DOIUrl":"https://doi.org/10.2478/bile-2020-0012","url":null,"abstract":"Summary We consider a new method of constructing non-orthogonal (incomplete) split-split-plot designs (SSPDs) for three (A, B, C) factor experiments. The final design is generated by some resolvable incomplete block design (for the factor A) and by square lattice designs for factors B and C using a modified Kronecker product of those designs (incidence matrices). Statistical properties of the constructed designs are investigated under a randomized-derived linear model. This model is strictly connected with a four-step randomization of units (blocks, whole plots, subplots, sub-subplots inside each block). The final SSPD has orthogonal block structure (OBS) and satisfies the general balance (GB) property. The statistical analysis of experiments performed in the SSPD is based on the analysis of variance often used for multistratum experiments. We characterize the SSPD with respect to the stratum efficiency factors for the basic estimable treatment contrasts. The structures of the vectors defining treatment contrasts are also given.","PeriodicalId":8933,"journal":{"name":"Biometrical Letters","volume":"17 1","pages":"177 - 194"},"PeriodicalIF":0.0,"publicationDate":"2020-12-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"86185941","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Summary This paper provides estimation and hypothesis testing procedures for experiments in split-plot designs. These experiments have been shown to have a convenient orthogonal block structure when properly randomized. Due to this property, the analysis of experimental data can be carried out in a relatively simple manner. Relevant simplification procedures are indicated. According to the adopted approach, the analysis of variance and hypothesis testing procedures can be performed directly, rather than by combining the results of analyses based on some stratum submodels. The practical application of the presented theory is illustrated by examples of real experiments in appropriate split-plot designs. The present paper is the fourth in the planned series of publications on the analysis of experiments with orthogonal block structure.
{"title":"On a new approach to the analysis of variance for experiments with orthogonal block structure. IV. Experiments in split-plot designs","authors":"T. Caliński, A. Łacka, I. Siatkowski","doi":"10.2478/bile-2020-0011","DOIUrl":"https://doi.org/10.2478/bile-2020-0011","url":null,"abstract":"Summary This paper provides estimation and hypothesis testing procedures for experiments in split-plot designs. These experiments have been shown to have a convenient orthogonal block structure when properly randomized. Due to this property, the analysis of experimental data can be carried out in a relatively simple manner. Relevant simplification procedures are indicated. According to the adopted approach, the analysis of variance and hypothesis testing procedures can be performed directly, rather than by combining the results of analyses based on some stratum submodels. The practical application of the presented theory is illustrated by examples of real experiments in appropriate split-plot designs. The present paper is the fourth in the planned series of publications on the analysis of experiments with orthogonal block structure.","PeriodicalId":8933,"journal":{"name":"Biometrical Letters","volume":"82 1","pages":"151 - 175"},"PeriodicalIF":0.0,"publicationDate":"2020-12-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"91247511","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Summary In Bayesian hypothesis testing, evidence for a statistical model is quantified by the Bayes factor, which represents the relative likelihood of observed data under that model compared to another competing model. In general, computing Bayes factors is difficult, as computing the marginal likelihood of data under a given model requires integrating over a prior distribution of model parameters. In this paper, I capitalize on a particular choice of prior distribution that allows the Bayes factor to be expressed without integral representation, and I develop a simple formula – the Pearson Bayes factor – that requires only minimal summary statistics as commonly reported in scientific papers, such as the t or F score and the degrees of freedom. In addition to presenting this new result, I provide several examples of its use and report a simulation study validating its performance. Importantly, the Pearson Bayes factor gives applied researchers the ability to compute exact Bayes factors from minimal summary data, and thus easily assess the evidential value of any data for which these summary statistics are provided, even when the original data is not available.
{"title":"The Pearson Bayes factor: An analytic formula for computing evidential value from minimal summary statistics","authors":"Thomas J. Faulkenberry","doi":"10.2478/bile-2021-0001","DOIUrl":"https://doi.org/10.2478/bile-2021-0001","url":null,"abstract":"Summary In Bayesian hypothesis testing, evidence for a statistical model is quantified by the Bayes factor, which represents the relative likelihood of observed data under that model compared to another competing model. In general, computing Bayes factors is difficult, as computing the marginal likelihood of data under a given model requires integrating over a prior distribution of model parameters. In this paper, I capitalize on a particular choice of prior distribution that allows the Bayes factor to be expressed without integral representation, and I develop a simple formula – the Pearson Bayes factor – that requires only minimal summary statistics as commonly reported in scientific papers, such as the t or F score and the degrees of freedom. In addition to presenting this new result, I provide several examples of its use and report a simulation study validating its performance. Importantly, the Pearson Bayes factor gives applied researchers the ability to compute exact Bayes factors from minimal summary data, and thus easily assess the evidential value of any data for which these summary statistics are provided, even when the original data is not available.","PeriodicalId":8933,"journal":{"name":"Biometrical Letters","volume":"54 1","pages":"1 - 26"},"PeriodicalIF":0.0,"publicationDate":"2020-11-18","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"85399020","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Anna Szczepaniak-Kozak, E. Bakinowska, Katerina Strani
Summary This paper focuses on pragmatic competence development in second or foreign language learners. In particular, it attempts to fill the significant research gap in measuring change in pragmatic competence and capturing pragmalinguistic and sociopragmatic development over time. The paper proposes an innovative approach by applying a logistic model with multinomial distribution for measuring change in InterLanguage Pragmatics Research (ILP). Applied in the field of pragmatics, this statistical tool offers a comprehensive and flexible approach to modelling relations between independent and dependent variables in ILP research. The model is tested in a longitudinal study of Polish undergraduate students learning English, and specifically in the way they formulate requests by means of requestive directness strategies. The paper concludes that, regardless of time elapsing, the factors P (power distance) and D (social distance) have a highly significant influence on the use of requestive directness strategies by Poles learning EFL. Furthermore, the analysis indicates that the pragmatic output of Poles learning EFL is dependent on one more independent variable: the estimation of future social distance (F).
{"title":"Measuring change in longitudinal research on pragmatic competence: A multinomial logistic model","authors":"Anna Szczepaniak-Kozak, E. Bakinowska, Katerina Strani","doi":"10.2478/bile-2020-0013","DOIUrl":"https://doi.org/10.2478/bile-2020-0013","url":null,"abstract":"Summary This paper focuses on pragmatic competence development in second or foreign language learners. In particular, it attempts to fill the significant research gap in measuring change in pragmatic competence and capturing pragmalinguistic and sociopragmatic development over time. The paper proposes an innovative approach by applying a logistic model with multinomial distribution for measuring change in InterLanguage Pragmatics Research (ILP). Applied in the field of pragmatics, this statistical tool offers a comprehensive and flexible approach to modelling relations between independent and dependent variables in ILP research. The model is tested in a longitudinal study of Polish undergraduate students learning English, and specifically in the way they formulate requests by means of requestive directness strategies. The paper concludes that, regardless of time elapsing, the factors P (power distance) and D (social distance) have a highly significant influence on the use of requestive directness strategies by Poles learning EFL. Furthermore, the analysis indicates that the pragmatic output of Poles learning EFL is dependent on one more independent variable: the estimation of future social distance (F).","PeriodicalId":8933,"journal":{"name":"Biometrical Letters","volume":"45 1","pages":"195 - 220"},"PeriodicalIF":0.0,"publicationDate":"2020-06-22","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"89358462","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
W. Mądry, Marcin Olik, B. Roszkowska-Mądra, M. Studnicki, D. Gozdowski, E. Wójcik-Gront
Summary High Nature Value farmlands in Europe are of greatest importance in the conservation of biodiversity. Their environmental importance has been recognized for some time, and has been studied mostly in Western Europe. This article describes the results of multivariate statistical analyses performed on data (13 variables) collected from the latest National Agricultural Census and the CORINE database to provide a typology of farmlands with respect to their nature value at municipality level (LAU 2, Local Administrative Units level 2) across Poland. All municipalities were grouped into eight categories (types). Some of the farmland categories were considered to be High Nature Value farmland (HNVf). The following interrelated variables mostly contributed to the identification of HNVf: share of protected areas and forest, grassland, arable land and fallow, farmland cover diversity, and rate of nitrogen fertilization. HNVf was identified in 958 out of 2173 municipalities, covering 44% of the territory of Poland. The identified HNVf also overlaps partially (61%) with LFAs (Less Favored Areas). Farmlands with the highest nature value are located mostly across mountain and hilly areas, close to forests, and protected areas on lowlands and river valleys. The identified HNV farmlands are characterized by low-input farming systems and a large share of semi-natural habitats with a high landscape mosaic.
{"title":"Identifying High Nature Value farmlands on a national scale based on multivariate typology at municipality (LAU 2) level","authors":"W. Mądry, Marcin Olik, B. Roszkowska-Mądra, M. Studnicki, D. Gozdowski, E. Wójcik-Gront","doi":"10.2478/bile-2020-0006","DOIUrl":"https://doi.org/10.2478/bile-2020-0006","url":null,"abstract":"Summary High Nature Value farmlands in Europe are of greatest importance in the conservation of biodiversity. Their environmental importance has been recognized for some time, and has been studied mostly in Western Europe. This article describes the results of multivariate statistical analyses performed on data (13 variables) collected from the latest National Agricultural Census and the CORINE database to provide a typology of farmlands with respect to their nature value at municipality level (LAU 2, Local Administrative Units level 2) across Poland. All municipalities were grouped into eight categories (types). Some of the farmland categories were considered to be High Nature Value farmland (HNVf). The following interrelated variables mostly contributed to the identification of HNVf: share of protected areas and forest, grassland, arable land and fallow, farmland cover diversity, and rate of nitrogen fertilization. HNVf was identified in 958 out of 2173 municipalities, covering 44% of the territory of Poland. The identified HNVf also overlaps partially (61%) with LFAs (Less Favored Areas). Farmlands with the highest nature value are located mostly across mountain and hilly areas, close to forests, and protected areas on lowlands and river valleys. The identified HNV farmlands are characterized by low-input farming systems and a large share of semi-natural habitats with a high landscape mosaic.","PeriodicalId":8933,"journal":{"name":"Biometrical Letters","volume":"32 1","pages":"63 - 84"},"PeriodicalIF":0.0,"publicationDate":"2020-06-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"80699283","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
京公网安备 11010802042870号
京ICP备2023020795号-1
扫码关注我们
求助内容：
应助结果提醒方式：