Characterization of the variable cow's age at last calving as a measurement of longevity by using the Kaplan-Meier estimator and the Cox model

In most studies on beef cattle longevity, only the cows reaching a given number of calvings by a specific age are considered in the analyses. With the aim of evaluating all cows with productive life in herds, taking into consideration the different forms of management on each farm, it was proposed to measure cow longevity from age at last calving (ALC), that is, the most recent calving registered in the files. The objective was to characterize this trait in order to study the longevity of Nellore cattle, using the Kaplan-Meier estimators and the Cox model. The covariables and class effects considered in the models were age at first calving (AFC), year and season of birth of the cow and farm. The variable studied (ALC) was classified as presenting complete information (uncensored = 1) or incomplete information (censored = 0), using the criterion of the difference between the date of each cow's last calving and the date of the latest calving at each farm. If this difference was >36 months, the cow was considered to have failed. If not, this cow was censored, thus indicating that future calving remained possible for this cow. The records of 11 791 animals from 22 farms within the Nellore Breed Genetic Improvement Program ('Nellore Brazil') were used. In the estimation process using the Kaplan-Meier model, the variable of AFC was classified into three age groups. In individual analyses, the log-rank test and the Wilcoxon test in the Kaplan-Meier model showed that all covariables and class effects had significant effects (P < 0.05) on ALC. In the analysis considering all covariables and class effects, using the Wald test in the Cox model, only the season of birth of the cow was not significant for ALC (P > 0.05). This analysis indicated that each month added to AFC diminished the risk of the cow's failure in the herd by 2%. Nonetheless, this does not imply that animals with younger AFC had less profitability. Cows with greater numbers of calvings were more precocious than those with fewer calvings. Copyright © The Animal Consortium 2012.

On the Accelerated Failure Time Model for Current Status and Interval Censored Data

This paper introduces a novel approach to making inference about the regression parameters in the accelerated failure time (AFT) model for current status and interval censored data. The estimator is constructed by inverting a Wald type test for testing a null proportional hazards model. A numerically efficient Markov chain Monte Carlo (MCMC) based resampling method is proposed to simultaneously obtain the point estimator and a consistent estimator of its variance-covariance matrix. We illustrate our approach with interval censored data sets from two clinical studies. Extensive numerical studies are conducted to evaluate the finite sample performance of the new estimators.

A SIMULATION STUDY OF THE SIZE AND POWER OF SOME TWO-SAMPLE TESTS FOR UNEQUALLY CENSORED SURVIVAL DATA (RANK TESTS, LOGRANK, PROPORTIONAL HAZARDS, WILCOXON TESTS, EXPONENTIAL DISTRIBUTION)

Sizes and power of selected two-sample tests of the equality of survival distributions are compared by simulation for small samples from unequally, randomly-censored exponential distributions. The tests investigated include parametric tests (F, Score, Likelihood, Asymptotic), logrank tests (Mantel, Peto-Peto), and Wilcoxon-Type tests (Gehan, Prentice). Equal sized samples, n = 18, 16, 32 with 1000 (size) and 500 (power) simulation trials, are compared for 16 combinations of the censoring proportions 0%, 20%, 40%, and 60%. For n = 8 and 16, the Asymptotic, Peto-Peto, and Wilcoxon tests perform at nominal 5% size expectations, but the F, Score and Mantel tests exceeded 5% size confidence limits for 1/3 of the censoring combinations. For n = 32, all tests showed proper size, with the Peto-Peto test most conservative in the presence of unequal censoring. Powers of all tests are compared for exponential hazard ratios of 1.4 and 2.0. There is little difference in power characteristics of the tests within the classes of tests considered. The Mantel test showed 90% to 95% power efficiency relative to parametric tests. Wilcoxon-type tests have the lowest relative power but are robust to differential censoring patterns. A modified Peto-Peto test shows power comparable to the Mantel test. For n = 32, a specific Weibull-exponential comparison of crossing survival curves suggests that the relative powers of logrank and Wilcoxon-type tests are dependent on the scale parameter of the Weibull distribution. Wilcoxon-type tests appear more powerful than logrank tests in the case of late-crossing and less powerful for early-crossing survival curves. Guidelines for the appropriate selection of two-sample tests are given. ^

The number and timing of time -dependent covariate measurements for the Cox regression model

The problem of analyzing data with updated measurements in the time-dependent proportional hazards model arises frequently in practice. One available option is to reduce the number of intervals (or updated measurements) to be included in the Cox regression model. We empirically investigated the bias of the estimator of the time-dependent covariate while varying the effect of failure rate, sample size, true values of the parameters and the number of intervals. We also evaluated how often a time-dependent covariate needs to be collected and assessed the effect of sample size and failure rate on the power of testing a time-dependent effect.^ A time-dependent proportional hazards model with two binary covariates was considered. The time axis was partitioned into k intervals. The baseline hazard was assumed to be 1 so that the failure times were exponentially distributed in the ith interval. A type II censoring model was adopted to characterize the failure rate. The factors of interest were sample size (500, 1000), type II censoring with failure rates of 0.05, 0.10, and 0.20, and three values for each of the non-time-dependent and time-dependent covariates (1/4,1/2,3/4).^ The mean of the bias of the estimator of the coefficient of the time-dependent covariate decreased as sample size and number of intervals increased whereas the mean of the bias increased as failure rate and true values of the covariates increased. The mean of the bias of the estimator of the coefficient was smallest when all of the updated measurements were used in the model compared with two models that used selected measurements of the time-dependent covariate. For the model that included all the measurements, the coverage rates of the estimator of the coefficient of the time-dependent covariate was in most cases 90% or more except when the failure rate was high (0.20). The power associated with testing a time-dependent effect was highest when all of the measurements of the time-dependent covariate were used. An example from the Systolic Hypertension in the Elderly Program Cooperative Research Group is presented. ^

The Value Of The 2005 International Society Of Urological Pathology (isup) Modified Gleason Grading System As A Predictor Of Biochemical Recurrence After Radical Prostatectomy.

To compare time and risk to biochemical recurrence (BR) after radical prostatectomy of two chronologically different groups of patients using the standard and the modified Gleason system (MGS). Cohort 1 comprised biopsies of 197 patients graded according to the standard Gleason system (SGS) in the period 1997/2004, and cohort 2, 176 biopsies graded according to the modified system in the period 2005/2011. Time to BR was analyzed with the Kaplan-Meier product-limit analysis and prediction of shorter time to recurrence using univariate and multivariate Cox proportional hazards model. Patients in cohort 2 reflected time-related changes: striking increase in clinical stage T1c, systematic use of extended biopsies, and lower percentage of total length of cancer in millimeter in all cores. The MGS used in cohort 2 showed fewer biopsies with Gleason score ≤ 6 and more biopsies of the intermediate Gleason score 7. Time to BR using the Kaplan-Meier curves showed statistical significance using the MGS in cohort 2, but not the SGS in cohort 1. Only the MGS predicted shorter time to BR on univariate analysis and on multivariate analysis was an independent predictor. The results favor that the 2005 International Society of Urological Pathology modified system is a refinement of the Gleason grading and valuable for contemporary clinical practice.

Cd8+ Tumour-infiltrating Lymphocytes And Cox2 Expression May Predict Relapse In Differentiated Thyroid Cancer.

There is an increasing rate of papillary thyroid carcinomas that may never progress to cause symptoms or death. Predicting outcome and determining tumour aggressiveness could help diminish the number of patients submitted to aggressive treatments. We aimed to evaluate whether markers of the immune system response and of tumour-associated inflammation could predict outcome of differentiated thyroid cancer (DTC) patients. Retrospective cohort study. We studied 399 consecutive patients, including 325 papillary and 74 follicular thyroid carcinomas. Immune cell markers were evaluated using immunohistochemistry, including tumour-associated macrophages (CD68) and subsets of tumour-infiltrating lymphocytes (TIL), such as CD3, CD4, CD8, CD16, CD20, CD45RO, GRANZYME B, CD69 and CD25. We also investigated the expression of cyclooxygenase 2 (COX2) in tumour cells and the presence of concurrent lymphocytic infiltration characterizing chronic thyroiditis. Concurrent lymphocytic infiltration characterizing chronic thyroiditis was observed in 29% of the cases. Among all the immunological parameters evaluated, only the enrichment of CD8+ lymphocytes (P = 0·001) and expression of COX2 (P =0·01) were associated with recurrence. A multivariate model analysis identified CD8+ TIL/COX2 as independent risk factor for recurrence. A multivariate analysis using Cox's proportional-hazards model adjusted for the presence of concurrent chronic thyroiditis demonstrated that the presence of concurrent chronic thyroiditis had no effect on prognostic prediction mediated by CD8+ TIL and COX2. In conclusion, we suggest the use of a relatively simple pathology tool to help select cases that may benefit of a more aggressive approach sparing the majority of patients from unnecessary procedures.

Analysis Of The Contribution Of Immunologically-detectable Her2, Steroid Receptors And Of The Triple-negative" Tumor Status To Disease-free And Overall Survival Of Women With Epithelial Ovarian Cancer."

We assessed associations between steroid receptors including: estrogen-alpha, estrogen-beta, androgen receptor, progesterone receptor, the HER2 status and triple-negative epithelial ovarian cancer (ERα-/PR-/HER2-; TNEOC) status and survival in women with epithelial ovarian cancer. The study included 152 women with primary epithelial ovarian cancer. The status of steroid receptor and HER2 was determined by immunohistochemistry. Disease-free and overall survival were calculated and compared with steroid receptor and HER2 status as well as clinicopathological features using the Cox Proportional Hazards model. A mean follow-up period of 43.6 months (interquartile range=41.4 months) was achieved where 44% of patients had serous tumor, followed by mucinous (23%), endometrioid (9%), mixed (9%), undifferentiated (8.5%) and clear cell tumors (5.3%). ER-alpha staining was associated with grade II-III tumors. Progesterone receptor staining was positively associated with a Body Mass Index≥25. Androgen receptor positivity was higher in serous tumors. In stand-alone analysis of receptor contribution to survival, estrogen-alpha positivity was associated with greater disease-free survival. However, there was no significant association between steroid receptor expression, HER2 status, or TNEOC status, and overall survival. Although estrogen-alpha, androgen receptor, progesterone receptor and the HER2 status were associated with key clinical features of the women and pathological characteristics of the tumors, these associations were not implicated in survival. Interestingly, women with TNEOC seem to fare the same way as their counterparts with non-TNEOC.

Menopausal status: a possible predictive factor for recurrence in women with cancer of the uterine cervix without pelvic lymph node metastasis

Abstract Objectives to evaluate risk factors for recurrence of carcinoma of the uterine cervix among women who had undergone radical hysterectomy without pelvic lymph node metastasis, while taking into consideration not only the classical histopathological factors but also sociodemographic, clinical and treatment-related factors. Study desin This was an exploratory analysis on 233 women with carcinoma of the uterine cervix (stages IB and IIA) who were treated by means of radical hysterectomy and pelvic lymphadenectomy, with free surgical margins and without lymph node metastases on conventional histopathological examination. Women with histologically normal lymph nodes but with micrometastases in the immunohistochemical analysis (AE1/AE3) were excluded. Disease-free survival for sociodemographic, clinical and histopathological variables was calculated using the Kaplan-Meier method. The Cox proportional hazards model was used to identify the independent risk factors for recurrence. Twenty-seven recurrences were recorded (11.6%), of which 18 were pelvic, four were distant, four were pelvic + distant and one was of unknown location. The five-year disease-free survival rate among the study population was 88.4%. The independent risk factors for recurrence in the multivariate analysis were: postmenopausal status (HR 14.1; 95% CI: 3.7-53.6; P < 0.001), absence of or slight inflammatory reaction (HR 7.9; 95% CI: 1.7-36.5; P = 0.008) and invasion of the deepest third of the cervix (HR 6.1; 95% CI: 1.3-29.1; P = 0.021). Postoperative radiotherapy was identified as a protective factor against recurrence (HR 0.02; 95% CI: 0.001-0.25; P = 0.003). (To continue) Postmenopausal status is a possible independent risk factor for recurrence even when adjusted for classical prognostic factors (such as tumour size, depth of tumour invasion, capillary embolisation) and treatment-related factors (period of treatment and postoperative radiotherapy status)

Different Approaches for Modeling Grouped Survival Data: A Mango Tree Study

Interval-censored survival data, in which the event of interest is not observed exactly but is only known to occur within some time interval, occur very frequently. In some situations, event times might be censored into different, possibly overlapping intervals of variable widths; however, in other situations, information is available for all units at the same observed visit time. In the latter cases, interval-censored data are termed grouped survival data. Here we present alternative approaches for analyzing interval-censored data. We illustrate these techniques using a survival data set involving mango tree lifetimes. This study is an example of grouped survival data.

Arrhythmias and mortality after myocardial infarction in diabetic patients - Relationship to diabetes treatment

OBJECTIVE- To assess the relationship between clinical course after acute myocardial infarction (AMI) and diabetes treatment. RESEARCH DESIGN AND METHODS- Retrospective analysis of data from all patients aged 25-64 years admitted to hospitals in Perth, Australia, between 1985 and 1993 with AMI diagnosed according to the International Classification of Diseases (9th revision) criteria was conducted. Short- (28-day) and long-term survival and complications in diabetic and nondiabetic patients were compared. For diabetic patients, 28-day survival, dysrhythmias, heart block, and pulmonary edema were treated as outcomes, and factors related to each were assessed using multiple logistic regression. Diabetes treatment was added to the model to assess its significance. Long-term survival was compared by means of a Cox proportional hazards model. RESULTS- Of 5,715 patients, 745 (12.9%) were diabetic. Mortality at 28 days was 12.0 and 28.1% for nondiabetic and diabetic patients, respectively (P < 0.001); there were no significant drug effects in the diabetic group. Ventricular fibrillation in diabetic patients taking glibenclamide (11.8%) was similar to that of nondiabetic patients (11.0%) but was lower than that for those patients taking either gliclazide (18.0%; 0.1 > P > 0.05) or insulin (22.8%; P < 0.05). There were no other treatment-related differences in acute complications. Long-term survival in diabetic patients was reduced in those taking digitalis and/or diuretics but type of diabetes treatment at discharge had no significant association with outcome. CONCLUSlONS- These results do not suggest that ischemic heart disease should influence the choice of diabetes treatment regimen in general or of sulfonylurea drug in particular.

A Randomized Trial of Sodium Fluoride (60 mg) ± Estrogen in Postmenopausal Osteoporotic Vertebral Fractures

Postmenopausal Caucasian women aged less than 80 years (n = 99) with one or more atraumatic vertebral fracture and no hip fractures, were treated by cyclical administration of enteric coated sodium fluoride (NaF) or no NaF for 27 months, with precautions to prevent excessive stimulation of bone turnover. In the first study 65 women, unexposed to estrogen (-E study), age 70.8 +/- 0.8 years (mean SEM) were all treated with calcium (Ca) 1.0-1.2 g daily and ergocalciferol (D) 0.25 mg per 25 kg once weekly and were randomly assigned to cyclical NaF (6 months on. 3 months off, initial dose 60 mg/day; group F CaD, n = 34) or no NaF (group CaD, n = 3 1). In the second study 34 patients. age 65.5 +/- 1.2 years, on hormone replacement therapy (E) at baseline, had this standardized, and were all treated with Ca and D and similarly randomized (FE CaD, n = 17, E CaD, n = 17) (+E study). The patients were stratified according to E status and subsequently assigned randomly to NaF. Seventy-five patients completed the trial. Both groups treated with NaF showed an increase in lumbar spinal density (by DXA) above baseline by 27 months: FE CaD + 16.2% and F CaD +9.3% (both p = 0.0001). In neither group CaD nor E CaD did lumbar spinal density increase. Peripheral bone loss occurred at most sites in the F CaD group at 27 months: tibia/fibula shaft -7.3% (p = 0.005); femoral shaft -7.1% (p = 0.004); distal forearm -4.0% (p = 0.004); total hip -4.1% (p = 0. 003); and femoral neck -3.5% (p = 0.006). No significant loss occurred in group FE CaD. Differences between the two NaF groups were greatest at the total hip at 27 months but were not significant [p < 0.05; in view of the multiple bone mineral density (BMD) sites, an alpha of 0.01 was employed to denote significance in BMD changes throughout this paper]. Using Cox's proportional hazards model, in the -E study there were significantly more patients with first fresh vertebral fractures in those treated with NaF than in those not so treated (RR = 24.2, p = 0.008, 95% CI 2.3-255). Patients developing first fresh fractures in the first 9 months were markedly different between groups: -23% of F CaD, 0 of CaD, 29% of FE CaD and 0 of E CaD. The incidence of incomplete (stress) fractures was similar in the two NaF-treated groups. Complete nonvertebral fractures did not occur in the two +E groups, there were no differences between groups F CaD and CaD. Baseline BMD (spine and femoral neck) was related to incident vertebral fractures in the control groups (no NaF), but not in the two NaF groups. Our results and a literature review indicate that fluoride salts. if used, should be at low dosage, with pretreatment and co-treatment with a bone resorption inhibitor.

Association between ADAMTS13 polymorphisms and risk of cardiovascular events in chronic coronary disease

Introduction: Association between ADAMTS13 levels and cardiovascular events has been described recently. However, no genetic study of ADAMTS13 in coronary patients has been described. Materials and Methods: Based on related populations frequencies and functional studies, we tested three ADAMTS13 polymorphisms: C1342G (Q448E), C1852G (P618A) and C2699T (A900V) in a group of 560 patients enrolled in the Medical, Angioplasty, or Surgery Study II (MASS II), a randomized trial comparing treatments for patients with coronary artery disease (CAD) and preserved left ventricular function. The incidence of the 5-year end-points of death and death from cardiac causes, myocardial infarction, refractory angina requiring revascularization and cerebrovascular accident was determined for each polymorphim`s allele, genotype and haplotype. Risk was assessed with the use of logistic regression and Cox proportional-hazards model and multivariable adjustment was employed for possible confounders. Results: Clinical characteristics and received treatment of each genotype group were similar at baseline. In an adjusted model for cardiovascular risk variables, we were able to observe a significant association between ADAMTS13 900V variant and an increased risk of death (OR: 1,92 CI: 1,14-3,23, p = 0,015) or death from cardiac cause (OR: 2,67, CI: 1,59-4,49, p = 0,0009). No association between events and ADAMTS13 Q448E or P618A was observed. Conclusions: This first report studying the association between ADAMTS13 genotypes and cardiovascular events provides evidence for the association between ADAMTS13 900V variant and an increased risk of death in a population with multi-vessel CAD. (C) 2009 Elsevier Ltd. All rights reserved.

Differential Persistence of Transmitted HIV-1 Drug Resistance Mutation Classes

Methods. We studied participants with acute and/or early HIV infection and TDR in 2 cohorts (San Francisco, California, and Sao Paulo, Brazil). We followed baseline mutations longitudinally and compared replacement rates between mutation classes with use of a parametric proportional hazards model. Results. Among 75 individuals with 195 TDR mutations, M184V/I became undetectable markedly faster than did nonnucleoside reverse-transcriptase inhibitor (NNRTI) mutations (hazard ratio, 77.5; 95% confidence interval [CI], 14.7-408.2; P < .0001), while protease inhibitor and NNRTI replacement rates were similar. Higher plasma HIV-1 RNA level predicted faster mutation replacement, but this was not statistically significant (hazard ratio, 1.71 log(10) copies/mL; 95% CI, .90-3.25 log(10) copies/mL; P = .11). We found substantial person-to-person variability in mutation replacement rates not accounted for by viral load or mutation class (P < .0001). Conclusions. The rapid replacement of M184V/I mutations is consistent with known fitness costs. The long-term persistence of NNRTI and protease inhibitor mutations suggests a risk for person-to-person propagation. Host and/or viral factors not accounted for by viral load or mutation class are likely influencing mutation replacement and warrant further study.

Prognostic value of NDRG1 and SPARC protein expression in breast cancer patients

An increasing number of studies have shown altered expression of secreted protein acidic and rich in cysteine (SPARC) and N-myc down-regulated gene (NDRG1) in several malignancies, including breast carcinoma; however, the role of these potential biomarkers in tumor development and progression is controversial. In this study, NDRG1 and SPARC protein expression was evaluated by immunohistochemistry on tissue microarrays containing breast tumor specimens from patients with 10 years of follow-up. NDRG1 and SPARC protein expression was determined in 596 patients along with other prognostic markers, such as ER, PR, and HER2. The status of NDRG1 and SPARC protein expression was correlated with prognostic variables and patient clinical outcome. Immunostaining revealed that 272 of the 596 cases (45.6%) were positive for NDRG1 and 431 (72.3%) were positive for SPARC. Statistically significant differences were found between the presence of SPARC and NDRG1 protein expression and standard clinicopathological variables. Kaplan-Meier analysis showed that NDRG1 positivity was directly associated with shorter disease-free survival (DFS, P < 0.001) and overall survival (OS, P < 0.001). In contrast, patients expressing low levels of SPARC protein had worse DFS (P = 0.001) and OS (P = 0.001) compared to those expressing high levels. Combined analysis of the two markers indicated that DFS (P < 0.001) and OS rates (P < 0.001) were lowest for patients with NDRG1-positive and SPARC-negative tumors. Furthermore, NDRG1 over-expression and SPARC down-regulation correlated with poor prognosis in patients with luminal A or triple-negative subtype breast cancer. On multivariate analysis using a Cox proportional hazards model, NDRG1 and SPARC protein expression were independent prognostic factors for both DFS and OS of breast cancer patients. These data indicate that NDRG1 over-expression and SPARC down-regulation could play important roles in breast cancer progression and serve as useful biomarkers to better define breast cancer prognosis.

Optimal patch-leaving behaviour: a case study using the parasitoid Cotesia rubecular

1. Parasitoids are predicted to spend longer in patches with more hosts, but previous work on Cotesia rubecula (Marshall) has not upheld this prediction, Tests of theoretical predictions may be affected by the definition of patch leaving behaviour, which is often ambiguous. 2. In this study whole plants were considered as patches and assumed that wasps move within patches by means of walking or flying. Within-patch and between-patch flights were distinguished based on flight distance. The quality of this classification was tested statistically by examination of log-survivor curves of flight times. 3. Wasps remained longer in patches with higher host densities, which is consistent with predictions of the marginal value theorem (Charnov 1976). tinder the assumption that each flight indicates a patch departure, there is no relationship between host density and leaving tendency. 4. Oviposition influences the patch leaving behaviour of wasps in a count down fashion (Driessen et al. 1995), as predicted by an optimal foraging model (Tenhumberg, Keller & Possingham 2001). 5. Wasps spend significantly longer in the first patch encountered following release, resulting in an increased rate of superparasitism.