r - t.test for all combinations of a grouping factor first by subsetting -

-

i trying t.test combinations of grouping factors, first choosing , subsetting data based on selection criteria column.

the structure of data: str(mydata)

'data.frame':   240 obs. of  6 variables:  $ group   : chr  "g1" "g1" "g1" "g1" ...  $ category: chr  "cat1" "cat1" "cat1" "cat1" ...  $ subgroup: chr  "sg1" "sg1" "sg1" "sg1" ...  $ score   : num  0.156 0.131 0.092 0.319 0.179 ...  $ sd      : num  0.0768 0.0768 0.0768 0.0768 0.0768 ...  $ se      : num  0.0172 0.0172 0.0172 0.0172 0.0172 ...

i have 3 groups in the group column: g1, g2 , g3
have 4 categories in category column: cat1, cat2, cat3 , cat4
have twelve subgroups: sg1, sg2, until sg12

currently, creating list of combinations subgroup names first subsetting data based of group id, in example g1 , g3:

combinations <- combn(unique(mydata[mydata$group %in% c("g1", "g3"),]$subgroup),2, simplify = false)

then doing t.test each of these combinations, referring answer here:

results <- lapply(seq_along(combinations), function (n) {   mydatatemp <- mydata[rownames(mydata$subgroup) %in% unlist(combinations[n]),]   result <- t.test(mydatatemp[,1], mydatatemp[,2], alternative="two.sided", var.equal=true)   return(result)}) results

the error following:

error in t.test.default(mydatatemp[, 1], mydatatemp[, 2], alternative = "two.sided",  :    not enough 'x' observations in addition: warning message: in mean.default(x) : argument not numeric or logical: returning na

is there more efficient way this? otherwise, how correct error?

update

actually problem how call values in score column in t.test formula?

your code should this:

results <- lapply(seq_along(combinations), function (n) {   mydatatemp <- mydata[with(mydata, subgroup %in% unlist(combinations2[n]) & group %in% c("g1", "g3")),]   result <- t.test(mydatatemp[                    mydatatemp$subgroup == sapply(combinations[n], "[",1),4],                    mydatatemp[mydatatemp$subgroup == sapply(combinations[n], "[", 2),4],                    alternative="two.sided", var.equal=true)              return(result)}) results

Comments

Post a Comment

Popular posts from this blog

ios - RestKit 0.20 — CoreData: error: Failed to call designated initializer on NSManagedObject class (again) -

-
this error has shown on many times, seems. prominent answer issue seems restkit 0.20 — coredata: error: failed call designated initializer on nsmanagedobject class however, solution posted doesn't work me. let's start beginning: i want upload file server. want parse response make sure uploaded. if not, try again later. model has boolean variable called nsnumber *senttoserver purpose (in coredata, bool saved nsnumber ) . because don't know whether sending going successful, save entity coredata model first before sending via restkit. when post request successful, server returns json string this: {"id":14,"created":true} notice that, of course, id not set on coredata model yet, both id , created new informations model. here's model: @interface recording : nsmanagedobject @property (nonatomic, retain) nsstring * audiofilepath; @property (nonatomic, retain) nsdate * createdat; @property (nonatomic, retain) nsstring * deviceid; @...
Read more

java - Digest auth with Spring Security using javaconfig -

-
so i'm trying create digest authentication spring following documentation trying translate xml "requirements" in java requirements. let's have xml in docs: <bean id="digestfilter" class= "org.springframework.security.web.authentication.www.digestauthenticationfilter"> <property name="userdetailsservice" ref="jdbcdaoimpl"/> <property name="authenticationentrypoint" ref="digestentrypoint"/> <property name="usercache" ref="usercache"/> </bean> <bean id="digestentrypoint" class= "org.springframework.security.web.authentication.www.digestauthenticationentrypoint"> <property name="realmname" value="contacts realm via digest authentication"/> <property name="key" value="acegi"/> <property name="noncevalidityseconds" value="10"/> </b...
Read more

laravel - PDOException in Connector.php line 55: SQLSTATE[HY000] [1045] Access denied for user 'root'@'localhost' (using password: YES) -

-
my .env: app_env=local app_debug=true app_key=myappkey app_url=http://localhost db_connection=mysql db_host=localhost db_port=3306 db_database=gestiontaller db_username=root db_password='root' cache_driver=file session_driver=file queue_driver=sync redis_host=127.0.0.1 redis_password=null redis_port=6379 mail_driver=smtp mail_host=mailtrap.io mail_port=2525 mail_username=null mail_password=null mail_encryption=null i've tried password ", ' , without, restarting server every time. i've tried changing localhost 172.0.0.1. yes, password correct, running mysql -u root -p introducing password root in prompt can go mysql. users table: +------------------+-----------+ | user | host | +------------------+-----------+ | debian-sys-maint | localhost | | mysql.sys | localhost | | root | localhost | | root1 | localhost | +------------------+-----------+ i've added homestead.yaml this: database_port: 3306 ...
Read more