通过XRJulia将公式从R传递到Julia

Question

library(XRJulia)
findJulia(test=T) # Works fine
juliaEval("using MixedModels")

jlmerj<-juliaEval("
 function(f,d)
   m=fit(LinearMixedModel,f,d)
   return(m)
 end
 ")
jlmer=JuliaFunction(jlmerj)

jdata<-juliaSend(data[,c('IRI_EC','learnid')])
jfrm<-juliaEval("@formula(IRI_EC ~  1 + (1|learnid))")

> jdata
Julia proxy object
Server Class: Dict{String,Any}; size: 7
> jfrm
Julia proxy object
Server Class: StatsModels.Formula; size: NA

res<-jlmer(jfrm,jdata)
Error: Evaluating Julia expression: R_1_7(R_1_9, R_1_29)

有人这样做过吗？怎么做？

Answer 1

你能给出一个示例数据集来运行你的代码吗？

再次检查后，我怀疑问题不在于公式，而在于xrjulia中Rdata.frame的转换，因为在您的消息中可以看到jdata是dict(dictionary)的代理，而不是julia中的dataframe或类似的东西。我还（简要地）签入了MixedModels，这个包似乎不能将dict作为数据集处理。

JuliaCall的用法有点类似于xrjulia，但它将Rdata.frame转换为juliadataframe。在您的示例中使用xrjulia对我来说似乎没问题。因此，JuliaCall可以使用，但XRJulia不能使用，这也意味着问题出在XRJulia中data.frame的转换上。

library(XRJulia)
findJulia(test = T) # Works fine
juliaEval("using MixedModels")

juliaEval("using DataFrames")

jlmerj <- juliaEval("
                  function(f,d)
                  m=fit(LinearMixedModel,f,d)
                  return(m)
                  end
                  ")
jlmer <- JuliaFunction(jlmerj)

mindata <- data.frame(IRI_EC = c(15, 14, 27, 0, 22, 16, 23, 17, 20, 26), learnid = factor(1:10))

jfrm <- juliaEval("@formula(IRI_EC ~  1 + (1|learnid))")

jIRI_EC <- juliaSend(mindata$IRI_EC) ## send columns one by one
## first convert factors to integers and convert back in julia
jlearnid <- juliaCall("CategoricalArrays.CategoricalArray", juliaSend(as.integer(mindata$learnid))

## combine columns to create dataframe in julia
jdata <- juliaEval(paste0('DataFrame(Dict([(:IRI_EC,', juliaName(jIRI_EC),
                          '), (:learnid, ', juliaName(jlearnid), ')]))'))

res <- jlmer(jfrm,jdata)

res