Posted on

How to subset a dataset for analysis in R (without creating a new dataset)

I got this advice from someone when  I needed to know how to apply a procedure on a subgroup of subjects within  the analysis dataset.  Thanks.

library(dplyr) library(magrittr)

ols_result <- data %>% dplyr::filter(year=1) %>% lm(y~x,.) summary(ols_result)

dplyr::filter(year==1)のところでデータセットを絞っています。 dply::を追加しているのは、filterがたまに他のパッケージに存在する同名の関数と競合するためです。

 

Posted on

Using R to run multilevel models

I'm learning how to run multilevel models in R.

I tried the analysis of variance model, AKA, the intercept-only model.

fit<lme(post_test~1,random=~1|school,data=mySASData,control=list(opt="optim"))
summary(fit)
anova(fit)
VarCorr(fit)
summary(fit)

 

I run this in SAS and get the same results.  I didn't get the same degree of freedom.

proc glimmix data=sashlm.core_2014_4_years;
class school;
model post_test=/solution ddfm=kr dist=normal link=identity;
random intercept /subject=school;
run;

Posted on

How to set a working directory in R

In R, you can set a working directory in this way (I use C:/temp as an example).

setwd("C:/temp")

If this generates an error message. suspect that you used a wrong quotation mark (this can happen when you copy and paste an example from the Internet).  Just retype using '.  You can read more about it here.

You also need to be careful about the slash.  Microsoft Windows uses backward slashes, "\", to indicate the folder structure.  R uses forward slash, "/" at least on my machine.  This may be OS-dependent, so to really determine which slash is used on your environment, try this command to know the currently working directory and see which is used:

getwd()

I got this on my personal PC ("/" is used).

"C:/Users/kaz/Documents"

Alternatively, go to File and click on Change dir... to specify the working directly.  (I don't see this option in R-Studio), so this must be specific to the regular version of R).