用R拟合市场单位销售额S曲线

我的任务是使用几个自变量，根据宏观经济数据拟合一些市场销售曲线；我有14年的历史。样本数据如下：

UnitSales,GDP,GDPPerCap,CPI,PropInvIndex,DispIncTopDecile,TransCommSecDecile,CivilVehOwn,Urban,AutoFin
" 2,833 ",1198243.4,949,81.62,4984,10643,618,1609,36.22,0
" 7,607 ",1324337.8,1042,81.38,6344,14219,782,1802,37.66,0
" 10,098 ",1453827.558,1135,80.8,7790.9223,18995.9,991.2,2053.17,39.08978381,0
" 18,649 ",1640958.735,1274,81.83,10153.8009,21837.3,1106,2382.93,40.53022975,0
" 14,525 ",1931644.33,1409,85.02,13158.2516,25377.2,1274.2,2693.71,41.76000862,0
" 19,149 ",2256902.591,1731,86.56,15909.2471,28773.1,1590.3,3160,42.98999663,0
" 22,474 ",2712950.885,2069,87.83,19422.9174,31967.3,1801,3697.3531,44.34301016,0
" 32,573 ",3494055.942,2651,92.02,25288.8373,36784.5,2467.7,4358.355,45.8892446,0.1
" 42,186 ",4521827.271,3414,97.45,31203.1942,43613.8,2632.9,5099.6094,46.98950317,0.12
" 52,765 ",4990233.519,3749,96.76,36241.808,46826.1,3181.9,6280.6086,48.34170101,0.14
" 140,855 ",5930502.27,4433,100,48259.403,51431.6,3630.6,7801.8259,49.94966105,0.16
" 88,873 ",7321891.955,5447,105.45,61796.8858,58841.9,3963,9356.3163,51.27027127,0.18
" 161,082 ",8229490.03,6093,108.22,71803.7869,63824.2,4304.1,10933.0912,52.57008656,0.22
" 149,197 ",9240270.452,6807,111.07,86013.3826,68874.37491,4748.669569,12670.1435,53.73,0.26

在我掌握的宏观经济变量中，如果自由度成为一个问题，我可以只使用一个子集。由于数据是遵循生命周期并具有市场饱和效应的产品，因此曲线应为S曲线或S形曲线。在我尝试R时，我只能拟合指数（J曲线）曲线；我已经尝试过tslm（来自fpp软件包）来适应这个模型

我现在的R代码是

# rm(list=ls())
ds <- read.csv(file='c:/RCode/Output/PremM.csv', header=T, sep=',')
tsdata <-ts(ds, s=2000)
results<- tslm(UnitSales~GDP+GPDPerCap+CPI+trend, data=tsdata)

write.csv(coef(results), file = "c:/Rcode/Output/CoefPremM.csv")

sink('C:/RCode/Output/PremMSummary2.txt',append=F,type="output")
summary(results)
sink(NULL)

#rm（list=ls（））
ds不发布
rm
语句。（我已经注释掉了。）示例代码中还有多个文件的读取语句，但只显示了一个数据文件。请把这个清理干净，使一切保持一致。此外，“不工作”也不是一个有用的描述。问题到底是什么还不清楚。谢谢你的反馈；我已经编辑了我的问题，删除了无关的部分。我的问题是我需要拟合一条S曲线，但tslm正在生成一条J曲线。R中的哪个包最适合拟合产品市场饱和的S曲线。请参见基准R中的
？nls
和
？SSlogis
？目前，结果变量可能被解读为一个因子，而它应该是一个计数。首先，请尝试将ds$UnitSales=As.numeric（gsub（“[，]”，“”），As.character（ds$UnitSales）