search for: lmat

...lors)) { #if (!is.matrix(ColSideColors)) #stop("'ColSideColors' must be a matrix") if (!is.character(ColSideColors) || dim(ColSideColors)[1] != nc) stop("'ColSideColors' dim()[2] must be of length ncol(x)") lmat <- rbind(lmat[1, ] + 1, c(NA, 1), lmat[2, ] + 1) lhei <- c(lhei[1], 0.2, lhei[2]) } if (!missing(RowSideColors)) { #if (!is.matrix(RowSideColors)) #stop("'RowSideColors' must be a matrix") if (!is.character(RowSideColors) || dim(Row...

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] = mat[1]; lmat[2] = mat[3]; F77_CALL(dppsv)("L", &d, &d, lmat, id, &d, &info); // id now contains L^(-T) F77_CALL(dpotri)("L", &d, mat, &d, &info); // mat c...

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 illustrates two possible > uses of a Cholesky decomp for inverting a matrix, equally va...

...feasible approach, but that's how I understood the MMC example) and a contrast vector (contrasting conditions 2 and 3): I have to admit that I don't quite understand what I'm doing here (not that you couldn't tell) > RT2.aov <- aov(terms(RT~Subj*Cond, WMU3C)) > Cond.lmat <- c(0,1,-1) > Tukey <- glht.mmc(RT2.aov, focus = "Cond", focus.lmat = Cond.lmat) yielding Error in mvt(lower = carg$lower, upper = carg$upper, df = df, corr = carg$corr, : NA/NaN/Inf in foreign function call (arg 6) In addition: Warning message: In cov2cor(covm) : diagon...

...end in the plot somehow. Again, I find it difficult to do without making your own myheatmap command and modifying the layout. And this seems easy to do spacewise, because there's a "0" in the layout matrix in the upper left corner over the row dendogram: e.g. if(!missing(legend.text)){lmat<-lmat+1} Clearly the plot could be really crowded if you choose all the options, but at least you'd have the options, and I think they're the devil to do yourself. And if you could choose the relative widths of the different parts in the layout (the widths and heights options in the lay...

Hi all, Does anyone know of a fairly easy way to "stretch" a heatmap vertically? I've got 42 arrays and would like to be able to see as many significant genes as possible (right now I can only get 50 genes with it still being readable). In some comparisons there are several hundred significant genes. I've fiddled with the "asp" argument, but that doesn't

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 ?????? ?????

Hello, I want to create a figure that consists of a collection of 16 graphs on 4 rows. I am using nf <- layout(matrix(seq(1,16), 4,4, byrow=TRUE), respect=TRUE) boxplot(... to create the layout of my 16 graphs. It works really well. However, I'd like to add sub-titles that would apply to each row of 4 graphs. More specifically, I'd like to have something like: subtitle 1

dear list, i am using a multifactorial design with two treatments (factor A: drugs, three levels; factor B: theraphy, two levels) and a time factor (three levels, different timepoint). hypothetically, i measured the same subjects for all treatements and timepoints, so its a repeated measurement design. now i ran an anova in R and also some Tukey post-hoc tests using glht. but what i am actually

...fun" "hcc" "hclustfun" [21] "hcr" "hline" "iy" "key" [25] "keysize" "labCol" "labRow" "lhei" [29] "linecol" "lmat" "lwid" "main" [33] "margins" "max.breaks" "max.raw" "max.scale" [37] "min.breaks" "min.raw" "min.scale" "mmat" [41] "na.color" &quot...

.... L and M > # > # Arguments: > # AMAT ... p by q matrix > # BMAT ... order p symmetric positive definite matrix > # CMAT ... order q symmetric positive definite matrix > # Returns: > # VALUES ... vector of length s = min(p,q) of eigenvalues > # LMAT ... p by s matrix L > # MMAT ... q by s matrix M > > It's not clear to me how it is used and exactly what it is doing and how that compares with Lapack. > > Berend >

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

...returns all buffers. winmm: Mark headers WHDR_DONE at the last possible time. Louis Lenders (6): dpnet: Fill the DPN_SP_CAPS structure in IDirectPlay8Peer_GetSPCaps with some sane values + tests. gdiplus: Return zero properties in stubbed GdipGetPropertyCount. include: Add lmat.h. netapi32: Add stub for NetScheduleJobEnum. wine.inf: Add key for InstallDate under CurrentVersionNT. dwmapi: Add stub for DwmRegisterThumbnail. Luca Bennati (1): po: Update Italian translation. Maarten Lankhorst (2): quartz: Fix support for mpeg layer 1 and 2....