similar to: question about matrix

Hello, I have a question about image command. When you draw a image plot, and let's say we use gray scale. then how can I generate a little bar next to the plot which shows the grey levels and the values which corresponds to the grey level? ( it doesn't matter whether I am using grey level or not. It could be just any colors...) The only way that I can do right now is overlaying

Hello, When we use optim and run into errors, such as generates NA/NaN/Inf, and the routine stops because of this error, is there a way to print the values of the argument of the functions where the error occurs? For example, I am doing minimizing negative log likelihood and when the optim stops because of those errors, I'd like to know what was the parameter values at that time. I heard

I'm sure this is probably simple, but I can't find an answer... I have a data frame (Galton's data on heights of parents and children), in grouped form, parent child frequency 73.5 72.2 1 73.5 73.2 3 72.5 68.2 1 72.5 69.2 2 72.5 70.2 1 72.5 71.2 2 72.5 72.2 7 72.5 73.2 2 72.5 74.2 4 71.5 65.2 1 71.5 66.2 3 71.5 67.2 4 71.5 68.2 3 71.5 69.2 5 71.5 70.2 10 ... and need the ungrouped

Full_Name: David Clifford Version: Version 1.5.1 (2002-06-17) OS: Red Hat 7.3 Submission from: (NULL) (128.135.149.55) For n values above 100 there appears to be a bug in contr.poly(n). The contrast matrix should have rank n-1. Running the code below gives output (ie errors) at n=98, 100 and every value greater than 102. for(n in 2:150) { K <- contr.poly(n) rnk <-

Hi, Thanks in advance for any advice you can give me, I am very stumped on this problem... I use R every day and consider myself a confident user, but this seems to be an elementary problem.. Outline of problem: I am analysing the results of a study on protein expression in cancer tissues. I have raw intensities from 2 different types of cancer and normal tissue, which can be taken from several

Hello list, Some days ago i sent an email to this ML about a problem im having with listing users on winXP/2k+3 computers. The thing is our samba PDC has +1600 users and, this time i counted, there only appear on the "security advanced settings" list box (where we select users that can read/write/execute on a windows directory share), about 500/550 users. There is no evident reason

Hi R users, I have a matrix of data similar to: > y=matrix(rnorm(55),ncol=5) I would like to know to which decile each number belongs compared to the numbers in its column. Say y[1,1] is the third decile among y[1:11,1] and y[2,1] is in the second decile I would like get a matrix that would return their ranks in decile, i.e., y[1,1] -> 3 y[2,1] -> 2 Your help is much appreciated!

I'm trying to do a sunflowerplot of Galton's data, with both regression lines and data ellipses, and I must be doing something wrong, because the lines do not intersect at \bar{x}, \bar{y}. The problem is likely in the line for x.y, but I don't know how to specify that correctly. The data is read in grouped form( galton), and then ungrouped (galton2): galton <-

Hello, Do you have any way to detect the client which is connecting to our IMAP server? I actually have an mail server which use dovecot but I just want to allow mobile device access to this server to get email, not desktop as Outlook, Thunderbird ... Thanks for your advice. -- Regards, Thu NGUYEN

Faraway's book titled "Practical Regression and Anova using R", with full text available online at: http://cran.r-project.org/doc/contrib/Faraway-PRA.pdf refers to a data set, stat500, which compares midterm and final grades. It can be used to illustrate similar concepts. A google search for faraway.zip will locate the actual data. --- Date: Sun, 15 Feb 2004 10:37:08 -0800

Dear R-Help, I'm trying to fit a model of the following form using gnls. I've fitted it using nlsList with the following syntax: nlsList(Y~log(exp(a0-a1*X)+exp(b0-b1*X))|K,start=list (a0=6,a1=0.2,b0=4.5,b1=0.001),data=data.frame(Y=y,X=X,K=k))) which works just fine: <snip> Coefficients: a0 a1 b0 b1 1 5.459381 0.5006811 5.137458 -0.0040548687

how do I check the rank of a matrix ? say A= 1 0 0 0 1 0 then rank(A)=2 what is this function? thanks I did try help.search("rank"), but all the returned help information seem irrelevant to what I want. I would like to know how people search for help information like this. rank(base) Sample Ranks SignRank(stats) Distribution of the

Hi! Maybe this is a silly question, but I need the column rank (http://en.wikipedia.org/wiki/Rank_matrix) of a matrix and R function 'rank()' only gives me the ordering of the elements of my matrix. How can I compute the column rank of a matrix? Is there not an R equivalent to Matlab's 'rank()'? I've been browsing for a time now and I can't find anything, so any

Dear R listers, Just a simple question. If we say an nxm matrix (n>m) is full rank of m, does this mean that this matrix has linearly independent columns? They are the same definition or needs some proof? Thanks for your answer. Fred [[alternative HTML version deleted]]

Dear All, a short and maybe simple question: I have to rank all values in a matrix from 0 to X, [1] [2] [3] [4] [1] 0.1 2 0 3 [2] 50 3 3 1 [3] 100 1 1 0 [4] 100 2 2 0 0->0 0.1->1 2->2 3->3 50->4 100->5 (X=5) is there any function for this? i have looked in several packages (vegan, labdsv etc.) because I am working with species by site tables,

I"ve got a data frame with a selection of columns I want to compute a rank-order correlation matrix from without disturbing the original data frame. foo[,c("a","b","d","f","g")] What I wanted to do, intuitively, was: > cor(rank(foo[,c("a","b","d","f","g")])) but rank in that context

Dear all, Rank of a matrix depends on which factors? Only on rows or coumns? or both ? If there is a collinearlity in the variables ( columns ) does it effects the rank? > X<-matrix((rnorm(10000)),50) > dim(X) [1] 50 200 > qr(X)$rank [1] 50 > X[,2]<-X[,30] > qr(X)$rank [1] 50 > X[10,]<-X[7,] > qr(X)$rank [1] 49 Thanks Alex [[alternative HTML

Hello R'users, I have a quick question. I wanted to know if there exist a function in R to compute the rank of a matrix. I could not find anything about it. Thank you, Raphael -------------- next part -------------- A non-text attachment was scrubbed... Name: raph.vcf Type: text/x-vcard Size: 303 bytes Desc: Card for Raphael Gottardo Url :

In summary.manova the qr decomposition of a NxN matrix is calculated and for some cases is giving me a rank < N. However, following suggestions of professor Ripley to calculate the rank of a Matrix On 7 Jun 2002, Brian Ripley wrote: > For a more reliable answer, look at the SVD > (function svd) and look at the > singular values. For example (from lda.default) X.s <-

? 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/20060728/957eb23c/attachment.pl