similar to: boxplot - bclust

Dear colleagues, What did I not understand ? ->my intention I want to create a new variable: In plain language: If someone is taking anithypertensive treatment (med.hyper==1) table(med.hyper) med.hyper 0 1 472 97 I want to subtract 5 mmHg (rr.dia.2m-5) from the measured diastolic blood pressure (rr.dia.2m) if not treated - the value of the measured diastolic blood pressure should

Hi everyone, Currently i'm trying to understand the bagged clustering algorithm, bclust {e1071}. When I run the given example in the help file (as below) data(iris) bc1 <- bclust(iris[,1:4], 3, base.centers=5) plot(bc1) and plot the bclust object, 2 graphs are produced. The first is a dendrogram, but what is the second plot? The axes are not labelled and what do the two

Hi All, Could you please help with the error I get from the following codes? > Library(cluster) > data(iris) > bc1 <- bclust(iris[,2:5], 3, base.method="clara", base.centers=5) Then I got: "Committee Member: 1(1)(2)(3)(4)(5)(6)(7)(8)(9)(10)(11)(12)(13)(14)(15)(16)(17)(18)(19)(20) Error in bclust(iris[, 2:5], 3, base.method = "clara",

Dear All: I'm a newbie to R and chemometrics. Now I'm trying apply bclust on fuzzy c-means like this: >bc1 <- bclust(iris[,1:4], 3, base.centers=20,iter.base=100, base.method="cmeans") Committee Member: 1(1)(2)(3)(4)(5)(6)(7)(8)(9)(10)(11)(12)(13)(14)(15)(16)(17)(18)(19)(20)Erro r in bclust(iris[, 1:4], 3, base.centers = 20, iter.base = 100, base.method =

Hi All, Just want to make sure, in function "bclust", do the following argument only have one option? argument "dist.method" has one option "Euclidian"; argument "hclust.method" has one option "average"; argument "base.method" has one option "kmeans". Thank you! [[alternative HTML version deleted]]

Dear all, I have a problem with the function "bagged clustering" of package e1071. When I try to run the example of bagged clustering with the iris data : data(iris) bc1 <- bclust(iris[,1:4], 3, base.centers=5) I got the following message error : Loading required package: class Committee Member: 1(1) 2(1) 3(1) 4(1) 5(1) 6(1) 7(1) 8(1) 9(1) 10(1)Error in bclust(iris[, 1:4], 3,

Hi, I am performing experiments in the field of visual perception where we often apply balanced designs. Within a group of normal subjects, we vary different stimulus conditions (like contrast, luminance, temporal frequency of stimulation) and derive some psychophysical or electrophysiological results from our subjects. Often, the main question is to test the effect of these parameters

Hi, I have a very large dataset of 3008 individuals and 800 numerical variables. In fact it is a table of 3008 36-monthes multivariated time series that I would like to classify with an unsupervised algorithm I had a look at the function kkmeans of e1071 package, which seems to be a kernel weighted version of the algorithm algorithm, and the bclust from the same package which does bootstrapping

Hi All, I run function "kmeans" to cluster a matrix. But when the matrix size is big enough, it keeps saying "did not converge in 10 iterations". Could you explain what it means and if the result is wrong? And the interesting thing is sometimes this warning happens when the sample size is around 51200 x 6, sometimes it happens around 30000 x 6. Does the warning related

Hi All, For the function "bclust"(e1071), the argument "base.method" is explained as "must be the name of a partitioning cluster function returning a list with the same components as the return value of 'kmeans'. In my understanding, there are three partitioning cluster functions in R, which are "clara, pam, fanny". Then I check each of them to

Hi all, I am trying to use in RStudio the latest code given in https://github.com/systematicinvestor/SIT/blob/master/R/bt.test.r, which seems to work fine but with the following warning for download limits (one for each of the tickers). I searched in options() something which could be related to this setting, w/o success. Any hint for me in order to raise or remove these limits? Where is this

Dear list I tried to read an xml file using the xml package. Unfortunately, some encoding problems occure. E.g. german Umlaut will be red correctly. I assume that the occurs due to (internal?) conversion to utf-8. To illustrate the problem, I have wrote to xml files. File Test 1 ----------- <?xml version="1.0" encoding="ISO-8859-1"?> <Daten> <ITEM>

Hello, I have just learnt some python, and would like to call R scripts from python, without using batch mode. I learnt how to do it from perl, in fact reading in the mail archive I found the following good example: *************************************** #!/usr/local/bin/perl open (FILE, ">test.txt"); print FILE "a,b,c,d,e\n1,2,3,4,5"; close FILE; ####### Start R

Hi all, I would like to fit a gamma pdf to my data using the method of RSS (Residual Sum-of-Squares). Here are the data: x <- c(86, 90, 94, 98, 102, 106, 110, 114, 118, 122, 126, 130, 134, 138, 142, 146, 150, 154, 158, 162, 166, 170, 174) y <- c(2, 5, 10, 17, 26, 60, 94, 128, 137, 128, 77, 68, 65, 60, 51, 26, 17, 9, 5, 2, 3, 7, 3) I have typed the following code, using nls method:

There are two misspellings in the german Error message for friedman test: Fehler in friedman.test.default(cont$score, group = cont$goup, blocks = cont$cont) : y, Gruppen und blöcke müssen die sekbe Länge haben The correct spelling would be: Fehler in friedman.test.default(cont$score, group = cont$goup, blocks = cont$cont) : y, Gruppen und *Blöcke* müssen die *selbe* Länge haben Instead of

Hello, I have written perl programs that extract data from a text file, process them, and create other text files, which I'd like to apply some statistics too (for example with R). I'd like to do it all in once , with a single script. I'm not familiar with R, I'd like to know if this task could be accomplished by creating a linux shells that launches the perl scripts and then

Hello together, this is ma first post, so please aplogize me if post this in the wrong section. I have problem concerning ma two matrizes. After a regressione and so on, I got two matrizes Matrixres contains the results of ma calculation. Matrixr contains my detiene, which where Aldo used for the regression. Please ser the following code: #Datei einlesen residual =

Hi, I want to calculate the t- and p-values for a linear model using the Newey West estimator. I tried this Code and it usually worked just fine: > oberlm <- lm(DYH ~ BIP + Infl + EOil, data=HU_H) > coeftest(oberlm, NeweyWest(oberlm, lag=2)) t test of coefficients: Estimate Std. Error t value Pr(>|t|) (Intercept) 0.1509950 0.0743832 2.0300 0.179486 BIP

Dear all, in one part of my code I want to sum two vectors element-wise the problem is that either the 1st vector or the 2nd vector always have one or two less elements example of my problem In TotalVector + (datalist2[[1]]$dataset$Results[[j]]$Results[[time]]$Sweep) :   Länge des längeren Objektes          ist kein Vielfaches der Länge des kürzeren Objektes Browse[1]> str(TotalVector)  int

Dear all, I am new R user and I am sure that this question has been asked quite often and I have also googled it and read about it! I understood that in order to read excel sheet into R you need to open it and saved it as csv or text, is this true? or you can use read.delim2 and read.csv2 to do this without the following error > dat <- read.csv2(file="C:\\Dokumente und