[R] Question Mixed-Design Anova in R

Hi Everyone, I have a question about Mixed-Design Anova in R. I want to obtain Mauchly�s test of Sphericity and the Greenhouse-Geisser correction. I have managed to do it in SPSS: GLM Measure1 Measure2 Measure3 Measure4 Measure5 Measure6 BY Grouping /WSFACTOR=Measure 6 Polynomial

Re: [R] Question about function scope

e becomes something like this: > >x <- y <- z <- 0 > ># here is my scope problem >result <- list(x = x, y = y, z = z) # c == 0 case >if (c==1) result <- bar1() >if (c==2) result <- bar2() >x <- result[["x"]] >y <-

Re: [R] Question about function scope

"y"]] z <- result[["z"]] Duncan Murdoch - Original Message - From: "Duncan Murdoch" To: "Sebastien Bihorel" , r-help@r-project.org Sent: Tuesday, October 30, 2018 4:13:05 PM Subject: Re: [R] Question about function scope On 30/10/2018 3:5

Re: [R] Question about function scope

Thanks a lot Eric, I think you are on the same page as Duncan (at least with his 2nd option). I will definitively explore this. From: "Eric Berger" To: "Duncan Murdoch" Cc: "Sebastien Bihorel" , "R mailing list" Sent: Tuesday, October 30, 20

Re: [R] Question about function scope

That's cool! I think this solution would fit better with what my intended setup. Thanks a lot - Original Message - From: "Duncan Murdoch" To: "Sebastien Bihorel" , r-help@r-project.org Sent: Tuesday, October 30, 2018 4:18:51 PM Subject: Re: [R] Question about fu

Re: [R] Question about function scope

. - Original Message - From: "Duncan Murdoch" To: "Sebastien Bihorel" , r-help@r-project.org Sent: Tuesday, October 30, 2018 4:13:05 PM Subject: Re: [R] Question about function scope On 30/10/2018 3:56 PM, Sebastien Bihorel wrote: > Hi, > > From the R user manual,

Re: [R] Question about function scope

Here's another modification to your code that also works. It's a lot uglier, but will allow bar1 and bar2 to be used in multiple functions, not just foo. bar1 <- function(env){ env$x <- 1 env$y <- 1 env$z <- 1 with(env, cat(sprintf('bar1: x=%d, y=%d, z=%d\n', x, y, z))) } bar2 <-

Re: [R] Question about function scope

Hi Sebastien, I like Duncan's response. An alternative approach is to pass around environments, as in the following: bar1 <- function(e) { e$x <- e$y <- e$z <- 1 cat(sprintf('bar1: x=%d, y=%d, z=%d\n', e$x, e$y, e$z)) } bar2 <- function(e) { e$x <- e$y <- e$z <- 2 cat(sprintf('bar2: x=%d,

Re: [R] Question about function scope

On 30/10/2018 3:56 PM, Sebastien Bihorel wrote: Hi, From the R user manual, I have a basic understanding of the scope of function evaluation but have a harder time understanding how to mess with environments. My problem can be summarized by the code shown at the bottom: - the foo function

[R] Question about function scope

Hi, >From the R user manual, I have a basic understanding of the scope of function >evaluation but have a harder time understanding how to mess with environments. My problem can be summarized by the code shown at the bottom: - the foo function performs some steps including the assignment of

[R] Question for help

How do I interpret p and q in the corARMA correlation structure? I understand that p is the autoregessive order and q is the moving average, but I don't know how to interpret what does it means when it is (2,3) or (2,2) for example. Can anyone give me a simple explanation? [[alternative

Re: [R] Question on Binom.Confint

t the code.) Göran -Original Message- From: Jim Lemon [mailto:drjimle...@gmail.com] Sent: Thursday, September 13, 2018 11:50 PM To: Guo, Fang (Associate) ; r-help mailing list Subject: Re: [R] Question on Binom.Confint Hi Fang, Let's assume that you are using the "binom.confint" f

Re: [R] Question on Binom.Confint

I did use library(binom). However, I was able to use the method "lrt" which is short for likelihood ratio test. -Original Message- From: Jim Lemon [mailto:drjimle...@gmail.com] Sent: Thursday, September 13, 2018 11:50 PM To: Guo, Fang (Associate) ; r-help mailing list Subje

Re: [R] Question on Binom.Confint

I used library(binom). -Original Message- From: Bert Gunter [mailto:bgunter.4...@gmail.com] Sent: Thursday, September 13, 2018 10:04 PM To: Guo, Fang (Associate) Cc: r-help-requ...@r-project.org; R-help Subject: Re: [R] Question on Binom.Confint In what package? Binomial confidence

Re: [R] Question on Binom.Confint

-- > From: Bert Gunter [mailto:bgunter.4...@gmail.com] > Sent: Thursday, September 13, 2018 10:04 PM > To: Guo, Fang (Associate) > Cc: r-help-requ...@r-project.org; R-help > Subject: Re: [R] Question on Binom.Confint > > In what package? > Binomial confidence inte

Re: [R] Question on Binom.Confint

Hi Fang, Let's assume that you are using the "binom.confint" function in the "binom" package and you have made a spelling mistake or two. This function employs nine methods for estimating the binomial confidence interval. Sadly, none of these is "lrt". The zero condition is discussed in the help

Re: [R] Question on Binom.Confint

On 2018-09-13 20:58, David Winsemius wrote: On Sep 13, 2018, at 1:15 PM, Guo, Fang (Associate) wrote: Hi, I have a question with the function Binom.Confint(x,n,"method"=lrt). For likelihood ratio test, I'd like to ask how you define the upper limit when the frequency of successes is

Re: [R] Question on Binom.Confint

In what package? Binomial confidence interval functions are in several. Bert Gunter "The trouble with having an open mind is that people keep coming along and sticking things into it." -- Opus (aka Berkeley Breathed in his "Bloom County" comic strip ) On Thu, Sep 13, 2018 at 6:38 PM Guo, Fang

Re: [R] Question on Binom.Confint

> On Sep 13, 2018, at 1:15 PM, Guo, Fang (Associate) > wrote: > > Hi, > > I have a question with the function Binom.Confint(x,n,"method"=lrt). For > likelihood ratio test, I'd like to ask how you define the upper limit when > the frequency of successes is zero. Thanks! First you need to

[R] Question on Binom.Confint

Hi, I have a question with the function Binom.Confint(x,n,"method"=lrt). For likelihood ratio test, I'd like to ask how you define the upper limit when the frequency of successes is zero. Thanks! Fang Guo Associate CORNERSTONE RESEARCH 699 Boylston Street, 5th Floor Boston, MA 02116-2836

Re: [R] Question

Your question is notable for what it is missing... any trace of R code. [1][2][3] Do read the Posting Guide. I don't see "Sargan" in base R, so your analysis likely used a contributed package... there seem to be a couple, so your example code would clarify. I don't see the number of IVs listed

[R] Question

Hi, My name is Luke and I come from Poland. I have one question, maybe very simple, but I can not resolve it. In dynamic panel data (GMM estimator) after running the model, I recieve a AR test and Sargan test, but the "number of instruments" are not displayed. In Stata and Gretl this

Re: [R] Question "breakdown point" --> robust statistics

> Olivier Crouzet > on Thu, 19 Apr 2018 18:13:30 +0200 writes: > Hi, > I think he's talking about how much a statistical estimator is influenced by extreme datapoints, e.g. >

Re: [R] Question

Hi, I think he's talking about how much a statistical estimator is influenced by extreme datapoints, e.g. https://en.m.wikipedia.org/wiki/Robust_statistics#Breakdown_point Olivier -- Olivier Crouzet Assistant Professor @LLING UMR6310 - Université de Nantes / CNRS Guest Scientist @UMCG -

Re: [R] Question

On 15/04/2018 17:26, Marc Girondot via R-help wrote: Le 15/04/2018 à 17:56, alireza daneshvar a écrit : break-down point Can you explain more what you plan to do and give an example of what you have tried to do until now to do a "break down point" in R. Perhaps a "break down point" is

Re: [R] Question

Le 15/04/2018 à 17:56, alireza daneshvar a écrit : break-down point Can you explain more what you plan to do and give an example of what you have tried to do until now to do a "break down point" in R. Perhaps a "break down point" is common in your field, but I have no idea about what it is

[R] Question

Hi, How to calculate break-down point in R? Thanks [[alternative HTML version deleted]] __ R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide

Re: [R] Question about subset

.org Sent: Monday, April 9, 2018 12:33:41 AM Subject: Re: [R] Question about subset Sent from my iPhone > On Apr 8, 2018, at 9:06 PM, Sebastien Bihorel > <sebastien.biho...@cognigencorp.com> wrote: > > Hi, > > The help page for subset states "subset: logical expr

Re: [R] Question about subset

Thanks. That works great! > df <- data.frame(x=c(1,1,NA,NA,2), y=c('a','a','a','b','b'), > z=c(TRUE,FALSE,TRUE,FALSE,TRUE)) > cond1 <- 'x==1' > cond2 <- 'x==1 & z' > df x y z 1 1 a TRUE 2 1 a FALSE 3 NA a TRUE 4 NA b FALSE 5 2 b TRUE > subset(df, subset =

Re: [R] Question about subset

Sent from my iPhone > On Apr 8, 2018, at 9:06 PM, Sebastien Bihorel > wrote: > > Hi, > > The help page for subset states "subset: logical expression indicating > elements or rows to keep: missing values are taken as false." > > Before I try to re-invent

Re: [R] Question about subset

Please read the "Details" section of the help file -- it essentially tells you what to do by directly using indexing; e.g. for a vector, index by: [is.na(subset) | subset] . Adjust as necessary for your data structure.Or look at the code of, e.g. subset.data.frame and create your own subset

[R] Question about subset

Hi, The help page for subset states "subset: logical expression indicating elements or rows to keep: missing values are taken as false." Before I try to re-invent the wheel, I would like to know if one of the base or recommended packages would contain a variant of the subset function that

Re: [R] question regarding the AICcmodavg package

Do not model average the regression coefficients. This makes no sense (see Cade 2015. Model averaging and muddled multimodel inferences). And, no there is no reasonable way to model average regression coefficients ala Burnham and Anderson approach when some of your models include interactions

[R] question regarding the AICcmodavg package

Dear moderator, If possible I would like to send in the following question for R-help: I am analyzing a small data set using PGLS with phylogenetic uncertainty taken into account and thereby including 100 potential phylogenetic tree scenarios. I've managed to run models on all of the different

Re: [R] question with integrate function

Oh ok. Thanks very much. I will have to restrict to a shorter interval. Hanna 2018-02-06 14:33 GMT-05:00 Göran Broström : > Hi Hanna, > > your function is essentially zero outside a short interval around 9. And > the help page states: "If the function is approximately

Re: [R] question with integrate function

Hi Hanna, your function is essentially zero outside a short interval around 9. And the help page states: "If the function is approximately constant (in particular, zero) over nearly all its range it is possible that the result and error estimate may be seriously wrong." You could try to

Re: [R] question with integrate function

Sorry. I meant in the previous email that the function h() is a monotone decreasing function. Thanks very much. 2018-02-06 13:32 GMT-05:00 li li : > Hi all, > The function h below is a function of c and it should be a monotone > increasing function since the integrand is

[R] question with integrate function

Hi all, The function h below is a function of c and it should be a monotone increasing function since the integrand is nonnegative and integral is taken from c to infinity. However, as we can see from the plot, it is not shown to be monotone. Something wrong with the usage of integrate function?

Re: [R] Question - TukeyHSD

I may be missing the point, but if you can do the calculations by hand could you not write a function in R to do the same and apply it to a data.frame of the values? On Saturday, November 4, 2017, 11:01:59 AM EDT, Paul Kent wrote: Hi, I've been

Re: [R] Question - TukeyHSD

Assuming you have as your model a one-way ANOVA, you can use the aovSufficient function in HH. ## install.packages("HH") ## if you don't have it yet library(HH) ?aovSufficient On Sat, Nov 4, 2017 at 9:30 AM, Paul Kent wrote: > Hi, > > > I've been performing some

Re: [R] Question - TukeyHSD

AFAICS, without knowing the structure of your data and what specifically you wish to do, how could one answer your question? -- which is probably yes, you can do it, but without further info, ??? Maybe someone with a better crystal ball can help -- or you could clarify. Cheers, Bert Bert

[R] Question - TukeyHSD

Hi, I've been performing some TukeyHSD tests in R and have come across papers that include all the necessary information (means, sample size and SE) for me to perform the piecewise comparison by hand. I can reach out to the authors to get the original dataset, but it raises the question, can

Re: [R] Question on Simultaneous Equations & Forecasting

gt; > Berend > >> Best, >> Bernhard >> >> -Ursprüngliche Nachricht- >> Von: Berend Hasselman [mailto:b...@xs4all.nl] >> Gesendet: Donnerstag, 13. Juli 2017 10:53 >> An: OseiBonsu, Frances >> Cc: Pfaff, Bernhard Dr.; r-help@r-project.org >

Re: [R] Question on Simultaneous Equations & Forecasting

g identities and/or >> technical equations - together in a format suitable for applying the >> Gauss-Seidel method. Hence, forget about 2SLS or 3SLS and Haavelmo-bias. >> >> Best wishes, >> Bernhard >> >> -Ursprüngliche Nachricht- >> Vo

Re: [R] Question on Simultaneous Equations & Forecasting

type. Best, Bernhard -Ursprüngliche Nachricht- Von: Berend Hasselman [mailto:b...@xs4all.nl] Gesendet: Donnerstag, 13. Juli 2017 10:53 An: OseiBonsu, Frances Cc: Pfaff, Bernhard Dr.; r-help@r-project.org Betreff: [EXT] Re: [R] Question on Simultaneous Equations & Forecasting Frances

Re: [R] Question on Simultaneous Equations & Forecasting

ether in a format suitable for applying the Gauss-Seidel > method. Hence, forget about 2SLS or 3SLS and Haavelmo-bias. > > Best wishes, > Bernhard > > -Ursprüngliche Nachricht- > Von: R-help [mailto:r-help-boun...@r-project.org] Im Auftrag von OseiBonsu, >

Re: [R] Question on Simultaneous Equations & Forecasting

2SLS or 3SLS and Haavelmo-bias. Best wishes, Bernhard -Ursprüngliche Nachricht- Von: R-help [mailto:r-help-boun...@r-project.org] Im Auftrag von OseiBonsu, Frances Gesendet: Mittwoch, 12. Juli 2017 22:36 An: r-help@r-project.org Betreff: [EXT] [R] Question on Simultaneous Equations

[R] Question on Simultaneous Equations & Forecasting

Hello, I have estimated a simultaneous equation model (similar to Klein's model) in R using the system.fit package. I have an identity equation, along with three other equations. Do you know how to explicitly identify the identity equation in R? I am also trying to forecast the dependent

Re: [R] Question about correlation

Please keep the conversation on the list: Others may be able to help you better than I can. On 2017-07-06 10:38, SEB140004 Student wrote: Ya. I had successfully got the result. Thank you very much. :) Isn't possible for me to obtain the network from the correlation matrix? I know nothing

Re: [R] Question about correlation

I think that Goran is right, I didn't take "cor(data)" literally. Jim On Thu, Jul 6, 2017 at 3:38 AM, Göran Broström wrote: > On 2017-07-05 11:56, Jim Lemon wrote: >> >> Hi Chin Yi, >> If you are trying to correlate "Health" with "Disease", i.e. >> >>

Re: [R] Question about correlation

On 2017-07-05 11:56, Jim Lemon wrote: Hi Chin Yi, If you are trying to correlate "Health" with "Disease", i.e. cydf<-read.table(text="OTU ID Health Disease Bacterial 1 0.29 0.34 Bacterial 2 0.25 0.07 Bacterial 3 0.06 0.06 Bacterial 4 0.07 0.09 Bacterial 5 0.02 0.05",

Re: [R] Question about correlation

Hi Chin Yi, If you are trying to correlate "Health" with "Disease", i.e. cydf<-read.table(text="OTU ID Health Disease Bacterial 1 0.29 0.34 Bacterial 2 0.25 0.07 Bacterial 3 0.06 0.06 Bacterial 4 0.07 0.09 Bacterial 5 0.02 0.05", header=TRUE) print(cor(cydf$Health,cydf$Disease)) [1]

[R] Question about correlation

Greeting. Dear Mr/Mrs/Miss, OTU ID Health Disease Bacterial 1 0.29 0.34 Bacterial 2 0.25 0.07 Bacterial 3 0.06 0.06 Bacterial 4 0.07 0.09 Bacterial 5 0.02 0.05 Above show the first 6 data sets, may I ask that the reason of R show the error like "Error in cor(data) : 'x' must be numeric" ? And

Re: [R] Question

This is not a Spark-help mailing list, either. -- Sent from my phone. Please excuse my brevity. On June 22, 2017 4:20:36 PM PDT, Amrith Deepak wrote: >This function won’t work with objects in spark as you can’t do a dfda$a >in spark as it’s not stored as a local variable. >

Re: [R] Question

Sorry, I don’t think you understand my problem. dta = spark_read_csv(sc,”data_tbl”,”extension/file",delimiter = "|") dta is just a pointer to spark where the data is stored. I am using sparklyr to run this on spark. I’m not running it locally so I can’t use the $. I know how to do this locally

Re: [R] Question

This function won’t work with objects in spark as you can’t do a dfda$a in spark as it’s not stored as a local variable. Thanks, Amrith > On Jun 22, 2017, at 4:15 PM, David Winsemius wrote: > > >> On Jun 22, 2017, at 11:22 AM, Amrith Deepak wrote:

Re: [R] Question

Rows are horizontal, columns are vertical. You really need to spend some time with an R tutorial. dta <- read.table( "yourfile", header=TRUE, as.is=TRUE ) dta2 <- dta dta2$D <- c( "awe", "abcd", "asdf", "xyz" ) dta2 <- dta2[ , c( "A", "D" ) ] -- Sent from my phone. Please excuse my brevity.

Re: [R] Question

> On Jun 22, 2017, at 11:22 AM, Amrith Deepak wrote: > > Hi, > > I am using Spark and the Sparklyr library in R. > > I have a file with several lines. For example > > A B C > awer.ttp.netCode554 > abcd.ttp.netCode747 > asdf.ttp.net

[R] Question

Hi, I am using Spark and the Sparklyr library in R. I have a file with several lines. For example A B C awer.ttp.netCode554 abcd.ttp.netCode747 asdf.ttp.netPart554 xyz.ttp.net Part747 I want to split just column A of the table and I want

Re: [R] An R question

Thank you! 2017-06-07 10:38 GMT-04:00 Ivan Calandra : > Hi, > > Check the FAQ 7.31 > https://cran.rstudio.com/doc/FAQ/R-FAQ.html#Why-doesn_0027t- > R-think-these-numbers-are-equal_003f > > And read the posting guide too... > https://www.r-project.org/posting-guide.html > > HTH,

Re: [R] An R question

Hi, Check the FAQ 7.31 https://cran.rstudio.com/doc/FAQ/R-FAQ.html#Why-doesn_0027t-R-think-these-numbers-are-equal_003f And read the posting guide too... https://www.r-project.org/posting-guide.html HTH, Ivan -- Dr. Ivan Calandra TraCEr, Laboratory for Traceology and Controlled Experiments

Re: [R] An R question

Hello, See FAQ 7.31 And try ?all.equal. Hope this helps, Rui Barradas Em 07-06-2017 15:32, li li escreveu: Hi all, In checking my R codes, I encountered the following problem. Is there a way to fix this? I tried to specify options(digits=). I did not fix the problem. Thanks so much

[R] An R question

Hi all, In checking my R codes, I encountered the following problem. Is there a way to fix this? I tried to specify options(digits=). I did not fix the problem. Thanks so much for your help! Hanna > cdf(pmass)[2,2]==pcum[2,2][1] FALSE> cdf(pmass)[2,2][1] 0.758> > pcum[2,2][1]

[R] Question on interpreting glmer() results

Hello, I originally posted this on the stats stack exchange site, but given its focus on R software, it was removed -- so I figured I'd post here. I'm having trouble interpreting a change in effect direction and significance when I add an interaction term to my glmer() model. *Part 1* I ran an

Re: [R] Question on function "scatterplot3d"

> On 1 Jun 2017, at 03:41, li li wrote: > > Hi all, > I have a question with regard to making plots using function > "scatterplot3d". > Please see the example below. It looks like, for y axis, the tickmark text > was cutoff. > The number "10" does not show up completely.

[R] Question on function "scatterplot3d"

Hi all, I have a question with regard to making plots using function "scatterplot3d". Please see the example below. It looks like, for y axis, the tickmark text was cutoff. The number "10" does not show up completely. I tried to work with par(mpg). It does not seem to work. Hope to get some

Re: [R] Question about change the length of a string.

Hi Jeff, Thank you a lot! Best, Yen -Original Message- From: Jeff Newmiller [mailto:jdnew...@dcn.davis.ca.us] Sent: Thursday, May 18, 2017 11:11 PM To: r-help@r-project.org; Yen Lee <b88207...@ntu.edu.tw> Subject: Re: [R] Question about change the length of a string.

Re: [R] Question about change the length of a string.

http://stackoverflow.com/questions/2261079/how-to-trim-leading-and-trailing-whitespace-in-r -- Sent from my phone. Please excuse my brevity. On May 18, 2017 8:57:51 PM PDT, Yen Lee wrote: >Hello everyone, > > > >I have a question and I need your precious kind help. > >

[R] Question about change the length of a string.

Hello everyone, I have a question and I need your precious kind help. I am working on matching two string. However, the length of the two strings is different. For example, one is "example" (nchar=7), the other one is "example " (nchar=10). The R considers them as different strings but

Re: [R] question: mediation results are not in line with compression of glmm consisted models

This is not a statistical help site, and your questions appear to be about statistics, not programming in R. I would suggest that you get local statistical help, but you might try posting on a stats.stackexchange.com for remote help. -- Bert Bert Gunter "The trouble with having an open mind is

[R] question: mediation results are not in line with compression of glmm consisted models

hi all, I'll begin with my two question and all the related information (description of the research and the data and full output) will follow. 1. When i execute model1 (glmm with random intercept only for subjects): predictor (suppBin) and outcome (DtlsBinUp) and pre-intervention variables, it

Re: [R] Question on accessing foreign files

I've attached data.restore4.txt, containing the function data.restore4(), which has the same argument list as foreign::data.restore() and is mean to be called by the latter if the first line of the file is "## Dump S Version 4 Dump". It can read version 4 of the 'S data dump' format, which for

[R] Question on accessing foreign files

Hi all, I have several data files provided in mtw format (Minitab) and sdd format (S-Plus) and I need to read them in R. I do not have access either to Minitab or to S-Plus. How can I accomplish this task ? Thank you, Daniel [[alternative HTML version deleted]]

Re: [R] question about the anova() function for deviance analysis

On 16/04/2017 9:46 AM, Sophie Dubois wrote: Dear Maintener, I have recently had a bad experience with the anova() function. Indeed, I wanted to process a deviance analysis between 2 mixed linear models and I was really surprise to see that depending on the ordre in which I gave my models, the

Re: [R] question about the anova() function for deviance analysis

You are sending your email to a whole mailing list of volunteers, not a specific "maintainer" (and I am not one). However, your assertions convey unfamiliarity with statistics rather than deficiencies in R, and this mailing list is not a stats tutoring list. I did a quick Google search and

Re: [R] question about the anova() function for deviance analysis

This list is about R programming; your question seems mostly about statistics, and is therefore off topic here. I suggest you consult a local statistical expert who *is* comfortable with such statistical analyses. In general, partitions of sums of squares in statistical models can depend on the

[R] question about the anova() function for deviance analysis

Dear Maintener, > I have recently had a bad experience with the anova() function. > Indeed, I wanted to process a deviance analysis between 2 mixed linear > models and I was really surprise to see that depending on the ordre in > which I gave my models, the function did not the same thing: once it

Re: [R] r question

Hello, Please keep this on the list. I'm cc-ing r-help and so should you always. I got no error, just warning messages. After a long output I got the following. NOTE: HARD MAXIMUM GENERATION LIMIT HIT Solution Fitness Value: -2.592094e+03 Parameters at the Solution (parameter, gradient):

Re: [R] r question

Hello, There's a paenthesis missing in > relativerisk<- matrix(log(c(1,1,2,2),ncol=4,byrow = TRUE) + beta_true<-relativerisk Error: unexpected symbol in: "relativerisk<- matrix(log(c(1,1,2,2),ncol=4,byrow = TRUE) beta_true" The correct instruction would be relativerisk<-

[R] r question

Hi ,I have some question about simulate, I don't know how to paste question to this website,so I paste below. I use genoud to find the maximum likelihood value, but when I use numcut=3 ,it will get error message,like this " coxph.wtest(fit$var[nabeta, nabeta], temp, control$toler.chol) :

[R] Question about logspline

I have a dataset with values of limited precision. I am trying to plot its density using package logspline. The logspline() call with default parameters causes oscillations in the density plot. My solution was to fuzzify the input values, using function: > fuzz <- function(x, prec) x +

Re: [R] Question about Cubist Model

> On Jan 12, 2017, at 5:37 PM, Lorenzo Isella wrote: > > Dear All, > I am fine tuning a Cubist model (see > https://cran.r-project.org/web/packages/Cubist/index.html). > I am a bit puzzled by its output. On a dataset which contains 275 > cases, I get non mutually

[R] Question about Cubist Model

Dear All, I am fine tuning a Cubist model (see https://cran.r-project.org/web/packages/Cubist/index.html). I am a bit puzzled by its output. On a dataset which contains 275 cases, I get non mutually exclusive rules. E.g., in the output below, rules 2 and 3 cover all the 275 cases of the data set

[R] Question about Cubist Model

Dear All, I am fine tuning a Cubist model (see https://cran.r-project.org/web/packages/Cubist/index.html). I am a bit puzzled by its output. On a dataset which contains 275 cases, I get non mutually exclusive rules. E.g., in the output below, rules 2 and 3 cover all the 275 cases of the data set

Re: [R] Question regarding package deSolve and R versions

My apologies. When I found the r-help email, I entirely missed the posting guide. I will try just update my version of R to be current. If I need to ask another qustion, I will be sure to look up the posting guide, first. Thank you very much!! Shirley On Dec 14, 2016 12:28 AM, "Jeff Newmiller"

Re: [R] Question regarding package deSolve and R versions

You are approaching this the wrong way... you should be upgrading your version of R to be current. The deSolve package is up to date... it is your R that is behind. You should also read the Posting Guide which would have told you this and also that you need to post using plain text and

Re: [R] Question regarding package deSolve and R versions

> On Dec 13, 2016, at 3:08 PM, Shirley R wrote: > > To whom it may concern, > > I've run into trouble installing the crqa package. R error messages tell me > it is because I don't have the deSolve package. I cannot install the > deSolve package, but according to

[R] Question regarding package deSolve and R versions

To whom it may concern, I've run into trouble installing the crqa package. R error messages tell me it is because I don't have the deSolve package. I cannot install the deSolve package, but according to everything I've found online, my version of R (3.2.1) should be fine to work with deSolve. Is

Re: [R] Question about proxy setting of R

.@googlemail.com> >To: jim holtman <jholt...@gmail.com> >Cc: qwertyui_per...@yahoo.co.jp; "r-help@r-project.org" <r-help@r-project.org> >Date: 2016/12/5, Mon 22:18 >Subject: Re: [R] Question about proxy setting of R > > >As far as I can see, the

Re: [R] Question about proxy setting of R

;http", not "https". What should I do ? J J - Original Message - >From: jim holtman <jholt...@gmail.com> >To: qwertyui_per...@yahoo.co.jp >Date: 2016/12/2, Fri 09:13 >Subject: Re: [R] Question about proxy setting of R > > >Try

Re: [R] Question about proxy setting of R

t "https". > > What should I do ? > > > > J J > > > > > > - Original Message - > > *From:* jim holtman <jholt...@gmail.com> > > *To:* qwertyui_per...@yahoo.co.jp > > *Date:* 2016/12/2, Fri 09:13 > > *Subject:* Re: [

Re: [R] Question about proxy setting of R

-- --- > > I assume the proxy server is only available for "http", not "https". > What should I do ? > > J J > > > - Original Message - > *From:* jim holtman <jholt...@gmail.com> > *To:* qwertyui_per...@ya

[R] Question about proxy setting of R

Hello, I use R 3.0.2 on Win 7 through proxy server using ".Rprofile" in home directory that includes "Sys.setenv(http_proxy=proxy_server:port)". There has been no problem to access the internet for some years. In this situation, I installed R 3.3.1 and then entered "update.packages ()", however,

[R] Question regarding Naive Bayes

Hello,   I am working with the naïve bayes function inlibrary(e1071).   The function calls are: transactions.train.nb = naiveBayes(as.factor(DealerID) ~    as.factor(Manufacturer)     + as.factor(RangeDesc)    

[R] Question on 'predict' in 'markovchain' package

Dear Sir/Madam, I want to report a problem of 'predict' function in the 'markovchain' package and I will use to examples to explain the problem. Problem: -It only follows the path with transition probability greater or equal than 0.5. Here are two examples. Example (1) I created an MC with

[R] Question about proxy setting of R

Hello, I use R 3.0.2 on Win 7 through proxy server using ".Rprofile" in home directory that includes "Sys.setenv(http_proxy=proxy_server:port)". There has been no problem to access the internet for some years. In this situation, I installed R 3.3.1 and then entered "update.packages ()", however,

Re: [R] Question about expression parser for "return" statement

Just to add on a bit, please note that the return is superfluous. If you write this: normalDensityFunction = function(x, Mean, Variance) { # no "return" value given at all (1/sqrt(2*pi*Variance))*exp(-(1/2)*((x - Mean)^2)/Variance) } normalDensityFunction(2,0,1) ...you get the

Re: [R] Question about expression parser for "return" statement

Sorry, I missed the operation-after-function call aspect of the OP question. However, I think my policy of avoiding the return function as much as possible serves as an effective antibugging strategy for this problem, in addition to its other benefits. -- Sent from my phone. Please excuse my

Re: [R] Question about expression parser for "return" statement

I stand corrected. I have been chided in the past for not explicitly returning my output by someone claiming it is not best practices. -Steve On Mon, 2016-11-14 at 12:22 -0500, Duncan Murdoch wrote: On 14/11/2016 11:26 AM, Wolf, Steven wrote: Just to add on a bit, please note that the

Re: [R] Question about expression parser for "return" statement

On 14/11/2016 11:26 AM, Wolf, Steven wrote: Just to add on a bit, please note that the return is superfluous. If you write this: normalDensityFunction = function(x, Mean, Variance) { # no "return" value given at all (1/sqrt(2*pi*Variance))*exp(-(1/2)*((x - Mean)^2)/Variance) }

<    1   2   3   4   5   6   7   8   9   10   >