Re: [R] pvals.fnc in lme4 and languageR

Daniela Reicheneder tz.agrar.tu-muenchen.de> writes: > I was using the function pvals.fnc from package 'languageR' until April. > I do not know which version. Yesterday I updated all my packages > and tried to run my loop again. Now I get the following error message: > > error in pvals.fnc(mm,

[R] An AIC model selection question

Dear R users, Assume I have three models with the following AIC values: model AIC df model1 -10 2 model2 -12 5 model3 -11 2 Obviously, model2 would be preferred, but it "wastes" 5 df compared to the other models. Would it be allowed to select model3 instead, simply becaus

Re: [R] Multiple hist(ograms) - One plot

Dear Michael, Try ggplot2. Use something like install.packages("ggplot2") library(ggplot2) recov <- 0:2 n <- 1000 all <- data.frame(bbED = rnorm(3 * n, mean = recov), recov = factor(rep(recov, n))) ggplot(data = all, aes(x = bbED)) + geom_histogram() + facet_grid(. ~ recov) ggplot(data = all, aes

Re: [R] question about working with workspaces

Duncan Murdoch wrote: Daniel Rabczenko wrote: Hello everybody, Two - I hope Simple questions about working with workspaces. Is there a way to force R to start in "clean" workspace / avoid "previously saved workspace restored"? Start with --no-restore. (Start with --help for the full li

Re: [R] An AIC model selection question

Christoph Scherber wrote: Dear R users, Assume I have three models with the following AIC values: modelAICdf model1 -102 model2-125 model3-112 Obviously, model2 would be preferred, but it "wastes" 5 df compared to the other models. Would it be allowed to se

Re: [R] Error: cannot allocate vector of size 117.3 Mb

ram basnet wrote: Dear R users, I am using RandomForest package. While using this package, i got "Error: cannot allocate vector of size 117.3 Mb" .message. I had this problem earlier too but could not manage. Is there any way to solve this problem or to increase vec

Re: [R] Bug or feature with finding a list element?

I think you will find that this follows the partial matching rules: * exact match * unique partial match * take evasive action if non-unique partial match Patrick Burns [EMAIL PROTECTED] +44 (0)20 8525 0696 http://www.burns-stat.com (home of S Poetry and "A Guide for the Unwilling S User") B Fo

Re: [R] lapply where each list object has multiple parts

The function in 'lapply' gets each component, so write a function that does what you want given a component of the list: lapply(BigList, function(x) mean(x$label1)) Patrick Burns [EMAIL PROTECTED] +44 (0)20 8525 0696 http://www.burns-stat.com (home of S Poetry and "A Guide for the Unwilling S Us

[R] Hexagon border in plotrix's math2D.plot

Hi there, I'm attempting to use color2D.matplot(x,cellcolors=cellcol, do.hex=TRUE) from the plotrix package, but do not want to have borders around the hexagons. I've tried various approaches such as fg="white", col="white", lty=0, lwd=0 and border=NULL/FALSE/"white", but none of them works.

[R] contingency table, several variables from dataframe

Hello R-Users! I need a little help to build up a contingency table out of several variables. A<-c("F","M","M","F","F","F","F","M","F","M","F","F") B<-c(0,0,0,0,0,0,1,1,1,1,0,1) C<-c(0,1,1,1,1,1,1,1,1,0,0,0) ABC<-as.data.frame(cbind(A,B,C)) ABC A B C 1 F 0 0 2 M 0 1 3 M 0 1 4

Re: [R] Lattice: don't draw unlabled tick marks

On Wed, Oct 01, 2008 at 01:19:04PM +, Dieter Menne wrote: > > I also would like to get rid of the unlabeled tick marks (top and > > right). How can I do that? > > > > Reminds me of a problem with a journal editor... See parameter tck under > scales > in xyplot documentation. Thanks! Aparent

Re: [R] Multiple hist(ograms) - One plot

Hello, If I use: p <- ggplot(dat, aes(x=bbContag, y=..density..)) + geom_histogram() and then: p + facet_grid(. ~ sc_recov) Its a little crazy because I have 48 different values in 'sc_recov'. Instead I want to select cases from 'sc_recov' and only use three at a time: I tried: > p + facet_gri

[R] R Training Courses

There are still places available on the following courses. For more detailed information and a registration form please contact [EMAIL PROTECTED], or visit our website at www.mango-solutions.com . R for Finance **

Re: [R] Multiple hist(ograms) - One plot

Dieter Menne: Thanks for the suggestion and link, it looks good. I think my trouble now is lack of basic R knowledge. Cheers, Michael For example, I tried: histogram(~bbContag | sc_recov, data=dat) This would work fine if I could select cases from sc_recov. How can I select more than one case?

[R] possible bug in function 'var' in R 2.7.2?

Dear R-Help, I have used R2.6.0 until I recently installed also R2.7.2 (see details below) In R 2.6.0, the following script using the function 'var' (cor(stats)): x.test <- c(NA, NA, NA, NA) var(x.test, na.rm=T) gives the output: NA In R2.7.2 the output of the same

[R] error message documentation: "Error: no function to return from, jumping to top level"

Dear R Users, Are error messages in R documented somewhere ? I am getting the following error message and would like to track down what it means and which specific routine triggers it: Error: no function to return from, jumping to top level Once I figure out which snippet of code triggers it,

Re: [R] error message documentation: "Error: no function to return from, jumping to top level"

I think you wrote a function, but pass only snippets of code from your script to your R-console. The part which return the error is the following: return() return can only be used to return from a function and not as a command in the R-console. Bart TU wrote: > > Dear R Users, > > Are err

Re: [R] Multiple hist(ograms) - One plot

Michael Just gmail.com> writes: > > sc_recov.21 <- dat[dat$sc_recov=21,] Be careful: it should be "==", not "=" when want a comparison. df = data.frame(grp = letters[1:5],val=rnorm(100)) dfab = df[df$grp=="a" | df$grp=="b",] # The simple way with "or" dfcde = df[df$grp %in% c("c","d","e"),] #

Re: [R] Multiple hist(ograms) - One plot

Michael, Use %in% to select multiple cases: dat[dat$sc_recov %in% c(21, 31, 41), ] ir. Thierry Onkelinx Instituut voor natuur- en bosonderzoek / Research Institute for Nature and Forest Cel biometrie, methodologie en

[R] Problem with read.table()

Hello everybody, I'm a new user and I'm trying to use read.table with a tab delimited file but the system tells me (i) that there are more columns than column names, and (ii) that headers and col.names are of different lengths. I have already checked my variable names and there's no blank

[R] RBloomberg :: blpConnect crash.

Hi there, How are you? I'm having some issues with the following set up: WINDOWS XP R-2.7.2 RBloomberg-0.1-10 RDCOMClient-0.92-0 zoo-1.5-4 chron-2.3-24 Bloomberg is installed in the machine. The following code crashes R: library(RBloomberg) .bbfields <- blpReadFields('C:\\path\\to\\blp\\API')

Re: [R] extracting diagonal from dist()

Stephen, I think this is what you want: > diag(as.matrix(d)[-1, ]) the '-1' above means that you ignore row one in the matrix. Given that it's a square matrix, I suppose you could also also do: > diag(as.matrix(d)[-1,-19]) the '-19' above mean you ignore the 19th column of hte matrix. Both give

Re: [R] Problem with read.table()

> I'm a new user and I'm trying to use read.table with a tab delimited > file but the system tells me (i) that there are more columns than > column names, and (ii) that headers and col.names are of different > lengths. I have already checked my variable names and there's no blank > spaces. Howe

Re: [R] contingency table, several variables from dataframe

First of all your construction of ABC leads to a structure with 3 factor variables due to the way cbind processes the input variables - which is not intended I think. You can do sth like ABC<-data.frame(A,B,C) aggregate(ABC[,2:3],by=list(A),sum) hth. Birgitle schrieb: Hello R-Users! I nee

Re: [R] lapply where each list object has multiple parts

Try this: rapply(BigList, mean)[c(1, 3)] On Wed, Oct 1, 2008 at 6:29 PM, eric lee <[EMAIL PROTECTED]> wrote: > Hi. I have a list where each object in the list has multiple parts. I'd > like to take the mean of just one part of each object. Is it possible to do > this with lapply? If not, can

Re: [R] contingency table, several variables from dataframe

Thanks for your answer. It is intended, that the variables are treated as class factor, because these are binary variables with, for example, the presence or the absence of a plant organ. As far as I understood, I have to treat them for other calculations as factor. Therefore I classified these v

Re: [R] function in R

you can see the different methods of mean by: methods(mean) Then you can type mean.default and you will see the complete code. Bart Alphonse Monkamg wrote: > > > > > > Dear ALL, > > Does anyone know how to get the complete code program for any build-in > function > in R, e.g. when I tape

Re: [R] function in R

On Thu, Oct 2, 2008 at 1:34 PM, Alphonse Monkamg <[EMAIL PROTECTED]> wrote: > > > > > Dear ALL, > > Does anyone know how to get the complete code program for any build-in > function > in R, e.g. when I tape mean in the R-console, I get the following: > > mean > > function (x, ...) > > UseMethod("

[R] Adding plane in a 3D scatterplot

I have drawn a 3D scatter plot : library(mnormt) library(scatterplot3d) dat = cbind(rmnorm(3, rep(0,2), diag(2)), 1:3) scatterplot3d(dat) Now I want to do 2 things : 1 : In the Z-axis (i.e. height), I want to see only numbers 1,2,3, etc NOT, 1,1.5,2,2.5. 2. I want to add two Horizo

Re: [R] Problem with read.table()

Joan-Josep Vallbé uab.cat> writes: > I'm a new user and I'm trying to use read.table with a tab delimited > file but the system tells me (i) that there are more columns than > column names, and (ii) that headers and col.names are of different > lengths. I have already checked my variable na

[R] nls with plinear and function on RHS

Dear R gurus, As part of finding initial values for a much more complicated fit I want to fit a function of the form y ~ a + bx + cx^d to fairly "noisy" data and have hit some problems. To demonstrate the specific R-related problem, here is an idealised data set, smaller and better fitting tha

Re: [R] contingency table, several variables from dataframe

Ok, then treat them as factors - but if they are really binary and coded 0 and 1, which kind of calculation would lead to different results for a "factor" instead of a numeric variable? Anyway, ABC<-as.data.frame(cbind(A,B,C)) aggregate(ABC[,2:3],by=list(A),FUN=function(x)sum(x=='1')) # '1' i

Re: [R] contingency table, several variables from dataframe

I think it makes a difference if I want to use a classification method like rpart () or if I use a modelling approach like glm(). Many thanks for the kind and fast help. I am still very untrained and it is difficult for me to create such codes. B. Eik Vettorazzi wrote: > > Ok, then treat them

[R] function in R

Dear ALL, Does anyone know how to get the complete code program for any build-in function in R, e.g. when I tape mean in the R-console, I get the following:  mean function (x, ...) UseMethod("mean") but I need the full mean function. Thank in advance, Alphonse.

Re: [R] Problem with read.table()

Hi, I've tried it with read.delim: ayuso <- read.delim("/Users/user/Rwork/ayusotriada.txt") and now the error message is "Error with make.names(col.names, unique = TRUE) : string multibyte 187 invalid" thanks! On 02/10/2008, at 12:21, Joan-Josep Vallbé wrote: Hello everybody, I'm a n

Re: [R] Problem with read.table()

Look at ?count.fields to diagnose this. On Thu, Oct 2, 2008 at 8:05 AM, Joan-Josep Vallbé <[EMAIL PROTECTED]> wrote: > Hi, I've tried it with read.delim: > > ayuso <- read.delim("/Users/user/Rwork/ayusotriada.txt") > > and now the error message is "Error with make.names(col.names, unique = > TRUE)

[R] constructing appropriate non-intercept formula

-BEGIN PGP SIGNED MESSAGE- Hash: SHA1 Trying to work out a model formula that will do what I want ... suppose I want to model y = b_i x + epsilon (i.e. a linear model with zero intercept and with slopes differing by groups), and I want to parameterize the slopes in the "usual" way of h

Re: [R] function in R

On Thu, Oct 2, 2008 at 6:34 AM, Alphonse Monkamg <[EMAIL PROTECTED]> wrote: > > > > > Dear ALL, > > Does anyone know how to get the complete code program for any build-in > function > in R, e.g. when I tape mean in the R-console, I get the following: > > mean > > function (x, ...) > > UseMethod("

Re: [R] Problem with read.table()

Hi, I've done ?count.fields and it helped me identify a blank space that I had not noticed (thanks!), but when I try to upload the table again with ayuso <- read.table("/Users/user/Rwork/ayusotriada.txt", header=TRUE, sep="\t", na.strings="NA", dec=".", strip.white=TRUE) or with read.deli

Re: [R] How do one suggest a patch for the Foreign package

In answer to your subject heading question: to suggest a patch, send an email to the package maintainer (as listed in the DESCRIPTION file). Your suggestion will get a better reception if you offer sample code/data to illustrate the problem, and a nicely written patch to the sources to fix it

Re: [R] possible bug in function 'var' in R 2.7.2?

risoe.dk> writes: > > Dear R-Help, > > I have used R2.6.0 until I recently installed also R2.7.2 (see details below) > > In R 2.6.0, the following script using the function 'var' (cor(stats)): > > x.test <- c(NA, NA, NA, NA) > > var(x.test, na.rm=T) > > gives the output: > > NA >

Re: [R] error message documentation: "Error: no function to return from, jumping to top level"

[EMAIL PROTECTED] wrote: Dear R Users, Are error messages in R documented somewhere ? I am getting the following error message and would like to track down what it means and which specific routine triggers it: Error: no function to return from, jumping to top level Once I figure out which s

Re: [R] An AIC model selection question

Uwe Ligges statistik.tu-dortmund.de> writes: > > > Christoph Scherber wrote: > > Dear R users, > > > > Assume I have three models with the following AIC values: > > > > modelAICdf > > model1 -102 > > model2-125 > > model3-112 > > > > Obviously, model2 would be

Re: [R] matrix with different type of column [SEC=UNCLASSIFIED]

?data.frame Lists and data.frames allow mixed components, matrices do not --- On Thu, 10/2/08, ZHU, Justin <[EMAIL PROTECTED]> wrote: > From: ZHU, Justin <[EMAIL PROTECTED]> > Subject: [R] matrix with different type of column [SEC=UNCLASSIFIED] > To: r-help@r-project.org > Received: Thursday, O

Re: [R] Aggregate issues with subset

I don't have R on this machine to check but are you sure that sc_recov still exists in s.all? Otherwise do str(s.all) and see what the object looks like. --- On Wed, 10/1/08, Michael Just <[EMAIL PROTECTED]> wrote: > From: Michael Just <[EMAIL PROTECTED]> > Subject: [R] Aggregate issues with

Re: [R] Problem with read.table()

Did you try skipping the header as one respondent mentioned? read.table(..., skip = 1, header = FALSE) On Thu, Oct 2, 2008 at 9:05 AM, Joan-Josep Vallbé <[EMAIL PROTECTED]> wrote: > Hi, I've done ?count.fields and it helped me identify a blank space that I > had not noticed (thanks!), but when I t

Re: [R] error message documentation: "Error: no function to return from, jumping to top level"

Duncan Murdoch wrote: > [EMAIL PROTECTED] wrote: >> Dear R Users, >> >> Are error messages in R documented somewhere ? I am getting the >> following error message and would like to track down what it means >> and which specific routine triggers it: >> >> Error: no function to return from, jumping t

Re: [R] error message documentation: "Error: no function to return from, jumping to top level"

Many thanks, that is very helpful Peter. Regards, Tolga Peter Dalgaard <[EMAIL PROTECTED]> 02/10/2008 14:24 To Duncan Murdoch <[EMAIL PROTECTED]> cc [EMAIL PROTECTED], r-help@r-project.org Subject Re: [R] error message documentation: "Error: no function to return from, jumping to top level"

[R] Help with aggregation

Hi All, I seem to be having a few troubles with aggregating data back onto the the dataframe, I want to take the max value of a user, and then apply this max value back against all id's that match (i.e a one to many matching) Can anyone offer any advice? is there a better way of doing this? Dummy

Re: [R] Problem with read.table()

Sorry, I didn't, and removing the first row it effectively works. Thank you very much for your help. On 02/10/2008, at 15:20, Gabor Grothendieck wrote: Did you try skipping the header as one respondent mentioned? read.table(..., skip = 1, header = FALSE) On Thu, Oct 2, 2008 at 9:05 AM, Joan-

Re: [R] Help with aggregation

See ?ave (in package stats) /hans -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf Of Michael Pearmain Sent: den 2 oktober 2008 15:28 To: r-help@r-project.org Subject: [R] Help with aggregation Hi All, I seem to be having a few troubles with aggregating data

[R] Time conversion

I'm trying to convert a variable that is imported from CSV into a datetime,I'm trying to use the strptime function but with no joy, can anyone offer any advise? i have a vector timestamp 07-07-2008-21:03:55 07-07-2008-19:52:55 07-07-2008-19:54:14 07-07-2008-15:24:01 07-07-2008-13:39:17 07-07-2008-

[R] non-mouse-sensitive pause between showing plots

Hello, I am using R version 2.7.2 on Windows. I would like to create a script that draws many plots and pauses before drawing a new plot, so that I could copy each plot consecutively to clipboard ("Copy as metafile" from context menu) or save it to a file ("Save as metafile" from context menu). I

Re: [R] Time conversion

> I'm trying to convert a variable that is imported from CSV into a datetime,I'm > trying to use the strptime function but with no joy, can anyone offer any > advise? > > i have a vector > timestamp > 07-07-2008-21:03:55 > 07-07-2008-19:52:55 > 07-07-2008-19:54:14 > 07-07-2008-15:24:01 > 07-07-2

Re: [R] non-mouse-sensitive pause between showing plots

> I am using R version 2.7.2 on Windows. I would like to create a script > that draws many plots and pauses before drawing a new plot, so that I > could copy each plot consecutively to clipboard ("Copy as metafile" from > context menu) or save it to a file ("Save as metafile" from context > menu).

Re: [R] nls with plinear and function on RHS

a) you should put values for Ca, Cb, Cc directly into the data list as data=list(Ca=1, b) you can simplify the call to # idealised data set aDF <- data.frame( x= c(1.80, 9.27, 6.48, 2.61, 9.86, 5.93, 6.76, 5.52, 6.06, 8.62), y= c(24.77, 2775.07, 895.15, 60.73, 3373.57, 677.82, 1021.92, 542.8

[R] specifying x-axis scale on random forest variable importance plot

i am new to R and using the random forest package. is there a way to specify the x-axis scale range for the variable importance plot? many thanks. -alison -- View this message in context: http://www.nabble.com/specifying-x-axis-scale-on-random-forest-variable-importance-plot-tp19780560p197805

Re: [R] fft inverse display help

My question is how does 'plot' know to implicilty call the plot.ts (in the case of the full "exact" spectrum being fed back into the inverse? Kevin Dieter Menne <[EMAIL PROTECTED]> wrote: > charter.net> writes: > > > > > ff <- complex(length(fs)) > > ff[9] <- fs[9] > > ff[5] <- fs[5] >

Re: [R] fft inverse display help

charter.net> writes: > > My question is how does 'plot' know to implicilty call the plot.ts (in the case of the full "exact" spectrum > being fed back into the inverse? So the title should be "How does the specific incarnation of object orientation in R work?" Try, for example, section 3.4 and

Re: [R] Running source() on a file in another directory

Thanks a lot! It seems to work fine now with what you suggested. I also tried another approach by running the following on the terminal R -f ~dir1/prog.R and it read in some lines in prog.R, but then the execution halted at one point for some reason I could not figure out. Any way I could get "R

[R] namespaces

I'd like to control my namespace thoroughly, separated by task. Is there a way, in R session, to introduce namespaces for tasks dynamically and switch them as needed? Or, is there a combination of load/save workspace steps which can facilitate this? Cheers, Alexy

[R] [solutions] "tapply versus by" in function with more than 1 arguments

Thanks to all. I summarized (in order to thank the list) the solutions to help future workers searching subjects like this at R help.   # Number of rows nr = 10 # Data set dataf = as.data.frame(matrix(c(rnorm(nr),rnorm(nr)*2,runif(nr),sort(c(1,1,2,2,3,3,sample(1:3,nr-6,replace=TRUE,ncol=4)) n

[R] SAS enterprise guide in R?

Hello, Our company has been looking at "SAS enterprise guide 4" (Insightful Miner is a similar product from Insightful, I believe) which seems to provide a nice graphical way of displaying/managing a process or project. It is also accessible to non-programmers. Does any of the R GUIs provide s

Re: [R] namespaces

Yes, I could prefix everything with taskN$, but that's boring! Can I now attach and detach tasks to have one prefix-less, since they look like dataframes? What does proto buy us? To respond to Uwe here as well, I want dynamic task switching and *also* an ability to load and save task from

Re: [R] namespaces

You could have an environment for each task and place your objects for each task in its environment. Note that this is getting close to object oriented ideas where each enviroment/task is an object containing your R variables and methods and the proto package can be used to facilitate that: libra

Re: [R] namespaces

What you mean is probably that you can have different workspaces. A Workspace is associated to the current R session and is saved/loaded in the directory from which R is started. Note that Namespaces (for packages) in R have a different meaning. Best, Uwe Ligges Alexy Khrabrov wrote: I'd like

Re: [R] namespaces

On 02/10/2008 11:03 AM, Alexy Khrabrov wrote: I'd like to control my namespace thoroughly, separated by task. Is there a way, in R session, to introduce namespaces for tasks dynamically and switch them as needed? Or, is there a combination of load/save workspace steps which can facilitate

Re: [R] namespaces

On Thu, Oct 2, 2008 at 11:26 AM, Alexy Khrabrov <[EMAIL PROTECTED]> wrote: > Yes, I could prefix everything with taskN$, but that's boring! Can I now > attach and detach tasks to have one prefix-less, since they look like > dataframes? Read the proto vignette and other info to get a better unders

Re: [R] nls with plinear and function on RHS

Wonderful, great! That solves my problem nicely (on to the next one!). Thanks a lot, Keith Jewell "Katharine Mullen" <[EMAIL PROTECTED]> wrote in message news:[EMAIL PROTECTED] > a) you should put values for Ca, Cb, Cc directly into the data list as > data=list(Ca=1, > > b) you can simplify

[R] combining bar and line plots with multiple axes

Hey list, I have a barplot with a line plot overlayed (on the 2nd y axis) and it looks fine except that the origin (0) on the 2nd y axis (side=4) is not lining up with the origin on the 1st y-axis (side=2, from the barplot). Is there a way to either shift the 2nd y-axis (line plot) down or readj

Re: [R] namespaces

Here is what I do: I have a different folder/directory for each of my different projects, then I start a new R session for each project that I am working on. If you start R with the folder/directory as the current working directory then the workspace for that project is loaded (and the script

Re: [R] combining bar and line plots with multiple axes

The updateusr function in the TeachingDemos package can be used to change the user coordinates of the current plot and may do what you want. -- Gregory (Greg) L. Snow Ph.D. Statistical Data Center Intermountain Healthcare [EMAIL PROTECTED] 801.408.8111 > -Original Message- > From: [EMAI

[R] spec.pgram help?

Hopefully this will not seem too ignorant of a question. I am having a hard time picking out the sources of the differences between: abs(fft(x))^2/length(x) and spec.pgram(x, taper=0, log="no", plot=FALSE) Also from the limited testing that I have done since the DC "frequency" is not return

Re: [R] namespaces

The trick about setting the windows title is quite nifty and I've just added that to my profile. Just wanted to add two additional ways to start up R in the current directory on Windows. They are provided in the batchfiles distribution. See home page at: http://batchfiles.googlecode.com 1. I

[R] StructTS hang?

I was wondering if I could get some advice on how best to diagnose a problem that I am having with StructTS. To illlustrate the problem I have generated a rigged time series: genseq <- function() { x <- numeric(4*365) s <- seq(as.Date("2005-01-01"), as.Date("2008-12-31"), by="month") ob <-

Re: [R] time segments intersection

help("%in%") may be of use. However I don't see any overlap in your example. --- On Wed, 10/1/08, [Ricardo Rodriguez] Your XEN ICT Team <[EMAIL PROTECTED]> wrote: > From: [Ricardo Rodriguez] Your XEN ICT Team <[EMAIL PROTECTED]> > Subject: [R] time segments intersection > To: r-help@r-project.o

[R] missing output in summary() and anova()

> y<-c(131.79, 131.79, 135.02, 135.55, 136.46, 136.83, 137.82, 138.00, 138.06, 138.04, 140.04, 142.44, 145.47, 144.34, 146.30, 147.54, 147.80) > x<-c(194.5, 194.3, 197.9, 198.4, 199.4, 199.9, 200.9, 201.1, 201.4, 201.3, 203.6, 204.6, 209.5,208.6, 210.7, 211.9, 212.2) > fitted.results<-lm(y~x) > sum

[R] back transforming output from negative binomial

Dear all, I used the glm.nb with the default values from the MASS package to run a negative binomial regression. Here is a simple example: set.seed(123) y <- c( rep(0, 30), rpois(70, lambda=2) ) smoke <- factor( sample( c("NO", "YES"), 100, replace=T ) ) height <- c( rnorm(30, me

Re: [R] fft inverse display help

Sorry, I must be looking at a different section but when I look at 3.4 in r-intro.pdf I see: 3.4 The class of an object All objects in R have a class, reported by the function class. For simple vectors this is just the mode, for example "numeric", "logical", "character" or "list", but "matrix",

[R] re moving a vector from a list.

Hi there I have 386 list of vectors.In that i want to delete 23 vector in the 53 list element. [[1]] "abc" "dfg" "dft" "err" [[2]] "elm" "erg"" trr" "rtt" similarly i have [[386]] now i want to delete or avoid trr in the [[2]]. i just tried something like this list[-list[[53]][3]] -- View

Re: [R] missing output in summary() and anova()

What version are you running? Under what OS. I don't seem to see the problem with 2.7.2 under Windows XP. --- On Thu, 10/2/08, Jordan Wong <[EMAIL PROTECTED]> wrote: > From: Jordan Wong <[EMAIL PROTECTED]> > Subject: [R] missing output in summary() and anova() > To: r-help@r-project.org > Rec

Re: [R] back transforming output from negative binomial

Adaikalavan Ramasamy imperial.ac.uk> writes: > > Dear all, > > I used the glm.nb with the default values from the MASS package to run a > negative binomial regression. Here is a simple example: [snip -- thanks for the example!] > The question now is how do I report the results, say, for height

Re: [R] back transforming output from negative binomial

Ben, fantastic. Thank you for confirming it. One more question. What do you call the back transformed variable? In my domain, people use something called the ratio of mean but I am not sure if it is the same. I am not what the "ratio" is between. Regards, Adai Ben Bolker wrote: Adaikalava

Re: [R] missing output in summary() and anova()

It is 2.7.2 under Windows XP SP3. I have tried 2.6.1,2.6.2 in another computer under Windows XP SP3 also and the problem still occur. 2008/10/3 John Kane <[EMAIL PROTECTED]> > > What version are you running? Under what OS. I don't seem to see the > problem with 2.7.2 under Windows XP. > > > ---

Re: [R] back transforming output from negative binomial

Adaikalavan Ramasamy wrote: > Ben, fantastic. Thank you for confirming it. > > One more question. What do you call the back transformed variable? > > In my domain, people use something called the ratio of mean but I am not > sure if it is the same. I am not what the "ratio" is between. It depe

Re: [R] time segments intersection

Try this: > library(zoo) > # create two time series to test with > z1 <- zooreg(0:10, Sys.Date()) > z2 <- lag(z1, 5) > > z1 2008-10-02 2008-10-03 2008-10-04 2008-10-05 2008-10-06 2008-10-07 2008-10-08 2008-10-09 2008-10-10 2008-10-11 2008-10-12 0 1 2 3

[R] 3 curves / 1 plot

I am a new user of R. Please does somebody knows how to plot 3 datasets (x1,a1),...,(xn,an), (x1,b1),...,(xn,bn), and (x1,c1),...,(xn,cn) on a single x,y plot, each of the three datasets being plotted with its own character pch() ? (three calls to plot() erase the two first datasets). Thank you ver

Re: [R] Adding plane in a 3D scatterplot

Megh Dal wrote: I have drawn a 3D scatter plot : library(mnormt) library(scatterplot3d) dat = cbind(rmnorm(3, rep(0,2), diag(2)), 1:3) scatterplot3d(dat) Now I want to do 2 things : 1 : In the Z-axis (i.e. height), I want to see only numbers 1,2,3, etc NOT, 1,1.5,2,2.5. 2. I wa

Re: [R] 3 curves / 1 plot

example(matplot) or ?points On Thu, Oct 2, 2008 at 11:43 AM, Michel PETITJEAN <[EMAIL PROTECTED]> wrote: > I am a new user of R. > Please does somebody knows how to plot 3 datasets > (x1,a1),...,(xn,an), (x1,b1),...,(xn,bn), and (x1,c1),...,(xn,cn) > on a single x,y plot, each of the three datas

Re: [R] 3 curves / 1 plot

imho better: ?xyplot using the groups argument. -- Bert -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf Of Kingsford Jones Sent: Thursday, October 02, 2008 11:17 AM To: Michel PETITJEAN Cc: r-help@r-project.org Subject: Re: [R] 3 curves / 1 plot example(

[R] Shadowed Plot

Hello R-users     Kindly assist me with the following plot problem in R- My data looks like this (shortened because I have about 5000 observations).   dat <- read.table(textConnection("Id  Time Y  1 0   194 1 5.22 179 1 5.97 190 2 1.61 265 2 2.1  234 2 16.4 300 2 29.5 345 3 0    212 3 0.36

Re: [R] 3 curves / 1 plot

Michel PETITJEAN-2 wrote: > > I am a new user of R. > Please does somebody knows how to plot 3 datasets > (x1,a1),...,(xn,an), (x1,b1),...,(xn,bn), and (x1,c1),...,(xn,cn) > on a single x,y plot, each of the three datasets being plotted with > its own character pch() ? > (three calls to plot()

[R] t.test() on a list

I have a list, myList, with each of its 9 components being a 15X15 matrix. I want to run a t-test across the list for each component in the matrix. For example, the first t-test is on myList[[1]][1, 1], myList[[2]][1, 1], ..., myList[[9]][1, 1]; and there are totally 15X15 t-tests. How can I run th

Re: [R] Adding plane in a 3D scatterplot

I got following error : > library(mnormt) > library(scatterplot3d) > dat = cbind(rmnorm(3, rep(0,2), diag(2)), 1:3) > s3d <- scatterplot3d(dat, lab.z=2, scale.y=0.7, angle=20) > s3d$plane3d(2,0,0, "solid", col="grey") Error in object$coefficients : $ operator is invalid for atomic vectors I am usi

[R] Load a program at the front end

I want to run a R program, prog.R, interactively. My question is, is there a way I can start prog.R on the shell terminal when invoking R, instead of using source() inside R? TIA, Gang __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/l

Re: [R] specifying x-axis scale on random forest variable importance plot

The default range is the range of the variable in the data. It wouldn't make much sense to extend beyond that. If you want to restrict the range, the varImpPlot() function does return the x-y pairs, so you can construct your own plot from that. Andy From: A Van Dyke > > i am new to R and usin

[R] Contour Plot Aspect Ratio

Hello there, I have a fairly simple request (I hope!) I have produced a filled contour plot like this: library(grDevices) library(gplots) library(plotrix) filled.contour(contour, axes=F, frame.plot=TRUE, color=terrain.colors, ylab= "Length Along Flume (m)", key.title = title(main="Velocity\n(m/

Re: [R] Bias in sample - Logistic Regression

Hi Shiva, Maybe you are interested in the following paper: Learning when Training Data are Costly: The Effect of Class Distribution on Tree Induction. G. Weiss and F. Provost. Journal of Artificial Intelligence Research 19 (2003) 315-354. For validating the models in those enviroments, Willia

Re: [R] t.test() on a list

I am sure you will get helpful answers. I am almost as sure that you shouldn't be doing this. I suggest you consult with your local statistician. -- Bert Gunter -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf Of Gang Chen Sent: Thursday, October 02, 2008 11

Re: [R] re moving a vector from a list.

Is this what you want to do - it was not clear from your example: > x [[1]] [1] 1 [[2]] [1] 2 [[3]] [1] 3 [[4]] [1] 4 > x <- x[-3] > x [[1]] [1] 1 [[2]] [1] 2 [[3]] [1] 4 > Or maybe this: > x <- list(c( 1,2 ,3,4),c(4,5,6,7),c(8,9)) > x [[1]] [1] 1 2 3 4 [[2]] [1] 4 5 6 7 [[3]] [1] 8 9

  1   2   >