similar to: A matrix problem

Is there some (much) more efficient way to do this? VECTOR=c(3,2,4,5,5,3,3,5,1,6,6); NEWVECTOR=VECTOR[1]; for(i in 1:(length(VECTOR)-1)) { if((identical(VECTOR[i], VECTOR[i+1]))==FALSE){ NEWVECTOR=c(NEWVECTOR,VECTOR[i+1])} } > VECTOR [1] 3 2 4 5 5 3 3 5 1 6 6 > NEWVECTOR [1] 3 2 4 5 3 5 1 6 _______________________________ Atte Tenkanen University of Turku, Finland

Hi, Can somebody tell me, which is the fastest way to make comparisons between all rows in a matrix (here A) and put the results to the new symmetric matrix? I have here used cosine distance as an example, but the comparison function can be any other, euclidean dist etc. A=rbind(c(2,3),c(4,5),c(-1,2),c(5,6)) M=matrix(nrow=length(A[,1]),ncol=length(A[,1])) for(i in 1:length(A[,1])) { for(j in

Is there some way to change the line widths of plot borders? I couldn't find any parameters for that purpose. Atte Tenkanen University of Turku, Finland

I have two numeric vectors each of length 17 and each is named the exact same way. so obsnum ppppp ppppm pppmp . dot dot dot...... temp1 is 1417 52 63 85 obsnum ppppp ppppm pppmp . dot dot dot...... temp2 is 1213 41 50 97 what i want to have is a resultant matrix with 2 rows and 16 columns where the 16

Is there any possibility to divide too long text in a plot to two or more lines, when using labels-parameter in the text()-command? Here is an example picture: http://users.utu.fi/attenka/253.jpeg My example script is something like this: text(1,0.7,labels=Chordnames[fnid(pcs%%12)]) # according to Larry Solomon's table http://solomonsmusic.net/pcsets.htm Chordnames is a long vector with

Hi, How do you get the most common row from a matrix? If I have a matrix like this array(1:3,dim=c(4,5)) [,1] [,2] [,3] [,4] [,5] [1,] 1 2 3 1 2 [2,] 2 3 1 2 3 [3,] 3 1 2 3 1 [4,] 1 2 3 1 2 in which rows 1 and 4 are similar, I want to find that vector c (1,2,3,1,2). Atte Tenkanen University of Turku, Finland

Hi, What's wrong here?: > v=c(`-`,`+`,1,`^`,`^`,NA,NA,"X",9,"X",2) > i2=16 > v[i2] [[1]] NULL > is.null(v[i2]) [1] FALSE Is it a bug or have I misunderstood something? Atte Tenkanen University of Turku, Finland

I'd like to pick every imbricated five character long subsets from a vector. I guess there is some efficient way to do this without loops... Here is a for-loop-version and a model for output: VECTOR=c(1,4,2,6,5,0,11,10,4,3,6,8,6); ADDRESSES=c(); for(i in 1:(length(VECTOR)-4)){ ADDRESSES[i]=paste(VECTOR[i:(i+4)],collapse="") } > ADDRESSES [1] "14265"

Could somebody program this kind of plot type to R, if none exists, based on mds or correlation tables or some more suitable method? What do you think about idea? Does it work? None similar or better exists? http://weightedassociationmap.blogspot.com/ Atte Tenkanen University of Turku, Finland

Hi, Is it possible to break using if-condition during the recursive function? Here is a function which almost works. It is for inorder-tree-walk. iotw<-function(v,i,Stack,Indexes) # input: a vector and the first index (1), Stack=c(), Indexes=c(). { print(Indexes) # if (sum(i)==0) break # Doesn't work... if (is.na(v[i])==FALSE & is.null(unlist(v[i]))==FALSE) {Stack=c(i,Stack);

Dear R-users, How do you save a big table or matrix as an independent object and attach it to your Rdata-project when needed? Atte Tenkanen University of Turku, Finland Department of Musicology +023335278

Dear R-users, I want to check if certain values are from random distribution, that includes values between 0-1. So, it is not really normal even though shapiro.test says it is highly normal... Can I do something like this and think that the values given are right. z.test is from package TeachingDemos. -------------------------------------------------------------------------------

Hello, Before I was able to interrupt R-calculations by pushing ESC, but now with R 2.6.2 (OSX Intel) it doesn't work. Is it bug or not? Regards, Atte Tenkanen University of Turku, Finland Department of Musicology +023335278

Dear R users, I wonder if it is possible to form a function from a character string. Here is an example: > x=3 > `-`(`+`(`^`(x,3),`^`(x,2)),1) # Here is my function evaluated. [1] 35 >

Dear R-users, I would like to test, whether a sample distribution differs significantly from a population distribution. They are not normally distributed. How should I proceed? Using somehow glm-models? How? The population and the sample data are here. They can be loaded using the load-command. http://users.utu.fi/attenka/D_Pop http://users.utu.fi/attenka/D_Samp Best regards, Atte Tenkanen

Hi, I would like to know whether there is some deeper rationale behind or is it just an established practice that the lenghts of principal components, giving for example by prcomp-function, are normalised to 1? Best regards, Atte Tenkanen University of Turku, Finland Department of Musicology +35823335278 http://users.utu.fi/attenka/

Hi, I'd like to know if it is possible to use Hershey vector fonts to create very primitive musical notation. If I can hang some whole notes on these lines X11() plot(0,0, xlim=c(0,10), ylim=c(0,10)) # Staves: for (i in c(seq(from=2,to=2.8,by=0.2),seq(from=4,to=4.8,by=0.2))) { abline(h=i) } it is enough. Best wishes, Atte Tenkanen University of Turku, Finland

Hello, I have tried to change the margin widths so that mtext (here "sd of consecutive pc intervals", look at the picture) and plot(...,xlab="bar") fits to the picture. Here is an example: http://users.utu.fi/attenka/margins.png This doesn't help: par(mar=c(5.1, 7.1, 4.1, 2.1)) And here are the commands: png(filename="/Users/kone/Vaitostutkimus/Pictures/

Hi, How can I test, whether all the values that is.element() -function returns are TRUE's Eg. > (is.element(c(1,4,2),c(1,2,3))) [1] TRUE FALSE TRUE This doesn't work: > (is.element(c(1,4,2),c(1,2,3)))==TRUE [1] TRUE FALSE TRUE Best, Atte Tenkanen, FT, MuM http://users.utu.fi/attenka/

Hi, How would you calculate the "mean colour" of several colours, for example c("#FF7C00","#00BF40","#FFFF00")? Yours, Atte Tenkanen