R'deki karışık model formüllerinde rastgele etkiler için (1 | id) gibi Wilkinson tarzı gösterimin kökeni


16

R'deki model formüller

y ~ x + a*b + c:d

Wilkinson notasyonuna dayanmaktadır : Wilkinson ve Rogers 1973, Varyans Analizi için Faktöriyel Modellerin Sembolik Açıklaması .

Bu makale karışık modeller (o zamanlar mevcut olmayabilir) için gösterimleri tartışmadı. Peki lme4R'de karışık paket formüller ve ilgili paketler nerede kullanıldı

y ~ x + a*b + c:d + (1|school) + (a*b||town)

dan geliyorum? Onları ilk kez kim ve ne zaman tanıttı? Onlar için "Wilkinson yazısı" gibi üzerinde anlaşılan herhangi bir terim var mı? Özellikle aşağıdaki gibi terimlere atıfta bulunuyorum

(model formula |  grouping variable)
(model formula || grouping variable)

Yanıtlar:


11

Gösterim 3.1-1 sürümünden beri dokümanlar |etrafında nlmeve muhtemelen 1999 sonlarında; CRAN nlme kod arşivinde kolayca kontrol edebiliriz . nlmebu gösterimi kullanır, örneğin deneyin library(nlme); formula(Orthodont); |çıkageldi - yani 2000'li kapalıdır. Öyleyse yapalım .... " Çoklu Yerleştirme Düzeylerine Sahip Veriler için Grafiksel Yöntemler " Pinheiro & Bates (1997) kurucunun groupedDatatanıtıldığı yer. Ve şöyle derler: " Gruplanmış bir veri nesnesindeki formül, S-PLUS'ta bir Trellis grafik işlevine yapılan çağrıda kullanılan xyplot gibi formülle aynı desene sahiptir. " Bu .... P&B'nin çalıştığı gibi mantıklıdır. .. Aslında operatörü kullanan Trellis grafik sistemini geliştiren Bell Labs (RIP)|zaten grupları belirtmek için. Muhtemelen ... Becker ve ark. Tarafından " Kafes Ekranının Görsel Tasarımı ve Kontrolü ". (1996) bununla bir ilgisi var. Gösterim bu makalede tanıtılmamıştır, ancak bulabildiğim ilk elektronik Trellis ekran referansıdır.

Aslında bu noktada görselleştirme literatürünü araştırmamız gerekiyor. Muhtemelen Cleveland'ın Görselleştirme Verileri (1993) kitabını ve Deepayan Sarkar'ın (geliştiren lattice) ilk eserlerini inceleyeceğim . Gerçek operatörün | (ve ||) operatörlerle ilişkilendirildikleri için gerçek ilkel operatörler olduğuna dikkat edin OR, bu nedenle birisi onları aşırı yükleyene kadar sadece bir zaman meselesiydi. Tam bir cevap olmasa da, P&B'nin meslektaşlarına havalı görselleştirme sistemini kontrol ettiğinden şüpheleniyorum (1996 tarihli makaledeki parseller 2010'un sonlarında standartlar için oldukça iyi) ve birisinin (Becker, Cleveland ve Shyu) zaten bu konuda bazı çalışmalar yaptığını fark ettim (belki hatta o zaman onlarla tartıştı) ve orada olanı takip ettim. Yani|operatör grafik gösteriminden kaynaklanır. Kafes neredeyse kesinlikle kullandı; Trellis'in potansiyel öncülleri de bunu yapmış olabilir, ancak e-ayak izlerini izlemek çok zor.

Genel olarak, bu sayfanın NLME: Bell Labs'ın karma efekt modelleri için daha tarihi bilgiler için yazılım isteyebileceğini düşünüyorum nlme.


1
Çok teşekkürler! nlmeKullanan doğru |ama ()rastgele etkileri göstermek için kullandığını sanmıyorum , değil mi? Rastgele efektler, işlev çağrısına ayrı bir bağımsız değişken olarak listelenir. O oldu lme4tanıttığı (x|id)bir parçası olarak aynı formülü?
amip diyor Reinstate Monica

2
Bence parantezin varlığını biraz fazla okuyorsunuz; Ben lme4tüm terimler için birleşik bir sözdizimi kullanır verilen ayrıştırma amaçları için var şüpheli . Örneğin fm1 <- lmer(Reaction ~ Days | Subject, sleepstudy)parantez olmadan sorunsuz çalışır.
usεr11852 diyor Reinstate Monic

Ah. Aslında. Asla bu şekilde düşünmedim :)
amip Reinstate Monica diyor

Bates sadece rastgele parça nlme yazarların (bunlardan o biridir) tarafından tanıtıldı doğruladı: twitter.com/BatesDmbates/status/1111283948615802881
Jonas Lindeløv

@ JonasLindeløv: Harika! Paylaştığın için teşekkürler, yarın gece cevapla bir bağlantı kuracağım.
usεr11852 diyor Reinstate Monic
Sitemizi kullandığınızda şunları okuyup anladığınızı kabul etmiş olursunuz: Çerez Politikası ve Gizlilik Politikası.
Licensed under cc by-sa 3.0 with attribution required.