[R] Determination of variables for optimizing one response using the desirability function

2018-12-25 Thread Thanh Tran
Dear all, I'm trying to use response surface methodology (rsm package) to . In my data, the response is KIC, and 4 factors are AC, AV, T, and Temp. A typical second-degree response modeling is as follows: > data<-read.csv("2.csv", header =T) > library(rsm) > f.quad <- rsm(KIC~SO(AC, AV, T,

Re: [R] Retrievable results in a procedure

2018-12-25 Thread Steven Yen
Thanks Sarah. Below, replacing "structure" with "invisible" does wonders--that serves my need. What I want is quite simple - I call a procedure and it does two things: (1) display results for all; (2) save retrievable results for use in further analysis, e.g., in knitr. Earlier, with "structure

Re: [R] fields package question

2018-12-25 Thread M P
Actually, let's set it grid_new.l <- list(abcissa=c(-15.0,-14.5),ordinate=y) to avoid out of bounds On Tue, Dec 25, 2018 at 4:41 PM M P wrote: > Thanks, Eric, for looking into that. > The values are below and since I subset the new abcissa is smaller range > grid_new.l <- list(abcissa=c(-15.0,-

Re: [R] fields package question

2018-12-25 Thread M P
Thanks, Eric, for looking into that. The values are below and since I subset the new abcissa is smaller range grid_new.l <- list(abcissa=c(-15.0,-14.),ordinate=y) I am emailing form gmail - don't know why is using html to format when all is in ascii x [1] -15.20180 -15.01948 -14.86533 -14.73180

Re: [R] Retrievable results in a procedure

2018-12-25 Thread Jeff Newmiller
You can use `capture.output`, but a far, far better solution is to remove the output statements from your computation functions entirely and let the caller decide whether to print the results. You can, for example, add a `debug` parameter to the function, and if true it can return a list of as

Re: [R] Problem with Kruskal–Wallis test

2018-12-25 Thread Bert Gunter
"So, I'm not an expert in R and statistics" So you need to seek local help from someone who is. Statistics is usually off-topic for this list -- it is about R programming primarily. And online is probably not a good venue for the sort of discussion you need anyway. Cheers, Bert Bert Gunter

Re: [R] Retrievable results in a procedure

2018-12-25 Thread Sarah Goslee
I'm a bit confused about what you actually want, but I think invisible() might be the answer. Note that there's already a base function try() so that's not a great name for test functions. Sarah On Tue, Dec 25, 2018 at 8:47 AM Steven Yen wrote: > I would like to suppressed printing of retrieva

[R] Retrievable results in a procedure

2018-12-25 Thread Steven Yen
I would like to suppressed printing of retrievable results in a procedure and to print only when retrieved. In line 10 below I call procedure "try" and get matrices A,B,C all printed upon a call to the procedure. I get around this unwanted printing by calling with v<-try(A,B) as in line 11. An

Re: [R] Problem with Kruskal–Wallis test

2018-12-25 Thread Giuseppe Cillis
Dear Michael, Thanks for your answer. So, I'm not an expert in R and statistics, how can I create this interval of confidence of groups? Thanks Gc Il giorno sab 22 dic 2018, 13:34 Michael Dewey ha scritto: > Dear Giuseppe > > If I understand you correctly you have a very large sample size so it

[R] Fwd: Problem with Kruskal–Wallis test

2018-12-25 Thread Giuseppe Cillis
Dear all, I am a beginner with R (and also with the statistics) for which I hope to be clear. I should do this non-parametric test on data I extracted from maps. In practice I have a column that represents the landscape Dynamics of a certain time period (there are 3 dynamics, each of them marked by