Evaluating different methods of microarray data normalization

Abstract Background With the development of DNA hybridization microarray technologies, nowadays it is possible to simultaneously assess the expression levels of thousands to tens of thousands of genes. Quantitative comparison of microarrays uncovers distinct patterns of gene expression, which define different cellular phenotypes or cellular responses to drugs. Due to technical biases, normalization of the intensity levels is a pre-requisite to performing further statistical analyses. Therefore, choosing a suitable approach for normalization can be critical, deserving judicious consideration. Results Here, we considered three commonly used normalization approaches, namely: Loess, Splines and Wavelets, and two non-parametric regression methods, which have yet to be used for normalization, namely, the Kernel smoothing and Support Vector Regression. The results obtained were compared using artificial microarray data and benchmark studies. The results indicate that the Support Vector Regression is the most robust to outliers and that Kernel is the worst normalization technique, while no practical differences were observed between Loess, Splines and Wavelets. Conclusion In face of our results, the Support Vector Regression is favored for microarray normalization due to its superiority when compared to the other methods for its robustness in estimating the normalization curve.

Procalcitonin (PCT) and C-reactive Protein (CRP) as severe systemic infection markers in febrile neutropenic adults

Abstract Background Procalcitonin (PCT) is an inflammatory marker that has been used as indicator of severe bacterial infection. We evaluated the concentrations of PCT as a marker for systemic infection compared to C-reactive protein (CRP) in patients neutropenic febrile. Methods 52 adult patients were enrolled in the study. Blood sample was collected in order to determine the serum concentrations of PCT, CRP and other hematological parameters at the onset of fever. The patients were divided into 2 groups, one with severe infection (n = 26) and the other in which the patients did not present such an infection (n = 26). Then PCT and CRP concentrations at the fever onset were compared between groups using non parametric statistical tests, ROC curve, sensitivity, specificity, likelihood ratio, and Spearman's correlation coefficient. Results The mean of PCT was significantly higher in the group with severe infection (6.7 ng/mL versus 0.6 ng/mL – p = 0.0075) comparing with CRP. Serum concentrations of 0.245 ng/mL of PCT displayed 100% de sensitivity and 69.2% specificity. PCT concentrations of 2,145 ng/mL presented a likelihood ratio of 13, which was not observed for any concentration of CRP. Conclusion PCT seems to be an useful marker for the diagnosis of systemic infection in febrile neutropenic patients, probably better than CRP.

Immunohistochemical evaluation of MMP-2, MMP-9 and CD31/microvascular density in squamous cell carcinomas of the floor of the mouth

The aim of this study was to evaluate the immunoexpression of MMP-2, MMP-9 and CD31/microvascular density in squamous cell carcinomas of the floor of the mouth and to correlate the results with demographic, survival, clinical (TNM staging) and histopathological variables (tumor grade, perineural invasion, embolization and bone invasion). Data from medical records and diagnoses of 41 patients were reviewed. Histological sections were subjected to immunostaining using primary antibodies for human MMP-2, MMP-9 and CD31 and streptavidin-biotin-immunoperoxidase system. Histomorphometric analyses quantified positivity for MMPs (20 fields per slide, 100?points grade, ×200) and for CD31 (microvessels <50?µm in the area of the highest vascularization, 5 fields per slide, 100?points grade, ×400). Statistical design was composed by non-parametric Mann-Whitney U test (investigating the association between numerical variables and immunostainings), chi-square frequency test (in contingency tables), Fisher's exact test (when at least one expected frequency was less than 5 in 2×2 tables), Kaplan-Meier method (estimated probabilities of overall survival) and Iogrank test (comparison of survival curves), all with a significance level of 5%. There was a statistically significant correlation between immunostaining for MMP-2 and lymph node metastasis. Factors associated negatively with survival were N stage, histopathological grade, perineural invasion and immunostaining for MMP-9. There was no significant association between immunoexpression of CD31 and the other variables. The intensity of immunostaining for MMP-2 can be indicative of metastasis in lymph nodes and for MMP-9 of a lower probability of survival

Temporomandibular disorders and parafunctional oral habits: an anamnestic study

OBJECTIVE: To assess the frequency and severity of the signs and symptoms of temporomandibular disorders (TMD), the frequency of parafunctional oral habits and the correlation between the variables by means of the patients' perception regarding their problem. METHODS: One hundred patients diagnosed with TMD, through a clinical examination of their masticatory system, answered the questions of a previously published protocol concerning the signs and symptoms most frequently reported in the literature. RESULTS: According to the results from the non parametric statistical analysis, the frequency for the following signs and symptoms was significant: Fatigue and muscle pain, joint sounds, tinnitus, ear fullness, headache, chewing impairment and difficulty to yawn (p<0.01) and otalgia (p<0.05). As to the parafunctional oral habits, there was a significant presence of teeth clenching during the day and night (p<0.01) and teeth grinding at night (p<0.05). The variable correlation analysis showed that there was a positive correlation between symptom frequency and severity; age was correlated with the presence of otalgia, cervical pain and teeth sensitivity, besides being correlated with muscle and joint pain severity. Habit frequency was negatively correlated with age. TMD duration was also positively correlated with the symptoms of tinnitus, ear fullness, muscle and joint pain. CONCLUSION: The study results showed that the anamnestic assessment using ProDTMMulti can predict the severity of the TMD case.

Comparative study of classic friction among different archwire ligation systems

OBJECTIVE: To describe and compare three alternative methods for controlling classical friction: Self-ligating brackets (SLB), special brackets (SB) and special elastomeric ligatures (SEB). METHODS: The study compared Damon MX, Smart Clip, In-Ovation and Easy Clip self-ligating bracket systems, the special Synergy brackets and Morelli's twin bracket with special 8-shaped elastomeric ligatures. New and used Morelli brackets with new and used elastomeric ligatures were used as control. All brackets had 0.022 x 0.028-in slots. 0.014-in nickel-titanium and stainless steel 0.019 x 0.025-in wires were tied to first premolar steel brackets using each archwire ligation method and pulled by an Instron machine at a speed of 0.5 mm/minute. Prior to the mechanical tests the absence of binding in the device was ruled out. Statistical analysis consisted of the Kruskal-Wallis test and multiple non-parametric analyses at a 1% significance level. RESULTS: When a 0.014-in archwire was employed, all ligation methods exhibited classical friction forces close to zero, except Morelli brackets with new and old elastomeric ligatures, which displayed 64 and 44 centiNewtons, respectively. When a 0.019 x 0.025-in archwire was employed, all ligation methods exhibited values close to zero, except the In-Ovation brackets, which yielded 45 cN, and the Morelli brackets with new and old elastomeric ligatures, which displayed 82 and 49 centiNewtons, respectively. CONCLUSIONS: Damon MX, Easy Clip, Smart Clip, Synergy bracket systems and 8-shaped ligatures proved to be equally effective alternatives for controlling classical friction using 0.014-in nickel-titanium archwires and 0.019 x 0.025-in steel archwires, while the In-Ovation was efficient with 0.014-in archwires but with 0.019 x 0.025-in archwires it exhibited friction that was similar to conventional brackets with used elastomeric ligatures.

A partial dependence factorial analysis to deal with selection bias in observational studis

This thesis presents a creative and practical approach to dealing with the problem of selection bias. Selection bias may be the most important vexing problem in program evaluation or in any line of research that attempts to assert causality. Some of the greatest minds in economics and statistics have scrutinized the problem of selection bias, with the resulting approaches – Rubin’s Potential Outcome Approach(Rosenbaum and Rubin,1983; Rubin, 1991,2001,2004) or Heckman’s Selection model (Heckman, 1979) – being widely accepted and used as the best fixes. These solutions to the bias that arises in particular from self selection are imperfect, and many researchers, when feasible, reserve their strongest causal inference for data from experimental rather than observational studies. The innovative aspect of this thesis is to propose a data transformation that allows measuring and testing in an automatic and multivariate way the presence of selection bias. The approach involves the construction of a multi-dimensional conditional space of the X matrix in which the bias associated with the treatment assignment has been eliminated. Specifically, we propose the use of a partial dependence analysis of the X-space as a tool for investigating the dependence relationship between a set of observable pre-treatment categorical covariates X and a treatment indicator variable T, in order to obtain a measure of bias according to their dependence structure. The measure of selection bias is then expressed in terms of inertia due to the dependence between X and T that has been eliminated. Given the measure of selection bias, we propose a multivariate test of imbalance in order to check if the detected bias is significant, by using the asymptotical distribution of inertia due to T (Estadella et al. 2005) , and by preserving the multivariate nature of data. Further, we propose the use of a clustering procedure as a tool to find groups of comparable units on which estimate local causal effects, and the use of the multivariate test of imbalance as a stopping rule in choosing the best cluster solution set. The method is non parametric, it does not call for modeling the data, based on some underlying theory or assumption about the selection process, but instead it calls for using the existing variability within the data and letting the data to speak. The idea of proposing this multivariate approach to measure selection bias and test balance comes from the consideration that in applied research all aspects of multivariate balance, not represented in the univariate variable- by-variable summaries, are ignored. The first part contains an introduction to evaluation methods as part of public and private decision process and a review of the literature of evaluation methods. The attention is focused on Rubin Potential Outcome Approach, matching methods, and briefly on Heckman’s Selection Model. The second part focuses on some resulting limitations of conventional methods, with particular attention to the problem of how testing in the correct way balancing. The third part contains the original contribution proposed , a simulation study that allows to check the performance of the method for a given dependence setting and an application to a real data set. Finally, we discuss, conclude and explain our future perspectives.

New insights into methodology and interpretation of osteoarthritis. The study of Frassetto identified skeletal collection.

Osteoarthritis (OA) or degenerative joint disease (DJD) is a pathology which affects the synovial joints and characterised by a focal loss of articular cartilage and subsequent bony reaction of the subcondral and marginal bone. Its etiology is best explained by a multifactorial model including: age, sex, genetic and systemic factors, other predisposing diseases and functional stress. In this study the results of the investigation of a modern identified skeletal collection will be presented. In particular, we will focus on the relationship between the presence of OA at various joints. The joint modifications have been analysed using a new methodology that allows the scoring of different degrees of expression of the features considered. Materials and Methods The sample examined comes from the Sassari identified skeletal collection (part of “Frassetto collections”). The individuals were born between 1828 and 1916 and died between 1918 and 1932. Information about sex and age is known for all the individuals. The occupation is known for 173 males and 125 females. Data concerning the occupation of the individuals indicate a preindustrial and rural society. OA has been diagnosed when eburnation (EB) or loss of morphology (LM) were present, or when at least two of the following: marginal lipping (ML), esostosis (EX) or erosion (ER), were present. For each articular surface affected a “mean score” was calculated, reflecting the “severity” of the alterations. A further “score” was calculated for each joint. In the analysis sexes and age classes were always kept separate. For the statistical analyses non parametric test were used. Results The results show there is an increase of OA with age in all the joints analyzed and in particular around 50 years and 60 years. The shoulder, the hip and the knee are the joints mainly affected with ageing while the ankle is the less affected; the correlation values confirm this result. The lesion which show the major correlation with age is the ML. In our sample males are more frequently and more severely affected by OA than females, particularly at the superior limbs, while hip and knee are similarly affected in the two sexes. Lateralization shows some positive results in particular in the right shoulder of males and in various articular surfaces especially of the superior limb of both males and females; articular surfaces and joints are quite always lateralized to the right. Occupational analyses did not show remarkable results probably because of the homogeneity of the sample; males although performing different activities are quite all employed in stressful works. No highest prevalence of knee and hip OA was found in farm-workers respect to the other males. Discussion and Conclusion In this work we propose a methodology to score the different features, necessary to diagnose OA, that allows the investigation of the severity of joint degeneration. This method is easier than the one proposed by Buikstra and Ubelaker (1994), but in the same time allows a quite detailed recording of the features. Epidemiological results can be interpreted quite simply and they are in accordance with other studies; more difficult is the interpretation of the occupational results because many questions concerning the activities performed by the individuals of the collection during their lifespan cannot be solved. Because of this, caution is suggested in the interpretation of bioarcheological specimens. With this work we hope to contribute to the discussion on the puzzling problem of the etiology of OA. The possibility of studying identified skeletons will add important data to the description of osseous features of OA, enriching the medical documentation, based on different criteria. Even if we are aware that the clinical diagnosis is different from the palaeopathological one we think our work will be useful in clarifying some epidemiological as well as pathological aspects of OA.

Identificazione dei fenotipi clinici dell’insufficienza renale cronica che predicono la progressione della malattia. Un’applicazione della tecnica degli alberi decisionali

The primary aim of this dissertation to identify subgroups of patients with chronic kidney disease (CKD) who have a differential risk of progression of illness and the secondary aim is compare 2 equations to estimate the glomerular filtration rate (GFR). To this purpose, the PIRP (Prevention of Progressive Kidney Disease) registry was linked with the dialysis and mortality registries. The outcome of interest is the mean annual variation of GFR, estimated using the Chronic Kidney Disease Epidemiology Collaboration (CKD-EPI) equation. A decision tree model was used to subtype CKD patients, based on the non-parametric procedure CHAID (Chi-squared Automatic Interaction Detector). The independent variables of the model include gender, age, diabetes, hypertension, cardiac diseases, body mass index, baseline serum creatinine, haemoglobin, proteinuria, LDL cholesterol, tryglycerides, serum phoshates, glycemia, parathyroid hormone and uricemia. The decision tree model classified patients into 10 terminal nodes using 6 variables (gender, age, proteinuria, diabetes, serum phosphates and ischemic cardiac disease) that predict a differential progression of kidney disease. Specifically, age <=53 year, male gender, proteinuria, diabetes and serum phosphates >3.70 mg/dl predict a faster decrease of GFR, while ischemic cardiac disease predicts a slower decrease. The comparison between GFR estimates obtained using MDRD4 and CKD-EPI equations shows a high percentage agreement (>90%), with modest discrepancies for high and low age and serum creatinine levels. The study results underscore the need for a tight follow-up schedule in patients with age <53, and of patients aged 54 to 67 with diabetes, to try to slow down the progression of the disease. The result also emphasize the effective management of patients aged>67, in whom the estimated decrease in glomerular filtration rate corresponds with the physiological decrease observed in the absence of kidney disease, except for the subgroup of patients with proteinuria, in whom the GFR decline is more pronounced.

Phylogeny and historical biogeography of the Australian Camphorosmeae (Chenopodiaceae)

Diese Studie befasst sich mit der Phylogenie und Biogeographie der australischen Camphorosmeae, die ein wichtiges Element der Flora arider Gebiete Australiens sind. Die molekularen Phylogenien wurden mit Hilfe Bayes’scher Statistik und „maximum likelihood”berechnet. Um das Alter der Gruppe und interner Linien abzuschätzen, wurden die Methoden „Nonparametric rate smoothing” und “penalized likelihood” benutzt. Morphologische Merkmale wurden nach Kriterien der Parsimonie auf den molekularen Baum aufgetragen. „Brooks parsimony analysis”, „cladistic analysis of distributions and endemism”, „dispersal-vicariance analysis”,„ancestral area analysis” und „weighted ancestral area analysis” wurden angewandt, um Abfolge und Richtungen der Ausbreitung der Gruppe in Australien zu analysieren.Von sieben getesteten Markern hatten nur die nukleären ETS und ITS genügend Variation für die phylogenetische Analyse der Camphorosmeae. Die plastidären Marker trnL-trnF spacer,trnP-psaJ spacer, rpS16 intron, rpL16 intron und trnS-trnG spacer zeigten kein ausreichendes phylogenetisches Signal. Die gefundenen phylogenetischen Hypothesen widersprechen der jetzigen Taxonomie der Gruppe. Neobassia, Threlkeldia, Osteocarpum und Enchylaena sollten den Gattungen Sclerolaena bzw. Maireana zugeordnet werden. Die kladistische Analyse der Fruchtanhängsel unterstützt die taxonomischen Ergebnisse der auf DNA basierenden Phylogenie. Allerdings hat die Behaarung, die bei anderen Gruppen der Chenopodiaceae als wichtiges taxonomisches Merkmal herangezogen wird, die Phylogenie nicht unterstützt. Vorfahren der heutigen Camphorosmeen sind im Miozän, vor ca. 8-14 Millionen Jahren, durch Fernausbreitung vermutlich aus Asien in Australien eingewandert. Anfängliche Diversifizierung fand während des späten Miozäns bis in das frühe Pliozän vor ca. 4-7 Millionen Jahren statt. Am Ende des Pliozäns existierten schon 45% - 72% der Abstammungslinien der jetzigen Camphorosmeen. Dies weist auf eine schnelle Ausbreitung hin. Das Alter stimmt mit dem Einsetzen der Aridisierung Australiens überein, und deutet darauf hin, dass die Ausbreitung der ariden Gebiete eine große Rolle bei der Diversifizierung der Gruppe spielte. Die Vorfahren der australischen Camphorosmeae scheinen die Südküste Australiens zuerst besiedeln zu haben. Dies geschah vor dem Einsetzen der Aridisierung des Kontinents. Die anschließende Ausbreitung erfolgte in verschiedene Richtungen und folgte der fortschreitenden Austrocknung im späten Tertiär und im ganzen Quartär. Durch ihre Anpassung an Trockenheit ist der Erfolg der Camphorosmeae in den ariden Gebieten zu erklären.Die Abwesenheit von klaren phylogenetischen und artspezifischen Signalen zwischen Arten der australischen Camphorosmeae ist auf das junge Alter und die schnelle Diversifizierung der Gruppe zurückzuführen, welche die Häufung von Mutationen und eine starke morphologische Differenzierung nicht zugelassen haben.

Analyzing convergence of Agricultural and Rural Areas in the Enlarged EU: The case of Hungary

The fall of the Berlin Wall opened the way for a reform path – the transition process – which accompanied ten former Socialist countries in Central and South Eastern Europe to knock at the EU doors. By the way, at the time of the EU membership several economic and structural weaknesses remained. A tendency towards convergence between the new Member States (NMS) and the EU average income level emerged, together with a spread of inequality at the sub-regional level, mainly driven by the backwardness of the agricultural and rural areas. Several progresses were made in evaluating the policies for rural areas, but a shared definition of rurality is still missing. Numerous indicators were calculated for assessing the effectiveness of the Common Agricultural Policy and Rural Development Policy. Previous analysis on the Central and Eastern European countries found that the characteristics of the most backward areas were insufficiently addressed by the policies enacted; the low data availability and accountability at a sub-regional level, and the deficiencies in institutional planning and implementation represented an obstacle for targeting policies and payments. The next pages aim at providing a basis for understanding the connections between the peculiarities of the transition process, the current development performance of NMS and the EU role, with particular attention to the agricultural and rural areas. Applying a mixed methodological approach (multivariate statistics, non-parametric methods, spatial econometrics), this study contributes to the identification of rural areas and to the analysis of the changes occurred during the EU membership in Hungary, assessing the effect of CAP introduction and its contribution to the convergence of the Hungarian agricultural and rural. The author believes that more targeted – and therefore efficient – policies for agricultural and rural areas require a deeper knowledge of their structural and dynamic characteristics.

Risikomaße in der Finanzmathematik

„Risikomaße in der Finanzmathematik“ Der Value-at -Risk (VaR) ist ein Risikomaß, dessen Verwendung von der Bankenaufsicht gefordert wird. Der Vorteil des VaR liegt – als Quantil der Ertrags- oder Verlustverteilung - vor allem in seiner einfachen Interpretierbarkeit. Nachteilig ist, dass der linke Rand der Wahrscheinlichkeitsverteilung nicht beachtet wird. Darüber hinaus ist die Berechnung des VaR schwierig, da Quantile nicht additiv sind. Der größte Nachteil des VaR ist in der fehlenden Subadditivität zu sehen. Deswegen werden Alternativen wie Expected Shortfall untersucht. In dieser Arbeit werden zunächst finanzielle Risikomaße eingeführt und einige ihre grundlegenden Eigenschaften festgehalten. Wir beschäftigen uns mit verschiedenen parametrischen und nichtparametrischen Methoden zur Ermittlung des VaR, unter anderen mit ihren Vorteilen und Nachteilen. Des Weiteren beschäftigen wir uns mit parametrischen und nichtparametrischen Schätzern vom VaR in diskreter Zeit. Wir stellen Portfoliooptimierungsprobleme im Black Scholes Modell mit beschränktem VaR und mit beschränkter Varianz vor. Der Vorteil des erstens Ansatzes gegenüber dem zweiten wird hier erläutert. Wir lösen Nutzenoptimierungsprobleme in Bezug auf das Endvermögen mit beschränktem VaR und mit beschränkter Varianz. VaR sagt nichts über den darüber hinausgehenden Verlust aus, während dieser von Expected Shortfall berücksichtigt wird. Deswegen verwenden wir hier den Expected Shortfall anstelle des von Emmer, Korn und Klüppelberg (2001) betrachteten Risikomaßes VaR für die Optimierung des Portfolios im Black Scholes Modell.

Multidimensional Measures of Firm Competitiveness: a Model-Based Approach

The concept of competitiveness, for a long time considered as strictly connected to economic and financial performances, evolved, above all in recent years, toward new, wider interpretations disclosing its multidimensional nature. The shift to a multidimensional view of the phenomenon has excited an intense debate involving theoretical reflections on the features characterizing it, as well as methodological considerations on its assessment and measurement. The present research has a twofold objective: going in depth with the study of tangible and intangible aspect characterizing multidimensional competitive phenomena by assuming a micro-level point of view, and measuring competitiveness through a model-based approach. Specifically, we propose a non-parametric approach to Structural Equation Models techniques for the computation of multidimensional composite measures. Structural Equation Models tools will be used for the development of the empirical application on the italian case: a model based micro-level competitiveness indicator for the measurement of the phenomenon on a large sample of Italian small and medium enterprises will be constructed.

Solutions to common issues in widefield microscopy: vignetting, mosaicing and depth of focus

In this thesis we have developed solutions to common issues regarding widefield microscopes, facing the problem of the intensity inhomogeneity of an image and dealing with two strong limitations: the impossibility of acquiring either high detailed images representative of whole samples or deep 3D objects. First, we cope with the problem of the non-uniform distribution of the light signal inside a single image, named vignetting. In particular we proposed, for both light and fluorescent microscopy, non-parametric multi-image based methods, where the vignetting function is estimated directly from the sample without requiring any prior information. After getting flat-field corrected images, we studied how to fix the problem related to the limitation of the field of view of the camera, so to be able to acquire large areas at high magnification. To this purpose, we developed mosaicing techniques capable to work on-line. Starting from a set of overlapping images manually acquired, we validated a fast registration approach to accurately stitch together the images. Finally, we worked to virtually extend the field of view of the camera in the third dimension, with the purpose of reconstructing a single image completely in focus, stemming from objects having a relevant depth or being displaced in different focus planes. After studying the existing approaches for extending the depth of focus of the microscope, we proposed a general method that does not require any prior information. In order to compare the outcome of existing methods, different standard metrics are commonly used in literature. However, no metric is available to compare different methods in real cases. First, we validated a metric able to rank the methods as the Universal Quality Index does, but without needing any reference ground truth. Second, we proved that the approach we developed performs better in both synthetic and real cases.

Clinicopathologic findings in thoroughbred horses during a high speed – short term competition

The aim of this study was to examine whether a real high speed-short term competition influences clinicopathological data focusing on muscle enzymes, iron profile and Acute Phase Proteins. 30 Thoroughbred racing horses (15 geldings and 15 females) aged between 4-12 years (mean 7 years), were used for the study. All the animals performed a high speed-short term competition for a total distance of 154 m in about 12 seconds, repeated 8 times, within approximately one hour (Niballo Horse Race). Blood samples were obtained 24 hours before and within 30 minutes after the end of the races. On all samples were performed a complete blood count (CBC), biochemical and haemostatic profiles. The post-race concentrations for the single parameter were corrected using an estimation of the plasma volume contraction according to the individual Alb concentration. Data were analysed with descriptive statistics and the percentage of variation from the baseline values were recorded. Pre- and post-race results were compared with non-parametric statistics (Mann Whitney U test). A difference was considered significant at p<0.05. A significant plasma volume contraction after the race was detected (Hct, Alb; p<0.01). Other relevant findings were increased concentrations of muscular enzymes (CK, LDH; p<0.01), Crt (p<0.01), significant increased uric acid (p<0.01), a significant decrease of haptoglobin (p<0.01) associated to an increase of ferritin concentrations (p<0.01), significant decrease of fibrinogen (p<0.05) accompanied by a non-significant increase of D-Dimers concentrations (p=0.08). This competition produced relevant abnormalities on clinical pathology in galloping horses. This study confirms a significant muscular damage, oxidative stress, intravascular haemolysis and subclinical hemostatic alterations. Further studies are needed to better understand the pathogenesis, the medical relevance and the impact on performance of these alterations in equine sport medicine.

The impact of grants in the agriculture sector in Albania: a counterfactual approach

The objective of this study is to measure the impact of the national subsidy scheme on the olive and fruit sector in two regions of Albania, Shkodra and Fier. From the methodological point of view, we use a non- parametric approach based on the propensity score matching. This method overcomes problem of the missing data, by creating a counterfactual scenario. In the first step, the conditional probability to participate in the program was computed. Afterwards, different matching estimators were applied to establish whether the subsidies have affected sectors performance. One of the strengths of this study stays in the data. Cross-sectional primary data was gathered through about 250 interviews.. We have not found empirical evidence of significant effects of government aid program on production. Differences in production found between beneficiaries and non-beneficiaries disappear after adjustment by the conditional probability of participating into the program. This suggests that subsidized farmers would have performed better than the subsidized households even in the absence of production grants, revealing program self-selection. On the other hand, the scheme has affected positively the farm structure increasing the area under cultivation, but yields has not increased for beneficiaries compared to non beneficiaries. These combined results shed light on the reason of the missed impact. It could be reasonable to believe that the new plantation, in particular in the case of olives, has not yet reached full production. Therefore, we have reasons to believe on positive impacts in the future. Concerning some qualitative results, the extension of area under cultivation is strongly conditioned by the small farm size. This together with a thin land market makes extremely difficult the expansion beyond farm boundaries.