COMMunication
RESearch.NET

User Tools

Site Tools

You are here: index » sequential_regression
Trace: • sequential_regression
sequential_regression

Differences

This shows you the differences between two versions of the page.

Link to this comparison view

Both sides previous revisionPrevious revision
Next revision		Previous revision
sequential_regression [2024/06/11 23:22] – [r] hkimscilsequential_regression [2026/06/14 23:47] (current) – hkimscil
Line 14: Line 14:
 | 9   | 320   | 1  |  | 9   | 320   | 1  | 
 | 9   | 270   | 3  |  | 9   | 270   | 3  | 
-<code>datavar <- read.csv("http://commres.net/wiki/_media/regression01-bankaccount.csv")</code>+<code>datavar <- read.csv("http://commres.net/_media/regression01-bankaccount.csv")</code>
  
 ====== Enter ====== ====== Enter ======
Line 109: Line 109:
 ====== r ====== ====== r ======
 <code> <code>
-datavar <- read.csv("http://commres.net/wiki/_media/regression01-bankaccount.csv")+datavar <- read.csv("http://commres.net/_media/regression01-bankaccount.csv")
 datavar  datavar 
 m1 <- lm(bankaccount~income+famnum, data=datavar) m1 <- lm(bankaccount~income+famnum, data=datavar)
Line 119: Line 119:
 </code> </code>
 <code> <code>
-> datavar <- read.csv("http://commres.net/wiki/_media/regression01-bankaccount.csv")+> datavar <- read.csv("http://commres.net/_media/regression01-bankaccount.csv")
 > datavar  > datavar 
    bankaccount income famnum    bankaccount income famnum
Line 279: Line 279:
  
 <code> <code>
-pcor.test(datavar$bankaccount, datavar$income, datavar$famnum) +pp.b.i <- pcor.test(datavar$bankaccount, datavar$income, datavar$famnum) 
-pcor.test(datavar$bankaccount, datavar$famnum, datavar$income)+p.b.i 
 +p.b.i$estimate
  
-spcor.test(datavar$bankaccount, datavar$income, datavar$famnum) +p.b.f <- pcor.test(datavar$bankaccount, datavar$famnum, datavar$income) 
-spcor.test(datavar$bankaccount, datavar$famnum, datavar$income)+p.b.f 
 +p.b.f$estimate 
 + 
 +sp.b.i <- spcor.test(datavar$bankaccount, datavar$income, datavar$famnum) 
 +sp.b.i 
 +sp.b.i$estimate 
 +sp.b.f <- spcor.test(datavar$bankaccount, datavar$famnum, datavar$income) 
 +sp.b.f 
 +sp.b.f$estimate 
 + 
 + 
 +zc.b.i <- cor(datavar$bankaccount, datavar$income) 
 +zc.b.i  
 +zc.b.f <- cor(datavar$bankaccount, datavar$famnum) 
 +zc.b.f 
 + 
 +zc.b.i^2 - (sp.b.i$estimate)^2 
 +zc.b.f^2 - (sp.b.f$estimate)^2
  
-cor(datavar$bankaccount, datavar$income) 
-cor(datavar$bankaccount, datavar$famnum) 
 </code> </code>
 . . .  . . . 
 <code> <code>
-> pcor.test(datavar$bankaccount, datavar$income, datavar$famnum) +> pp.b.i <- pcor.test(datavar$bankaccount, datavar$income, datavar$famnum) 
-   estimate    p.value statistic  n gp  Method +> p.b.i 
-1 0.7825112 0.01267595  3.325102 10  1 pearson +  estimate p.value statistic  n gp  Method 
-> pcor.test(datavar$bankaccount, datavar$famnum, datavar$income) +1   0.7825 0.01268     3.325 10  1 pearson 
-   estimate    p.value statistic  n gp  Method +> p.b.i$estimate 
-1 -0.672856 0.04702022 -2.406425 10  1 pearson +[1] 0.7825
-> +
-> +
-> spcor.test(datavar$bankaccount, datavar$income, datavar$famnum) +
-   estimate  p.value statistic  n gp  Method +
-1 0.5646726 0.113182  1.810198 10  1 pearson +
-> spcor.test(datavar$bankaccount, datavar$famnum, datavar$income) +
-    estimate   p.value statistic  n gp  Method +
-1 -0.4086619 0.2748117 -1.184655 10  1 pearson+
 >  > 
-> +> p.b.f <- pcor.test(datavar$bankaccount, datavar$famnum, datavar$income) 
-> cor(datavar$bankaccount, datavar$income)+> p.b.f 
 +  estimate p.value statistic  n gp  Method 
 +1  -0.6729 0.04702    -2.406 10  1 pearson 
 +> p.b.f$estimate 
 +[1] -0.6729 
 +>  
 +> sp.b.i <- spcor.test(datavar$bankaccount, datavar$income, datavar$famnum) 
 +> sp.b.i 
 +  estimate p.value statistic  n gp  Method 
 +1   0.5647  0.1132      1.81 10  1 pearson 
 +> sp.b.i$estimate 
 +[1] 0.5647 
 +> sp.b.f <- spcor.test(datavar$bankaccount, datavar$famnum, datavar$income) 
 +> sp.b.f 
 +  estimate p.value statistic  n gp  Method 
 +1  -0.4087  0.2748    -1.185 10  1 pearson 
 +> sp.b.f$estimate 
 +[1] -0.4087 
 +>  
 +>  
 +> zc.b.i <- cor(datavar$bankaccount, datavar$income) 
 +> zc.b.i 
 [1] 0.7944 [1] 0.7944
-> cor(datavar$bankaccount, datavar$famnum)+> zc.b.f <- cor(datavar$bankaccount, datavar$famnum) 
 +> zc.b.f
 [1] -0.6923 [1] -0.6923
->   +>  
- +> zc.b.i^2 - (sp.b.i$estimate)^2 
 +[1] 0.3123 
 +> zc.b.f^2 - (sp.b.f$estimate)^2 
 +[1] 0.3123 
 +>  
 +>  
 +> 
 </code> </code>
  
 ====== e.g. 3. College enrollment in New Mexico University ====== ====== e.g. 3. College enrollment in New Mexico University ======
 <code> <code>
-> datavar <- read.csv("http://commres.net/wiki/_media/r/dataset_hlr.csv")+> datavar <- read.csv("http://commres.net/_media/r/dataset_hlr.csv")
 > str(datavar) > str(datavar)
 'data.frame': 29 obs. of  5 variables: 'data.frame': 29 obs. of  5 variables:
Line 427: Line 463:
 # Import data (simulated data for this example) # Import data (simulated data for this example)
 # myData <- read.csv('http://static.lib.virginia.edu/statlab/materials/data/hierarchicalRegressionData.csv')   # myData <- read.csv('http://static.lib.virginia.edu/statlab/materials/data/hierarchicalRegressionData.csv')  
-myData <- read.csv("http://commres.net/wiki/_media/hierarchical.regression.data.csv")+myData <- read.csv("http://commres.net/_media/hierarchical.regression.data.csv")
 </code> </code>
  
sequential_regression.1718148122.txt.gz · Last modified: by hkimscil
Donate Powered by PHP Valid HTML5 Valid CSS Driven by DokuWiki