formula公式中的因变量如何递增？？

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

公式formula=log(X1)~1，这里的1需要递增，请问as.formula()的表达式该如何表示？？？我知道formula=X1~1，这种递增的方法，formula=as.formula(paste0('X',i,'~1'))即可。
可是多了个log不会了：我试了formula=as.formula(paste0('log(X',i,')~1'))，报错！！！
请问该如何表示呢？？？

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享0 收藏0 回帖

关键词：Formula form ORM Mul For

可以如下实现方式

1.   x.i <- get(paste0("x",i))
   
2.   fit <- lm(log(xi)~1)

复制代码

你好，感谢您的回复，我试了依然不行。我用的包是
library(raster)
library(rgdal)
library(dismo)
library(gstat)
我想做克里金插值，可是插值单次运行就可以，循环起来就出错，运行代码(出错部分）：library(gstat)
v<- variogram(log(rain) ~ 1, data =dsp)
出错代码：
Error in if (rat > 1e+06 || rat < 1e-06) print("a possible solution MIGHT be to scale semivariances and/or distances") :   需要TRUE/FALSE值的地方不可以用缺少值

也就是建半方差函数的过程中，log(rain)这里我改成我的列名log(X1)~1，因为我有很多数据要算，依次X1,X2...X3652这么多列。这个式子表示就成问题。
我看了帮助，log(X1)~1的部分它是variogram.formula，可能跟一般formula的定义不一样~！！
请用过这个的帮我看看解决方法~谢谢！