Ошибка метода: метод не соответствует isless(::Symbol,::Int64)

Я пытаюсь построить в Джулии образец Гиббса процесса Дирихле и представить каждый кластер как объект. Я сделал следующее:

using BayesianNonparametrics
using DataFrames
using Distances

type cluster
    m::Vector
    Sigma::LinAlg.Cholesky
    nu::Int
    kappa::Int
    nk::Int
  end

function logPredPdf(f::cluster,x::Array{Float64,1}):

  kappa <- f.kappa
  S <- f.Sigma
  mu <- f.m
  nu <- f.nu

  d <- size(x)[1]
  v = nu-d+1
  U = sqrt((1+1/kappa)/v) * LinAlg.lowrankdowndate!(S, sqrt(f.kappa)*mu)[:U]

  x = x - m
  Q = \(transpose(U),x)
  q= vecdot(Q,Q)
  o = -log(1+q/v)*((v+d)/2)
  c = lgamma((v+d)/2)-lgamma(v/2)-(d*log(v*pi)+2*sum(log(diag(U))))/2
  y = c + o

  return y
end

data = readtable("PCA_transformed_data_gt1000.csv",header= true);
data = delete!(data, :1);
n,d = size(data);

s = 6.6172
S0 = s*eye(d)
kappa_0 = 1
nu_0 = d
mu_0 = zeros(d)
S1 = LinAlg.cholfact(S0+kappa_0*(mu_0*mu_0'))

X= DataFrame(Matrix(data)');# transpose data Matrix
prior = cluster(mu_0, S1,nu_0, kappa_0, 0)

x = X[:,1]
kt = logPredPdf(prior,x)

Когда я запускаю строки внутри logPredPdf вне функции, она работает отлично! Но если я пытаюсь запустить приведенный выше пример, когда он идет в конце, он выдает мне следующую ошибку:

MethodError: no method matching isless(::Symbol, ::Int64)
Closest candidates are:
  isless(!Matched::Char, ::Integer) at deprecated.jl:49
  isless(::Symbol, !Matched::Symbol) at strings/basic.jl:137
  isless(!Matched::DataArrays.NAtype, ::Any) at /Users/u1560476/.julia/v0.5/DataArrays/src/operators.jl:510
  ...
 in logPredPdf(::cluster, ::Array{Float64,1}) at Dirichlet_Process_Gibbs_Sampler.jl:33
 in include_string(::String, ::String) at loading.jl:441
 in include_string(::String, ::String, ::Int64) at eval.jl:30
 in include_string(::Module, ::String, ::String, ::Int64, ::Vararg{Int64,N}) at eval.jl:34
 in (::Atom.##53#56{String,Int64,String})() at eval.jl:50
 in withpath(::Atom.##53#56{String,Int64,String}, ::String) at utils.jl:30
 in withpath(::Function, ::String) at eval.jl:38
 in macro expansion at eval.jl:49 [inlined]
 in (::Atom.##52#55{Dict{String,Any}})() at task.jl:60

Типы "prior" и "x" - это cluster и Array{Float64,1}, как требует функция. Есть идеи, что мне здесь не хватает?

Источник

2 ответа

Решение

<- не оператор в Юлии. использование = для назначения в юлию. Ни один пакет (или даже макрос) никогда не позволит вам написать код Julia, который содержит функционал <- оператор. x <- y всегда значит x меньше чем -y:

julia> expand(:(x <- y))
:(x < -y)

julia> Meta.show_sexpr(ans)
(:call, :<, :x, (:call, :-, :y))

Я полагаю, что в вашей глобальной рабочей области эти переменные определены так, как вы ожидали, но внутри этого файла это выглядит так kappa вместо этого связан с символом.

Изменить: Вы также должны удалить : от function logPredPdf(f::cluster,x::Array{Float64,1}):, Определения функций (и блоки в целом) не используют : в юли.

Источник

Меняя все <- в = не решил проблему, но добавил

  let
     global k

после

function logPredPdf(f::cluster,x::Array{Float64,1}):

кажется, чтобы решить проблему! Все еще не знаю почему, хотя!

Источник
Другие вопросы по тегам