2014-04-25 2 views
1

Я пытаюсь взять набор независимых переменных и проверить, если они (статистически значимо) по-разному коррелированы с двумя группами данных.Итерация через переменные в регрессии в SAS/JMP

Я информирована о том, что способ сделать это в СОМ, чтобы сделать серию линейных регрессий, как в следующем,

result = group + varA + group*varA 

, а затем изучить значимость эффекта взаимодействия, например, " Prob> F "в этом примере« Country * Displacement »: http://i.stack.imgur.com/EcCdd.png (у меня нет репутации публикации изображения.)

Теперь мне нужно иметь возможность отключить одну из этих переменных; то есть, для списка ~ 350 переменных, скажем varA, varB и т.д., мне нужно выполнить следующие регрессии,

result = group + varA + group*varA 
result = group + varB + group*varB 
result = group + varC + group*varC 
... 

и получить значение этого эффекта взаимодействия. Предыдущие попытки создания сценариев привели к появлению ~ 350 окон результатов или диалогов моделей ~ 350. , , Любой совет будет принят во внимание.

Редактировать:

Например, при использовании Авиакомпании Задержки набора данных СПМ образца, это результат от одного из шагов: http://i.stack.imgur.com/HVFL8.png. Мне нужно извлечь значение эффекта взаимодействия (0.1397 в тестах эффектов) для каждого из множества переменных; например, переставляя переменную «Расстояние» с «Истекшее время». Но мне нужно изменить эту переменную для каждого в наборе ~ 350.

+0

Чтобы быть понятным, вы НЕ хотите делать это в базе SAS, только в JMP? – Joe

+0

Я отлично справляюсь с этим в SAS; Я просто не знаю, как подойти к нему. –

+1

Перед вызовом процедуры вы можете использовать несколько операторов модели в «PROC REG» и «ODS SELECT ANOVA». –

ответ

0

Предполагая, что вы знаете, как это сделать. Это даст вам эффект P Values.

fit = Fit Model(
    Y(:Arrival Delay), 
    Effects(:Distance, :Day of Week, :Distance * :Day of Week), 
    Personality(Standard Least Squares), 
    Emphasis(Minimal Report), 
    Run(
     :Arrival Delay << {Lack of Fit(0), Plot Actual by Predicted(0), 
     Plot Residual by Predicted(0), Plot Effect Leverage(0)} 
    ) 
); 

hash = associative array(fit<<Get Effect Names, fit<<Get Effect PValues); 

value = hash["Distance*Day of Week"]; 

затем просто плотное прилегание < < Закрыть окно; и перейдите к следующему параметру.

Смежные вопросы