similar to: Tukey HSD (or other post hoc tests) following repeated measures ANOVA

Dear R experts, I would like to get an R^2 - like value for a multilevel regression using lme. I followed an archived suggestion by José Pinheiro to use the squared correlation between fitted and observed values, i.e., (cor(fitted(model.lme), getResponse(model.lme))^2 but getResponse returns the error message Error in val[, level] : incorrect number of dimensions The same happens with

Dear R contributors, Considering the following sample C code, that illustrates two possible uses of a Cholesky decomp for inverting a matrix, equally valid at least in theory: SEXP test() { int d = 2; int info = 0; double mat[4] = {2.5, 0.4, 0.4, 1.7}; double id[4] = {1.0, 0.0, 0.0, 1.0}; double lmat[3]; F77_CALL(dpotrf)("L", &d, mat, &d, &info); lmat[0] = mat[0]; lmat[1]

Full_Name: Elizabeth Purdom Version: 1.9.1 OS: Windows XP Submission from: (NULL) (171.64.102.199) It would be great if heatmap and/or image had the option of printing a legend bar on it somewhere that would indicate the ranges given by the colors in the heatmap (i.e. a very small image rectangle with the same colors as in the heatmap). Because heatmap is using layout, it seems pretty

Hello everyone, I'm trying to test the accurracy of R on the Eckerle4 dataset from NIST and I don't understand how the control option of the nls function works. I tought nls(...) was equivalent to nls(...control=nls.control()) i.e nls.control() was the default value of control, but here is the error I get : > n2=nls(V1~(b1/b2) *

> > Message: 6 > Date: Mon, 30 Sep 2013 20:23:16 +0200 > From: David Demelier <demelier.david at gmail.com> > To: freebsd-stable at freebsd.org > Subject: Re: 9.2-PRE: switch off that stupid "Nakatomi Socrates" > Message-ID: <5249C194.8030504 at gmail.com> > Content-Type: text/plain; charset=ISO-8859-1 > > On 30.09.2013 15:50, Matthieu Volat

Dear R-experts, my goal is to visualize the following polynomial regression as a 3D-surface: Z = b0 + b1*X + b2*Y + b3*XY + b4*X^2 + b5*Y^2 I believe that a solution to this problem may be of interest to a wider range of scientists because the problem is a derivative of a more general problem, i.e.: how to describe the relationship between one dependent variable and the DIFFERENCE between two

Hello again, lately I have been playing a lot with heatmap functions, in particular I complemented heatmap.2 with a nice feature from heatmap.plus, i.e. the possibility of assigning a matrix to ColSideColors so to take into account more annotations. Below is the code I introduced in heatmap.2, everything else is maintained equal to heatmap.2: if (!missing(ColSideColors)) { #if

Does anybody know how to deffine in the test tukeyHSD the pairs of comparison that you want to get? It is throw the lmat commnad, but I don't know the correct usage... Thanks a lot! Jos? -- MSc Jos? Alberto F. Monteiro Botanisches Institut Universit?t Basel ?????? ?????

Hy there, I have a problem with Samba 2.2.5 (on AIX 4.3) concerning a printer share. I set up a few printer share which work well. Now I want to restrict the access to one of them. In order to achieve this I put in a hosts deny=all hosts allow=<IP-Adresses of hosts, which should be allowed to print> in the corresponding printer-share-section. What happens is this: non

This isn't the help list for LAPACK, but as far as I can tell, dppsv expects a symmetric matrix input compacted as triangular, not a Choleski decomposed one. So try assigning lmat before the call to dpotrf. -pd > On 20 Dec 2014, at 22:06 , Pierrick Bruneau <pbruneau at gmail.com> wrote: > > Dear R contributors, > > Considering the following sample C code, that

Dear *, I would like to obtain for each factor of my anova model the "response variable vs factor" plot with means and 95% Tukey HSD intervals. I would appreciate any information on how to do that. Cheers -------------------------------------------------------------------- Max MANFRIN Tel.: +32 (0)2 650 3168 IRIDIA - CoDE, CP 194/6

Hi together, with our first implementation of Samba running under Linux on s/390 we are running in problems because some of our users authentificate to the windows domain with their full name, i.e. "maenz ullrich". Logon to Samba services failed because user (and password) verification seems to stop at that whitespace without ignoring it - by the way, it's impossible to create

Dear all, I would be glad if someone could tell me if Tukey HSD (function of stats library) accept the NA values and could correct the unbalanced design (different number of sample in each group) Thanks for help Alice Le Bars

I've finally took the time to figure the last few bugs (that I know of). This patch will be submit to be included in a few weeks. This patch should be generic enough for portable without too much hassle. This patch mimics OpenBSD's ftp behavior. I'm not sure like that (e.g. it doesn't put / at the end of directories by default), but that is more a question for the community

Hi, I need to determinate HSD value from a matrix like that Thesis Days A1 Cx 0 66.07 Cx 0 60.24 Cx 0 42.86 Tw 0 66.07 Tw 0 60.24 Tw 0 42.86 Aa 0 66.07 Aa 0 60.24 Aa 0 42.86 Qe 0 66.07 Qe 0 60.24 Qe 0 42.86 Cx 56 310.86 Cx 56 223.17 Cx 56 186.77 Tw 56 149.42 Tw 56 127.75 Tw 56 138.59 Aa 56 130.24 Aa 56 214.83 Aa 56 137.95 Qe 56 186.64 Qe 56 189.09 Qe 56 187.87 this is my script

Dear all, I am doing a one way between subjects anova in an unbalanced data set. Suppose we have "a" levels of the one factor. I want to merge the not so significantly different levels into the same cluster. Can I do a Tukey Kramer HSD and then use the following algorithm: For i in 2 : "a" For j in 1 : i-1 if mean of level i is not significantly different

Dear all Is there any way of doing a Tukey HSD post-hoc test after a Kruskal- Wallis or Friedman rank sum test (in the ctest package)? Thanks in advance, Albertus Dr. Albertus J. Smit Department of Botany University of Cape Town Private Bag Rondebosch 7700 South Africa Tel. +27 21 689 3032

Dear list members, I'm running some tests looking at differences between means for various levels of a factor, using Tukey's HSD method. I would like to plot the data as boxplots or dotplots, with horizontal significance lines indicating which groups are statistically significantly different, according to Tukey HSD. Here's a nice image showing an example of such a graphical

Hi, I was wondering if you could help me: The following are the first few lines of my data set: subject group condition depvar s1 c ver 114.87 s1 c feet 114.87 s1 c body 114.87 s2 c ver 73.54 s2 c feet 64.32 s2 c body 61.39 s3 a ver 114.87 s3 a feet 97.21 s3 a body 103.31 etc. I entered the following ANOVA command: >dat <- read.table("mydata.txt", header=T)

dear all effect of a, b and c on d, total 48 comparisons, got one anova result in model1=aov(d~A*B*C). can we get all the result in one command? or can we interpret the whole comparisons from this result? how it work in Tukey HSD? jose ____________________________________________________________________________________ Access over 1 million songs.