similar to: clustering in R

Hi , I want to assign names to the rows in my matrix so that when I use the 'agnes' function from R cluster package , the dendogram that is produced represents the rows of the matrix. This way I would know what elements in the matrix are clustered together. But when I do the following, rownames(distF) <- paste("A",dim(distF)[1]) it gives the following error; Error in

Hi, I'm trying to calculate euclidean distance for my matrix rx1. Its a big matrix with 947 elements. However when i do the following, for(k in 1:947) for(m in (k+1):947){ A =rbind(A,c(k,m,ED(rx1[,k], rx1[,m]))) } I get the following error: ED(rx1[,k], rx1[,m]) : subscript out of bounds I checked for k in 1:100 and it runs fine; its a huge matrix and takes atleast an hour to run before I

Hi, I am very new to R. Just started yesterday. I have to generate a sequence of 40 random treatments using permuted block randomization with a block size of 8. Then I have to plot moving averages for the resulting sequence. I have tried the sample function but I dont know if wht i did is right or not > sample(c(1,1,1,1,2,2,2,2),8) Any ideas please? -- Ayesha Khan MS Bioinformatics Dept. of

Hi, I'm trying to do a non-metric multidimensional scaling using isoMDS. However, I have some '0' distances in my data, and I'm not sure how to deal with them. I'd rather not drop rows from the original data, as I am comparing several datasets (morphology and molecular data) for the same individuals, and it's interesting to see how much morphological variation can be

Hi All I am tying to create a user in SAMBA 4 AD with ldapadd, but its not allow me to add multiple mail attributes, here is my ldif of user. dn: CN=ayesha,CN=Users,DC=samba4pdc,DC=net objectClass: top objectClass: person objectClass: organizationalPerson objectClass: user cn: ayesha name: ayesha sAMAccountName: ayesha userPrincipalName: Ayesha Umar objectCategory:

Ho Mourik Thanks for your reply, any other attribute which we can duplicate? Br. Umar On Fri, Jan 23, 2015 at 1:47 PM, mourik jan heupink - merit < heupink at merit.unu.edu> wrote: > Hi, > > In AD, the attribute mail can only exist once. > > MJ > > On 01/23/2015 05:27 AM, Umar Draz wrote: > >> Hi All >> >> I am tying to create a user in SAMBA 4

Dear all, this might be not only an R-question but also a statistical. When I do a logistic regression analysis (species distribution modeling) with function lrm (Design package) I get the follwoing error message: > tadl1<-lrm(triad~fd+dista+fd2+dista2+fd:dista+dista:geo2, x=T, y=T) Error in if (!length(fname) || !any(fname == zname)) { : missing value where TRUE/FALSE needed The

Hi, I have tried looking at the archives but havent found any answer that works till now (Sorry if i have missed anything) I am a newbie to R and i am trying to carry out hierarchical clustering using hclust -> as.dendrogram and then plotting the results as a dendrogram using the plot function plot(object). My question is : In the function "plot", can one decrease the leaf label

Hi the list, I need to include some C code in R, but the behavior of the C code is strange : Here is my code : --- 8< --- Rprintf("\n XXXX mTraj=%f mClus=%f",mTraj[i+nbId*c],mClustersCenter[j+nbClusters*c]); Rprintf("\nDistA=%d Tmp=%d",dist,tmp); tmp = mTraj[i+nbId* c] - mClustersCenter [j+nbClusters* c]; Rprintf("\nDistB=%d

Hi. I am getting an "unexpected symbol" error, which I haven't gotten before. I have checked my data and searched the FAQ. I hoping someone has a suggestion on what I might be doing wrong? Thanks! A. Gray Fisheries Consultant Here's the code I'm using and the error: > 0708smolts <- read.csv("C:\Documents and Settings\ayesha\My

Hi folks, I'm running several samba servers that use "security= server" and "password server" to point authentications to a WinNT pdc. Snooping at the network traffic, I find that the samba server initially tries to connect to the WinNT server through port 445 (which fails), then reverts to 139 (which succeeds). I know it's possible to use the "samba

Hey Everyone! I wanted to carry out Hierarchical clustering using distance matrices i have calculated ( instead of euclidean distance etc.) I understand as.dist is the function for this, but the distances in the dendrogram i got by using the following script(1) were not the distances defined in my distance matrices. script: var<-read.table("the distance matrix i calculated",

On 23/01/15 09:44, Umar Draz wrote: > Ho Mourik > > Thanks for your reply, any other attribute which we can duplicate? > > Br. > > Umar > > On Fri, Jan 23, 2015 at 1:47 PM, mourik jan heupink - merit < > heupink at merit.unu.edu> wrote: > >> Hi, >> >> In AD, the attribute mail can only exist once. >> >> MJ >> >> On

Dear all, below you find enclosed my message from January 9th and my program (attention: beginner). Thanks for both answers! a. However, as far as I know the cosine is not the same as the Pearson correlation (only in special cases). b. Reid Huntsinger's hint was very useful, however I had to transpose the matrix first, for I want to calculate the distance of the _rows_. Regards, Petra

Metrum Institute (http://metruminstitute.org/) is offering a postdoctoral fellowship in pharmacometrics. It will focus on the development of computational methods and open-source software tools for Bayesian modeling and simulation relevant to pharmacometric applications. An ideal candidate will have a post-graduate degree in statistics, pharmacokinetics, bioengineering or other related

I am not sure whether there is a bug. When I tested the example given for "glht" in the help, I entered the following error: Running commands: amod <- aov(minutes ~ blanket, data = recovery) rht <- glht(amod, linfct = mcp(blanket = "Dunnett"), alternative = "less") Errors are: Error in try(coef.(model)) : could not find function

Dear R-users, For the purpose of model selection I am looking for a way to exhaustively (and efficiently) search for best subsets of predictor variables for a logistic regression model. I am looking for something like leaps() but that works with glm. Any feedback highly appreciated. -- Harald von Waldow <hvwaldow at chem.ethz.ch> Safety and Environmental Technology Group Institute for

Dear Peter, I am running your g.test() with the william's correction but I have a question about the input numbers. These are my data: "Our data are consistent with those obtained using microarray comparative genome hybridization in that we found significantly fewer variants per Mb on the X compared to the autosomal chromosomes (152 versus 336 respectively, G = 93.4, P < 2e-16, df =

Two Postdoctoral fellows (fulltime) Leiden University, Leiden, The Netherlands Job description ________________________________________ TOP Institute Pharma (TI Pharma) has granted our proposal to set up a mechanism-based PK-PD modelling platform. This platform focuses on the transfer of knowledge from academia to the pharmaceutical industry and is a collaborative effort of

________________________________ Principal Post-Doctoral Fellow PKPD modelling platform (full time) Leiden University, Leiden, The Netherlands <javascript:doredirect();> Job description ________________________________ TOP Institute Pharma (TI Pharma) has granted our proposal to set up a mechanism-based PK-PD modelling platform. This platform focuses on the transfer of