Extracting coefficients for a particular factor

[Asaf Weinstein]

Hi,

I used lm() on ANOVA data. 

Is there an easy way to extract coefficients corresponding to a particular factor?

Thanks a lot,

Asaf

[Jonathan Rosenblatt]

you can always do coef(lm.object)

and then extract from that vector by index or by name.

are you looking for something else?

--
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-group+unsub...@googlegroups.com.
For more options, visit https://groups.google.com/d/optout.

--

--

Jonathan Rosenblatt

Dept. of Industrial Engineering and Management

Ben Gurion University of the Negev

www.john-ros.com

[Diklah Geva]

האם אתה מתכוון ל
Coeff()
Lmfit$coefficients

--

[Asaf Weinstein]

תודה, 

הכוונה היתה אם אפשר לחלץ את המקדמים לפי משתנה קטגוריאלי (לצורך המחשה, אם התאמתי

(lm(yield ~ fertilizer + variety

ו-variety משתנה עם 3 רמות ו-fertilizer עם 5 רמות, הייתי רוצה לקבל את ה-coeff רק של המקדמים שמתאימים ל-variety. 

תודה

On 20 September 2016 at 03:34, Diklah Geva <dikla...@gmail.com> wrote:

האם אתה מתכוון ל
Coeff()
Lmfit$coefficients

On 31 Aug 2016 22:02, "Asaf Weinstein" <assafwe...@gmail.com> wrote:

Hi,

I used lm() on ANOVA data. 

Is there an easy way to extract coefficients corresponding to a particular factor?

Thanks a lot,

Asaf

--
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-group+unsubscribe...@googlegroups.com.

For more options, visit https://groups.google.com/d/optout.

[Tal Galili]

שלום אסף,

אתה רוצה רק את המקדמים, או את שאר הטבלה?

כך או אחרת, נראה הגיוני שמה שתעשה פה זה תשתמש ב coeff, ותצמצמם אותו רק לערכים שרלוונטיים למשתנה שמעניין אותך (על ידי שימוש ב grepl על ה names של הפלט של coeff).

----------------Contact Details:-------------------------------------------------------
Contact me: Tal.G...@gmail.com | 
Read me: www.talgalili.com (Hebrew) | www.biostatistics.co.il (Hebrew) | www.r-statistics.com (English)
----------------------------------------------------------------------------------------------

[amit gal]

לדעתי הבעיה חורגת מהתחום הטכני גרידא.

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

אם יש רק משתנה אחד, אחד הפתרונות הטכניים הוא להעיף מהמודל את המקדם החופשי - כלומר לרשום משהו כמו y~0+x, אבל גם זה מתבחבש אם יש יותר ממשתנה אחד.

אז אם הפתרון הוא לא טכני גרידא, אז נשאלת השאלה לשם מה צריכים את המקדמים (שאלה, לאסף). הרבה פעמים בכלל רוצים להשוות בין הקבוצות, ואז המשתנים עצמם לא מעניינים אלא ההפרשים ביניהם. ואז הפתרון הנכון הוא להסתכל על קונטרסטים בין מקדמים ולא המקדמים עצמם. לדעתי אנליזת פוסט-הוק של טוקי או משהו דומה יכולה לעזור.

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

[Tal Galili]

אם המטרה היא לדעת את המובהקות של כל המשתנה הקטגוריאלי, אז הגיוני להשתמש בפונקציה anova על אובייקט הרגרסיה (מה שיתן את המובהקות של הסרת כל המשתנה הקטגוריאלי מהמודל). הפונקציה Anova בחבילה car מכילה אופציות נוספות לשליטה בסוג ה type of errors שרוצים לעבוד איתו.

----------------Contact Details:-------------------------------------------------------
Contact me: Tal.G...@gmail.com | 
Read me: www.talgalili.com (Hebrew) | www.biostatistics.co.il (Hebrew) | www.r-statistics.com (English)
----------------------------------------------------------------------------------------------