최근에 R에서 sum(a*b)를 해야 하는 상황에 직면했다. 문제는 a와 b가 매우 크거나 작은 수라는 것!Recently, I came across the situation in which I had to sum(a*b) in R. The problem was that a and b are vectors of very large or small numbers. 다음의 R implementation은 log a b = log c b/log c a라는 것을 안다면 쉽게 이해할 수 있을 것이다.The implementation of log-sum-exp below must be easily understandable if you know log a b = log c b/log c a # log-sum-exp..
매우 작은 양수(0과 가까운 수)나 매우 큰 수를 나타내는 한 가지 방법은 log를 활용하는 것이다. 이 때 두 수의 곱셈은 자명하다.Very small or large numbers can be easily stored in computers by taking logs. And mutiplying two numbers expressed in log-scale is straight-forward. 하지만 덧셈은 약간의 trick(소위 말해 log-sum-exp trick)이 필요하다.But the addition of two extreme numbers need some trick(so called, log-sum-exp trick). 두 큰 수 exp(800)과 exp(900)을 더하고자 한다. 문제는..
수학자의 무한히 정확한 세계에서는 underflow 따위는 없지만, 컴퓨터에는 있다.In the world of Mathematics where infinite precision is guaranteed, there's no room for Underflow, but computers has it. 최근에 IRT 모형의 일종인 GRM 모형의 log Likelihood 최대화를 하는 과정에서 underflow의 의한 심각한 문제가 발생했다.I recently programmed a routine to maximize the log likelihood of GRM, which is one kind of IRT models and found some problems related to UNDERFLOW! 다음..
R에서 응답자료를 gpcm 모형으로 분석하고자 할 때 ltm 또는 gpcm 패키지를 쓸 수 있다.You can use the package or to fit your data to gpcm. 두 패키지의 가장 큰 차이점은 모수를 다른 방식으로 사용한다는 점이다.The difference between two functions, gpcm and ltm is that they use different parameterization. 이제 두 함수 ltm::gpcm과 mirt::mirt를 통해 주어진 자료를 gpcm에 적합시키자.Now, let's see how we can fit the data to gpcm via and . library("mirt")library("ltm") apply(R.gpcm, 2..
Polychoric, tetrachoric, polyserial, and biserial correlation는 정규 분포를 띄는 잠재변수를 가정합니다.Polychoric, tetrachoric, polyserial, and biserial correlations assume latent variables which are normally distributed. 그리고 잠재변수가 0/1 또는 0/1/2/3/.../n 등의 이분/ordinal 변수로 관찰될 때, 잠재변수 간의 상관계수를 구합니다.And when the observations are either 0/1 or ordinal which is manifestation of continuous latent variables, polychoric, ..
보호되어 있는 글입니다.
http://stackoverflow.com/questions/7019912/using-the-rjava-package-on-win7-64-bit-with-r JAVA_HOME 환경 변수 때문에 생기는 경향? jvm.dll이 있는 폴더 예) C:\Program Files\Java\jdk1.6.0_25\jre\bin\server 또는 C:\Program Files\Java\jre6\jre\bin\client 를 PATH에 포함시킨다. 내 컴퓨터 -(우클릭)- 속성 - 좌: 고급시스템설정 - (고급 탭 아래) 환경설정
보호되어 있는 글입니다.
보호되어 있는 글입니다.