Fitting model-based psychometric functions to simultaneity and temporal-order judgment data: MATLAB and R routines

Research on temporal-order perception uses temporal-order judgment (TOJ) tasks or synchrony judgment (SJ) tasks in their binary SJ2 or ternary SJ3 variants. In all cases, two stimuli are presented with some temporal delay, and observers judge the order of presentation. Arbitrary psychometric functions are typically fitted to obtain performance measures such as sensitivity or the point of subjective simultaneity, but the parameters of these functions are uninterpretable. We describe routines in MATLAB and R that fit model-based functions whose parameters are interpretable in terms of the processes underlying temporal-order and simultaneity judgments and responses. These functions arise from an independent-channels model assuming arrival latencies with exponential distributions and a trichotomous decision space. Different routines fit data separately for SJ2, SJ3, and TOJ tasks, jointly for any two tasks, or also jointly for the three tasks (for common cases in which two or even the three tasks were used with the same stimuli and participants). Additional routines provide bootstrap p-values and confidence intervals for estimated parameters. A further routine is included that obtains performance measures from the fitted functions. An R package for Windows and source code of the MATLAB and R routines are available as Supplementary Files.

A genetic two-factor model of the covariation among a subset of Multidimensional Aptitude Battery and Wechsler Adult Intelligence Scale - Revised subtests

The phenotypic and genetic factor structure of performance on five Multidimensional Aptitude Battery (MAB) subtests and one Wechsler Adult Intelligence Scale-Revised (WAIS-R) subtest was explored in 390 adolescent twin pairs (184 monozygotic [MZ]; 206 dizygotic (DZ)). The temporal stability of these measures was derived from a subsample of 49 twin pairs, with test-retest correlations ranging from .67 to .85. A phenotypic factor model, in which performance and verbal factors were correlated, provided a good fit to the data. Genetic modeling was based on the phenotypic factor structure, but also took into account the additive genetic (A), common environmental (C), and unique environmental (E) parameters derived from a fully saturated ACE model. The best fitting model was characterized by a genetic correlated two-factor structure with specific effects, a general common environmental factor, and overlapping unique environmental effects. Results are compared to multivariate genetic models reported in children and adults, with the most notable difference being the growing importance of common genes influencing diverse abilities in adolescence. (C) 2003 Elsevier Inc. All rights reserved.

Direction of Causation Modeling Between Cross-Sectional Measures of Parenting and Psychological Distress in Female Twins

Under certain conditions, cross-sectional analysis of cross-twin intertrait correlations can provide important information about the direction of causation (DOC) between two variables. A community-based sample of Australian female twins aged 18 to 45 years was mailed an extensive Health and Lifestyle Questionnaire (HLQ) that covered a wide range of personality and behavioral measures. Included were self-report measures of recent psychological distress and perceived childhood environment (PBI). Factor analysis of the PBI yielded three interpretable dimensions: Coldness, Overprotection, and Autonomy. Univariate analysis revealed that parental Overprotection and Autonomy were best explained by additive genetic, shared, and nonshared environmental effects (ACE), whereas the best-fitting model for PBI Coldness and the three measures of psychological distress (Depression, Phobic Anxiety, and Somatic Distress) included only additive genetic and nonshared environmental effects (AE). A common pathway model best explained the covariation between (1) the three PBI dimensions and (2) the three measures of psychological distress. DOC modeling between latent constructs of parenting and psychological distress revealed that a model which specified recollected parental behavior as the cause of psychological distress provided a better fit than a model which specified psychological distress as the cause of recollected parental behavior. Power analyses and limitations of the findings are discussed.

Effect of surfactant on PDLC films with and without permanent memory effect

Dissertação para obtenção do Grau de Mestre em Engenharia Química e Bioquímica

The influence of cationic lipid type on in-vitro release kinetic profiles of antisense oligonucleotide from cationic nanoemulsions.

Novel formulations of cationic nanoemulsions based on three different lipids were developed to strengthen the attraction of the polyanionic oligonucleotide (ODN) macromolecules to the cationic moieties on the oil nanodroplets. These formulations were developed to prolong the release of the ODN from the nanoemulsion under appropriate physiological dilutions as encountered in the eye following topical application. Increasing the concentration of the new cationic lipid exhibiting two cationic amine groups (AOA) in the emulsion from 0.05% to 0.4% did not alter markedly the particle size or zeta potential value of the blank cationic nanoemulsion. The extent of ODN association did not vary significantly when the initial concentration of ODN remained constant at 10 microM irrespective of the cationic lipid nature. However, the zeta potential value dropped consistently with the low concentrations of 0.05% and 0.1% of AOA in the emulsions suggesting that an electrostatic attraction occurred between the cationic lipids and the polyanionic ODN molecules at the o/w interface. Only the nanoemulsion prepared with N-[1-(2,3-dioleoyloxy)propyl]-N,N,N-trimethylammonium salts (DOTAP) remained physically stable over time. DOTAP cationic lipid nanoemulsion was the most efficient formulation capable of retaining the ODN despite the high dilution of 1:100 with simulated tear solution (STS). Less than 10% of the ODN was exchanged in contrast to 40-50% with the other cationic nanoemulsions. The in-vitro release kinetic behavior of ODN exchange with physiological anions present in the STS appears to be complex and difficult to characterize using mathematical fitting model equations. Further pharmacokinetic studies are needed to verify our kinetic assumptions and confirm the in-vitro ODN release profile from DOTAP cationic nanoemulsions.

The average inter-crossing number of equilateral random walks and polygons

In this paper, we study the average inter-crossing number between two random walks and two random polygons in the three-dimensional space. The random walks and polygons in this paper are the so-called equilateral random walks and polygons in which each segment of the walk or polygon is of unit length. We show that the mean average inter-crossing number ICN between two equilateral random walks of the same length n is approximately linear in terms of n and we were able to determine the prefactor of the linear term, which is a = (3 In 2)/(8) approximate to 0.2599. In the case of two random polygons of length n, the mean average inter-crossing number ICN is also linear, but the prefactor of the linear term is different from that of the random walks. These approximations apply when the starting points of the random walks and polygons are of a distance p apart and p is small compared to n. We propose a fitting model that would capture the theoretical asymptotic behaviour of the mean average ICN for large values of p. Our simulation result shows that the model in fact works very well for the entire range of p. We also study the mean ICN between two equilateral random walks and polygons of different lengths. An interesting result is that even if one random walk (polygon) has a fixed length, the mean average ICN between the two random walks (polygons) would still approach infinity if the length of the other random walk (polygon) approached infinity. The data provided by our simulations match our theoretical predictions very well.

Exploring the proposed DSM-5 diagnostic threshold and severity scale for the alcohol use disorder

Aims :¦Several studies have questioned the validity of separating the diagnosis of alcohol abuse from that of alcohol dependence, and the DSM-5 task force has proposed combining the criteria from these two diagnoses to assess a single category of alcohol use disorders (AUD). Furthermore, the DSM-5 task force has proposed including a new 2-symptom threshold and a severity scale based on symptom counts for the AUD diagnosis. The current study aimed to examine these modifications in a large population-based sample.¦Method :¦Data stemmed from an adult sample (N=2588 ; mean age 51.3 years (s.d.: 0.2), 44.9% female) of current and lifetime drinkers from the PsyCoLaus study, conducted in the Lausanne area in Switzerland. AUDs and validating variables were assessed using a semi-structured diagnostic interview for the assessment of alcohol¦and other major psychiatric disorders. First, the adequacy of the proposed 2- symptom threshold was tested by comparing threshold models at each possible cutoff and a linear model, in relation to different validating variables. The model with the smallest Akaike Criterion Information (AIC) value was established as the best¦model for each validating variable. Second, models with varying subsets of individual AUD symptoms were created to assess the associations between each symptom and the validating variables. The subset of symptoms with the smallest AIC value was established as the best subset for each validator.¦Results :¦1) For the majority of validating variables, the linear model was found to be the best fitting model. 2) Among the various subsets of symptoms, the symptoms most frequently associated with the validating variables were : a) drinking despite having knowledge of a physical or psychological problem, b) having had a persistent desire or unsuccessful efforts to cut down or control drinking and c) craving. The¦least frequent symptoms were : d) drinking in larger amounts or over a longer period than was intended, e) spending a great deal of time in obtaining, using or recovering from alcohol use and f) failing to fulfill major role obligations.¦Conclusions :¦The proposed DSM-5 2-symptom threshold did not receive support in our data. Instead, a linear AUD diagnosis was supported with individuals receiving an increasingly severe AUD diagnosis. Moreover, certain symptoms were more frequently associated with the validating variables, which suggests that these¦symptoms should be considered as more severe.

La compétence en emploi peut-elle prédire la santé psychologique des enseignants?

La présente thèse de doctorat porte sur la relation entre la santé psychologique au travail, les besoins fondamentaux et la compétence en emploi. En plus de valider un questionnaire de compétences en emploi (QCE), cette thèse propose que les compétences en emploi puissent prédire la santé psychologique au travail par la satisfaction des besoins fondamentaux dans le domaine de l’enseignement. Le premier article a pour objectif de présenter la conception et la validation de l’instrument de mesure des compétences en emploi dans le domaine de l’enseignement. Une première étape dans notre étude a permis l’identification de trois solutions factorielles possibles : structure à un seul facteur, structure à deux facteurs et structure à un facteur de deuxième ordre. Les matrices des corrélations des énoncés étaient fiables, se prêtaient bien à l’analyse factorielle exploratoire et présentaient de bonnes cohérences internes. La deuxième étape de notre étude a été de type confirmatif. Chacune des trois solutions factorielles proposées a été analysée, ce qui a permis l’identification du modèle le mieux ajusté compte tenu des seuils empiriques des indices retenus pour l’analyse par équations structurelles. L’étude réalisée sur un échantillon d’enseignants québécois démontre que nous pouvons conclure aux bonnes qualités psychométriques de l’instrument analysé. Les limites et les apports de cette étude seront aussi présentés. Le deuxième article examine les liens entre les compétences en emploi dans le domaine de l’enseignement, la santé psychologique au travail et la satisfaction des trois besoins fondamentaux (autonomie, compétence et affiliation). Un modèle de médiation est testé par des analyses de modélisation par équations structurelles. Dans ce modèle, on considère que la satisfaction des trois besoins fondamentaux (autonomie, compétence et affiliation sociale) agisse comme médiateur dans la relation entre les compétences en emploi et la santé psychologique au travail. Une relation de médiation partielle (Baron & Kenny, 1986) a été trouvée dans notre échantillon, en utilisant l’analyse de rééchantillonnage par « bootstrap », dans le cadre de la modélisation par équations structurelles. Les limites de la recherche, ainsi que des suggestions de recherches futures seront proposées.

The relationship between multi-species catch and effort: among fishery comparisons

The control of fishing mortality via fishing effort remains fundamental to most fisheries management strategies even at the local community or co-management level. Decisions to support such strategies require knowledge of the underlying response of the catch to changes in effort. Even under adaptive management strategies, imprecise knowledge of the response is likely to help accelerate the adaptive learning process. Data and institutional capacity requirements to employ multi-species biomass dynamics and age-structured models invariably render their use impractical particularly in less developed regions of the world. Surplus production models fitted to catch and effort data aggregated across all species offer viable alternatives. The current paper seeks models of this type that best describe the multi-species catch–effort responses in floodplain-rivers, lakes and reservoirs and reef-based fisheries based upon among fishery comparisons, building on earlier work. Three alternative surplus production models were fitted to estimates of catch per unit area (CPUA) and fisher density for 258 fisheries in Africa, Asia and South America. In all cases examined, the best or equal best fitting model was the Fox type, explaining up to 90% of the variation in CPUA. For lake and reservoir fisheries in Africa and Asia, the Schaefer and an asymptotic model fitted equally well. The Fox model estimates of fisher density (fishers km−2) at maximum yield (iMY) for floodplain-rivers, African lakes and reservoirs and reef-based fisheries are 13.7 (95% CI [11.8, 16.4]); 27.8 (95% CI [17.5, 66.7]) and 643 (95% CI [459,1075]), respectively and compare well with earlier estimates. Corresponding estimates of maximum yield are also given. The significantly higher value of iMY for reef-based fisheries compared to estimates for rivers and lakes reflects the use of a different measure of fisher density based upon human population size estimates. The models predict that maximum yield is achieved at a higher fishing intensity in Asian lakes compared to those in Africa. This may reflect the common practice in Asia of stocking lakes to augment natural recruitment. Because of the equilibrium assumptions underlying the models, all the estimates of maximum yield and corresponding levels of effort should be treated with caution.

Estimation of genetic variance for macro- and micro-environmental sensitivity using double hierarchical generalized linear models

Background: Genetic variation for environmental sensitivity indicates that animals are genetically different in their response to environmental factors. Environmental factors are either identifiable (e.g. temperature) and called macro-environmental or unknown and called micro-environmental. The objectives of this study were to develop a statistical method to estimate genetic parameters for macro- and micro-environmental sensitivities simultaneously, to investigate bias and precision of resulting estimates of genetic parameters and to develop and evaluate use of Akaike’s information criterion using h-likelihood to select the best fitting model. Methods: We assumed that genetic variation in macro- and micro-environmental sensitivities is expressed as genetic variance in the slope of a linear reaction norm and environmental variance, respectively. A reaction norm model to estimate genetic variance for macro-environmental sensitivity was combined with a structural model for residual variance to estimate genetic variance for micro-environmental sensitivity using a double hierarchical generalized linear model in ASReml. Akaike’s information criterion was constructed as model selection criterion using approximated h-likelihood. Populations of sires with large half-sib offspring groups were simulated to investigate bias and precision of estimated genetic parameters. Results: Designs with 100 sires, each with at least 100 offspring, are required to have standard deviations of estimated variances lower than 50% of the true value. When the number of offspring increased, standard deviations of estimates across replicates decreased substantially, especially for genetic variances of macro- and micro-environmental sensitivities. Standard deviations of estimated genetic correlations across replicates were quite large (between 0.1 and 0.4), especially when sires had few offspring. Practically, no bias was observed for estimates of any of the parameters. Using Akaike’s information criterion the true genetic model was selected as the best statistical model in at least 90% of 100 replicates when the number of offspring per sire was 100. Application of the model to lactation milk yield in dairy cattle showed that genetic variance for micro- and macro-environmental sensitivities existed. Conclusion: The algorithm and model selection criterion presented here can contribute to better understand genetic control of macro- and micro-environmental sensitivities. Designs or datasets should have at least 100 sires each with 100 offspring.

Fertilidade do caranguejo de mangue Ucides cordatus (Linnaeus) (Crustacea, Brachyura, Ocypodidae), em Iguape (São Paulo, Brasil)

A total of 58 Ucides cordatus (Linnaeus, 1763) ovigerous females in final embryonic stage were collected in mangrove forest at Iguape (São Paulo). Each female were kept in laboratory until the larvae hatching and had their carapace width (CW) measured with caliper (0.05 mm) and the larval number (LN) calculated by volumetric method. The hatching rates (HR) were calculated for each animal's size class (5 mm). The values of LN and CW obtained were subjected to regression analyses and the best fitting model was chosen by comparing the obtained determination coefficients (R²). The fertility varied from 71,200 to 220,800 larvae (147,169 ± 32,070 larvae), corresponding to 41.7 to 76.8 mm CW (63.7±7.9 mm), respectively. The relationship LNxCW was expressed by the equation LN = 284.1CW1.5º¹ (N = 58; R² = 0.74; p < 0.01), that showed a positive correlation between the variables. Most part of ovigerous females (89.7%) showed their hatching during the nightly period. A negative correlation between HR and CW was observed in these animals (p < 0.05), probably due to reduction of spermatophores in the seminal receptacle (multiparous females) or due to senility in the females with large size. The biological studies of U. cordatus are very important for practices adoption that preserve and make the rational use of this recourse in Brazilian's mangrove forests.

Seleção para peso pós-desmama em um rebanho Gir. 2. Estimativas de variâncias e parâmetros genéticos dos efeitos direto e materno para características de crescimento

Componentes de variância e parâmetros genéticos para características de crescimento foram estimados usando diferentes modelos em um rebanho da raça Gir. Utilizou-se o método da máxima verossimilhança restrita sob modelo animal univariado. Os modelos de análise incluíram os efeitos fixos de mês de nascimento, grupo contemporâneo e idade da vaca. Cinco modelos diferindo quanto aos efeitos aleatórios foram testados. Para todas as características da fase pré-desmama, o teste de razão de verossimilhança (LRT) indicou o modelo com efeito genético aditivo direto e efeitos maternos (genético e de ambiente permanente) como o de melhor ajuste. As estimativas de herdabilidade direta para peso ao nascer (PN), peso aos quatro meses corrigido para 120 dias (P120), peso à desmama corrigido para 210 dias (P210) e ganho diário na fase pré-desmama (GPRE) foram, respectivamente, 0,31± 0,07; 0,14 ± 0,06; 0,23 ± 0,07 e 0,22 ± 0,07. Para as características da fase pós-desmama, o modelo que forneceu o melhor ajuste aos dados incluiu apenas o efeito genético aditivo direto. As estimativas de herdabilidade direta para peso de machos ao final da prova de ganho de peso (P378), peso de fêmeas corrigido para 550 dias (P550), ganho diário na prova de ganho de peso (G112), altura aos 378 dias em machos (AM) e altura aos 550 dias, em fêmeas (AF) foram, respectivamente: 0,45 ± 0,11; 0,29 ± 0,11; 0,37 ± 0,11; 0,79 ± 0,13 e 0,36 ± 0,0. Os efeitos maternos, tanto o genético quanto o de ambiente permanente, foram fontes de variação importantes para as características da fase pré-desmama, não sendo verificada influência desses efeitos sobre as características da fase pós-desmama.

Componentes de variância para o período de gestação em bovinos de corte

Estimaram-se os componentes de variância do período de gestação (PG) considerando-se o efeito direto do bezerro e os efeitos direto da vaca e aleatório do touro (pai do bezerro). Além dos efeitos aleatórios, os modelos estatísticos incluíram os efeitos fixos de grupo de contemporâneos, composto pela data juliana de inseminação, ano de nascimento e sexo do bezerro, e a idade da vaca ao parto (covariável linear e quadrática). As médias de PG para os animais da raça Nelore foram 294,55 dias (machos) e 293,34 dias (fêmeas) e para os animais cruzados, 292,49 dias (machos) e 292,55 dias (fêmeas). Os componentes de variância observados no Nelore, utilizando-se o modelo 1, que considerou o PG como característica do bezerro, foram 14,47; 72,78 e 57,31 para os componentes aditivo direto (sigma2a), fenotípico total (sigma2p) e residual (sigma2e), respectivamente, e a herdabilidade foi 0,21. Para os animais cruzados, pelo mesmo modelo, os componentes de variância foram 90,40 (sigma2a), 127,35 (sigma2p ) e 36,95 (sigma2e), e a herdabilidade, 0,71. Os componentes de variância do PG do Nelore sob o modelo 2, que considerou o PG como característica da vaca, foram 12,78; 5,01; 74,84 e 57,05, para sigma2a; sigma2pb (variância do pai do bezerro); sigma2p e sigma2e, respectivamente. A fração fenotípica atribuída ao touro (c²) foi 0,07 e a repetibilidade, 0,17. Para os cruzados, estimaram-se 22,11; 22,97; 127,70 e 82,61 para sigma2a; sigma2pb; sigma2p e sigma2e, respectivamente, enquanto o c² foi 0,18 e o coeficiente de repetibilidade, 0,17. Sugere-se, para fins de seleção, que o coeficiente de herdabilidade utilizado seja o obtido na análise em que PG foi considerado como característica do bezerro.

Estimação de parâmetros genéticos para produção de leite de cabras da raça Alpina

Foram utilizados 9.374 registros semanais de produção de leite de 302 primeiras lactações de cabras da raça Alpina. A produção de leite no dia do controle foi analisada por meio de um modelo animal, unicarater, de regressão aleatória, em que as funções de covariâncias para os componentes genéticos aditivos e de ambiente permanente foram modeladas por meio das funções de Wilmink, Ali e Schaeffer e por polinômios ortogonais, em uma escala de Legendre de ordens cúbica e quíntica. Assumiu-se, ainda, variância residual homogênea durante toda a lactação e heterogênea com três e quatro classes de variância residual. Os modelos foram comparados pelo critério de informação de Akaike (AIC), pelo critério de informação Bayesiano de Schwar (BIC), pela função de verossimilhança (Ln L), pela visualização das estimativas de variâncias genéticas, de ambiente permanente, fenotípicas e residuais e pelas herdabilidades. O polinômio de Legendre de ordem quíntica, com quatro e três classes de variâncias residuais, e a função de Ali e Schaeffer, com quatro classes de variâncias residuais, foram indicados como os mais adequados pelo AIC, BIC e Ln L. Estes modelos diferiram na partição da variância fenotípica para as variâncias de ambiente permanente, genética e residual apenas no início e no final da lactação. Contudo, a função de Ali e Schaeffer resultou em estimativas negativas de correlação genética entre os controles mais distantes. O polinômio de Legendre de ordem quíntica, assumindo variância residual heterogênea, mostrou-se mais adequado para ajustar a produção de leite no dia do controle de cabras da raça Alpina.

Avaliação de modelos estatísticos considerando a incerteza da paternidade

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES)