请注意以下代码:
Xij <- scan(n=45)
6398400 6273897 6038777 5810740 5673521 5688332 5669445 5682840 5679432
5723561 5555929 5345696 5321179 5199592 5165409 5130744 5132372
4717909 4925673 4999103 4960733 4840036 4824080 4821902
7115151 7114401 7039423 6967723 6967513 6901684
8203359 8286980 8222974 8323470 8067521
5930080 5862383 5994123 6017566
5558436 5754304 5613530
4595506 5074887
3443322
n <- length(Xij); TT <- trunc(sqrt(2*n))
i <- rep(1:TT,TT:1); j <- sequence(TT:1)
i <- as.factor(i); j <- as.factor(j)
如果我现在尝试运行以下命令:
Xij.1 <- xtabs(Xij˜i+j)
我收到错误'错误:“Xij.1&lt; - xtabs(Xij〜”
中的意外输入然而,这个练习类似于“使用R的现代精算风险理论”一书中的一个例子。
有人知道可能出现的问题吗?
答案 0 :(得分:3)
工作正常:
xtabs(Xij~i+j)
请注意,在R公式中,您必须使用波浪号~
而不是˜
个字符。这是两个不同的角色。