Блог

Главная - Вопросы по программированию - Попытка перевести SAS GLM в R sasLM ::GLM

Попытка перевести SAS GLM в R sasLM ::GLM

#r #sas #glm

#r #sas #glm

Вопрос:

У меня есть такой набор данных:

 X1 Record  Plot   Row Column    Cp   Csp Entry  Year Location Genotype Trait Value Whole_plot
 <dbl>  <dbl> <dbl> <dbl>  <dbl> <dbl> <dbl> <dbl> <dbl> <chr>    <chr>    <chr> <dbl> <chr>     
 1  3256    717   566     6      7     0     2   717  2019 Preston  Novelty  STD       5 6   7     
 2  3263    716   567     6      7     0     1   716  2019 Preston  Flanders STD       4 6   7     
 3  3893    716   657     7      8     0     1   716  2019 Preston  Flanders STD       2 7   8     
 4  3900    717   658     7      8     0     2   717  2019 Preston  Novelty  STD       2 7   8     
 5  4698    716   772     9      3     0     1   716  2019 Preston  Flanders STD       3 9   3     
 6  4712    717   774     9      3     0     2   717  2019 Preston  Novelty  STD       3 9   3     
 7  3257    717   566     6      7     0     2   717  2019 Preston  Novelty  V1        5 6   7     
 8  3264    716   567     6      7     0     1   716  2019 Preston  Flanders V1        4 6   7     
 9  3894    716   657     7      8     0     1   716  2019 Preston  Flanders V1        3 7   8     
10  3901    717   658     7      8     0     2   717  2019 Preston  Novelty  V1        3 7   8

код SAS является

 PROC GLM DATA=MAD.MDA_plot_subplot_controls outstat=MAD.subplot_control_anova_stat;
    BY Year Location Trait;
    CLASS  Whole_plot Genotype;
    MODEL Value = Whole_plot Genotype/SS3;
    RANDOM Whole_plot; 
    LSMEANS Whole_plot;
RUN;
QUIT;

что я попробовал, так это

   require(sasLM)
  glm(Value ~ factor(Row)   factor(Column),MDA_plot_controls)

но я получаю сообщение об ошибке

 Error in glm(Value ~ factor(Row)   factor(Column), MDA_plot_controls) : 
  'family' not recognized
In addition: Warning message:
Unknown or uninitialised column: `family`.

я не понимаю столбец ‘family’.

Комментарии:

1. рассмотрите возможность использования data = MDA_plot_controls, family = binomial , например

2. Что??? использовать , family = binomial в команде R ? Я пытаюсь преобразовать это в R? Я не понимаю, о чем вы говорите

3. glm(Value ~ factor(Row) factor(Column),binomial(), MDA_plot_controls) это то, что вы ищете

4. Я не знаю пакет sasLM. Что не так с хорошим ole lm()

5. Ага, я пытаюсь воссоздать функциональность SAS PROC GLM… вот почему я пытаюсь использовать sasLM