[R] histogram scaling

Hi, I'm drawing two histograms in the same plot.However, my point of comparison is the difference in their x coordinates.But my problem is one histogram is much taller than the other.How can I get them both to the same height? -- Rajesh.J [[alternative HTML version deleted]] __

Re: [R] ggplot 2 semi-transparency error

I really don't know anything much about this but I tried "device=pdf" (no-quotes) and I got a pdf file with transparent plots. ?device.print gives help with some links to other devices hth rajesh j wrote: > > Hi, > I used the command ggplot as follows... > p<-ggplot(a,aes(x=a$V1,colour=a$V2,f

Re: [R] Plot(x,y)

Hi malcolm, sorry but your post continue with pour quality, because it is not reproducible, it is!? By the way, gavin.simp...@ucl.ac.uk wrote: . You really shouldn't call your data, 'data' as it will get very confusing when you start using modelling functions and using the data argument, or u

Re: [R] trouble with Rcmdr

Hello, On 8/18/09, Erin Hodgess wrote: > I used it on Saturday, and it worked fine. > > When I used it today, the window did not appear. > > This is version 2.9.1 on Ubuntu Jaunty Jackalope. > Did you make any system updates between Sat and today? Liviu ___

[R] Off-topic: New website using R

Dear list, Sorry for the off-topic post. I've recently built a website using R and Rapache to process survey and reaction time data. I'm using ggplot2 to generate graphical reports for participants, and thought some of you might be interested in taking a look. I'd appreciate any feedback on the gra

[R] ggplot 2 semi-transparency error

Hi, I used the command ggplot as follows... p<-ggplot(a,aes(x=a$V1,colour=a$V2,fill=a$V2)) p + geom_density(alpha = 0.2,xlim=c(-10,10),ylim=c(0,0.5)) when I say, dev.print(device=postscript,file="/alpha/dct.pdf") I get Warning message: In grid.Call.graphics("L_polygon", x$x, x$y, list(as.integer

Re: [R] Plot(x,y)

On Tue, Aug 18, 2009 at 6:59 PM, David Winsemius wrote: > ITYM: > > plot(data$V6, data$V5, col="red") Yup! My mistake. -- Michael Knudsen micknud...@gmail.com http://lifeofknudsen.blogspot.com/ __ R-help@r-project.org mailing list https://stat.ethz.c

Re: [R] Sweave output from print.summary.glm is too wide

Hi Peter Sorry I must have been reading too much into your email. I had thought that either latex or xtable packages handled summaries effectively. I have not really got into summaries in Sweave. I do not know whether SweaveListingUtils may help. I had a look at it but I did not have the time t

Re: [R] ggplot2 legend problem - SOLVED

Trawling the online and pdf manuals for density plots further I found a reference to position and tried it with histograms. It worked! So here is an example that gives me overlapping histograms, alpha transparency so they can both be seen in the area of overlap, and with properly labelled and col

[R] radial.plot() in plotrix module

Hello, Is it possible to plot differently colored points with the radial.plot() function in the plotrix package? For example, the code: radial.plot( c(.5,.6, .7), c(1,2,3), rp.type = "s", point.col=c("red", "green", "blue") ) produces a plot whose three points are all red, instead of one red, one

Re: [R] aggregating using a non-summary function?

On Tue, Aug 18, 2009 at 7:09 PM, Sherri Heck wrote: > Dear all- > > I have a data set that looks as follows (data are taken every 5 minutes): > >      LST   in     mph    Deg   DegF  DegF2    %    volts   Deg    mph2   w/m2 > 090501   -999    7.5  259.3   33.0   -999   41.5   -999   -999    9.1

Re: [R] ggplot2 legend problem

Still struggling with this. A further example using a slightly different organisation of the data. The factors "A" and "B" are included in the dataframe in an attempt to get ggplot to generate a legend automatically. x <- data.frame(value=rnorm(5000, mean=0), case="A") y <- data.frame(value=rnorm

Re: [R] How to use R to perform prediction based on history data

Thank you very much, it works cool! On Tue, Aug 18, 2009 at 9:38 AM, Gabor Grothendieck wrote: > Please read the last line to every message on r-help and note the > request to provide reproducible code. > > Anyways, try this: > > Lines <- "Field1 Field2 Field3     Field4  ResultField > 1231    CA

Re: [R] jpeg() does not work

I figured out already thanks for your help. Edward Chen -Original Message- From: Benilton Carvalho Sent: 8/19/2009 12:33:12 AM To: Edward Chen Cc: r-help@r-project.org Subject: Re: [R] jpeg() does not work what do you get with: capabilities()[["jpeg"]] ? b On Aug 18, 2009, at 7:27 PM,

Re: [R] Rcompression on Omegahat

Duncan, Please see below: ## > in.file [1] "~/Downloads/Rcompression_0.6-0.tar.gz" > install.packages( + pkgs = in.file, + repos = NULL, + contriburl = NULL + ) > > > library("Rcompression") Error in library("Rcompression") : 'Rcompression' is not a

Re: [R] Sweave output from print.summary.glm is too wide

Oops I forgot to mention ?xtable and ?latex in Frank Harrell's Hmisc or Design packages (I cannot remember) Duncan Hi Peter If all else fails you can make the space between the columns smaller by changing the column spacing in latex by: \setlength\tabcolsep{3pt} The default is 7pt. You m

[R] ggplot2 legend problem

I'm trying to overlay two histograms using transparency to enable viewing of multiple distributions on a single scale. So far ggplot2 seems to do what I want. However I'm having a problem generating the legend coloring appropriate to each distribution in the plot. Here is a test case to show my b

Re: [R] Sweave output from print.summary.glm is too wide

Hi Peter If all else fails you can make the space between the columns smaller by changing the column spacing in latex by: \setlength\tabcolsep{3pt} The default is 7pt. You may have to go a little higher if it looks cramped. I have gone down as far as 2pt using a small font size Also changi

[R] lmer with random slopes for 2 or more first-level factors?

I have data from a design in which items are completely nested within subjects. Subject is the only second-level factor, but I have multiple first-level factors (IVs). Say there are 2 such independent variables that I am interested in. What is the proper syntax to fit a mixed-effects mod

Re: [R] RGoogleDocs/RCurl through proxy

This might be a little difficult in the current version. But it is a good point. I have a few ideas about this could be done conveniently but with a reasonably clear semantic model. It may take a couple of days before I have time to implement that. D. Remko Duursma wrote: Dear list, I am t

[R] Problem with predict.coxph

We occasionally utilize the coxph function in the survival library to fit multinomial logit models. (The breslow method produces the same likelihood function as the multinomial logit). We then utilize the predict function to create summary results for various combinations of covariates. For exa

[R] RGoogleDocs/RCurl through proxy

Dear list, I am trying to use RGoogleDocs, but I am connecting through a proxy server. I know RCurl is used for the connection, which should be able to deal with proxies and such. How do I set this up for RCurl? And can I use those settings with RGoogleDocs as well? I have the name of the proxy

[R] moving color key in heatmap

Dear list, I have a question on moving color keys when side color bars are added to a heatmap. The R code below produces the color key in the upper left corner. Notice I have added side bars to the heatmap, but how could I move the color key below the image? -- library(gplots) data(mtcar

Re: [R] Rcompression on Omegahat

Hi Stanley Yes, I install it on my macs and need it as part of working with .docx, .xlsx files and downloading content via RCurl. So it gets some exercise and does work on Intel Macs. (Although my systems are always non-standard!) You will have to give us more information about what the pro

Re: [R] jpeg() does not work

what do you get with: capabilities()[["jpeg"]] ? b On Aug 18, 2009, at 7:27 PM, Edward Chen wrote: My script: raw = read.table("c:\\Documents and Settings\\protein\\My Documents\ \My Fragments\\file1.txt", header = TRUE) #Normal Average Signal VS Pool jpeg("Normal Average Signal vs Pool_t

[R] Sweave output from print.summary.glm is too wide

Hi all I am preparing a document using Sweave; a really useful tool. But I am having a problem. Consider this toy example Sweave file: \documentclass{article} \begin{document} <>= options(width=40) # Set width to 40 characters hide <- capture.output(example(glm)) # Create an example

[R] need help for t.test using "by"

I am trying to do ttest for each plate which has equal number of disease and controls. by searching this forum I found one posting suggesting OP to use "by(eo,PlateID, function(.sub) t.test(IL1Ra~Group,data=.sub))". when i modified this for my use I used to get the pvalues for each plate, recently

Re: [R] combined value in forestplot function

carol white wrote: Hi, How is the combined value (displayed by a diamond on the forest plot) of a meta-analysis illustrated by a forest plot calculated? The applied method and all calculation details are welcome. Thanks Carol This is not a very well posed question. First of all there are t

[R] aggregating using a non-summary function?

Dear all- I have a data set that looks as follows (data are taken every 5 minutes): LST in mphDeg DegF DegF2%volts Degmph2 w/m2 090501 -9997.5 259.3 33.0 -999 41.5 -999 -9999.10.2 0905010005 -999 11.0 261.2 33.1 -999 42.6

Re: [R] Simulation function

You need to store your results in a list and then access the information in the list to get the values: > boot<-function(a,b,c){ + media<-(a+b+c)/3 + var_app<-var(a) + list(media=media,var_app=var_app) + } > boot(2,4,10) $media [1] 5.33 $var_app [1] NA > > simul<-function(S){ +

[R] adding points to a wireframe

A quick question. I'm trying to plot a surface from a fitted model along with the original points, as in the following example: df<-data.frame(expand.grid(100*runif(1:100), 100*runif(1:100))) df$Var3<-rnorm(length(df$Var1), mean=df$Var1*df$Var2, sd=10) my.lm<-lm(Var3 ~ Var1*Var2, data=df) m

[R] jpeg() does not work

My script: raw = read.table("c:\\Documents and Settings\\protein\\My Documents\\My Fragments\\file1.txt", header = TRUE) #Normal Average Signal VS Pool jpeg("Normal Average Signal vs Pool_test.jpg", width = 800, height = 600) xx_norm = raw[,3] yy_norm = rowMeans(raw[,4:13]) plot(xx_norm, yy_norm)

Re: [R] Is there a construct for conditional comment?

On 18/08/2009 5:54 PM, Peng Yu wrote: Hi, In C++, I can use the following construct to choice either of the two blocks the comment but not both. Depending on whether the number after "#if" is zero or not, the commented block can be chose. I'm wondering if such thing is possible in R? No. The

Re: [R] second y-axis in lattice xyplot

At 07:50 19/08/2009, you wrote: Dear R-help, I recently discovered lattice xyplot, but have been struggling with plotting a second y-axis. I am able to do this with plot() but can't figure it out in xyplot. Any help would be appreciated. I have found help on adding a second Y-axis when the data

Re: [R] How to Dedup a Spatial Points Data Set

On 19/08/2009, at 9:10 AM, Paul Heinrich Dietrich wrote: I'm new to spatial analysis and am exploring numerous packages, mostly enjoying sp, gstat, and spBayes. Is there a function that allows the user to dedup a data set with multiple values at the same coordinates and replace those duplic

[R] Is there a construct for conditional comment?

Hi, In C++, I can use the following construct to choice either of the two blocks the comment but not both. Depending on whether the number after "#if" is zero or not, the commented block can be chose. I'm wondering if such thing is possible in R? #if 0 commented with 0 #else commented with 1 #end

[R] second y-axis in lattice xyplot

Dear R-help, I recently discovered lattice xyplot, but have been struggling with plotting a second y-axis. I am able to do this with plot() but can't figure it out in xyplot. Any help would be appreciated. I have found help on adding a second Y-axis when the data is all part of the same data fra

[R] kernel density estimation for univariate data using splancs

Hi, I previously received help in extract data from a shapefile and now my question is about kernel density estimation. My objective is to have 3 kernel density plots; 2 for the each set of cases and the 3rd is the difference in kernel densities between the 2 sets of cases. Previously, I used

Re: [R] Rcompression on Omegahat

Hi, Has anyone successfully installed Rcompression on Intel-based Mac? I used the command: install.packages( pkgs = in.file, repos = NULL, contriburl = NULL ) where in.file is the path to Rcompression_0.6-0.tar.gz on local drive. It did not return error mess

[R] heatmap and automatic box sizes

Dear R, I have a list of X,Y coordinates and a ratio associated with each coordinate. The X and Y coordinates are continuous but random from 50-500, I would like to make a continuous heatmap of the ratios at each coordinate. One caveat is that the coordinates are clustered together do some bixes m

Re: [R] Applying Logical statement to DateTime string as factor

As mentioned, read R News 4/1 and its references for info on dates and times. On Tue, Aug 18, 2009 at 5:06 PM, Pat Schmitz wrote: > Fantastic, that works fantastically.  Is there a optimal way to search R > News topics, so that I can catch references like this? For all my searching > of list serve

Re: [R] aggregating values at discreet irregular time intervals into hourly values

Try this: library(tseries) # create 3 hours of test data set.seed(123) st <- as.POSIXct("2008-12-19 19:11:03") tt <- seq(from = st, to = st + 60 * 60 * 3, by = "min") x <- irts(tt, sample(0:1, length(tt), replace = TRUE)) library(zoo) # convert to zoo and chron z <- zoo(coredata(x), as.chron(fo

[R] How to Dedup a Spatial Points Data Set

I'm new to spatial analysis and am exploring numerous packages, mostly enjoying sp, gstat, and spBayes. Is there a function that allows the user to dedup a data set with multiple values at the same coordinates and replace those duplicated values with the mean at those coordinates? I've written s

Re: [R] how to draw pentagon?

Look at the my.symbols and ms.polygon functions in the TeachingDemos package. -- Gregory (Greg) L. Snow Ph.D. Statistical Data Center Intermountain Healthcare greg.s...@imail.org 801.408.8111 > -Original Message- > From: r-help-boun...@r-project.org [mailto:r-help-boun...@r- > project.or

Re: [R] Applying Logical statement to DateTime string as factor

Fantastic, that works fantastically. Is there a optimal way to search R News topics, so that I can catch references like this? For all my searching of list serves, and documentation, I found very little which could help me. I have however, had excellent help with R-help, Gabor, I appreciate your

[R] aggregating values at discreet irregular time intervals into hourly values

Hello R users, I'm a newby to R (and programming software at large) and I would need some help to sum up event data at discreet time and irregular time interval into a hourly frequency. Here is an example of my time series frame (irregular time-serie object - irts in the tseries package): ti

Re: [R] Transpose a dataset

## if x is a matrix this should work library(reshape) x <- structure(c(36.41099, 73.60079, 171.94, 67.48221, 131.917, 85.17079, 0.4015699, 9.4656186, 9.201167, 11.7657645, 14.4986667, 17.3150434, 35.1, 50, 153.5, 40, 97, 57.3, 36.2, 67, 166.5, 60, 122, 70.9, 36.5, 73, 173, 68, 132, 83.1, 36.666

Re: [R] Transpose a dataset

You mean ***transform*** a dataset, not ``transpose''. Transposition is something done to a matrix. If you do not get your terminology straight, no one can help you. You may get some mileage out of the reshape() function; possibly you may need the more elaborate facilities in the reshape packa

[R] trouble with Rcmdr

Dear R People: I am having trouble with Rcmdr. I used it on Saturday, and it worked fine. When I used it today, the window did not appear. This is version 2.9.1 on Ubuntu Jaunty Jackalope. Has anyone run into this, please? Any suggestions much appreciated. Thanks, Sincerely, Erin -- Erin

[R] combined value in forestplot function

Hi, How is the combined value (displayed by a diamond on the forest plot) of a meta-analysis illustrated by a forest plot calculated? The applied method and all calculation details are welcome. Thanks Carol [[alternative HTML version deleted]]

Re: [R] R CMD BATCH question under Ubuntu 9.04

Did you try just running 'R CMD BATCH test.R' from within the /var/www directory? or if you're calling from a php script, perhaps something funky is going on with what it's interpreting as directories? I've seen some weird stuff happen where scripts will take web_home (/var/www) and do someth

Re: [R] R CMD BATCH question under Ubuntu 9.04

On Aug 18, 2009, at 1:54 PM, julien cuisinier wrote: Dear list, I could not find a mailing list of R under ubuntu, so I give it a try here...any direction/suggestion welcome I am under the impression that the Ubuntu's are Debian ports: https://stat.ethz.ch/mailman/listinfo/r-sig-debian

[R] Error with acf/pacf functions

I'm trying to apply a acf/pacf function to a xts object and get the following error: Error in if (frequency > 1 && abs(frequency - round(frequency)) < ts.eps) frequency <- round(frequency) : missing value where TRUE/FALSE needed Any ideas? Thanks. -- View this message in context: http://ww

[R] R CMD BATCH question under Ubuntu 9.04

Dear list, I could not find a mailing list of R under ubuntu, so I give it a try here...any direction/suggestion welcome I run R 2.8.1 under Ubuntu 9.04 >From the terminal I am able to run a R CMD BATCH command when I am on my home >folder /home/user j...@jc-laptop:~$ R CMD BATCH /var/www/t

Re: [R] Installing quantreg package under Ubuntu

> -Original Message- > From: Frank E Harrell Jr [mailto:f.harr...@vanderbilt.edu] > Sent: Tuesday, August 18, 2009 5:20 AM > To: Ben Bolker > Cc: r-help@r-project.org; William Dunlap; roger koenker > Subject: Re: [R] Installing quantreg package under Ubuntu > > Ben Bolker wrote: > > > >

Re: [R] memory limit in R

do 'object.size' on all the objects in 'ls()'; also show the output of 'gc()' Sent from my iPhone On Aug 18, 2009, at 13:35, Hongwei Dong wrote: > The size of my .Rdata workspace is about 9.2 M and the data I'm > using is the only one object in this workspace. Is it a large one? > > Thanks.

Re: [R] three dimensions barchart

On Aug 18, 2009, at 7:18 AM, Yichih Hsieh wrote: Dear R community, I have one problem with figures. I draw the Relative Distribution graph, it looks like barchart(X,Y plot), but I have ten(year) Relative Distribution grapgs, have any command am can combine ten barcharts(X,Y plot) to become a

Re: [R] memory limit in R

The size of my .Rdata workspace is about 9.2 M and the data I'm using is the only one object in this workspace. Is it a large one? Thanks. Harry On Tue, Aug 18, 2009 at 4:21 AM, jim holtman wrote: > About 2GB is the limit of the address space on 32-bit windows (you can > get upto 3GB with a spec

Re: [R] libname version in R

Michael Yutzi wrote: > > Hello friends, > in SAS there is the 'libname' source. You associate a way.. ex: ab .. > "C:\My Paste\Works" and when u do > > data ab.example; > (...) > run; > > you are saving the archive "example" in a SAS format, and you can see > C:\My Paste\Works\example "gro

[R] Odd results with Chi-square test. (Not an R problem, but general statistics, I think)

I'm far from an expert on stats but what I think you are saying is if you try and compare Baseline with Version 3 you don't think your p-value is as good as version 1 and 2. I'm not 100% sure you are meant to do that with p-values but I'll let someone else comment on that!. tot

Re: [R] Plot(x,y)

On Aug 18, 2009, at 1:40 AM, Michael Knudsen wrote: On Sun, Aug 16, 2009 at 9:19 PM, malcolm Crouch wrote: plot(V6,V5, col="red") or plot(V6,V5) It seems that V5 and V6 are column names in your data matrix. If your matrix is called data, you should use plot(x$V6,x$V5,col="red") ITYM:

Re: [R] graph label greek symbol failure

On Tue, 2009-08-18 at 13:06 +0100, e-letter wrote: > On 17/08/2009, Michael Knudsen wrote: > > On Mon, Aug 17, 2009 at 12:51 PM, e-letter wrote: > > > >> I have tried to add the delta (δ) symbol to the y axis label and the > >> result is &D, using the command: > >> > >> ...ylab="δt"... > > > > Tr

Re: [R] Simulation Function - Save results

On Aug 17, 2009, at 1:40 PM, MarcioRibeiro wrote: Ok, the LIST function I understood... I didn't see how you got a "random" input to that function. Would seem to need one of the r functions as input. What I would like now is to simulate this Function A many times (S) in order to get S

Re: [R] value of nth percentile

Inline below. Bert Gunter Genentech Nonclinical Biostatisics -Original Message- From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org] On Behalf Of Steve Lianoglou Sent: Tuesday, August 18, 2009 8:07 AM To: Ajay Singh Cc: r-help@r-project.org Subject: Re: [R] value of nt

[R] Simulation function

Hi listers, I've been looking for a procedure, but I am not succeding... I have a function that give multiple results... Then, I would like to simulate this function n times, so I need to save/keep the n multiple results, in order to calculate my desired statistics... I have already tried with the

Re: [R] Newbie that don't understand R code

Yap -- View this message in context: http://www.nabble.com/Newbie-that-don%27t-understand-R-code-tp25009034p25028474.html Sent from the R help mailing list archive at Nabble.com. __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listin

[R] Odd results with Chi-square test. (Not an R problem, but general statistics, I think.)

Hi, I am working on a system which automatically answers user questions (such systems are commonly called "Question Answering systems"). I evaluated different versions of the same system on a publicly available test sets. Naturally, there is a fixed number of questions in the test set, and the sy

Re: [R] Replacing NA values in one column of a data.frame

Clearly the poor user needs to export something to one of those archaic stats systems like SAS or SPSS :) --- On Tue, 8/18/09, Bert Gunter wrote: > From: Bert Gunter > Subject: RE: [R] Replacing NA values in one column of a data.frame > To: "'Steve Lianoglou'" , "'John Kane'" > > Cc: r-help@

Re: [R] reading in MS Word files

Thanks guys, as I wanted to do a little preprocessing before importing into tm (the files have all sorts of stuff in them that I don't need), I used a "system" to invoke Abiword and do the batch conversions. Mark Mark W. Kimpel MD ** Neu

Re: [R] postscript, options

On Tue, 18 Aug 2009 16:08:36 +0100 e-letter wrote: EL> > postscript("/pathto/filename.eps",horizontal=FALSE,onefile=FALSE,paper="special") EL> > EL> This command created a small blank eps file. Its hard to say why without any example code. Does postscript("/pathto/filename.eps",width=8,height=

Re: [R] Specify Database location from R

I believe you will have to go into your ODBC configuration software outside of R, find where "SOUTHAMPTON" is defined, and change it so that it connects to that "different part of the server". Note that "different part of the server" is not standard database terminology, so it's kind of hard to

Re: [R] transformation/link glm in R

1. It sounds as if you say "left skewed" when you mean "right skewed" and v.v. Right skewed distributions have a long tail to the right. 2. Your 3rd variable may be censored, as you indicate it is bounded at 8 hr and that is also the mode. If so, this will complicate your analysis. 3. Your variab

[R] Posterior use of iterations output in a for() loop

Dear all, Please apologyse if this is too simple. I am a beginner. I am trying to extract the resulting output of each iteration in a for() loop in order to concatenate them and make an histogram out of that. Example: a<-c(1,2,3,3,3,4,6,7,7,8) b<-c(5,8) # function to obtain 2 (length of "a") ve

Re: [R] Tr : create a table in the console!!

Hi, On Aug 18, 2009, at 10:52 AM, Inchallah Yarab wrote: - Message transféré De : Inchallah Yarab À : r-help@r-project.org Envoyé le : Mardi, 18 Août 2009, 16h26mn 20s Objet : create a table in the console!! HI I want to do a table with R (in the console) GWP_Max NumberOfPolicies N

Re: [R] open txt

Stela Valenti Raupp wrote: Não consigo abrir a pasta txt no R, dá a mensagem: Warning message: In file(file, "r") : cannot open file 'plantula.txt': No such file or directory O arquivo está na mesma página do Scrip. Não sei qual é o problema [[alternative HTML version deleted]]

Re: [R] Remove columns

Hi, This line of code does the trick: a[,which(apply(a, 2, sum) != 0)] cheers, Paul Alberto Lora M wrote: Hi Everbody Could somebody help me.? I need to remove the columns where the sum of it components is equal to zero. For example a<-matrix(c(0,0,0,0,0,0,0,1,0,1,0,0,0,0,0,0,0,0,1,1,

Re: [R] create a table in the console!!

> # not sure what you want your data to look like, but > # one alternative is to create a dataframe: > > x <- data.frame("GWP_Max"=20, + No_GWPMax= 8, + "[0-1000]"= 4, + "[1000-3000]"= 3, + "[> 3000]"= 5, check.names=FALSE) > t(x) # transpose for output [,1] GWP_Max

Re: [R] open txt

On Aug 18, 2009, at 9:41 AM, Stela Valenti Raupp wrote: Não consigo abrir a pasta txt no R, dá a mensagem: Warning message: In file(file, "r") : cannot open file 'plantula.txt': No such file or directory O arquivo está na mesma página do Scrip. Não sei qual é o problema Tu file "plantul

Re: [R] Replacing NA values in one column of a data.frame

-Original Message- From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org] On Behalf Of Steve Lianoglou Sent: Tuesday, August 18, 2009 7:25 AM To: John Kane Cc: r-help@r-project.org; Steve Murray Subject: Re: [R] Replacing NA values in one column of a data.frame ...** B

Re: [R] value of nth percentile

Hi Ajay, On Aug 18, 2009, at 8:16 AM, Ajay Singh wrote: Dear All, I have to get the value of say 90th percentile of precipitation time series.. The series is of daily precipitation value of 96 years, I have to to get 90the percentile value of daily precipitation each year. If you know the

Re: [R] Newbie that don't understand R code

Hi, Comments inline and at end: On Aug 17, 2009, at 11:36 AM, kfcnhl wrote: I got some R code that I don't understand. Question as comment in code //where is t comming from, what is phi inverse rAC <- function(name, n, d, theta){ #generic function for Archimedean copula simulation illegalpar

Re: [R] printing a dataframe summary to disk

Philip A. Viton wrote: I'd like to write the summary of a dataframe to disk, so that it looks essentially the same as what you'd see on screen; but I can't seem to do it. Can someone tell me how? Thanks! Philip A. Viton City Planning, Ohio State University 275 West W

Re: [R] postscript, options

On 18/08/2009, Stefan Grosse wrote: > On Tue, 18 Aug 2009 13:52:58 +0100 e-letter wrote: > > EL> On 18/08/2009, Rodrigo Aluizio wrote: > EL> > I'm not shure but I guess that you miss a " and putted it in the > EL> > wrong place! > EL> > Try this: > EL> > postscript("/pathto/filename.eps",horizon

Re: [R] reading in MS Word files

On Tue, Aug 18, 2009 at 12:00:07PM +0200, Mark Kimpel wrote: > I am familiar with packages that read and write Excel files on both Windows > and Linux platforms. > > Do any packages provide similar functionality for MS Word files? I have a > lot of text processing to do and the text is embedded in

Re: [R] lm.fit algo

Hi, On Aug 17, 2009, at 5:09 PM, Pavlo Kononenko wrote: Hi, everyone, This is a little silly, but I cant figure out the algorithm behind lm.fit function used in the context of promax rotation algorithm: The promax function is: promax <- function(x, m = 4) { if(ncol(x) < 2) return(x) dn

[R] Tr : create a table in the console!!

- Message transféré De : Inchallah Yarab À : r-help@r-project.org Envoyé le : Mardi, 18 Août 2009, 16h26mn 20s Objet : create a table in the console!! HI I want to do a table with R (in the console) GWP_Max NumberOfPolicies No_GWPMax 8 [0-1000] 4 [1000-3000] 3 [> 3000] 5 i

Re: [R] Remove columns

Try this, a[ ,as.logical(colSums(a))] mind an unfortunate logical vs integer indexing trap: isTRUE(all.equal( a[ ,!!colSums(a)] , a[ ,colSums(a)] )) [1] FALSE HTH, baptiste 2009/8/18 Alberto Lora M : > Hi Everbody > > Could somebody help me.? > > I need to remove the columns where the sum o

Re: [R] postscript, options

On Tue, 18 Aug 2009 13:52:58 +0100 e-letter wrote: EL> On 18/08/2009, Rodrigo Aluizio wrote: EL> > I'm not shure but I guess that you miss a " and putted it in the EL> > wrong place! EL> > Try this: EL> > postscript("/pathto/filename.eps",horizontal=FALSE,onefile=FALSE) EL> > EL> > Hope It helps

Re: [R] printing a dataframe summary to disk

Hi Philip Try: ?sink. Best, Tal On Mon, Aug 17, 2009 at 5:54 PM, Philip A. Viton wrote: > > I'd like to write the summary of a dataframe to disk, so that it looks > essentially the same as what you'd see on screen; but I can't seem to do > it. Can someone tell me how? Thanks! > > > ---

Re: [R] Remove columns

Hi Alberto, On Aug 18, 2009, at 4:14 AM, Alberto Lora M wrote: Hi Everbody Could somebody help me.? I need to remove the columns where the sum of it components is equal to zero. For example a<-matrix(c(0,0,0,0,0,0,0,1,0,1,0,0,0,0,0,0,0,0,1,1,0,0,1,0), ncol=4) a [,1] [,2] [,3] [,4]

Re: [R] printing a dataframe summary to disk

Hi, On Aug 17, 2009, at 10:54 AM, Philip A. Viton wrote: I'd like to write the summary of a dataframe to disk, so that it looks essentially the same as what you'd see on screen; but I can't seem to do it. Can someone tell me how? Thanks! Look at: ?write.table and friends You can write the

Re: [R] Replacing NA values in one column of a data.frame

OOPs clumsy cut and paste. :( --- On Tue, 8/18/09, Uwe Ligges wrote: > From: Uwe Ligges > Subject: Re: [R] Replacing NA values in one column of a data.frame > To: "John Kane" > Cc: r-help@r-project.org, "Steve Murray" > Received: Tuesday, August 18, 2009, 10:25 AM > > > John Kane wrote: >

Re: [R] Replacing NA values in one column of a data.frame

Hi, On Aug 18, 2009, at 10:19 AM, John Kane wrote: Perhaps testdata$onecolumn[testdata$onecolumn==NA] <- I don't think this would work -- the is.na function is the way to go: R> a <- c(1,2,3,NA,NA,10,20,NA,NA) R> a [1] 1 2 3 NA NA 10 20 NA NA R> a[a == NA] <-999 R> a [1] 1 2 3 NA

[R] create a table in the console!!

HI I want to do a table with R (in the console) GWP_Max NumberOfPolicies No_GWPMax 8 [0-1000] 4 [1000-3000] 3 [> 3000] 5 i begin by calculate the number of policies in each class :  Data1 <- read.csv2("c:/Total1.csv", sep=",") > Data2 <- read.csv2("c:/GWPMax1.csv",sep=",")[1:20,1:2] > M <

[R] open txt

Não consigo abrir a pasta txt no R, dá a mensagem: Warning message: In file(file, "r") : cannot open file 'plantula.txt': No such file or directory O arquivo está na mesma página do Scrip. Não sei qual é o problema [[alternative HTML version deleted]] ___

[R] value of nth percentile

Dear All, I have to get the value of say 90th percentile of precipitation time series.. The series is of daily precipitation value of 96 years, I have to to get 90the percentile value of daily precipitation each year. If you know the R code or command for this please let me know. I would appreci

Re: [R] Replacing NA values in one column of a data.frame

John Kane wrote: Perhaps testdata$onecolumn[testdata$onecolumn==NA] <- Well, with is.na() there is a much better chance to make it work Uwe Ligges --- On Mon, 8/17/09, Steve Murray wrote: From: Steve Murray Subject: [R] Replacing NA values in one column of a data.frame To: r-h

Re: [R] Log logistic Distribution - Parameter estimation

Maithili Shiva wrote: Dear R helpers, I have following data c(19679.50,5975.72,17555.52,20078.08,1956.97,8383.25,7899.27,4881.75,2719.29,11701.96,13104.25,10081.01,9417.72,14976.89,7119.62,8373.56,24580.52,8330..70,9110.07,7848.08,12910.46,10310.53,20390.34,15824.71,21703.33) How do I fit 3 p

Re: [R] Multiple comparison on lme model with 2 fixed factors

Sorry to come back again, but with this and other (real) data, I'm not able to obtain those comparisons. I'm using the following command and obtain the following result (lmer model). glht(lmelac.ex, linfct = mcp("ath1:tim20 - ath1:tim40 = 0"="Tukey"))->comp Erreur dans mcp2matrix(model, linfct =

Re: [R] Differing rows in two data.frames?

Mark Knecht wrote: On Mon, Aug 17, 2009 at 2:25 PM, Mark Knecht wrote: Hi, What commands would I look at to compare row-by-row two data frames of the same size and print out any differences between the two? Thanks, Mark So I got an answer using the following code: XXX = cbind(Syst

  1   2   3   >