[R] Strange behavior with poisosn and glm

Hi, I'm just learning about poison links for the glm function. One of the data sets I'm playing with has several of the variables as factors (i.e. month, group, etc.) When I call the glm function with a formula that has a factor variable, R automatically converts the variable to a series of

Re: [R] two questions for R beginners

On Tue, 2 Mar 2010 08:58:25 +1300 Peter Alspach peter.alsp...@plantandfood.co.nz wrote: This brings up another confusion for new users. Simply typing the object name at the command line gives just one view of the object (that provided by print()). Good point. Any good introduction to R

Re: [R] file reading /problems with encoding

Quoting Uwe Ligges lig...@statistik.tu-dortmund.de: R is not able to re-encode the file to the native encoding. But if you keep it in UTF-8, what is the problem to grep for the specific characters (as grep and friends support the argument useBytes these days)? The Problem with UTF-8 is that

[R] embedded nuls in 2.10 versus 2.11

I have been reading binary files, and parsing the output, for some time now. I have tried to develop a technique that is as robust as possible to all the strange things that appear in text fields, not to mention different global/regional encodings. I have no control over the data generated by

Re: [R] Strange behavior with poisosn and glm

On 02-Mar-10 08:02:27, Noah Silverman wrote: Hi, I'm just learning about poison links for the glm function. One of the data sets I'm playing with has several of the variables as factors (i.e. month, group, etc.) When I call the glm function with a formula that has a factor variable, R

Re: [R] Simple Linear Autoregressive Model with R Language

Emil Davtyan wrote: Hello - I need to do simple linear autoregressive model with R software for my thesis. I looked into all your documentation and I am not able to find anything too helpful. Can someone help me with the codes? Thanks Emil [[alternative HTML version deleted]]

Re: [R] Simple Linear Autoregressive Model with R Language

On Mon, 1 Mar 2010, Emil Davtyan wrote: Hello - I need to do simple linear autoregressive model with R software for my thesis. I looked into all your documentation and I am not able to find anything too helpful. Can someone help me with the codes? By all documentation you mean that you

Re: [R] the predict.lda function

On Mon, 2010-03-01 at 16:55 -0500, Diana Connett wrote: Hello. I just downloaded R onto a new computer, and after entering library(MASS), I still get the message Error: could not find function predict.lda when I try to use the predict.lda function (even just predict.lda()) Can anyone help

Re: [R] Changepoints estimation in a data series

On Mon, 1 Mar 2010, FMH wrote: Dear All, I'm trying to find changepoints in a data series  which only consist of  11 measurements of altitude(x) and temperature(y), respectively, in which the data are as followed: y = 16.3, 16.2, 16.1, 15.6, 14.2, 10, 8.2,  8.0, 7.5, 7.3, 7.2 x = 1,  2,  5,

Re: [R] capturing errors in Sweave

Thanks, Berwin. That works just great! --sundar On Tue, Mar 2, 2010 at 12:57 AM, Berwin A Turlach ber...@maths.uwa.edu.auwrote: G'day Sundar, On Mon, 1 Mar 2010 23:46:55 -0800 Sundar Dorai-Raj sdorai...@gmail.com wrote: Thanks for the input, but I don't want try in the Sweave output. I

Re: [R] repeated measures anova, car package

Hello John, As you said, I could also take a means model and test linear hypothesis for the desired effects - would this also be the case for the repeated measure i did in the first place. I copied the model from the car model where you first call: modx-lm(cbind(div_h, div_l) ~ site, divrep)

Re: [R] Thougt I understood factors but??

On Mon, 1 Mar 2010 14:23:04 -0500 Liaw, Andy andy_l...@merck.com wrote: Indeed this is one of the (few, I believe) traps of R, Oh, no; there are many more: http://www.burns-stat.com/pages/Tutor/R_inferno.pdf :-) -- Karl Ove Hufthammer __

Re: [R] why a text editor?

On Mon, 01 Mar 2010 16:26:37 - (GMT) ted.hard...@manchester.ac.uk ted.hard...@manchester.ac.uk wrote: In vim (to which I'm wedded for life) it will pick up matching (), {} and []. You can also easily move between matching delimeter by typing '%'. A similar feature should be available in

Re: [R] two questions for R beginners

On Mon, 01 Mar 2010 10:00:07 -0500 Duncan Murdoch murd...@stats.uwo.ca wrote: Suppose X is a dataframe or a matrix. What would you expect to get from X[1]? What about as.vector(X), or as.numeric(X)? All this of course depends on type of object one is speaking of. There are plenty of

[R] lm.influence on glm objects

Dear R users Today I discovered that function lm.influence() stops when applied to glm objects with the following error message Error in if (NROW(e) != n) stop(non-NA residual length does not match cases used in fitting) : argument is of length zero After inspecting lm.influence.R (both into

Re: [R] Strange behavior with poisosn and glm

Ted, Brilliant explanation (as usual) I'm back in school, just starting on a post-graduate degree in stats so the help is really appreciated. Now, I have a slightly trickier question about the same model. I've seen more than one way to get values out of the glm model. i.e. If we're

Re: [R] Reading sas7bdat files directly

The dsread output is little-endian, as that's the native format for floats on the Wintel platform. The byte order should stay the same if converting directly to a float, using a data structure like (C/C++): union { char bytes[8]; double value; } If reading the values with a SAS HEX

Re: [R] capturing errors in Sweave

G'day Sundar, On Mon, 1 Mar 2010 23:46:55 -0800 Sundar Dorai-Raj sdorai...@gmail.com wrote: Thanks for the input, but I don't want try in the Sweave output. I want the output to look just like it does in the console, as if an uncaptured error really did occur. I don't think that you will get

Re: [R] Install R 2.10.1 on Windows XP Errors

From CRAN: 2.8 What's the best way to upgrade? That's a matter of taste. For most people the best thing to do is to uninstall R (see the previous Q), install the new version, copy any installed packages to the library folder in the new installation, run update.packages(checkBuilt=TRUE,

[R] Problem with package fpc

I am trying to load package fpc in order to use the 'plotcluster' function however everytime I attempt to do so I get the following warning message: library(fpc) Loading required package: MASS Error: package 'MASS' could not be loaded In addition: Warning messages: 1:

Re: [R] two questions for R beginners

On Mon, Mar 1, 2010 at 11:49 PM, Liviu Andronic landronim...@gmail.com wrote: On 3/1/10, Keo Ormsby keo.orms...@gmail.com wrote:  Perhaps my biggest problem was that I couldn't (and still haven't) seen *absolute beginners* documents. there was once a link posted on r-sig-teaching that would

Re: [R] bwplot with pch = |

On Mon, Mar 1, 2010 at 10:49 AM, Duncan Mackay mac...@northnet.com.au wrote: Dear All Below is a toy example of a modified standard bwplot. require(lattice) DF - data.frame(site = rep(1:5, each = 20),           height = rnorm(100)) bwplot(site ~ height,DF, pch = |, par.settings =

[R] question to define a matrix with some vectors with different lengths

Hi, I have some vector v1,v2,...,vk, with different lengths. I want to consider these vectors as a matrix with k rows. Can you please guide me how I can do it? Regards khazaei __ R-help@r-project.org mailing list

Re: [R] two questions for R beginner

What were your biggest misconceptions or stumbling blocks to getting up and running with R? Easy. I terms of materials I have been unable to find good books that introduce users to R from the perspective of someone familiar only with packages like SPSS or STATA, or not familiar with

Re: [R] two questions for R beginners

Please take what follows not as an ad hominem statement, but rather as an attempt to improve what is already an excellent program, that has been built as a result of many, many hours of dedicated work by many, many unpaid, unsung volunteers. It troubles me a bit that when a confusing aspect of

[R] code for empirical copula

Hi,   I hope somebody can give me an idea where can I can find the code for empirical copula. I have a bivariate data.   Thank you so much for your help. [[alternative HTML version deleted]] __ R-help@r-project.org mailing list

Re: [R] two questions for R beginner

Brandon Zicha wrote: What were your biggest misconceptions or stumbling blocks to getting up and running with R? Easy. I terms of materials I have been unable to find good books that introduce users to R from the perspective of someone familiar only with packages like SPSS or STATA, or not

[R] Random real numbers

Hi, How could i generate random real numbers between 0 en 2*pi? Thanks, Frederik [[alternative HTML version deleted]] __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide

[R] Problem With Pasting (Mac OS X)

Greetings, fellow travelers. When I paste a series of commands into R, they execute serially, and when I go back through commands by hitting the up key, they show up as a block, rather than as individual lines. Is there any way to change this behavior? I'm running the 32-bit build of R 2.10.1 on

Re: [R] two questions for R beginner

On Tue, 2 Mar 2010 12:31:45 +0100 Brandon Zicha brandon.zi...@ua.ac.be wrote: Easy. I terms of materials I have been unable to find good books that introduce users to R from the perspective of someone familiar only with packages like SPSS or STATA, Have you read these books: R for SAS

Re: [R] question to define a matrix with some vectors with different lengths

On Tue, 2 Mar 2010 12:29:28 +0100 (CET) khaz...@ceremade.dauphine.fr khaz...@ceremade.dauphine.fr wrote: I have some vector v1,v2,...,vk, with different lengths. I want to consider these vectors as a matrix with k rows. Can you please guide me how I can do it? What do you want to do with the

Re: [R] two questions for R beginners

John Sorkin wrote: Please take what follows not as an ad hominem statement, but rather as an attempt to improve what is already an excellent program, that has been built as a result of many, many hours of dedicated work by many, many unpaid, unsung volunteers. It troubles me a bit that when

Re: [R] two questions for R beginner

Hi Brandon,   I just read this book, which I am sure you will be interested in: http://www.amazon.com/SAS-SPSS-Users-Statistics-Computing/dp/0387094172 Cheers!! Albert-Jan ~~ In the face of ambiguity, refuse the temptation to

[R] simple data transformation question

Hi all, I have a (hopefully) simple newbie-level question.   # I have data like this: dtf - data.frame(read.table(textConnection(var  value   company  9887.1   company  91117.0   blaah  91.1   etc  11   etc  97111), header=TRUE))   # I would like to have output like this (the index number may

Re: [R] two questions for R beginners

On Tue, Mar 2, 2010 at 7:27 AM, Duncan Murdoch murd...@stats.uwo.ca wrote: John Sorkin wrote: Please take what follows not as an ad hominem statement, but rather as an attempt to improve what is already an excellent program, that has been built as a result of many, many hours of dedicated

Re: [R] Random real numbers

On Tue, 2 Mar 2010 11:51:39 +0100 frederik vanhaelst frederik.vanhae...@gmail.com wrote: How could i generate random real numbers between 0 en 2*pi? Ten such numbers from the uniform distribution: 2*pi*runif(10) -- Karl Ove Hufthammer __

Re: [R] Random real numbers

runif(20,0,2*pi) [1] 1.29417642 1.10933879 4.31669186 2.41339484 4.83705630 3.12713657 4.50893007 6.23232980 2.38783146 4.88483239 5.87292617 [12] 1.33293077 4.09458703 0.7593 1.67899698 2.42602639 0.08413394 2.40261439 5.46442874 2.13847582 On Tue, Mar 2, 2010 at 5:51 AM, frederik

[R] Row-wisely converting a data frame into a list

Hello, is there an elegant way, how I can convert each row of a data frame into distinct elements of a list? In essence, what I'm looking for is something like rows.to.lists - function( df ) { ll - NULL for( i in 1:nrow(df) ) ll - append( ll, list(df[i,]) )

Re: [R] question to define a matrix with some vectors with different lengths

This is what you would use a 'list' for. On Tue, Mar 2, 2010 at 6:29 AM, khaz...@ceremade.dauphine.fr wrote: Hi, I have some vector v1,v2,...,vk, with different  lengths. I want to consider these vectors as a matrix with k rows. Can you please guide me how I can do it? Regards khazaei

Re: [R] Random real numbers

From what distribution? If the uniform, runif(100,0,2*pi) If another, install package Runuran, and do this ?Runuran Vignette(Runuran) HTH Dr. Rubén Roa-Ureta AZTI - Tecnalia / Marine Research Unit

[R] Double Colors in Main

Dear All, Consider the following trivial code snippet rm(list=ls()) name_vec - c(color1, color2) pdf(test_color.pdf) plot(seq(5), seq(5), main=paste(name_vec[1], and ,name_vec[2], sep=)) dev.off() What I would like to achieve is rather simple to explain, but it is giving me a headache: how

Re: [R] Row-wisely converting a data frame into a list

as.data.frame(t(df)) For example x - as.data.frame(t(mtcars)) typeof(x) [1] list -Original Message- From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org] On Behalf Of Sebastian Bauer Sent: Tuesday, March 02, 2010 8:12 AM To: r-help@r-project.org Subject: [R]

Re: [R] Double Colors in Main

See this thread : http://finzi.psych.upenn.edu/Rhelp10/2009-January/185693.html On 03/02/2010 02:18 PM, Lorenzo Isella wrote: Dear All, Consider the following trivial code snippet rm(list=ls()) name_vec - c(color1, color2) pdf(test_color.pdf) plot(seq(5), seq(5), main=paste(name_vec[1], and

Re: [R] Row-wisely converting a data frame into a list

Hi! On 03/02/2010 02:22 PM, Nutter, Benjamin wrote: as.data.frame(t(df)) For example x- as.data.frame(t(mtcars)) typeof(x) [1] list Thanks for the quick reply! I would never have guessed that as.data.frame() works that way! BTW This one seems also to do the trick: rows.to.list -

Re: [R] ANOVA Types and Regression models: the same?

If memory serves, Bill Venables said in the paper cited several times here, that there's only one type of sums of squares. So there's only one type of ANOVA (if I understand what you mean by ANOVA). Just forget about the different types of tests, and simply ask yourself this (hopefully simple

Re: [R] Stack type

Here is an example using proto based on converting Duncan's example: library(proto) Stack - proto(new = function(.) proto(Stack, stack = NULL, push = function(., el) .$stack - c(list(el), .$stack), pop = function(.) { stopifnot(length(.$stack) 0) out -

Re: [R] simple data transformation question

Try this: reshape(cbind(id = as.numeric(dtf$var), dtf, time = with(dtf, ave(value, var, FUN = seq))), timevar=time, direction=wide) Or: xtabs(value ~ var + ave(value, var, FUN = seq), data = dtf) On Tue, Mar 2, 2010 at 9:40 AM, Albert-Jan Roskam fo...@yahoo.com wrote: Hi all, I have a

Re: [R] Row-wisely converting a data frame into a list

On Mar 2, 2010, at 8:11 AM, Sebastian Bauer wrote: Hello, is there an elegant way, how I can convert each row of a data frame into distinct elements of a list? split(dfrm, rownames(dfrm)) In essence, what I'm looking for is something like rows.to.lists - function( df ) { ll -

Re: [R] ANOVA Types and Regression models: the same?

Ravi Kulkarni wrote: Hello, I think I am beginning to understand what is involved in the so-called Type-I, II, ... ANOVAS (thanks to all the replies I got for yesterday's post). I have a question that will help me (and others?) understand it better (or remove a misunderstanding): I know that

Re: [R] two questions for R beginner

Brandon Zicha wrote: Hey Paul, Hey Brandon, (adding R-help in the cc) I agree with you that the documentation of R could be better, especially with more examples in code showing not only the common cases, but also more esoteric cases. It would be great if everyone invested a lot of time to

Re: [R] ANOVA Types and Regression models: the same?

Sorry there were 2 typos in my note: John Fox's posting provided a lot of information. Briefly, the Types refer to whether effects are adjusted for all other effects in the model (Types II, III, IV) or no (Type I or sequential tests only adjust for no - not EARLIER terms in the model).

Re: [R] simple data transformation question

Hi Henrique,   *Thank you!* The reshape code does precisely what I want. Cheers!! Albert-Jan ~~ In the face of ambiguity, refuse the temptation to guess. ~~

[R] nu-SVM crashes in e1071

Hello ! I`m using SVMs for multi-class classification problems. Therefore I`m using the svm() function in the package e1071. If I use svm(...type=C-classification) everything works fine. But if I want to use nu-SVM with svm(..., type=nu-classification, nu=0.5) R crashes immediately. No error

Re: [R] Strange behavior with poisosn and glm

On Tue, 2010-03-02 at 00:58 -0800, Noah Silverman wrote: Ted, Brilliant explanation (as usual) I'm back in school, just starting on a post-graduate degree in stats so the help is really appreciated. Now, I have a slightly trickier question about the same model. I've seen more than

Re: [R] nu-SVM crashes in e1071

On 02.03.2010 15:41, Häring, Tim (LWF) wrote: Hello ! I`m using SVMs for multi-class classification problems. Therefore I`m using the svm() function in the package e1071. If I use svm(...type=C-classification) everything works fine. But if I want to use nu-SVM with svm(...,

[R] Fwd: Building R packages in Windows 7

-- Forwarded message -- From: Duncan Murdoch murd...@stats.uwo.ca Date: 25 February 2010 22:35 Subject: Re: [R] Building R packages in Windows 7 To: RHelp r-help@r-project.org, Eric Ferreira ericbferre...@gmail.com On 25/02/2010 7:56 PM, Eric Ferreira wrote: Thank you, Sir, but

Re: [R] nu-SVM crashes in e1071

I`m using SVMs for multi-class classification problems. Therefore I`m using the svm() function in the package e1071. If I use svm(...type=C-classification) everything works fine. But if I want to use nu-SVM with svm(..., type=nu-classification, nu=0.5) R crashes immediately. No error

Re: [R] Random real numbers

frederik vanhaelst wrote: Hi, How could i generate random real numbers between 0 en 2*pi? Thanks, Frederik [[alternative HTML version deleted]] __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do

[R] Reading data file with both fixed and tab-delimited fields

Hello R wizards, What is the best way to read a data file containing both fixed-width and tab-delimited files? (More detail follows.) _*Details:*_ The U.S. Bureau of Labor Statistics provides local area unemployment statistics at ftp://ftp.bls.gov/pub/time.series/la/, and the data are

Re: [R] plotting a subset of a time series

The indexing in xts is very nice; it may do what you want. library(xts) x.xts - as.xts(x) plot(x.xts) plot(x.xts['2005::2006-10']) HTH, David Reiner -Original Message- From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org] On Behalf Of Erin Hodgess Sent: Tuesday, March

Re: [R] nu-SVM crashes in e1071

On 02.03.2010 17:05, Häring, Tim (LWF) wrote: I`m using SVMs for multi-class classification problems. Therefore I`m using the svm() function in the package e1071. If I use svm(...type=C-classification) everything works fine. But if I want to use nu-SVM with svm(...,

Re: [R] plotting a subset of a time series

Try this: plot(window(x, end = c(2006, 10))) On Tue, Mar 2, 2010 at 2:19 AM, Erin Hodgess erinm.hodg...@gmail.com wrote: Dear R People: I have the following time series and plot: x - ts(rnorm(50),start=2005,freq=12) plot(x) which works fine. I would like to plot a subset of that time

Re: [R] Random real numbers

On Tue, Mar 2, 2010 at 4:05 PM, Paul Hiemstra p.hiems...@geo.uu.nl wrote: frederik vanhaelst wrote: Hi, How could i generate random real numbers between 0 en 2*pi? Thanks, Googeling for R generate random number gave this as a second hit (on my machine):

Re: [R] kriging with geoR package

hi all. If someone have the same problems this is the answer: -vertical legend : legend.krige(x.leg=c(X,X),y.leg=c(X,X),kr$pred,vert=TRUE,col=gray(seq(.7,0,l=10))) -sample's positions on the map: ###coords.dat=table$coords### like in

Re: [R] capturing errors in Sweave

G'day Sundar, On Tue, 2 Mar 2010 01:03:54 -0800 Sundar Dorai-Raj sdorai...@gmail.com wrote: Thanks, Berwin. That works just great! You are welcome. I noticed by now that cat(tmp) is sufficient; the tmp[1] in cat(tmp[1]) was a left over from earlier attempts to get the output to look correct.

Re: [R] nu-SVM crashes in e1071

Hi, On Tue, Mar 2, 2010 at 11:05 AM, Häring, Tim (LWF) tim.haer...@lwf.bayern.de wrote: I`m using SVMs for multi-class classification problems. Therefore I`m using the svm() function in the package e1071. If I use svm(...type=C-classification) everything works fine. But if I want to use

[R] problem with choose.files

I have recently upgraded to R 2.10.1 on Windows XP and am using scripts that I've used in previous versions successfully. I'm having a problem with choose.files. The lines read: fura_scan_file-choose.files(caption=Select log file (*.log) for fura-2 scans) PI_scan_file-choose.files(caption=Select

Re: [R] Reading data file with both fixed and tab-delimited fields

I tried to shoehorn the read.* functions and match both the fixed width and the variable width fields in the data but it doesn't seem evident to me. (read.fwf reads fixed width data properly but the rest of the fields must be processed separately -- maybe insert NULL stubs in the remaining fields

[R] add a header to a forest plot (metafor)

Dear R-community, I'm currently trying to assemble a forest plot using the forest function from package metaphor. Works well. Even the regular main-argument works for adding a title to the graph. However, I would like to add one top row which explains the nature of the columns. Very much

Re: [R] Reading data file with both fixed and tab-delimited fields

Ah, I should have mentioned this. Personally I work on Macs (Leopard) and PC's (XP Pro and XP Pro x64). Even though the PC's do have Cygwin, I'm trying to make this code portable. So I want to avoid such things as sed, perl, etc. I want to do this in R, even if processing is a bit slower.

[R] adding row ID numbers by group

Hello R community, I am hoping for some help with the following problem. I have a data frame containing various groups. These groups are identified by a grouping variable. I would like to add a sequential ID number to each group to later sort these individuals within each group by this ID

Re: [R] two questions for R beginners

On 02/03/2010 11:53 AM, William Dunlap wrote: -Original Message- From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org] On Behalf Of John Sorkin Sent: Tuesday, March 02, 2010 3:46 AM To: Karl Ove Hufthammer; r-h...@stat.math.ethz.ch Subject: Re: [R] two

Re: [R] adding row ID numbers by group

Try this: data$ID - with(data, ave(group, group, FUN = seq)) On Tue, Mar 2, 2010 at 2:53 PM, Alexander Schwall alexander.schw...@gmail.com wrote: Hello R community, I am hoping for some help with the following problem. I have a data frame containing various groups. These groups are

Re: [R] two questions for R beginners

William, I agree that changing syntax can lead to problems. I don't, however think extending the language will break existing code. Providing a common syntax for accessing matrices and dataframes will not change the way things have been done to date, but rather how things will be done in the

Re: [R] adding row ID numbers by group

Like this? group- c(1,1,1,1,2,2,2,2,2,3,3,3,3,3,3) var2- c(1,2,3,4,5,6,7,8,9,10,11,12,13,14,15) data-data.frame(group, var2) data ddply(data,group,transform,ID=1:length(group))   Felipe D. Carrillo Supervisory Fishery Biologist Department of the Interior US Fish Wildlife Service California, USA

[R] how to import map data (maptools?) from a html set of 'coords'

Dear R users, I would like to draw map and import it in maptools/spatstat packages. The 'raw data' I have come from a web page (map.../map) and are basically a list of coordinates of a polygon. I would like to know how to import them in R; I checked the maptools packages, but all the examples

[R] Fwd: Building R packages in Windows 7

-- Forwarded message -- From: Duncan Murdoch murd...@stats.uwo.ca Date: 25 February 2010 22:30 Subject: Re: [R] Building R packages in Windows 7 To: Eric Ferreira ericbferre...@gmail.com On 25/02/2010 7:56 PM, Eric Ferreira wrote: Thank you, Sir, but how can I demand it to

[R] Creating matrix from long table in database (pivoting)

Hi all, I have a table in database that is very long and when simplified it has only two columns in it (id, text). id is the row, and text is the column. Technically the text is a term and and id is the document. If simplifying this and assuming there is only one occurrence of the term per the

Re: [R] two questions for R beginners

Liviu Andronic escribió: On Mon, Mar 1, 2010 at 11:49 PM, Liviu Andronic landronim...@gmail.com wrote: On 3/1/10, Keo Ormsby keo.orms...@gmail.com wrote: Perhaps my biggest problem was that I couldn't (and still haven't) seen *absolute beginners* documents. there was once a

Re: [R] adding row ID numbers by group

Thank you gentlemen, all three solutions are working and very insightful. Your help and time is very much appreciated. Alexander On Tue, Mar 2, 2010 at 1:08 PM, Felipe Carrillo mazatlanmex...@yahoo.comwrote: Like this? group- c(1,1,1,1,2,2,2,2,2,3,3,3,3,3,3) var2-

[R] plotting fitted lme values as a smooth line

I am trying to plot fitted lme values as a smooth line of a graph showing the exponential relationship between temperature and soil respiration. In the plot, the x-axis has temperature, and the y-axis has soil respiration. When I try to add a line showing temperature versus the fitted

[R] Compare parameter estimates of a nlsList object

Hello together, Is there a tool to test the statistical differences between parameter estimates of an nlsList fit, with more than two groups? I am able to complete the nlme function for two groups after getting starting paramaters in nlsList, as seen below. fit.nlme - nlme(rate ~

Re: [R] Creating matrix from long table in database (pivoting)

Try this: DF - read.table(textConnection(1 this 1 is 1 the 1 first 1 row 2 this 2 is 2 the 2 send 2 row)) reshape(DF, v.names = 'V2', idvar = 'V1', timevar = 'V2', direction = 'wide') On Tue, Mar 2, 2010 at 3:35 PM, Jan Hornych jh.horn...@gmail.com wrote: Hi all, I have a table in database

Re: [R] Creating matrix from long table in database (pivoting)

Jan - Here's one way: tbl = data.frame(id=c(1,1,1,1,1,2,2,2,2,2), text=c('this','is','the','first','row','this','is','the','second','row')) xtabs(~id+text,tbl) text id first is row second the this 1 1 1 1 0 11 2 0 1 1 1 11

Re: [R] Creating matrix from long table in database (pivoting)

Or better: reshape(cbind(DF, value = 1), v.names = 'value', idvar = 'V1', timevar = 'V2', direction = 'wide') On Tue, Mar 2, 2010 at 3:49 PM, Henrique Dallazuanna www...@gmail.com wrote: Try this: DF - read.table(textConnection(1 this 1 is 1 the 1 first 1 row 2 this 2 is 2 the 2 send

Re: [R] capturing errors in Sweave

What I ended up using was: cat(unclass(tmp)) --sundar On Tue, Mar 2, 2010 at 8:58 AM, Berwin A Turlach ber...@maths.uwa.edu.auwrote: G'day Sundar, On Tue, 2 Mar 2010 01:03:54 -0800 Sundar Dorai-Raj sdorai...@gmail.com wrote: Thanks, Berwin. That works just great! You are welcome. I

Re: [R] Gradient Boosting Trees with correlated predictors in gbm

On Mon, 01-Mar-2010 at 12:01PM -0500, Max Kuhn wrote: | In theory, the choice between two perfectly correlated predictors is | random. Therefore, the importance should be diluted by half. | However, this is implementation dependent. | | For example, run this: | | set.seed(1) | n - 100 | p

Re: [R] Gradient Boosting Trees with correlated predictors in gbm

In most implementations of boosting, and for that matter, single tree, the first variable wins when there are ties. In randomForest the variables are sampled, and thus not tested in the same order from one node to the next, thus the variables are more likely to share the glory. Best, Andy

Re: [R] Strange behavior with poisosn and glm

On 2/03/2010, at 9:02 PM, Noah Silverman wrote: Hi, I'm just learning about poison links for the glm function. One of the data sets I'm playing with has several of the variables as factors (i.e. month, group, etc.) When I call the glm function with a formula that has a factor

[R] distribution for random effects?

Hi R users, I am using the following model to analyze data from a factorial experiment (randomized complete block design with no replication within blocks): model - glm(survival ~ density * vegetation + (1|block), data=sal2005, family=binomial) Does R use a binomial distribution in this

[R] Three most useful R package

Hi R-fans, I would like put out a question to all R users on this list and hope it will create some feedback and discussion. 1) What are your 3 most useful R package? and 2) What R package do you still miss and why do you think it would make a useful addition? Pulling answers together for

Re: [R] Gradient Boosting Trees with correlated predictors in gbm

On Tue, Mar 2, 2010 at 2:43 PM, Liaw, Andy andy_l...@merck.com wrote: In most implementations of boosting, and for that matter, single tree, the first variable wins when there are ties. They must be in a union :-) What happens if there's a third? If they were P perfectly correlated

Re: [R] add a header to a forest plot (metafor)

Hi Sebastian, Here is an example showing a forest plot with some column headings: library(metafor) data(dat.bcg) dat - dat.bcg res - rma(ai=tpos, bi=tneg, ci=cpos, di=cneg, data=dat, measure=RR) windows(width=6.5,

Re: [R] problem with choose.files

Try this; specify where you want the second one to start: files.a - choose.files() # now change to the directory of the first file name to continue search files.b - choose.files(paste(dirname(files.a[1L]), *, sep='/')) On Tue, Mar 2, 2010 at 12:17 PM, Caleb Rounds caleb.rou...@gmail.com wrote:

Re: [R] variable substitution in for loops

Friends Seems I've run into another snag. More of the nitty-gritty r-details I don't understand. So, as I mentioned below, dataset[[var_sub]] seems to be understood well by the functions I previously used and I was able to run my loop successfully with the [[var_sub]] as a

[R] R / R+ Webminar *** R-PLUS Rocks: Interactive, Comprehensible and Highly Visual. March 12th @ 12PM ET (USA Time)

Welcome to R/ R-PLUS Webminar Series. R-PLUS 3.3 Rocks: Interactive, Comprehensible and Highly Visual. http://www.xlsolutions-corp.com/webminar.asp. March 12th @ 12PM ET (USA Time) Increase your productivity with R-PLUS 3.3 by attending the webminar and learning how to: 1. Interactively

[R] Creating a timeSeries Data Frame

Hello I have 2000 univariate timeSeries of about 20 observations each, as the following, I would like to store all of them in one object, sort of a data frame, and to be able to recall each by its column name, which by the way is the same as the first date. Do you know how can I do this. Thank you

[R] selection simulation

Hi, I was having trouble developing a looping function in a selection simulation that I'm trying to develop and I was hoping if someone could help. Basically, I have a matrix with a random generated numbers representing scores on the variables. The rows represent applicants and columns represent

[R] Binding a matrix to a matrix

Hello I have a 2x10x200 matrix and I would like to bind to it another 2x10 matrix in order to end up with an 2x10x2001 matrix, which command should i use in order to do this? Thank you Felipe Parra [[alternative HTML version deleted]] __

Re: [R] Binding a matrix to a matrix

If it is 3 dimensional then it is an array, not a matrix. The abind function in the abind package is probably what you want. -- Gregory (Greg) L. Snow Ph.D. Statistical Data Center Intermountain Healthcare greg.s...@imail.org 801.408.8111 -Original Message- From:

  1   2   >