similar to: error message from gls(), unstructured in lower triangle, identical diagonal, help.

Hi, I got a computational problem with lme (nlme library R 1.9.1) using corSymm(). Here is the data: [,1] [,2] [,3] [,4] [,5] [,6] [1,] 0.19639793 0.09127954 0.11733288 0.07598273 0.06545106 0.06211532 [2,] 0.22773467 0.10981912 0.16052847 0.38101187 0.18353474 0.24072918 [3,] 0.46743388 0.45733836 0.32191178 0.43356107 0.39159746 0.53984221 [4,]

Hi, lme() users, Can some one tell me how to do this. I model Orthodont with the same G for random variables, but different R{i}'s for boys and girls, so that I can get sigma1_square_hat for boys and sigma2_square_hat for girls. The model is Y{i}=X{i}beta + Z{i}b + e{i} b ~ iid N(0,G) and e{i} ~ iid N(0,R{i}) i=1,2 orth.lme <- lme(distance ~ Sex * age, data=Orthodont, random=~age|Subject,

Hi, R or S+ users, I want to make a simple transformation for the model, but for the whole design matrix. The model is distance ~ age * Sex, where Sex is a factor. So the design matrix may look like the following: (Intercept) age SexFemale age:SexFemale 1 1 8 0 0 2 1 10 0 0 3 1 12 0 0 4

? stato filtrato un testo allegato il cui set di caratteri non era indicato... Nome: non disponibile Url: https://stat.ethz.ch/pipermail/r-help/attachments/20070905/6d1002c1/attachment.pl

Dear R community, I have the following problem I hoped you could help me with. My data is save in thousand of files with a weird extension containing for numbers and a z. For example *.1405z. With list.files I managed to load this data into R. It looks like this (the row numbers are not in the original file): 35 :LATEST STAGE 3.60 FT AT 730 AM CST ON 0102 36

Two new packages have been added to CRAN. The first is a simple package which provides the standard functions for a triangle distribution (rtriangle, dtriangle, ptriangle, qtriangle). The second package generates and augments Latin Hypercube Samples. Functions for generating random and optimized designs are provided. I look forward hearing your feedback and improving the documentation!

Two new packages have been added to CRAN. The first is a simple package which provides the standard functions for a triangle distribution (rtriangle, dtriangle, ptriangle, qtriangle). The second package generates and augments Latin Hypercube Samples. Functions for generating random and optimized designs are provided. I look forward hearing your feedback and improving the documentation!

Dear all, mydata<-data.frame(x1=c(1,4,6),x2=c(3,1,2),x3=c(2,1,3)) cor(mydata) x1 x2 x3 x1 1.0000000 -0.5960396 0.3973597 x2 -0.5960396 1.0000000 0.5000000 x3 0.3973597 0.5000000 1.0000000 I wonder if it is possible to fill only lower triangle of this correlation matrix? Using 'dist' doesn't seem to be useful as it doesnt allow to convert this table

Hello, Looking to generate Pascal's triangle in R. How do I get started? Thanks, -- Matilda Gogos matildaelizabethv@gmail.com [[alternative HTML version deleted]]

so, is it safe to put canreinvite=yes on a 7960? on a 1750? on a spa-x000? an xten? how the heck do i find out other than the hard way? randy -- ps: pun intended

Hi, I want to compute the lower triangle of a square matrix (optionally, sans diagonal). With for() loops I can do something like this: ## 5 by 5 matrix rtn for (j in 1:5) { for (k in 1:j) { if (j != k) { ## optional rtn[j, k] <- my.func(j, k) } } } I'd like to do this with apply(). Is there some way I can do this kind of 'short-circuit'? Thanks, Mark Wilkinson

Hi.. I a have a symmetric matrix to plot . I would like to plot only the Upper triangle but with the diagonal as the Base of the rectangle. Is there an easy way to do it. Thanks. Harsh --------------------------------- [[alternative HTML version deleted]]

Sorry if this answer was already given, or if I miss the point, but did you have a look to "lowerTriangle" and "upperTriangle" functions in the "gdata" package ? # example # A<-matrix(rnorm(9),3,3) # B<-B<-matrix(NA,dim(A)[1],dim(A)[2]) # lowerTriangle(B)<-lowerTriangle(A) # upperTriangle(B)<-lowerTriangle(A) # diag(B)<-diag(A) Hope this helps,

Hi there, I'd like to plot a triangle with each point having a certain color defined and fill it with the interpolated values. The following code shall represent somehow, three points with x,y, and the "amount" of red for example... point1 <- c(1.1, 1.7, 255) point2 <- c(2.2, 1.5, 180) point3 <- c(1.8, 2.2, 60) can anybody give me a hint which methods to use to realize

Hello, I am using triangle.class in ade4 package and a factor. I would like to use this graphic parameter (pch) to set a different kind of point for each group in my factor (lakes). I have tried pch=1:4, I guess it does what it is supposed to do but not what I want. I can't use colors. Thank you in advance, -- View this message in context:

Hi, I am trying to use function pairs to plot the scatterplot, but I only want to keep the upper triangle, what's the argument to do this? Thanks, -Jack [[alternative HTML version deleted]]

Hello R users... how to align this Zmat (triangle image) in X axis? I would like that the triangle's base become in the X axis and the triangle's height become in the Y axis Is there some trick for make this? Thanks. Cleber ######################## my test and try f <- function(x,y){ z=1-x-y z[ z < (-1e-15) ] <- NA return( -100*x+0*y+100*z ) } x = seq( 1, 0,

Hi, I have made some nice triangle plots showing the projected airborne, landborne and oceanborne fraction of anthropogenic emitted carbon dioxide at 10 year intervals for the next century. Currently all ten points on my plot are identical so I can't distinguish between them. I would like to be able to track the trend as the century progresses. How can I show the time progression on the

Hi, does anyone have any idea how I might make a plot in R similar to this: http://assets.sbnation.com/assets/272612/BATTER-1B-PENA.png Specifically looking to plot 3 numbers in a triangle like that. [[alternative HTML version deleted]]

hello, I am trying to develop a triangle plot but am having difficultly assigning the row.names to the 3 columns in the data.frame Here is what I've done, attach(SoilVegHydro) dim(SoilVegHydro) 1292 39 # now take 3 variables from main data.frame for plotting dat <- cbind.data.frame(TP, meanAnnualDepthAve, BulkDensity) # These are variables held in the data frame SoilVegHydro