I created custom function according to my requirement which is given below:
*selection.index<- function(ID, phen_mat, gen_mat, weight_mat, GAY){ ID = toString(ID) p<- as.matrix(phen_mat) g<- as.matrix(gen_mat) w<- as.matrix(weight_mat) bmat<- solve(phen_mat) %*% gen_mat %*% weight_mat GA<- 2.063 * t(bmat) %*% g %*% w / (t(bmat) %*% p %*% bmat)^0.5 PRE<- if(missing(GAY)){ (GA/GA) * 100 } else { (GA/GAY) * 100 } result<- list("ID" = ID, "b" = matrix(round(bmat,4), nrow = 1), "GA" = round(GA,4), "PRE" = round(PRE,4)) return(data.frame(result))}* *sc<- list()* *sc[[1]]<- selection.index(ID = 12, pmat[1:2,1:2], gmat[1:2,1:2], wmat[1:2,1])* *sc[[2]]<- selection.index(ID = 13, pmat[c(1,3),c(1,3)], gmat[c(1,3),c(1,3)], wmat[c(1,3),1])* for more detail about question follow stack overflow link: https://stackoverflow.com/questions/66734928/how-to-solve-this-through-loop *"Healthy soil, Healthy life."* *"A war based on Satyagraha is always of two kinds. One is the war we wage against injustice, and the other we fight our won weaknesses."* - *Sardar Patel* *"You have to dream before your dreams can come true."* - *A. P. J.* *Abdul Kalam* *"Think before you print and save a tree."* *ZANKRUT GOYANI* *B.Sc. (Hons.) Agriculture* [[alternative HTML version deleted]] ______________________________________________ [hidden email] mailing list -- To UNSUBSCRIBE and more, see https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code. |
Hi Goyani,
In its present form, the function stalls because you haven't defined pmat before trying to pass it to the function. gmat and wmat suffered the same fate. Even if I define these matrices as I think you have, "solve" fails because at least one is singular. First, put the function in order as below. I think this is what you sent made readable. selection.index<- function(ID, phen_mat, gen_mat, weight_mat, GAY){ ID<-toString(ID) p<-as.matrix(phen_mat) g<-as.matrix(gen_mat) w<-as.matrix(weight_mat) bmat<- solve(phen_mat) %*% gen_mat %*% weight_mat GA<- 2.063 * t(bmat) %*% g %*% w / (t(bmat) %*% p %*% bmat)^0.5 if(missing(GAY)) PRE<-(GA/GA) * 100 else PRE<-(GA/GAY) * 100 result<-list(ID=ID,b=matrix(round(bmat,4),nrow=1), GA=round(GA,4),PRE=round(PRE,4)) return(data.frame(result)) } Next, what sort of matrices do you want to pass? Then an answer may emerge. Jim On Mon, Mar 22, 2021 at 6:03 AM Goyani Zankrut <[hidden email]> wrote: > > I created custom function according to my requirement which is given below: > > > > > > > > > > > > > > > *selection.index<- function(ID, phen_mat, gen_mat, weight_mat, GAY){ ID = > toString(ID) p<- as.matrix(phen_mat) g<- as.matrix(gen_mat) w<- > as.matrix(weight_mat) bmat<- solve(phen_mat) %*% gen_mat %*% weight_mat > GA<- 2.063 * t(bmat) %*% g %*% w / (t(bmat) %*% p %*% bmat)^0.5 PRE<- > if(missing(GAY)){ (GA/GA) * 100 } else { (GA/GAY) * 100 } result<- > list("ID" = ID, "b" = matrix(round(bmat,4), nrow = 1), "GA" = round(GA,4), > "PRE" = round(PRE,4)) return(data.frame(result))}* > > *sc<- list()* > *sc[[1]]<- selection.index(ID = 12, pmat[1:2,1:2], gmat[1:2,1:2], > wmat[1:2,1])* > *sc[[2]]<- selection.index(ID = 13, pmat[c(1,3),c(1,3)], > gmat[c(1,3),c(1,3)], wmat[c(1,3),1])* > > for more detail about question follow stack overflow link: > https://stackoverflow.com/questions/66734928/how-to-solve-this-through-loop > *"Healthy soil, Healthy life."* > *"A war based on Satyagraha is always of two kinds. One is the war we wage > against injustice, and the other we fight our won weaknesses."* - *Sardar > Patel* > *"You have to dream before your dreams can come true."* - *A. P. J.* *Abdul > Kalam* > *"Think before you print and save a tree."* > > *ZANKRUT GOYANI* > *B.Sc. (Hons.) Agriculture* > > [[alternative HTML version deleted]] > > ______________________________________________ > [hidden email] mailing list -- To UNSUBSCRIBE and more, see > https://stat.ethz.ch/mailman/listinfo/r-help > PLEASE do read the posting guide http://www.R-project.org/posting-guide.html > and provide commented, minimal, self-contained, reproducible code. ______________________________________________ [hidden email] mailing list -- To UNSUBSCRIBE and more, see https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code. |
Hi Goyani,
You are setting "PRE" to the return value of "if" which is one of TRUE (1), FALSE(0) or NULL. Because GAY is always missing in your example, "PRE" is always set to 1. Then you always want to pass 1 in the sample list, and that will not assign anything to PRE. By correcting the "if" clause and defining matrices that are unlikely to be singular, I can run a "for" loop as follows: selection.index<- function(ID, phen_mat, gen_mat, weight_mat, GAY){ p<-as.matrix(phen_mat) g<-as.matrix(gen_mat) w<-as.matrix(weight_mat) bmat<- solve(phen_mat) %*% gen_mat %*% weight_mat GA<- 2.063 * t(bmat) %*% g %*% w / (t(bmat) %*% p %*% bmat)^0.5 if(missing(GAY)) PRE<-(GA/GA) * 100 else PRE<-(GA/GAY) * 100 result<-list(ID=ID,b=matrix(round(bmat,4),nrow=1), GA=round(GA,4),PRE=round(PRE,4)) return(data.frame(result)) } pmat<-matrix(sample(1:16,16),4) gmat<-matrix(sample(17:32),16,4) wmat<-matrix(sample(1:4,4),4) mi<-combn(1:4,2) sc<-list() for(i in 1:ncol(matindx)) { as.numeric(ID<-paste0(mi[,i])) sc[[i]]<-selection.index(ID,pmat[mi[,i],mi[,i]],gmat[mi[,i],mi[,i]], wmat[mi[,i]],1) } This produces output for me. Good luck with whatever you are doing with this. Jim On Mon, Mar 22, 2021 at 2:51 PM Goyani Zankrut <[hidden email]> wrote: > > Greetings of the day, > Thank you for your response, Sir. > The full problem statement is given below: > > In our case, I'm taking 4 traits. > library(arrangements) > a<- combinations(4,2) # gives 6 pairwise combinations > class(a) # it's a "matrix" "array" > > now hypothetical data of three matrix for further calculation: > pmat<- matrix(1:16, nrow = 4) > gmat<- matrix(17:32, nrow = 4) > wmat<- matrix(1:4, nrow = 4) > > My custom function for further calculations: > selection.index<- function(ID, phen_mat, gen_mat, weight_mat, GAY){ > ID = toString(ID) > p<- as.matrix(phen_mat) > g<- as.matrix(gen_mat) > w<- as.matrix(weight_mat) > bmat<- solve(phen_mat) %*% gen_mat %*% weight_mat > GA<- 2.063 * t(bmat) %*% g %*% w / (t(bmat) %*% p %*% bmat)^0.5 > PRE<- if(missing(GAY)){ > (GA/GA) * 100 > } else { > (GA/GAY) * 100 > } > result<- list("ID" = ID, "b" = matrix(round(bmat,4), nrow = 1), "GA" = round(GA,4), "PRE" = round(PRE,4)) > return(data.frame(result)) > } > > Now I want to store this data into a list for further calculation: > sc<- list() > sc[[1]]<- selection.index(ID = 12, phen_mat = pmat[c(1,2),c(1,2)], gen_mat = gmat[c(1,2),c(1,2)], weight_mat = wmat[c(1,2),1]) > sc[[2]]<- selection.index(ID = 13, phen_mat = pmat[c(1,3),c(1,3)], gen_mat = gmat[c(1,3),c(1,3)], weight_mat = wmat[c(1,3),1]) > sc[[3]]<- selection.index(ID = 14, phen_mat = pmat[c(1,4),c(1,4)], gen_mat = gmat[c(1,4),c(1,4)], weight_mat = wmat[c(1,4),1]) > sc[[4]]<- selection.index(ID = 23, phen_mat = pmat[c(2,3),c(2,3)], gen_mat = gmat[c(2,3),c(2,3)], weight_mat = wmat[c(2,3),1]) > sc[[5]]<- selection.index(ID = 24, phen_mat = pmat[c(2,4),c(2,4)], gen_mat = gmat[c(2,4),c(2,4)], weight_mat = wmat[c(2,4),1]) > sc[[6]]<- selection.index(ID = 34, phen_mat = pmat[c(3,4),c(3,4)], gen_mat = gmat[c(3,4),c(3,4)], weight_mat = wmat[c(3,4),1]) > above list code is monotonous and time consuming for large data combination cycles like (7,2) = 21 combinations, (10,2) = 45 combinations. So I want to use the matrix a's each row as a vector in the selection.index function and result stores in a list. > > I hope now you will understand the full problem. I have checked the selection.index which has no issues and works well. > Thank you. > ______________________________________________ [hidden email] mailing list -- To UNSUBSCRIBE and more, see https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code. |
On 22/03/2021 1:59 a.m., Jim Lemon wrote:
> Hi Goyani, > You are setting "PRE" to the return value of "if" which is one of TRUE > (1), FALSE(0) or NULL. That's not true at all. The statement was PRE<- if(missing(GAY)){ (GA/GA) * 100 } else { (GA/GAY) * 100 } so the result is (GA/GA) * 100 or (GA/GAY)*100. > Because GAY is always missing in your example, If that's true and GA isn't missing, the result will always be 100. Duncan Murdoch > "PRE" is always set to 1. Then you always want to pass 1 in the sample > list, and that will not assign anything to PRE. By correcting the "if" > clause and defining matrices that are unlikely to be singular, I can > run a "for" loop as follows: > > selection.index<- function(ID, phen_mat, gen_mat, weight_mat, GAY){ > p<-as.matrix(phen_mat) > g<-as.matrix(gen_mat) > w<-as.matrix(weight_mat) > bmat<- solve(phen_mat) %*% gen_mat %*% weight_mat > GA<- 2.063 * t(bmat) %*% g %*% w / (t(bmat) %*% p %*% bmat)^0.5 > if(missing(GAY)) PRE<-(GA/GA) * 100 > else PRE<-(GA/GAY) * 100 > result<-list(ID=ID,b=matrix(round(bmat,4),nrow=1), > GA=round(GA,4),PRE=round(PRE,4)) > return(data.frame(result)) > } > > pmat<-matrix(sample(1:16,16),4) > gmat<-matrix(sample(17:32),16,4) > wmat<-matrix(sample(1:4,4),4) > > mi<-combn(1:4,2) > sc<-list() > for(i in 1:ncol(matindx)) { > as.numeric(ID<-paste0(mi[,i])) > sc[[i]]<-selection.index(ID,pmat[mi[,i],mi[,i]],gmat[mi[,i],mi[,i]], > wmat[mi[,i]],1) > } > > This produces output for me. Good luck with whatever you are doing with this. > > Jim > > > > > > On Mon, Mar 22, 2021 at 2:51 PM Goyani Zankrut <[hidden email]> wrote: >> >> Greetings of the day, >> Thank you for your response, Sir. >> The full problem statement is given below: >> >> In our case, I'm taking 4 traits. >> library(arrangements) >> a<- combinations(4,2) # gives 6 pairwise combinations >> class(a) # it's a "matrix" "array" >> >> now hypothetical data of three matrix for further calculation: >> pmat<- matrix(1:16, nrow = 4) >> gmat<- matrix(17:32, nrow = 4) >> wmat<- matrix(1:4, nrow = 4) >> >> My custom function for further calculations: >> selection.index<- function(ID, phen_mat, gen_mat, weight_mat, GAY){ >> ID = toString(ID) >> p<- as.matrix(phen_mat) >> g<- as.matrix(gen_mat) >> w<- as.matrix(weight_mat) >> bmat<- solve(phen_mat) %*% gen_mat %*% weight_mat >> GA<- 2.063 * t(bmat) %*% g %*% w / (t(bmat) %*% p %*% bmat)^0.5 >> PRE<- if(missing(GAY)){ >> (GA/GA) * 100 >> } else { >> (GA/GAY) * 100 >> } >> result<- list("ID" = ID, "b" = matrix(round(bmat,4), nrow = 1), "GA" = round(GA,4), "PRE" = round(PRE,4)) >> return(data.frame(result)) >> } >> >> Now I want to store this data into a list for further calculation: >> sc<- list() >> sc[[1]]<- selection.index(ID = 12, phen_mat = pmat[c(1,2),c(1,2)], gen_mat = gmat[c(1,2),c(1,2)], weight_mat = wmat[c(1,2),1]) >> sc[[2]]<- selection.index(ID = 13, phen_mat = pmat[c(1,3),c(1,3)], gen_mat = gmat[c(1,3),c(1,3)], weight_mat = wmat[c(1,3),1]) >> sc[[3]]<- selection.index(ID = 14, phen_mat = pmat[c(1,4),c(1,4)], gen_mat = gmat[c(1,4),c(1,4)], weight_mat = wmat[c(1,4),1]) >> sc[[4]]<- selection.index(ID = 23, phen_mat = pmat[c(2,3),c(2,3)], gen_mat = gmat[c(2,3),c(2,3)], weight_mat = wmat[c(2,3),1]) >> sc[[5]]<- selection.index(ID = 24, phen_mat = pmat[c(2,4),c(2,4)], gen_mat = gmat[c(2,4),c(2,4)], weight_mat = wmat[c(2,4),1]) >> sc[[6]]<- selection.index(ID = 34, phen_mat = pmat[c(3,4),c(3,4)], gen_mat = gmat[c(3,4),c(3,4)], weight_mat = wmat[c(3,4),1]) >> above list code is monotonous and time consuming for large data combination cycles like (7,2) = 21 combinations, (10,2) = 45 combinations. So I want to use the matrix a's each row as a vector in the selection.index function and result stores in a list. >> >> I hope now you will understand the full problem. I have checked the selection.index which has no issues and works well. >> Thank you. >> > > ______________________________________________ > [hidden email] mailing list -- To UNSUBSCRIBE and more, see > https://stat.ethz.ch/mailman/listinfo/r-help > PLEASE do read the posting guide http://www.R-project.org/posting-guide.html > and provide commented, minimal, self-contained, reproducible code. > ______________________________________________ [hidden email] mailing list -- To UNSUBSCRIBE and more, see https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code. |
If he's setting PRE to the return value of "if", that is the logical
value of the expression in the if statement as far as I know. I think that the expression within the else clause would be evaluated but not assigned to anything and since it is within the loop, would just be lost. PRE<-ifelse(missing(GAY),(GA/GA)*100,(GA/GAY)*100) would return either value depending upon whether GAY was missing. That's what I get from the help pages. Jim On Mon, Mar 22, 2021 at 8:34 PM Duncan Murdoch <[hidden email]> wrote: > > On 22/03/2021 1:59 a.m., Jim Lemon wrote: > > Hi Goyani, > > You are setting "PRE" to the return value of "if" which is one of TRUE > > (1), FALSE(0) or NULL. > > That's not true at all. The statement was > > PRE<- if(missing(GAY)){ > (GA/GA) * 100 > } else { > (GA/GAY) * 100 > } > > so the result is (GA/GA) * 100 or (GA/GAY)*100. > > > Because GAY is always missing in your example, > > If that's true and GA isn't missing, the result will always be 100. > > Duncan Murdoch > > > "PRE" is always set to 1. Then you always want to pass 1 in the sample > > list, and that will not assign anything to PRE. By correcting the "if" > > clause and defining matrices that are unlikely to be singular, I can > > run a "for" loop as follows: > > > > selection.index<- function(ID, phen_mat, gen_mat, weight_mat, GAY){ > > p<-as.matrix(phen_mat) > > g<-as.matrix(gen_mat) > > w<-as.matrix(weight_mat) > > bmat<- solve(phen_mat) %*% gen_mat %*% weight_mat > > GA<- 2.063 * t(bmat) %*% g %*% w / (t(bmat) %*% p %*% bmat)^0.5 > > if(missing(GAY)) PRE<-(GA/GA) * 100 > > else PRE<-(GA/GAY) * 100 > > result<-list(ID=ID,b=matrix(round(bmat,4),nrow=1), > > GA=round(GA,4),PRE=round(PRE,4)) > > return(data.frame(result)) > > } > > > > pmat<-matrix(sample(1:16,16),4) > > gmat<-matrix(sample(17:32),16,4) > > wmat<-matrix(sample(1:4,4),4) > > > > mi<-combn(1:4,2) > > sc<-list() > > for(i in 1:ncol(matindx)) { > > as.numeric(ID<-paste0(mi[,i])) > > sc[[i]]<-selection.index(ID,pmat[mi[,i],mi[,i]],gmat[mi[,i],mi[,i]], > > wmat[mi[,i]],1) > > } > > > > This produces output for me. Good luck with whatever you are doing with this. > > > > Jim > > > > > > > > > > > > On Mon, Mar 22, 2021 at 2:51 PM Goyani Zankrut <[hidden email]> wrote: > >> > >> Greetings of the day, > >> Thank you for your response, Sir. > >> The full problem statement is given below: > >> > >> In our case, I'm taking 4 traits. > >> library(arrangements) > >> a<- combinations(4,2) # gives 6 pairwise combinations > >> class(a) # it's a "matrix" "array" > >> > >> now hypothetical data of three matrix for further calculation: > >> pmat<- matrix(1:16, nrow = 4) > >> gmat<- matrix(17:32, nrow = 4) > >> wmat<- matrix(1:4, nrow = 4) > >> > >> My custom function for further calculations: > >> selection.index<- function(ID, phen_mat, gen_mat, weight_mat, GAY){ > >> ID = toString(ID) > >> p<- as.matrix(phen_mat) > >> g<- as.matrix(gen_mat) > >> w<- as.matrix(weight_mat) > >> bmat<- solve(phen_mat) %*% gen_mat %*% weight_mat > >> GA<- 2.063 * t(bmat) %*% g %*% w / (t(bmat) %*% p %*% bmat)^0.5 > >> PRE<- if(missing(GAY)){ > >> (GA/GA) * 100 > >> } else { > >> (GA/GAY) * 100 > >> } > >> result<- list("ID" = ID, "b" = matrix(round(bmat,4), nrow = 1), "GA" = round(GA,4), "PRE" = round(PRE,4)) > >> return(data.frame(result)) > >> } > >> > >> Now I want to store this data into a list for further calculation: > >> sc<- list() > >> sc[[1]]<- selection.index(ID = 12, phen_mat = pmat[c(1,2),c(1,2)], gen_mat = gmat[c(1,2),c(1,2)], weight_mat = wmat[c(1,2),1]) > >> sc[[2]]<- selection.index(ID = 13, phen_mat = pmat[c(1,3),c(1,3)], gen_mat = gmat[c(1,3),c(1,3)], weight_mat = wmat[c(1,3),1]) > >> sc[[3]]<- selection.index(ID = 14, phen_mat = pmat[c(1,4),c(1,4)], gen_mat = gmat[c(1,4),c(1,4)], weight_mat = wmat[c(1,4),1]) > >> sc[[4]]<- selection.index(ID = 23, phen_mat = pmat[c(2,3),c(2,3)], gen_mat = gmat[c(2,3),c(2,3)], weight_mat = wmat[c(2,3),1]) > >> sc[[5]]<- selection.index(ID = 24, phen_mat = pmat[c(2,4),c(2,4)], gen_mat = gmat[c(2,4),c(2,4)], weight_mat = wmat[c(2,4),1]) > >> sc[[6]]<- selection.index(ID = 34, phen_mat = pmat[c(3,4),c(3,4)], gen_mat = gmat[c(3,4),c(3,4)], weight_mat = wmat[c(3,4),1]) > >> above list code is monotonous and time consuming for large data combination cycles like (7,2) = 21 combinations, (10,2) = 45 combinations. So I want to use the matrix a's each row as a vector in the selection.index function and result stores in a list. > >> > >> I hope now you will understand the full problem. I have checked the selection.index which has no issues and works well. > >> Thank you. > >> > > > > ______________________________________________ > > [hidden email] mailing list -- To UNSUBSCRIBE and more, see > > https://stat.ethz.ch/mailman/listinfo/r-help > > PLEASE do read the posting guide http://www.R-project.org/posting-guide.html > > and provide commented, minimal, self-contained, reproducible code. > > > ______________________________________________ [hidden email] mailing list -- To UNSUBSCRIBE and more, see https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code. |
No, I am confounded, It does return the value of the expressions
within the respective braces, just like ifelse. Learn something every day. Jim On Mon, Mar 22, 2021 at 9:35 PM Jim Lemon <[hidden email]> wrote: > > If he's setting PRE to the return value of "if", that is the logical > value of the expression in the if statement as far as I know. I think > that the expression within the else clause would be evaluated but not > assigned to anything and since it is within the loop, would just be > lost. > > PRE<-ifelse(missing(GAY),(GA/GA)*100,(GA/GAY)*100) > > would return either value depending upon whether GAY was missing. > That's what I get from the help pages. > > Jim > > On Mon, Mar 22, 2021 at 8:34 PM Duncan Murdoch <[hidden email]> wrote: > > > > On 22/03/2021 1:59 a.m., Jim Lemon wrote: > > > Hi Goyani, > > > You are setting "PRE" to the return value of "if" which is one of TRUE > > > (1), FALSE(0) or NULL. > > > > That's not true at all. The statement was > > > > PRE<- if(missing(GAY)){ > > (GA/GA) * 100 > > } else { > > (GA/GAY) * 100 > > } > > > > so the result is (GA/GA) * 100 or (GA/GAY)*100. > > > > > Because GAY is always missing in your example, > > > > If that's true and GA isn't missing, the result will always be 100. > > > > Duncan Murdoch > > > > > "PRE" is always set to 1. Then you always want to pass 1 in the sample > > > list, and that will not assign anything to PRE. By correcting the "if" > > > clause and defining matrices that are unlikely to be singular, I can > > > run a "for" loop as follows: > > > > > > selection.index<- function(ID, phen_mat, gen_mat, weight_mat, GAY){ > > > p<-as.matrix(phen_mat) > > > g<-as.matrix(gen_mat) > > > w<-as.matrix(weight_mat) > > > bmat<- solve(phen_mat) %*% gen_mat %*% weight_mat > > > GA<- 2.063 * t(bmat) %*% g %*% w / (t(bmat) %*% p %*% bmat)^0.5 > > > if(missing(GAY)) PRE<-(GA/GA) * 100 > > > else PRE<-(GA/GAY) * 100 > > > result<-list(ID=ID,b=matrix(round(bmat,4),nrow=1), > > > GA=round(GA,4),PRE=round(PRE,4)) > > > return(data.frame(result)) > > > } > > > > > > pmat<-matrix(sample(1:16,16),4) > > > gmat<-matrix(sample(17:32),16,4) > > > wmat<-matrix(sample(1:4,4),4) > > > > > > mi<-combn(1:4,2) > > > sc<-list() > > > for(i in 1:ncol(matindx)) { > > > as.numeric(ID<-paste0(mi[,i])) > > > sc[[i]]<-selection.index(ID,pmat[mi[,i],mi[,i]],gmat[mi[,i],mi[,i]], > > > wmat[mi[,i]],1) > > > } > > > > > > This produces output for me. Good luck with whatever you are doing with this. > > > > > > Jim > > > > > > > > > > > > > > > > > > On Mon, Mar 22, 2021 at 2:51 PM Goyani Zankrut <[hidden email]> wrote: > > >> > > >> Greetings of the day, > > >> Thank you for your response, Sir. > > >> The full problem statement is given below: > > >> > > >> In our case, I'm taking 4 traits. > > >> library(arrangements) > > >> a<- combinations(4,2) # gives 6 pairwise combinations > > >> class(a) # it's a "matrix" "array" > > >> > > >> now hypothetical data of three matrix for further calculation: > > >> pmat<- matrix(1:16, nrow = 4) > > >> gmat<- matrix(17:32, nrow = 4) > > >> wmat<- matrix(1:4, nrow = 4) > > >> > > >> My custom function for further calculations: > > >> selection.index<- function(ID, phen_mat, gen_mat, weight_mat, GAY){ > > >> ID = toString(ID) > > >> p<- as.matrix(phen_mat) > > >> g<- as.matrix(gen_mat) > > >> w<- as.matrix(weight_mat) > > >> bmat<- solve(phen_mat) %*% gen_mat %*% weight_mat > > >> GA<- 2.063 * t(bmat) %*% g %*% w / (t(bmat) %*% p %*% bmat)^0.5 > > >> PRE<- if(missing(GAY)){ > > >> (GA/GA) * 100 > > >> } else { > > >> (GA/GAY) * 100 > > >> } > > >> result<- list("ID" = ID, "b" = matrix(round(bmat,4), nrow = 1), "GA" = round(GA,4), "PRE" = round(PRE,4)) > > >> return(data.frame(result)) > > >> } > > >> > > >> Now I want to store this data into a list for further calculation: > > >> sc<- list() > > >> sc[[1]]<- selection.index(ID = 12, phen_mat = pmat[c(1,2),c(1,2)], gen_mat = gmat[c(1,2),c(1,2)], weight_mat = wmat[c(1,2),1]) > > >> sc[[2]]<- selection.index(ID = 13, phen_mat = pmat[c(1,3),c(1,3)], gen_mat = gmat[c(1,3),c(1,3)], weight_mat = wmat[c(1,3),1]) > > >> sc[[3]]<- selection.index(ID = 14, phen_mat = pmat[c(1,4),c(1,4)], gen_mat = gmat[c(1,4),c(1,4)], weight_mat = wmat[c(1,4),1]) > > >> sc[[4]]<- selection.index(ID = 23, phen_mat = pmat[c(2,3),c(2,3)], gen_mat = gmat[c(2,3),c(2,3)], weight_mat = wmat[c(2,3),1]) > > >> sc[[5]]<- selection.index(ID = 24, phen_mat = pmat[c(2,4),c(2,4)], gen_mat = gmat[c(2,4),c(2,4)], weight_mat = wmat[c(2,4),1]) > > >> sc[[6]]<- selection.index(ID = 34, phen_mat = pmat[c(3,4),c(3,4)], gen_mat = gmat[c(3,4),c(3,4)], weight_mat = wmat[c(3,4),1]) > > >> above list code is monotonous and time consuming for large data combination cycles like (7,2) = 21 combinations, (10,2) = 45 combinations. So I want to use the matrix a's each row as a vector in the selection.index function and result stores in a list. > > >> > > >> I hope now you will understand the full problem. I have checked the selection.index which has no issues and works well. > > >> Thank you. > > >> > > > > > > ______________________________________________ > > > [hidden email] mailing list -- To UNSUBSCRIBE and more, see > > > https://stat.ethz.ch/mailman/listinfo/r-help > > > PLEASE do read the posting guide http://www.R-project.org/posting-guide.html > > > and provide commented, minimal, self-contained, reproducible code. > > > > > ______________________________________________ [hidden email] mailing list -- To UNSUBSCRIBE and more, see https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code. |
Free forum by Nabble | Edit this page |