Re: [R] Quantile Regression without intercept

2015-10-06 Thread Lorenz, David
: > > > On Oct 6, 2015, at 7:58 AM, Lorenz, David wrote: > > > > Did you verify that the correct percentages were above/below the > regression > > lines? I did a quick check and for example did not consistently get 50% > of > > the observed response values gre

Re: [R] Quantile Regression without intercept

2015-10-06 Thread Lorenz, David
Did you verify that the correct percentages were above/below the regression lines? I did a quick check and for example did not consistently get 50% of the observed response values greater than the tau=.5 line. I did when I included the nonzero intercept term. > Date: Mon, 5 Oct 2015 21:14:04 +05

Re: [R] wilcox.test - difference between p-values of R and online calculators

2014-09-04 Thread Lorenz, David
I think that the issue, at least with the online calculator that I looked at, is that it does not adjust the standard deviation of the test statistic for ties, so the standard deviation is larger and hence larger p-value. I was able to reproduce the reported z-score using the equation for the st

Re: [R] Simulative data production

2014-04-30 Thread Lorenz, David
Merve, I'm not 100 percent sure I understand everything that you want. But start with the simulated likert scale data. The code that you have is not very efficient and it has at least one typo. I do not know if columns or rows represent the persons, so I'll set up as NROW and NCOL. An efficient

Re: [R] predefined area under the curve

2014-01-17 Thread Lorenz, David
Elisa, Part of the issue is that there is no unique solution. You could increase each of the 6 values by a fixed amount to make up the 0.242003 difference in area, or you could increase them proportionally, among many other schemes. I'll outline an approach and you can decide how you want to pr

Re: [R] R-help Digest, Vol 128, Issue 30

2013-10-28 Thread Lorenz, David
Pavlos, There are several ways to evaluate how well new data fit an old regression.Part of the answer depends on what you are concerned about. For example, if you are concerned about bias, you can test whether the mean of the new data is within the expected range of the mean of that many new valu

Re: [R] Aggregate rows with same fields, within factors

2013-09-17 Thread Lorenz, David
Andrea, The argument na.action controls how missings are treated. is this what you wanted? aggregate(cbind(IND, DUM)~.,data=net1,sum, na.action=na.pass) Dave >Date: Mon, 16 Sep 2013 11:42:07 -0400 >From: Andrea Goijman >To: arun >Cc: R help >Subject: Re: [R] Aggregate rows with same fields,