fitsaemodel

[עמית קורנגוט‎]

שלום,

 

אני מנסה להריץ fitsaemodel והוא כותב לי את השגיאה הבאה:

Error in .initmethod(model, control$init) :

  BLAS/LAPACK routine 'DGELS ' gave error code -8

לא ממש ברור לי למה כי אני פשוט מכניסה לו:

p=data.frame(index=seq(1,n_domain),y=y+u+e,x=x)

  fhmodel=saemodel(formula=y ~ x, area =~index, data = p)

  mlfit=fitsaemodel("ml", fhmodel)

 

 

אשמח לעזרתכם..

תודה,

עמית 

[Jonathan Rosenblatt]

BLAS/LAPACK are the linear algebra routines.

R typically ships with these are it cannot assume their are present in the hosting operating system.

Assuming you are running Windows, you have several options (yes- Linux and Mac users would have had an easier life):

(a) Play with `fitsaemodel()` parameters. Namely- try changing the "method" argument to method="huberm".

Alternatively, look at the control arguments (fitsaemodel.control) and play with the "dec" parameter by changing from the default dec=0 (SVD) to dec=1 (Cholesky). 

(b) You can replace R's BLAS/LAPACK by the more advanced openBLAS and LAPACK, but I am really unsure how to do it in Windows. Google should know better, but only go down this path if you know a thing or two about system administration.

Please let us know if any of these helped. 

--
You received this message because you are subscribed to the Google Groups "Israel R User Group" group.
To unsubscribe from this group and stop receiving emails from it, send an email to israel-r-user-g...@googlegroups.com.
For more options, visit https://groups.google.com/d/optout.

--

Jonathan Rosenblatt
www.john-ros.com

[Yotam Hechtlinger]

תנסי לבדוק אם הבעיה היא עם השורה הראשונה של הקוד שלך או השנייה על ידי זה שתריצי את השורה השנייה עם דטא פריים סטנדרטי פחות או יותר. אם הבעיה עם השורה השנייה, אז תעשי את מה שיונתן אמר, אחרת תפשטי את הדרך בה את מגדירה את הדטא פריים.

--

[עמית קורנגוט‎]

תודה על העזרה! ממה שהבנתי (אחרי מספר משחקים עם הנתונים) אני חושבת שהבעיה נובעת מכך שכדי להשתמש בfitsaemodel הוא צריך data.frame שיש בה לפחות אזור אחד שמופיע פעמיים וזה כנראה בגלל שהוא משתמש במודל B:

basic unit-level SAE model.

אבל אני רוצה להשתמש במודל A- מודל FH ולהפעיל עליו MLE. השאלה היא האם יש אופציה להגיד לו באיזה מודל אני רוצה שישתמש?

תודה,

עמית

בתאריך 31 במרץ 2015 בשעה 21:06, מאת Yotam Hechtlinger <hecht...@gmail.com>:

You received this message because you are subscribed to a topic in the Google Groups "Israel R User Group" group.
To unsubscribe from this topic, visit https://groups.google.com/d/topic/israel-r-user-group/ITfnAwDsdHo/unsubscribe.
To unsubscribe from this group and all its topics, send an email to israel-r-user-g...@googlegroups.com.