Non-muscle invasive bladder cancer is a heterogenous disease whose management is dependent upon the risk of progression to muscle invasion. Although the recurrence rate is high, the majority of tumors are indolent and can be managed by endoscopic means alone. The prognosis of muscle invasion is poor and radical treatment is required if cure is to be obtained. Progression risk in non-invasive tumors is hard to determine at tumor diagnosis using current clinicopathological means. To improve the accuracy of progression prediction various biomarkers have been evaluated. To discover novel biomarkers several authors have used gene expression microarrays. Various statistical methods have been described to interpret array data, but to date no biomarkers have entered clinical practice. Here, we describe a new method of microarray analysis using neurofuzzy modeling (NFM), a form of artificial intelligence, and integrate it with artificial neural networks (ANN) to investigate non-muscle invasive bladder cancer array data (n=66 tumors). We develop a predictive panel of 11 genes, from 2800 expressed genes, that can significantly identify tumor progression (average Logrank p = 0.0288) in the analyzed cancers. In comparison, this panel appears superior to those genes chosen using traditional analyses (average Logrank p = 0.3455) and tumor grade (Logrank, p = 0.2475) in this non-muscle invasive cohort. We then analyze panel members in a new non-muscle invasive bladder cancer cohort (n=199) using immunohistochemistry with six commercially available antibodies. The combination of 6 genes (LIG3, TNFRSF6, KRT18, ICAM1, DSG2 and BRCA2) significantly stratifies tumor progression (Logrank p = 0.0096) in the new cohort. We discuss the benefits of the transparent NFM approach with respect to other reported methods.