Re: [R] Problem with dummy variables controlling for years

Hi Johannes. My guess is that you will get no Std. Error, if you have only 1 observation for that year. Could that be the case ? On Mon, Mar 2, 2009 at 8:41 PM, Johannes Habel wrote: > > Hello everybody, > > I'm running a panel regression with dummy variables controlling for years. > While thi

Re: [R] R-help

So, it seems you are using Windows. One possibilty is that your Windows security settings are stopping you displaying CHM help files: if you have control over these you could try to resolve that. Another alternative is to change the default type of help used: the settings for that are in the

Re: [R] R - need more memory, or rejection sampling algorithm doesn't work?

> b12=log(.4/.6) > b13=log(.85/.15) Can you point out the syntax errors in the first one? I have been using a couple of books for the code but I am no expert. The choice of the runif(1,0,2) density and choice of 2 as the constant in the envelope function was made to envelope the whole variable;

Re: [R] vertically aligned X axis labels disappear off R Graphics window

On Wed, Feb 25, 2009 at 12:51 PM, Uwe Ligges wrote: > > > R User R User wrote: >> >> Hi guys, >> I'm evaluating R for basic data exploration. I produce a bar plot of the >> data, with the x axis labels aligned vertically. However, the start of >> labels longer than about 10 characters are cut off

Re: [R] regular expression question

G'day Mark, On Tue, 03 Mar 2009 00:16:34 -0600 (CST) markle...@verizon.net wrote: > can someone show me how to use a regular expression to break the > string at the bottom up into its three components : > > (-0.791,-0.263] > (-38,-1.24] > (0.96,2.43] > > I tried to use strplit because of my reg

Re: [R] Large data set in R

Thanks Kjetil. This is exactly what I wanted. Hardi From: Kjetil Halvorsen Cc: r-help Sent: Monday, March 2, 2009 9:45:43 PM Subject: Re: [R] Large data set in R install.packages("biglm", dep=TRUE) library(help=biglm) kjetil Hello, I'm trying to us

Re: [R] [R Question]: regarding exponential axis

Sun, Larry amgen.com> writes: > How to plot exponential axis (i.e. equally spaced .1, 1, 10, 100..)? > I want to use it with xyplot; it would be great if there is a detailed > example. > Thanks a lot! > Sound like what is usually called a logarithmic axis. See parameter log in scales of xyplo

Re: [R] Writing R package and do.call

Thanks Rolf, I used your solution in my original, actual function and it worked very fine and a much more sophisticated solution than what I originally came up with! Thanks to the other people who replied as well -- I know that the message said it was a "Note" rather than a "Warning", but it sound

Re: [R] R - need more memory, or rejection sampling algorithm doesn't work?

> -Original Message- > From: r-help-boun...@r-project.org > [mailto:r-help-boun...@r-project.org] On Behalf Of ekwaters > Sent: Monday, March 02, 2009 9:39 PM > To: r-help@r-project.org > Subject: [R] R - need more memory, or rejection sampling > algorithm doesn't work? > > > Hi all, >

[R] regular expression question

can someone show me how to use a regular expression to break the string at the bottom up into its three components : (-0.791,-0.263] (-38,-1.24] (0.96,2.43] I tried to use strplit because of my regexpitis ( it's not curable. i've been to many doctors all over NYC. they tell me there's no cure

[R] R - need more memory, or rejection sampling algorithm doesn't work?

Hi all, I am trying to run rejection sampling for the quantity z11 in the function below. Unfortunately I can't simplify the function further so that z11 only appears once. Whenever I run the algorithm, R looks as if it is running it (no error messages or anything), but then nothing happens fo

[R] Error loading packages at sprintf()

Hello all, I'm running R2.8.1 on a Mac OS 10.4.11. While trying to install the package gdata, I was presented with the following (error at end of report): R > install.packages("gdata") also installing the dependency ‘gtools’ trying URL 'http://cran.stat.ucla.edu/bin/macosx/universal/contrib/2

Re: [R] Inefficiency of SAS Programming

for an " inefficient " language , it sure has dominated the predictive analytics world for 3 plus decades. I referred once to intellectual jealousy between newton and liebnitz. i am going ahead and creating the R package called "Anne". It basically is meant only for SAS users who want to learn R

Re: [R] rounding problem

> -Original Message- > From: r-help-boun...@r-project.org [mailto:r-help-boun...@r- > project.org] On Behalf Of Prof Brian Ripley > Sent: Monday, March 02, 2009 12:38 AM > To: tedzzx > Cc: r-help@r-project.org > Subject: Re: [R] rounding problem > I nominate the following for the fortune

Re: [R] Bold Face in Plot

Here is one way: thestring <- "the-actual-string" plot(0:1,0:1,type="n") text(x=0.5, y=0.5, labels=bquote(bold(.(thestring There is an example in the plotmath help page that shows using bquote (but that page has so much information that I did not find it on my first scan through even knowin

[R] SPSS data import: problems & work arounds for GSS surveys

I'm using R 2.8.1 on Ubuntu 8.10. I'm writing partly to ask what's wrong, partly to tell other users who search that there is a work around. The General Social Survey is a long standing series of surveys provided by NORC (National Opinion Research Center). I have downloaded some years of the sur

Re: [R] Inefficiency of SAS Programming

This does not really address my point. Yes, if the few nerds who want to do funny stuff are the ones making the purchase, then there is a good chance (but still not guaranteed) that they will get IML, but do all companies that buy SAS actually think about that, or do they just see the extra pri

Re: [R] portable R editor

Thanks for the plug on NppToR. Yes it is portable, but a few of the features don't work. I have it on my plan to have a launcher for the portable apps menu. Also the website for NppToR is now https://sourceforge.net/projects/npptor/ On my personal website I also keep, semi-up-to-date a portablea

Re: [R] how to pass a command variable in DOS to R program in R CMD BATCH

Create a file test.R containing this single line: print(commandArgs(TRUE)) and then try this: Rscript test.R a b c On Mon, Mar 2, 2009 at 6:22 PM, Xianming Wei wrote: > Hi all, > > > > I need to run a program (asreml) thousands of times and each time I have > to provide a slightly diffe

[R] [solved, please disregard] FW: how to pass a command variable in DOS to R program in R CMD BATCH

Regards, Xianming From: Xianming Wei Sent: Tuesday, 3 March 2009 9:22 AM To: 'r-help@r-project.org' Subject: how to pass a command variable in DOS to R program in R CMD BATCH Hi all, I need to run a program (asreml) thousands of times and each time I have to provide a slightly d

Re: [R] Writing R package and do.call

Try this: temp.fn.2 <- function(mean = 0, sd = 1, n = 101) { x <- seq(mean - sd, mean + sd, length = n) plot(x, dnorm(x), type = "l") } # test temp.fn.2() On Mon, Mar 2, 2009 at 8:46 PM, Rob Foxall wrote: > Dear R-Help, > > I have written a function, which in simplified format can be repr

Re: [R] Writing R package and do.call

On Mon, Mar 2, 2009 at 7:46 PM, Rob Foxall wrote: > Dear R-Help, > > I have written a function, which in simplified format can be represented as: > > temp.fn <- function(my.mean,my.sd){ >   Parameters <- list(mean = my.mean, sd = my.sd) >   curve(do.call(dnorm,c(list(x), Parameters)), from = my.me

Re: [R] Writing R package and do.call

On 3/03/2009, at 2:46 PM, Rob Foxall wrote: Dear R-Help, I have written a function, which in simplified format can be represented as: temp.fn <- function(my.mean,my.sd){ Parameters <- list(mean = my.mean, sd = my.sd) curve(do.call(dnorm,c(list(x), Parameters)), from = my.mean-my.sd,

Re: [R] Reading from Google Docs

Download and install each of R, Rtools and MiKTeX. Each of these have automated installers so its just a matter of pressing Enter repeatedly. You don't have to change your path if you place Rcmd.bat from http://batchfiles.googlecode.com anywhere on your path (so it can be found). To see your path

Re: [R] portable R editor

Werner Wernersen yahoo.de> writes: > > > Hi, > > I have been dreaming about a complete R environment on my USB stick for a long time. Now I finally want to > realize it but what I am missing is a good, portable editor for R which has tabs and syntax highlighting, can > execute code, has book

Re: [R] Reading from Google Docs

I was able to extract RGoogleDocs with 7-zip. That was when I began to appreciate the gravity of what Prof Ripley was saying. I am on a windows xp machine and the RGoogleDocs download does not have a windows binary version. Instead it has only a source package. I started using the R-admin manual to

[R] Writing R package and do.call

Dear R-Help, I have written a function, which in simplified format can be represented as: temp.fn <- function(my.mean,my.sd){ Parameters <- list(mean = my.mean, sd = my.sd) curve(do.call(dnorm,c(list(x), Parameters)), from = my.mean-my.sd, to = my.mean+my.sd) } This works as I want it do.

Re: [R] Cross Tables with odfTable in odfweave

on 03/02/2009 03:11 PM Max Kuhn wrote: > PJ, > >> Hi, I've been trying to prepare some crosstables for some survey questions >> for a client. I have been using the CrossTable function in the gmodels >> package. However, this command only seems to be able to create CrossTables >> in text documents.

[R] Two Column-group names using Latex(Hmisc)

Hi, Using the following code I can generate the attachment "Original.pdf": latex(test, file="C:/Documents and Settings/darej/My Documents/James/READY.tex", title = "", rowname = NULL, cgroup = c("Fishstock", stocknames), n.cgroup = c(1, rep(2,nostocks)), colhea

Re: [R] (no subject)

On Mar 2, 2009, at 9:57 AM, Brajkovic J. wrote: spec=garchSpec(model=list(ar= 0.440270860, omega=0.000374365,alpha=0.475446583 , mu=0, beta=0)) sim<-garchSim(spec, length(dp)-1,extended=T) I get: > spec=garchSpec(model=list(ar= 0.440270860, omega=0.000374365,alpha=0.475446583 , mu=0, bet

Re: [R] latex output of regressions with standardized regression coefficients and t-statistics based on Huber-White

Jan Schulz wrote: Hello, first of all: I'm new to R and have only used SPSS befor this (which can't do this at all...). I'm trying to output some regression results to latex. The regressions are normal OLS and I'm trying to output the results with standardized regression coefficients and t-

[R] latex output of regressions with standardized regression coefficients and t-statistics based on Huber-White

Hello, first of all: I'm new to R and have only used SPSS befor this (which can't do this at all...). I'm trying to output some regression results to latex. The regressions are normal OLS and I'm trying to output the results with standardized regression coefficients and t-statistics based on

[R] Unrealistic dispersion parameter for quasibinomial

I am running a binomial glm with response variable the no of mites of two species y->cbind(mitea,miteb) against two continuous variables (temperature and predatory mites) - see below. My model shows overdispersion as the residual deviance is 48.81 on 5 degrees of freedom. If I use quasibinomial t

Re: [R] power analysis for prop trend test

Breslow & Day has a nice three page discussion in volume 2 of their "Statistical Methods in Cancer Research". See pages 285-7. Most of the gain in power comes from the decrease in degrees of freedom and only if the trend is approximately linear. Alternatives that are quadratic are not well

Re: [R] plotting on map

On Mon, 2009-03-02 at 23:15 +, Dr. Alireza Zolfaghari wrote: > Hi list, > I want to plot state Texas using map function. I wrote the following code to > plot the grid boundary for Texas: > > library(maps) > require("mapproj") > longlatLimit<-c(-106.65,

Re: [R] Error setting rowname if rowname currently NULL

Hi, At the 1st step, you only assign a name for the 1st row. However, fred has 4 rows which mean you need to assign 4 rownames for these 4 rows. At the 2nd step, you only "CHANGE" the 1st rowname from 1 to APPLE If I am wrong please correct me. Thanks -Original Message- From: r-help-

Re: [R] Error setting rowname if rowname currently NULL

THe problem is that 'rownames(fred)' is NULL and therefore you will have to define the names of all the rows: > fred<-matrix(,4,2) > fred [,1] [,2] [1,] NA NA [2,] NA NA [3,] NA NA [4,] NA NA > rownames(fred) NULL > rownames(fred) <- c("Apple",'','','') > fred [,1] [,2]

[R] how to pass a command variable in DOS to R program in R CMD BATCH

Hi all, I need to run a program (asreml) thousands of times and each time I have to provide a slightly different dataset. Because I have to run asreml under Windows (DOS or scripts) environment, I have trouble to pass a command variable (or pointer variable or %counter in the following example)

[R] R-help

Hi list, When I type my question in R console using ? sign (example: ?par, when I want to query for par), the following error message popps up: Error in print.help_files_with_topic("C:/PROGRA~1/R/R-24~1.1/library/maps/chm/map") : CHM file could not be displayed I appreciate if any body c

[R] Error setting rowname if rowname currently NULL

Hi, My first post here and new to R so please bear with me (long time programmer though, helping a friend with some scripts). I've noticed a behaviour when using rownames() that I think is odd, wondering if I'm doing something wrong. To illustrate, say I create a very simple matrix (called fred

[R] plotting on map

Hi list, I want to plot state Texas using map function. I wrote the following code to plot the grid boundary for Texas: library(maps) require("mapproj") longlatLimit<-c(-106.65, -93.53 , 25.93 , 36.49) par(plt=c(0,1,0,1),cex=1,cex.main=1) #Set plotti

[R] power analysis for prop trend test

Hi all, Is there any way we can to power analysis for prop trend test? Many thanks! __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provi

Re: [R] Modifying a built-in R function

Ok, thank you for your guidelines. Additionally, I have had to put the biplotes.princomp function in the first place of the script, just after biplotes <- function(x, ...) UseMethod("biplotes"). By doing this it works, althrough I don´t really know why. I hope to acquire higher skills in R program

Re: [R] density > 1?

On Mon, 2009-03-02 at 13:27 +0100, Johannes Elias wrote: > Dear R-Gurus, > > I wonder why 'density' values as shown in hist or plot(density(x)) are > sometimes over 1. How can that be? > > Example > > >hist(rnorm(1000,sd=.5),freq=FALSE) > > The resulting plot shows density values below 1 on the

Re: [R] trying to get rid of for loop in simple scenario

Dear Juan Carlos, > x = c(1,1,3,3,6,7,8,8,9) > table(x) x 1 3 6 7 8 9 2 2 1 1 2 1 > z<-table(x) > z x 1 3 6 7 8 9 2 2 1 1 2 1 > as.vector(z) [1] 2 2 1 1 2 1 Take a look at ?table and ?vector for more information. HTH, Jorge On Mon, Mar 2, 2009 at 4:41 PM, Juancarlos Laguardia wrote: > Hi, > >

Re: [R] portable R editor

Hello Werner, Did you already try to use the new version of Tinn-R (2.2.0.2) and R in USB stick? I use it without problem! A good idea is to have a backup of your system configuration also in the USB stick. Unzip it (in any place) and replace the folder Tinn-R where Tinn-R stores all ini files i

Re: [R] TinnR & Philips Webcam

Tinn-R was programed (intentionally) to detect any opened windows with the word 'Console' inside of the caption as a R instance. So, for while, it is necessary to close others process that has the word 'Console' in the caption. It is because the necessity of Tinn-R detects others GUIs (JCR, SciVi

Re: [R] TinnR & Philips Webcam

Tinn-R was programed (intentionally) to detect any opened windows with the word 'Console' inside of the caption as a R instance. So, for while, it is necessary to close others process that has the word 'Console' in the caption. It is because the necessity of Tinn-R detects others GUIs (JCR, SciVi

Re: [R] Modifying a built-in R function

Alternatively, gives a looked in the package bpca and biplotGUI, both available on CRAN. HTH, JCFaria japal wrote: > > Hello, > > Something incredible (at least for me) has happen. Yesterday night I > downloaded biplot.R to edit this function and add new features I wished. > Namely I wanted

[R] trying to get rid of for loop in simple scenario

Hi, I have a problem in which i am having trouble vectorizing a snippet of my code to avoid a for loop. A basic example of what I'm trying to do is below. #simple example of problem x = c(1,1,3,3,6,7,8,8,9) # I create a vector that contains the unique values of x y = unique(x) y = c(1,3

Re: [R] rgl persp3d bounding box color problem.

On 3/2/2009 2:48 PM, chaogai wrote: Oops, sorry R version 2.8.1 (2008-12-22) i686-pc-linux-gnu locale: LC_CTYPE=en_US.UTF-8;LC_NUMERIC=C;LC_TIME=en_US.UTF-8;LC_COLLATE=en_US.UTF-8;LC_MONETARY=C;LC_MESSAGES=en_US.UTF-8;LC_PAPER=en_US.UTF-8;LC_NAME=C;LC_ADDRESS=C;LC_TELEPHONE=C;LC_MEASUREMENT=en_

Re: [R] logistic regression model validation through bootstrapping

Install the Design and Hmisc package then do: library(Design) ?validate.lrm ?calibrate Resampling methods such as the bootstrap and cross-validation assume you have done no model/variable selection outside the bootstrap look, i.e., you let the bootstrap repeat all modeling steps that used Y. v

[R] [R Question]: regarding exponential axis

Hi there, How to plot exponential axis (i.e. equally spaced .1, 1, 10, 100..)? I want to use it with xyplot; it would be great if there is a detailed example. Thanks a lot! Larry __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listin

Re: [R] Cross Tables with odfTable in odfweave

PJ, > Hi, I've been trying to prepare some crosstables for some survey questions > for a client. I have been using the CrossTable function in the gmodels > package. However, this command only seems to be able to create CrossTables > in text documents. > > I've been trying to use odfTable in odfwea

Re: [R] xy-plot alb font size

Try this with the built in CO2: xyplot(uptake ~ conc | Treatment, CO2, xlab = list("A1C", cex = 2), ylab = list("Fasting Plasma Glucose", cex = 2)) On Mon, Mar 2, 2009 at 3:47 PM, Veerappa Chetty wrote: > I would appreciate any help for the following question. How to increase the > font size

Re: [R] using par funtions in bargraph.CI()

On Mon, 2009-03-02 at 11:23 -0600, maxa0...@umn.edu wrote: > I'm trying to create a bargraph of means with standard error bars using the > function bargraph.CI (in the sciplot package). Like this: > > bargraph.CI(x.factor, response,data,xlab, ylab, par(family="serif"),font=11) > > However, an er

[R] xy-plot alb font size

I would appreciate any help for the following question. How to increase the font size of labels for the x and y axes in the following code? - xyplot(lbxglu~lbxgh|eth,reg.dat.5, panel = p, col = 1:2, par.strip.text = list(cex = 1.25,col="brown"), xlab="A1C",ylab="Fa

Re: [R] initial gradient and vmmin not finite

Hi June: As don't have your data (yogurtnp.csv), I cannot replicate your problem. What does fr( c(0,0,0,1,0,-25,-0.2) ) return? Arne On Monday, 2. March 2009 11:55:56, June Wong wrote: > Dear Rhelpers > > I have the problem with initial values, could you please tell me how to > solve it? Tha

Re: [R] ways to put multiple graphs on single page (using ggplot2)

I usually have a function like this: vplayout <- function(x, y) viewport(layout.pos.row=x, layout.pos.col=y) draw4 <- function(pngname, a,b,c,d,w,h) { png(pngname,width=w, height=h) grid.newpage() pushViewport(viewport(layout=grid.layout(2,2) ) ) print(a, vp=vplayout(1,1))

Re: [R] rgl persp3d bounding box color problem.

Oops, sorry R version 2.8.1 (2008-12-22) i686-pc-linux-gnu locale: LC_CTYPE=en_US.UTF-8;LC_NUMERIC=C;LC_TIME=en_US.UTF-8;LC_COLLATE=en_US.UTF-8;LC_MONETARY=C;LC_MESSAGES=en_US.UTF-8;LC_PAPER=en_US.UTF-8;LC_NAME=C;LC_ADDRESS=C;LC_TELEPHONE=C;LC_MEASUREMENT=en_US.UTF-8;LC_IDENTIFICATION=C attached

Re: [R] rgl persp3d bounding box color problem.

On 3/2/2009 2:00 PM, chaogai wrote: Whatever it is, it is also happening on my Acer, Suse 11.1, 32 bits, self build R. Versions of R and rgl? Duncan Murdoch Kees Bo Zhou wrote: Thanks Duncan. I haven't got a chance to try it in windows yet. But I won't be surprised if it's my driver's

[R] Cross Tables with odfTable in odfweave

Hi, I've been trying to prepare some crosstables for some survey questions for a client. I have been using the CrossTable function in the gmodels package. However, this command only seems to be able to create CrossTables in text documents. I've been trying to use odfTable in odfweave to create

Re: [R] ways to put multiple graphs on single page (using ggplot2)

Juliet Hannah gmail.com> writes: > > Hi, Here are three plots: > (How to layout on one page) library(ggplot2) data(diamonds) randind <- sample(nrow(diamonds),1000,replace=FALSE) dsmall <- diamonds[randind,] a= qplot(carat, data=dsmall, geom="histogram",binwidth=1) b= qplot(carat, data=dsmall,

Re: [R] Finding Lambda in Poisson distribution

Hi, assuming the assumption of Poisson distribution is true (or at least approximately true), you can use MLE as mentioned by Roland using (fitdistr of MASS, mle of stats4, MLEstimator of distrMod, probably more). Using package distrMod you can compute minimum distance estimators (cf. functi

Re: [R] ways to put multiple graphs on single page (using ggplot2)

Hi Juliet, Have a look at the last section in http://had.co.nz/ggplot2/book/polishing.pdf Hadley On Mon, Mar 2, 2009 at 12:35 PM, Juliet Hannah wrote: > Hi, Here are three plots: > > library(ggplot2) > data(diamonds) > randind <- sample(nrow(diamonds),1000,replace=FALSE) > dsmall <- diamonds[ra

[R] Problem with dummy variables controlling for years

Hello everybody, I'm running a panel regression with dummy variables controlling for years. While this works fine for most of my models, in some cases I get the error message "In sqrt(diag(object$vcov)) : NaNs wurden erzeugt" and missing results for some or all of my year dummies:

Re: [R] rgl persp3d bounding box color problem.

Whatever it is, it is also happening on my Acer, Suse 11.1, 32 bits, self build R. Kees Bo Zhou wrote: > Thanks Duncan. > > I haven't got a chance to try it in windows yet. But I won't be surprised if > it's my driver's fault. I had OpenGL problems on this dell laptop before. > > >> Date

Re: [R] logistic regression model validation through bootstrapping

Vivienne Ozohili wrote: >> I am currently trying to implement a boostrapping approach to validation >> but don't know where to start. Help please. Take a look at Prof. Harrell's ?validate.lrm in package Design. Regards, Mark. Vivienne_O.Ozohili wrote: > > > Hi, > > I was wondering whether

Re: [R] Help with Wilcoxon Test

Hi, type is(sampA) is(sampV) in the R-prompt. It will show you that at least one of the vectors is not a vector of numerical values in R's sense. Then you have to convert using the as.numeric() function, but you will have to make sure that this really gives you the original vectors. Since your

Re: [R] xyplot color question

Using data frame CO2, built into R, try this: library(lattice) p <- function(..., col) { panel.xyplot(..., type = "p", col = col[1]) panel.xyplot(..., type = "smooth", col = col[2]) } xyplot(conc ~ uptake | Treatment, CO2, panel = p, col = 1:2) On Mon, Mar 2, 2009 at 9:13 AM, Pau

Re: [R] FW: partial residual plots

Dwayne, Sorry, I forgot about this: also look at ?ceres.plots, also in package car. >From Prof. Fox's documentation: "Ceres plots are a generalization of component+residual (partial residual) plots that are less prone to leakage of nonlinearity among the predictors." Regards, Mark. Culp, Dway

Re: [R] FW: partial residual plots

Dwayne Culp wrote: >> Is it possible to produce partial residual plots in R? See ?cr.plots in package car. This does "Component+Residual (Partial Residual) Plots." Regards, Mark. Culp, Dwayne wrote: > > I was advised someone might be able to help me with this. > > Very truly yours, > > D

Re: [R] Finding Lambda in Poisson distribution

Hi, > -Original Message- > From: r-help-boun...@r-project.org > [mailto:r-help-boun...@r-project.org] On Behalf Of Saeed Ahmadi > Sent: Monday, March 02, 2009 3:16 PM > To: r-help@r-project.org > Subject: [R] Finding Lambda in Poisson distribution > > > Hi, > > I have a dataset. First

[R] ways to put multiple graphs on single page (using ggplot2)

Hi, Here are three plots: library(ggplot2) data(diamonds) randind <- sample(nrow(diamonds),1000,replace=FALSE) dsmall <- diamonds[randind,] qplot(carat, data=dsmall, geom="histogram",binwidth=1) qplot(carat, data=dsmall, geom="histogram",binwidth=.1) qplot(carat, data=dsmall, geom="histogram",bin

[R] using par funtions in bargraph.CI()

I'm trying to create a bargraph of means with standard error bars using the function bargraph.CI (in the sciplot package). Like this: bargraph.CI(x.factor, response,data,xlab, ylab, par(family="serif"),font=11) However, an error message comes up when I try to use the par funtion. Does the char

Re: [R] Fwd: Converting R to Sweave (Rnw)

On 3/2/2009 10:12 AM, Duncan Murdoch wrote: On 3/2/2009 9:06 AM, Rainer M Krug wrote: Hi I am thinking about using Sweave more frequently, especially for documenting code. But the syntax is slightly awkward for me (<>= ... @), and I was thinking if there would be a way of importing the type of

Re: [R] TinnR & Philips Webcam

On Mon, 2 Mar 2009 12:48:16 -0500 (EST) wrote: WE> Now when I start up TinnR to use with R, it opens the webcam. If I WE> manually start up R, Tinn R will no longer synch with it. More detail would help to answer those questions properly. Which Tinn-R are you using, and how are you starting a

[R] initial gradient and vmmin not finite

Dear Rhelpers I have the problem with initial values, could you please tell me how to solve it? Thank you June > p = summary(maxLik(fr,start=c(0,0,0,1,0,-25,-0.2))) Error in maxRoutine(fn = logLik, grad = grad, hess = hess, start = start, : NA in the initial gradient > p = summary(maxL

[R] RWeka dataset running

Hi to all list members, I'm a newbie to R and will work on it from now. I'm preparing my thesis based on rpart and RWeka. I've managed to run already existing test data but couldn't import a data set from out of R. (By the way I'm an ubuntu user.) Thanks for your help Ahmet KOCYIGIT

[R] (no subject)

Greetings, I am using fGarch package to estimate and simulate GARCH models. What I would like to do is to perform Monte Carlo simulation. Unfortunately I cannot figure how to modify the code to achieve this. I use the following code to run a single simulation: spec=garchSpec(model=list(ar= 0.44

Re: [R] logistic regression model validation through bootstrapping

Hi, I was wondering whether this query was addressed on how to perform validation through boostrapping. I am currently trying to implement a boostrapping approach to validation but don't know where to start. Help please. Thank you and Regards, Vivienne Ozohili Risk Model Validation Manager Gro

[R] barplot with specific order of x axis labels

Hi all, I'd appreciate your help with this problem. I need to plot a barplot with the categories in a specific order. My data might be: hours Freq AN 10 MO 14 LU 30 I need the categories to be in the order: MO LU AN Is there some way to pass sort the dataset into this order so that

[R] FW: partial residual plots

I was advised someone might be able to help me with this. Very truly yours, Dwayne E. Culp, P.E., CFM Effective February 8, 2009, my contact information becomes: Dwayne E. Culp, P.E., CFM | JACOBS ENGINEERING GROUP : North American Infrastructure | Manager: Hydrology & Hydraulics Section

[R] xyplot color question

Hi, I am plotting scatterplots of horsepower by torque, conditional on brand (I'm just making up the variables for this example), and the goal is to see both the scatterplot points as well as the smoothed line. When I do the following code, I get the same color for the points and line, and would

[R] QQplot

Hi All, this might be an easy question but I do not know how to find R-squared for a QQplot when doing linear regression? -- View this message in context: http://www.nabble.com/QQplot-tp22288218p22288218.html Sent from the R help mailing list archive at Nabble.com. _

[R] TinnR & Philips Webcam

Hi All, I have a Philips Webcam, model SPC110NC. I had to download a diver in order for the webcam to work with Skype from the philips support page (www.philips.com/support). Now when I start up TinnR to use with R, it opens the webcam. If I manually start up R, Tinn R will no longer synch w

Re: [R] Newton-Raphson method and Quasi Newton

Check out: ?optim ?nlminb Ravi. Ravi Varadhan, Ph.D. Assistant Professor, Division of Geriatric Medicine and Gerontology School of Medicine Johns Hopkins University Ph. (410) 502-2619 email: rvarad...@jhmi.edu - Origina

[R] Goldbach partitions code

Folks, I put up a brief note describing my naive attempts to compute Goldbach partitions, starting with a brute-force approach and refining progressively. http://jostamon.blogspot.com/2009/02/goldbachs-comet.html I'd welcome your suggestions on improvements, alternatives, other optimisations

Re: [R] Bash script that uses an R command

stephen sefick wrote: > echo 'hal-device | grep battery.remaining_time | awk '{print$3/3600}' > | awk '{print int($1)}":"int(60*($1-int($1)))'` > > here is the final shell script is anyone is interested - this is > written and working in debian linux > it can't work -- you have unmatched quotes

Re: [R] Trouble - Installing R on RedHat el5

on 03/02/2009 11:07 AM steve_fried...@nps.gov wrote: > I am lucky and am now working with a new Redhat Linux 64-bit OS > > But I am getting the following error can anyone provide some graciously > needed assistance: > > > [r...@bluebird system-files]# rpm -i R-2.8.1-1.rh5.x86_64.rpm > warning: R

Re: [R] Need help extracting info from XML file using XML package

Duncan Temple Lang wrote: > >> not really: the xpath pattern '//coordinates' does say 'find all >> coordinates nodes searching from the root', but the root here is not the >> original root of the whole document, but each polygon node in turn. >> try: >> >> root = xmlInternalTreeParse(' >>

[R] Trouble - Installing R on RedHat el5

I am lucky and am now working with a new Redhat Linux 64-bit OS But I am getting the following error can anyone provide some graciously needed assistance: [r...@bluebird system-files]# rpm -i R-2.8.1-1.rh5.x86_64.rpm warning: R-2.8.1-1.rh5.x86_64.rpm: Header V3 DSA signature: NOKEY, key ID 99b6

Re: [R] Help with Wilcoxon Test

What do you get with: str(SampA) str(SampV) The error message suggests you are not giving it numeric vectors. -- David Winsemius On Mar 2, 2009, at 10:59 AM, Amit Patel wrote: Hi I have 2 sets of data that I want to do a Wilcoxon test on. They are of the same dimension. One has 4 zero va

Re: [R] comment on this book "A Handbook of Statistical Analyses Using R by Brian S. Everitt (Author), Torsten Hothorn (Author)"

Yes I think it is. I used it. On Mon, Mar 2, 2009 at 11:46 AM, choonhong ang wrote: > Is this book a good reference to learn R for statistical analysis ? > > A Handbook of Statistical Analyses Using R by Brian S. > Everitt

Re: [R] portable R editor

I feel emacs is portable enough for me. On 3/2/09, Werner Wernersen wrote: > > Hi, > > I have been dreaming about a complete R environment on my USB stick for a > long time. Now I finally want to realize it but what I am missing is a good, > portable editor for R which has tabs and syntax highlig

Re: [R] handle graph size in eps

See the help for argument 'paper' in ?postscript On Mon, 2 Mar 2009, Benoit Boulinguiez wrote: Sorry for the lack of information. I'm indeed under Windows. I indeed used the menu "save as" in the graph window. The matter with the eps obtained was the width of the graph which is lower than wha

[R] comment on this book "A Handbook of Statistical Analyses Using R by Brian S. Everitt (Author), Torsten Hothorn (Author)"

Is this book a good reference to learn R for statistical analysis ? A Handbook of Statistical Analyses Using R by Brian S. Everitt(Author), Torste

Re: [R] combining identify() and locator()

I am new to R and programming in general, so although your framework wasn't exactly what I needed, it was a tremendous help in getting what I needed. Here is what I ended up doing: idloc <- function(xy,tol=0.10,hitcol="red",mispch=19,miscol="red"){ if (type=="flag") {col="red"}else if

Re: [R] handle graph size in eps

Try adding paper = "special" to your postscript arguments. Clint BowmanINTERNET: cl...@ecy.wa.gov Air Dispersion Modeler INTERNET: cl...@math.utah.edu Air Quality Program VOICE: (360) 407-6815 Department of Ecology FAX:

Re: [R] handle graph size in eps

Sorry for the lack of information. I'm indeed under Windows. I indeed used the menu "save as" in the graph window. The matter with the eps obtained was the width of the graph which is lower than what I had on the screen or what I got when I saved it as a JPEG file. I tried the "postscript" comm

  1   2   >