A comparison between joint coordinate system and attitude vector for multi-segment foot kinematics.

The joint angles of multi-segment foot models have been primarily described using two mathematical methods: the joint coordinate system and the attitude vector. This study aimed to determine whether the angles obtained through these two descriptors are comparable, and whether these descriptors have similar sensitivity to experimental errors. Six subjects walked eight times on an instrumented walkway while the joint angles among shank, hindfoot, medial forefoot, and lateral forefoot were measured. The angles obtained using both descriptors and their sensitivity to experimental errors were compared. There was no overall significant difference between the ranges of motion obtained using both descriptors. However, median differences of more than 6° were noticed for the medial-lateral forefoot joint. For all joints and rotation planes, both descriptors provided highly similar angle patterns (median correlation coefficient: R>0.90), except for the medial-lateral forefoot angle in the transverse plane (median R=0.77). The joint coordinate system was significantly more sensitive to anatomical landmarks misplacement errors. However, the absolute differences of sensitivity were small relative to the joints ranges of motion. In conclusion, the angles obtained using these two descriptors were not identical, but were similar for at least the shank-hindfoot and hindfoot-medial forefoot joints. Therefore, the angle comparison across descriptors is possible for these two joints. Comparison should be done more carefully for the medial-lateral forefoot joint. Moreover, despite different sensitivities to experimental errors, the effects of the experimental errors on the angles were small for both descriptors suggesting that both descriptors can be considered for multi-segment foot models.

Decision-support tools to manage drug incompatibilities: evaluation by nurses

Introduction Preventing drug incompatibilities has a high impact onthe safety of drug therapy. Although there are no internationalguidelines to manage drug incompatibilities, different decision-supporttools such as handbooks, cross-tables and databases are available.In a previous study, two decision-support tools have been pre-selectedby pharmacists as fitting nurses' needs on the wards1. The objective ofthis study was to have these both tools evaluated by nurses todetermine which would be the most suitable for their daily practice.Materials & Methods Evaluated tools were:1. Cross-table of drug pairs (http://files.chuv.ch/internet-docs/pha/medicaments/pha_phatab_compatibilitessip.pdf)2. Colour-table (a colour for each drug according to the pH: red =acid; blue = basic; yellow = neutral; black = to be infused alone)2Tools were assessed by 48 nurses in 5 units (PICU, adult andgeriatric intensive care, surgery, onco-hematology) using a standardizedform1. The scientific accuracy of the tools was evaluated bydetermining the compatibility of five drugs pairs (rate of correctanswers according to the Trissel's Handbook on Injectable Drugs,chi-square test). Their ergonomics, design, reliability and applicabilitywere estimated using visual analogue scales (VAS 0-10; 0 =null, 10 = excellent). Results are expressed as the median and interquartilerange (IQR) for 25% and 75% (Wilcoxon rank sum test).Results The rate of correct answers was above 90% for both tools(cross-table 96.2% vs colour-table 92.5%, p[0.05).The ergonomics and the applicability were higher for the crosstable[7.1 (IQR25 4.0, IQR75 8.0) vs 5.0 (IQR25 2.7, IQR75 7.0), p =0.025 resp. 8.3 (IQR25 7.4, IQR75 9.2) vs 7.6 (IQR25 5.9, IQR75 8.8)p = 0.047].The design of the colour-table was judged better [4.6 (IQR25 2.9,IQR75 7.1) vs 7.1 (IQR25 5.4, IQR75 8.4) p = 0.002].No difference was observed in terms of reliability [7.3 (IQR25 6.5,IQR75 8.4) vs 6.7 (IQR25 5.0, IQR758.6) p[0.05].The cross-table was globally preferred by 65% of the nurses (27%colour-table, 8% undetermined) and 68% would like to have thisdecision-support tool available for their daily practice.Discussion & Conclusion Both tools showed the same accuracy toassess drug compatibility. In terms of ergonomics and applicabilitythe cross-table was better than the colour-table, and was preferred bythe nurses for their daily practice. The cross-table will be implementedin our hospital as decision-support tool to help nurses tomanage drug incompatibilities.

Experimental support for the make-up hypothesis in nestling tawny owls (Strix aluco)

Body condition can affect coloration of traits used in sexual selection and parent-offspring communication by inducing rapid internal changes in pigment concentration or aggregation, thickness of collagen arrays, or blood flux. The recent "makeup hypothesis" proposes an alternative honesty-reinforcing mechanism, with behaviorally mediated deposition of substances on body surfaces ("cosmetics") generating covariation between body condition and coloration. In birds, the uropygial gland wax is actively spread on feathers using the bill and changes in its deposition rate may cause rapid changes in bill and plumage coloration. Using tawny owl nestlings, we tested 3 predictions of the makeup hypothesis, namely that 1) quantity of preen wax deposited accounts for variation in bill coloration, 2) an immune stimulation (induced by injection of a lipopolysaccharide [LPS]) impairs uropygial gland wax production, and 3) different intensities of immune stimulations (strong vs. weak stimulations induced by injections of either LPS or phytohemagglutinin [PHA], respectively) and high versus low food availabilities result in different bill colorations. We found that 1) preen wax reduced bill brightness, 2) a challenge with LPS impaired uropygial gland development, and 3) nestlings challenged with LPS had a brighter bill than PHA-injected nestlings, whereas diet manipulation had no significant effect. Altogether, these results suggest that a strong immune challenge may decrease preen wax deposition rate on the bill of nestling birds, at least by impairing gland wax production, which causes a change in bill coloration. Our study therefore highlights that cosmetic colors might signal short-term variation in immunological status.

Are the light-harvesting I complexes from Rhodospirillum rubrum arranged around the reaction centre in a square geometry?

The basic photosynthetic unit containing the reaction centre and the light-harvesting I complex (RC-LHI) of the purple non-sulphur bacterium Rhodospirillum rubrum was purified and reconstituted into two-dimensional (2D) membrane crystals. Transmission electron microscopy using conventional techniques and cryoelectron microscopy of the purified single particles and of 2D crystals yielded a projection of the RC-LHI complex at a resolution of at least 1.6 nm. In this projection the LHI ring appears to have a square symmetry and packs in a square crystal lattice. The square geometry of the LHI ring was observed also in images of single isolated particles of the RC-LHI complex. However, although the LHI units are packed identically within the crystal lattice, a new rotational analysis developed here showed that the reaction centres take up one of four possible orientations within the ring. This fourfold disorder supports our interpretation of a square ring symmetry and suggests that a hitherto undetected component may be present within the photosynthetic unit.

Optimal priming of poxvirus vector (NYVAC)-based HIV vaccine regimens requires 3 DNA injections. Results of the EV03/ANRS Vac20 Phase I/II Trial

Background: Two or three DNA primes have been used in previous smaller clinical trials, but the number required for optimal priming of viral vectors has never been assessed in adequately powered clinical trials. The EV03/ANRS Vac20 phase I/II trial investigated this issue using the DNA prime/poxvirus NYVAC boost combination, both expressing a common HIV-1 clade C immunogen consisting of Env and Gag-Pol-Nef polypeptide. Methods: 147 healthy volunteers were randomly allocated through 8 European centres to either 3xDNA plus 1xNYVAC (weeks 0, 4, 8 plus 24; n¼74) or to 2xDNA plus 2xNYVAC (weeks 0, 4 plus 20, 24; n¼73), stratified by geographical region and sex. T cell responses were quantified using the interferon g Elispot assay and 8 peptide pools; samples from weeks 0, 26 and 28 (time points for primary immunogenicity endpoint), 48 and 72 were considered for this analysis. Results: 140 of 147 participants were evaluable at weeks 26 and/ or 28. 64/70 (91%) in the 3xDNA arm compared to 56/70 (80%) in the 2xDNA arm developed a T cell response (P¼0.053). 26 (37%) participants of the 3xDNA arm developed a broader T cell response (Env plus at least to one of the Gag, Pol, Nef peptide pools) versus 15 (22%) in the 2xDNA arm (P¼0.047). At week 26, the overall magnitude of responses was also higher in the 3xDNA than in the 2xDNA arm (similar at week 28), with a median of 545 versus 328 SFUs/106 cells at week 26 (P<0.001). Preliminary overall evaluation showed that participants still developed T-cell response at weeks 48 (78%, n¼67) and 72 (70%, n¼66). Conclusion: This large clinical trial demonstrates that optimal priming of poxvirus-based vaccine regimens requires 3 DNA regimens and further confirms that the DNA/NYVAC prime boost vaccine combination is highly immunogenic and induced durable T-cell responses.

Value added logistical support service: logistics cost structure and performance in the new concept. Part 3

The aim of this study is to present an Activity-Based Costing spreadsheet tool for analyzing the logistics costs. The tool can be used both by customer-companies and logistics service providers. The study discusses the influence of different activity models on costs. Additionally this paper discusses about the logistical performance across the total supply chain This study is carried out using ananalytical research approach and literature material has been used for supplementing the concerned research approach. Cost structure analysis was based on the theory of activity-based management. This study was outlined to spare part logistics in machine-shop industry. The outlines of logistics services and logisticalperformance discussed in this report are based on the new logistics business concept (LMS-concept), which has been presented earlier in the Valssi-project. Oneof the aims of this study is to increase awareness of different activity modelson logistics costs. The report paints an overall picture about the business environment and requirements for the new logistics concept.

Torque vibration model of axial-flux surface-mounted permanent magnet synchronous machine

In order that the radius and thus ununiform structure of the teeth and otherelectrical and magnetic parts of the machine may be taken into consideration the calculation of an axial flux permanent magnet machine is, conventionally, doneby means of 3D FEM-methods. This calculation procedure, however, requires a lotof time and computer recourses. This study proves that also analytical methods can be applied to perform the calculation successfully. The procedure of the analytical calculation can be summarized into following steps: first the magnet is divided into slices, which makes the calculation for each section individually, and then the parts are submitted to calculation of the final results. It is obvious that using this method can save a lot of designing and calculating time. Thecalculation program is designed to model the magnetic and electrical circuits of surface mounted axial flux permanent magnet synchronous machines in such a way, that it takes into account possible magnetic saturation of the iron parts. Theresult of the calculation is the torque of the motor including the vibrations. The motor geometry and the materials and either the torque or pole angle are defined and the motor can be fed with an arbitrary shape and amplitude of three-phase currents. There are no limits for the size and number of the pole pairs nor for many other factors. The calculation steps and the number of different sections of the magnet are selectable, but the calculation time is strongly depending on this. The results are compared to the measurements of real prototypes. The permanent magnet creates part of the flux in the magnetic circuit. The form and amplitude of the flux density in the air-gap depends on the geometry and material of the magnetic circuit, on the length of the air-gap and remanence flux density of the magnet. Slotting is taken into account by using the Carter factor in the slot opening area. The calculation is simple and fast if the shape of the magnetis a square and has no skew in relation to the stator slots. With a more complicated magnet shape the calculation has to be done in several sections. It is clear that according to the increasing number of sections also the result will become more accurate. In a radial flux motor all sections of the magnets create force with a same radius. In the case of an axial flux motor, each radial section creates force with a different radius and the torque is the sum of these. The magnetic circuit of the motor, consisting of the stator iron, rotor iron, air-gap, magnet and the slot, is modelled with a reluctance net, which considers the saturation of the iron. This means, that several iterations, in which the permeability is updated, has to be done in order to get final results. The motor torque is calculated using the instantaneous linkage flux and stator currents. Flux linkage is called the part of the flux that is created by the permanent magnets and the stator currents passing through the coils in stator teeth. The angle between this flux and the phase currents define the torque created by the magnetic circuit. Due to the winding structure of the stator and in order to limit the leakage flux the slot openings of the stator are normally not made of ferromagnetic material even though, in some cases, semimagnetic slot wedges are used. In the slot opening faces the flux enters the iron almost normally (tangentially with respect to the rotor flux) creating tangential forces in the rotor. This phenomenon iscalled cogging. The flux in the slot opening area on the different sides of theopening and in the different slot openings is not equal and so these forces do not compensate each other. In the calculation it is assumed that the flux entering the left side of the opening is the component left from the geometrical centre of the slot. This torque component together with the torque component calculated using the Lorenz force make the total torque of the motor. It is easy to assume that when all the magnet edges, where the derivative component of the magnet flux density is at its highest, enter the slot openings at the same time, this will have as a result a considerable cogging torque. To reduce the cogging torquethe magnet edges can be shaped so that they are not parallel to the stator slots, which is the common way to solve the problem. In doing so, the edge may be spread along the whole slot pitch and thus also the high derivative component willbe spread to occur equally along the rotation. Besides forming the magnets theymay also be placed somewhat asymmetric on the rotor surface. The asymmetric distribution can be made in many different ways. All the magnets may have a different deflection of the symmetrical centre point or they can be for example shiftedin pairs. There are some factors that limit the deflection. The first is that the magnets cannot overlap. The magnet shape and the relative width compared to the pole define the deflection in this case. The other factor is that a shifting of the poles limits the maximum torque of the motor. If the edges of adjacent magnets are very close to each other the leakage flux from one pole to the other increases reducing thus the air-gap magnetization. The asymmetric model needs some assumptions and simplifications in order to limit the size of the model and calculation time. The reluctance net is made for symmetric distribution. If the magnets are distributed asymmetrically the flux in the different pole pairs will not be exactly the same. Therefore, the assumption that the flux flows from the edges of the model to the next pole pairs, in the calculation model from one edgeto the other, is not correct. If it were wished for that this fact should be considered in multi-pole pair machines, this would mean that all the poles, in other words the whole machine, should be modelled in reluctance net. The error resulting from this wrong assumption is, nevertheless, irrelevant.

Direct torque control of permanent magnet synchronous machines - analysis and implementation

Thedirect torque control (DTC) has become an accepted vector control method besidethe current vector control. The DTC was first applied to asynchronous machines,and has later been applied also to synchronous machines. This thesis analyses the application of the DTC to permanent magnet synchronous machines (PMSM). In order to take the full advantage of the DTC, the PMSM has to be properly dimensioned. Therefore the effect of the motor parameters is analysed taking the control principle into account. Based on the analysis, a parameter selection procedure is presented. The analysis and the selection procedure utilize nonlinear optimization methods. The key element of a direct torque controlled drive is the estimation of the stator flux linkage. Different estimation methods - a combination of current and voltage models and improved integration methods - are analysed. The effect of an incorrect measured rotor angle in the current model is analysed andan error detection and compensation method is presented. The dynamic performance of an earlier presented sensorless flux estimation method is made better by improving the dynamic performance of the low-pass filter used and by adapting the correction of the flux linkage to torque changes. A method for the estimation ofthe initial angle of the rotor is presented. The method is based on measuring the inductance of the machine in several directions and fitting the measurements into a model. The model is nonlinear with respect to the rotor angle and therefore a nonlinear least squares optimization method is needed in the procedure. A commonly used current vector control scheme is the minimum current control. In the DTC the stator flux linkage reference is usually kept constant. Achieving the minimum current requires the control of the reference. An on-line method to perform the minimization of the current by controlling the stator flux linkage reference is presented. Also, the control of the reference above the base speed is considered. A new estimation flux linkage is introduced for the estimation of the parameters of the machine model. In order to utilize the flux linkage estimates in off-line parameter estimation, the integration methods are improved. An adaptive correction is used in the same way as in the estimation of the controller stator flux linkage. The presented parameter estimation methods are then used in aself-commissioning scheme. The proposed methods are tested with a laboratory drive, which consists of a commercial inverter hardware with a modified software and several prototype PMSMs.

Lack of support for the association between facial shape and aggression: a reappraisal based on a worldwide population genetics perspective

Antisocial and criminal behaviors are multifactorial traits whose interpretation relies on multiple disciplines. Since these interpretations may have social, moral and legal implications, a constant review of the evidence is necessary before any scientific claim is considered as truth. A recent study proposed that men with wider faces relative to facial height (fWHR) are more likely to develop unethical behaviour mediated by a psychological sense of power. This research was based on reports suggesting that sexual dimorphism and selection would be responsible for a correlation between fWHR and aggression. Here we show that 4,960 individuals from 94 modern human populations belonging to a vast array of genetic and cultural contexts do not display significant amounts of fWHR sexual dimorphism. Further analyses using populations with associated ethnographical records as well as samples of male prisoners of the Mexico City Federal Penitentiary condemned by crimes of variable level of inter-personal aggression (homicide, robbery, and minor faults) did not show significant evidence, suggesting that populations/individuals with higher levels of bellicosity, aggressive behaviour, or power-mediated behaviour display greater fWHR. Finally, a regression analysis of fWHR on individual"s fitness showed no significant correlation between this facial trait and reproductive success. Overall, our results suggest that facial attributes are poor predictors of aggressive behaviour, or at least, that sexual selection was weak enough to leave a signal on patterns of between- and within-sex and population facial variation.

Preferred Spatial Frequencies for Human Face Processing Are Associated with Optimal Class Discrimination in the Machine

Psychophysical studies suggest that humans preferentially use a narrow band of low spatial frequencies for face recognition. Here we asked whether artificial face recognition systems have an improved recognition performance at the same spatial frequencies as humans. To this end, we estimated recognition performance over a large database of face images by computing three discriminability measures: Fisher Linear Discriminant Analysis, Non-Parametric Discriminant Analysis, and Mutual Information. In order to address frequency dependence, discriminabilities were measured as a function of (filtered) image size. All three measures revealed a maximum at the same image sizes, where the spatial frequency content corresponds to the psychophysical found frequencies. Our results therefore support the notion that the critical band of spatial frequencies for face recognition in humans and machines follows from inherent properties of face images, and that the use of these frequencies is associated with optimal face recognition performance.

Development of support and drive systems for defiberization drum

Diplomityön tarkoituksena oli kuidutusrumpulaitteiston käytön- ja kannatuksen kehittä-minen. Työ rajattiin laajuutensa vuoksi koskemaan tuotesarjan viittä pienintä kokoa. Työn alkuosassa käsitellään kuidutuksen teoriaa ja siihen soveltuvia laitteistoja. Käytön suunnittelun kannalta olennaista käynnistystehon tarvetta on tarkasteltu lähtökohdaisesti fysiikan avulla. Perustietoja teorialle on haettu aiemmista tutkimuksista sekä kirjallisuu-desta. Tarkastelun tuloksena teoriaa on kehitty ja se on saatu vastaamaan todellisuutta aiempaa paremmin. Kannatuksen ja käytön toteuttamisvaihtoja etsittäessä on käytetty systemaattisen koneen-suunnittelun keinoja. Saatuja ideoita on arvioitu teknillis-taloudellisin perustein ja näistä on valittu parhaat vaihtoehdot jatkokehitykseen. Jatkokehitysvaiheessa ratkaisuvaihto-ehtoja on tarkasteltu komponenttitasolla ja näistä on tehty yksityiskohtaiset kustannus-laskelmat. Työn tuloksena on esitetty kannatuksen ja käytön toteutusvaihtoehto, jonka avulla voidaan saavuttaa merkittäviä kustannussäästöjä. Korkea, 30 prosentin kustannussäästö-tavoite saavutettiin.

Designing a new type of support for rotary kilns

Lime kiln is used as a part of the modern kraft pulp process in order to produce burnt lime from lime mud. This rotating kiln is supported by support rollers, which are traditionally supported by journal bearings. Since the continuous growth in the production of pulp mills requires larger lime kilns, the traditional bearing construction has become unreliable. The main problem especially involves the running-in phase of the bearings. In the present thesis, a new type of support roller was developed by using the systematic approach of machine design. Structural analysis was conducted on the critical parts of the selected structure by the finite element method. The operation of hydrodynamic bearings was examined by analytical methods. As a result of this work, a new type of support for rotating kilns was designed, which is more reliable and easier to service. A new support roller geometry is described, which pro¬vides for significant cost savings.

Salient Pole Synchronous Machine Modelling in an Industrial Direct Torque Controlled Drive Application

Synchronous motors are used mainly in large drives, for example in ship propulsion systems and in steel factories' rolling mills because of their high efficiency, high overload capacity and good performance in the field weakening range. This, however, requires an extremely good torque control system. A fast torque response and a torque accuracy are basic requirements for such a drive. For large power, high dynamic performance drives the commonly known principle of field oriented vector control has been used solely hitherto, but nowadays it is not the only way to implement such a drive. A new control method Direct Torque Control (DTC) has also emerged. The performance of such a high quality torque control as DTC in dynamically demanding industrial applications is mainly based on the accurate estimate of the various flux linkages' space vectors. Nowadays industrial motor control systems are real time applications with restricted calculation capacity. At the same time the control system requires a simple, fast calculable and reasonably accurate motor model. In this work a method to handle these problems in a Direct Torque Controlled (DTC) salient pole synchronous motor drive is proposed. A motor model which combines the induction law based "voltage model" and motor inductance parameters based "current model" is presented. The voltage model operates as a main model and is calculated at a very fast sampling rate (for example 40 kHz). The stator flux linkage calculated via integration from the stator voltages is corrected using the stator flux linkage computed from the current model. The current model acts as a supervisor that prevents only the motor stator flux linkage from drifting erroneous during longer time intervals. At very low speeds the role of the current model is emphasised but, nevertheless, the voltage model always stays the main model. At higher speeds the function of the current model correction is to act as a stabiliser of the control system. The current model contains a set of inductance parameters which must be known. The validation of the current model in steady state is not self evident. It depends on the accuracy of the saturated value of the inductances. Parameter measurement of the motor model where the supply inverter is used as a measurement signal generator is presented. This so called identification run can be performed prior to delivery or during drive commissioning. A derivation method for the inductance models used for the representation of the saturation effects is proposed. The performance of the electrically excited synchronous motor supplied with the DTC inverter is proven with experimental results. It is shown that it is possible to obtain a good static accuracy of the DTC's torque controller for an electrically excited synchronous motor. The dynamic response is fast and a new operation point is achieved without oscillation. The operation is stable throughout the speed range. The modelling of the magnetising inductance saturation is essential and cross saturation has to be considered as well. The effect of cross saturation is very significant. A DTC inverter can be used as a measuring equipment and the parameters needed for the motor model can be defined by the inverter itself. The main advantage is that the parameters defined are measured in similar magnetic operation conditions and no disagreement between the parameters will exist. The inductance models generated are adequate to meet the requirements of dynamically demanding drives.

Machine Learning and Clinical Text. Supporting Health Information Flow

Fluent health information flow is critical for clinical decision-making. However, a considerable part of this information is free-form text and inabilities to utilize it create risks to patient safety and cost-­effective hospital administration. Methods for automated processing of clinical text are emerging. The aim in this doctoral dissertation is to study machine learning and clinical text in order to support health information flow.First, by analyzing the content of authentic patient records, the aim is to specify clinical needs in order to guide the development of machine learning applications.The contributions are a model of the ideal information flow,a model of the problems and challenges in reality, and a road map for the technology development. Second, by developing applications for practical cases,the aim is to concretize ways to support health information flow. Altogether five machine learning applications for three practical cases are described: The first two applications are binary classification and regression related to the practical case of topic labeling and relevance ranking.The third and fourth application are supervised and unsupervised multi-class classification for the practical case of topic segmentation and labeling.These four applications are tested with Finnish intensive care patient records.The fifth application is multi-label classification for the practical task of diagnosis coding. It is tested with English radiology reports.The performance of all these applications is promising. Third, the aim is to study how the quality of machine learning applications can be reliably evaluated.The associations between performance evaluation measures and methods are addressed,and a new hold-out method is introduced.This method contributes not only to processing time but also to the evaluation diversity and quality. The main conclusion is that developing machine learning applications for text requires interdisciplinary, international collaboration. Practical cases are very different, and hence the development must begin from genuine user needs and domain expertise. The technological expertise must cover linguistics,machine learning, and information systems. Finally, the methods must be evaluated both statistically and through authentic user-feedback.

PREDICTING THE BOILING POINT OF PCDD/Fs BY THE QSPR METHOD BASED ON THE MOLECULAR DISTANCE-EDGE VECTOR INDEX

The quantitative structure property relationship (QSPR) for the boiling point (Tb) of polychlorinated dibenzo-p-dioxins and polychlorinated dibenzofurans (PCDD/Fs) was investigated. The molecular distance-edge vector (MDEV) index was used as the structural descriptor. The quantitative relationship between the MDEV index and Tb was modeled by using multivariate linear regression (MLR) and artificial neural network (ANN), respectively. Leave-one-out cross validation and external validation were carried out to assess the prediction performance of the models developed. For the MLR method, the prediction root mean square relative error (RMSRE) of leave-one-out cross validation and external validation was 1.77 and 1.23, respectively. For the ANN method, the prediction RMSRE of leave-one-out cross validation and external validation was 1.65 and 1.16, respectively. A quantitative relationship between the MDEV index and Tb of PCDD/Fs was demonstrated. Both MLR and ANN are practicable for modeling this relationship. The MLR model and ANN model developed can be used to predict the Tb of PCDD/Fs. Thus, the Tb of each PCDD/F was predicted by the developed models.