Last chance! 50% off unlimited learning
Sale ends in
Number of patents, R&D spending and sales in millions of dollar for 70 pharmaceutical and biomedical companies in 1976.
data("patent")A data frame with 70 observations on the following 4 variables.
Name of company.
Number of patents.
R&D spending per sales.
Logarithmized R&D spendings (in millions of dollars).
The data is taken from the National Bureau of Economic Research R&D Masterfile.
B.H. Hall, C. Cummins, E. Laderman and J. Mundy. The R&D Master File Documentation. Technical Working Paper 72, National Bureau of Economic Research, 1988. Cambridge, MA.
data("patent", package = "flexmix")
patentMix <- initFlexmix(Patents ~ lgRD, k = 3,
model = FLXMRglm(family = "poisson"),
concomitant = FLXPmultinom(~RDS),
nrep = 5, data = patent)
plot(Patents ~ lgRD, data = patent,
pch = as.character(clusters(patentMix)))
ordering <- order(patent$lgRD)
apply(fitted(patentMix), 2, function(y)
lines(sort(patent$lgRD), y[ordering]))
Run the code above in your browser using DataLab