Data analysis for flow-based in-vitro receptomics array, like a tongue-on-a-chip, is complicated by the relatively large variability within and between arrays, transfected DNA types, spots, and cells within spots. Simply averaging responses of spots of the same type would lead to high variances and low statistical power. This paper presents an approach based on linear mixed models, allowing a quantitative and robust comparison of complex samples and indicating which receptors are responsible for any differences. These models are easily extended to take into account additional effects such as the build-up of cell stress and to combine data from replicated experiments. The increased analytical power this brings to receptomics research is discussed.