search for: statmethods

...ke to do, but don't know how to do it in R. The site should be particularly useful for users of statistical packages like SAS, SPSS, Stata, and Systat who would like to get up and running in R quickly. However, I think that it will also be useful for even veteran R users. The URL is http://www.statmethods.net Topics are organized around these headings: * The R Interface * Data Input * Data Management * Basic Statistics * Advanced Statistics * Basic Graphs * Advanced Graphs Topic range from GUIs to creating publication quality output, from database access to date values, from t-tests to generalized...

...kyphosis) plot(fit) text(fit, use.n=TRUE) The text labels represent the count of each class at the leaf node. Unfortunately, the numbers are rounded and in scientific notation rather than the exact number of examples sorted by that node in each class. The plot is supposed to look like http://www.statmethods.net/advstats/images/ctree.png as per http://www.statmethods.net/advstats/cart.html. I'm running 2.14.1 on a mac. Can anyone verify or point out if I am doing something obviously wrong for displaying the counts rounded and in scientific notation rather than the true counts in each class at eac...

Dear All, In a package, I want to use PCA function. The structure I used follow this page: http://www.statmethods.net/advstats/factor.html. fit<-principle(mydata, nfactors=9, rotation=TRUE) or: result<-PCA(mydata) But I don't known why R language in my computer noticed: "not found principle", "not found PCA". I download and installed R-2.6.2-win32.exe. Thanks...

Hola, escribo porque tengo una gran duda como gráficas tan fáciles en Excel son tan difíciles en R? No consigo hacer una gráfica en R de estas características. Osea como puedo dividir por sitios por tamaños. En una planilla enorme tengo la primera columna son las estaciones (pero solo quiero las tres ultimas) y después muchos parámetros, en las columnas 46, 53 y 60 los promedios ponderados. Y no

I ran this code (several times) from the Quick-R web page ( http://www.statmethods.net/advstats/cart.html) but my cross-validation errors increase instead of decrease (same thing happens with an unrelated data set). Why does this happen? Am I doing something wrong? # Classification Tree with rpart library(rpart) # grow tree fit <- rpart(Kyphosis ~ Age + Number + Start, m...

...562706518953 0.0151515151515152, 0.368991215730364 0.0263157894736842, 0.256406702156839 0.0263157894736842, 0.344003157058329 0.0303030303030303, 0.200307950418176 0.0303030303030303, 0.558093143666938 0.04, 0.363786583569412 0.04, 0.462668853163967 and I'd like to have a plot like http://www.statmethods.net/stats/images/meanplot.jpg Position shall be plotted to the x-axis and Relevance to y-axis. Foreach value that occurs in the column Position the of all corresponding values in the second column and their variance shall be plotted. Furthermore all means shall be connected by a curve. I have no i...

Hello all, I would like to compare the similarity of two cluster solutions using a validation criteria (such as Hubert's gamma coefficient, the Dunn index the corrected rand index and so on) I see (from here:http://www.statmethods.net/advstats/cluster.html) that the function cluster.stats() in the fpc package provides a mechanism for comparing 2 cluster solutions - *BUT* - it requires me to give the the distance matrix among objects. *My question *is: What ways can you suggest for comparing two cluster solutions, while usin...

Hi, I would like to display inter-parameter scatter plots like those with the corrgram package (see upper triangle here: http://www.statmethods.net/advgraphs/images/corrgram2.png ), just that I would like to plot two datasets instead of one. Say one with black and one with red dots. Or a merged dataset where an indicator column is used to assign different colors to particular dots - with still the same layout. Is there a method to generat...

I am interested in applying a "classification tree" analysis where the response variable is a censored variable (survival data). I've discovered the package 'party' through this page: http://www.statmethods.net/advstats/cart.html. However, as my sample is not very big I would like to apply 'bootstrap' and use 'random forests', but with my censored response variable. Are there any packages for that?? Looking forward to your answer, -- vicent @vginer_upv [[alternative HTML version...

Hello, I am trying to write a script with the end goal of graphing power (y) as a result of sample size (x) at a variety of effects sizes. I am new to loops, and I think my problem is there. Here's the script, which is modified from the script found at the bootom of http://www.statmethods.net/stats/power.html. ANy help would be much appreciated!mtesche at ualberta.ca library(pwr) library(lattice) # range of sample sizes r <- seq(5,300,5) nr <- length(r) # effect sizes p <- seq(.4,1.2,.2) np <- length(p) # obtain power power <- array(numeric(nr*np), dim=c(nr,np))...

Hi, This is my first time so bear with me if I do anything silly. I'll learn! I'm a relatively new user of R as well so will appreciate any help. I have a data set in the form Trial.Group MeanHeart Rate Upper confidence level Lower confidence level 333subj: 0-5 M 80 60 120etcetc I want to plot a graph that has Trial.Group on the x axis(names, not just an

Hi! as my subject says I am struggling with the different of a two-way ANOVA and a (two-way) ANCOVA. I found the following examples from this webpage: http://www.statmethods.net/stats/anova.html # One Way Anova (Completely Randomized Design) fit <- aov(y ~ A, data=mydataframe) # Randomized Block Design (B is the blocking factor) fit <- aov(y ~ A + B, data=mydataframe) # Two Way Factorial Design fit <- aov(y ~ A + B + A:B, data=mydataframe) fit <- aov(y...

...OM 0.2317915 0.1610068 1.6941637 2 AAPL 0.6735488 0.7433611 0.1594102 3 GE 1.2554160 0.9237384 1.6767711 4 IBM 1.6296938 0.3730387 0.5858115 5 CVX 0.9194169 0.4785705 0.1803601 6 PG 0.7768241 1.7622060 0.7640163 . . . I would like to produce something similar to what is found at: http://www.statmethods.net/graphs/bar.html # the grouped barplot example or http://had.co.nz/ggplot2/geom_bar.html # the Dodged bar charts example Across the X-axis, for each set(row) of 3 data points(V2, V3, V4) associated with a symbol(V1), I would like to create a group of 3 bars reflecting th...

MY REQUIREMENT: I have a database and i want to connect to the database through R Project and using the tables from the database i need to generate the reports. can anyone plese tell me whether this is possible with R project or not. Thanks in advance Chinna. -- View this message in context: http://n4.nabble.com/Generating-reports-from-database-tp1564238p1564238.html Sent from the R help

Hi, excuse me for my english, i am using R on windows and i have to do several graphs with axis labels and the axis text thicks has a specified font type, (Arial) and a specified font size. How can i do these? Thank you in advance -- View this message in context: http://www.nabble.com/graph%3A-axis-label-font-tp24463974p24463974.html Sent from the R help mailing list archive at Nabble.com.

...classification tree using either tree or rpart functions but when it comes to plotting the results the formatting I get is different than what I see in all the tutorials (like http://www.youtube.com/watch?v=9XNhqO1bu0A or http://www.youtube.com/watch?v=m3mLNpeke0I&feature=related or http://www.statmethods.net/advstats/cart.html "tree for kyphosis"). I am trying to take a large demographic population and create a tree which systematically and accurately divides them into 2 pre-defined classifications using multiple predictor variables. What I would like to see is what I have seen in the t...

Hola: Estoy intentando dibujar una espiral a la que le quería colocar unas marcas equidistantes. Estoy utilizando: T = seq(0, 20*pi, length.out=1000) X = T * cos(T) Y = T * sin(T) plot(X, Y) que dibuja puntos siguiendo una espiral. Como podría convertirlo en una espiral continua con marcas equidistantes? He buscado en r-cran, y no he encontrado nada. Muchas gracias y saludos!! Griera

Buenas tardes: Quería lograr el efecto "expansivo" que se consigue con la opción paper="a4" de la función pdf, pero sin volcar en pdf; es decir, yo quiero enviar mi gráfico (luego definiré sus propiedades) a la salida gráfica (RGui) de modo que quede con verticalidad, como en un A4, y luego el usuario final es libre de guardar en pdf o no, con el propio RGui. Ahora voy a

Hi, I'm stuck with parsing data into R for heatmap representation. The data looks like: 1 id1 x1 x2 x3 .... x20 2 id1 x1 x2 x3 .... x20 3 id1 x1 x2 x3 .... x20 4 id1 x1 x2 x3 .... x20 ......... 348 id2 x1 x2 x3 .... x20 349 id2 x1 x2 x3 .... x20 350 id2 x1 x2 x3 .... x20 351 id2 x1 x2 x3 .... x20 ......... The data is sorted for the IDs (id1,id2 .....id40) and I like to

Hello dear R help list, I wish to compile a check-list for diagnostic measures for OLS regression. My question: Can you offer more (or newer) tests/measures for the validity of a linear model then what is given here: http://www.statmethods.net/stats/rdiagnostics.html This resource gives a list of measures to test for: OUTLIERS, INFLUENTIAL OBSERVATIONS, NON-NORMALITY, NON-CONSTANT ERROR VARIANCE, MULTI-COLLINEARITY, NONLINEARITY, NON-INDEPENDENCE OF ERRORS and some global validation. I came across it after searching online for ways...