Re: [R] Replace NAs in a data frame with elements randomly generated from the beta distribution

Thank you so much for your assistance. A.K's suggestion is exactly what I was looking for. Thanks so much... On Wed, 10 Apr 2013 23:12:25 +0200 peter dalgaard pda...@gmail.com wrote On Apr 10, 2013, at 13:09 , Sibusiso Ndzukuma wrote: Hi there! Please help me, I am trying to replace

Re: [R] Plot two separate curves in R Graphics and R Lattice package

You need to apply the same transformation to the x argument of panel.lines as you are setting up for the panel.xyplot. -- David Sent from my iPhone On Apr 10, 2013, at 8:10 PM, jpm miao miao...@gmail.com wrote: Hi David, Thanks. How can I make the graph in log scale? Thanks

[R] Encircle the closest points

Dear R-user, For the matplot of the following data drawn by matplot(t, pch=1) structure(c(2372L, 2494L, 2170L, 2238L, 1918L, 715L, 1573L, 2437L, 2481L, 1881L, 440L, 2165L, 343L, 121L, 597L, 1163L, 1459L, 2093L, 1454L, 1761L, 1872L, 174L, 2040L, 552L, 1403L, 6381L, 1649L, 4136L, 4415L, 10650L,

[R] Encircle the closest points

Dear R-user, For the matplot of the following data drawn by matplot(t, pch=1) structure(c(2372L, 2494L, 2170L, 2238L, 1918L, 715L, 1573L, 2437L, 2481L, 1881L, 440L, 2165L, 343L, 121L, 597L, 1163L, 1459L, 2093L, 1454L, 1761L, 1872L, 174L, 2040L, 552L, 1403L, 6381L, 1649L, 4136L, 4415L, 10650L,

[R] Cannot find ldfortran (R on Cygwin)

Hi I am new to Cygwin and Linux. I installed R under Cygwin as part of the setup I chose All during installation, for all packages. So I have the FULL installlation of cygwin up and running, including gfortran. *Under Cygwin, how do I check and configure the path to the various libraries?* I am

Re: [R] Calculating std errors of marginal effects in interactions

Hi Richard, Have a look at this blog: http://www.r-bloggers.com/recovering-marginal-effects-and-standard-errors-from-interaction-terms-in-r/ Hope it helps! José José Iparraguirre Chief Economist Age UK T 020 303 31482 E jose.iparragui...@ageuk.org.uk Twitter @jose.iparraguirre@ageuk Tavis

Re: [R] Encircle the closest points

You might wish to use the nndist() and nnwhich() functions from the spatstat package. E.g., using your data, Y - structure(c(2372L, 2494L, 2170L, 2238L, 1918L, 715L, 1573L, 2437L, 2481L, 1881L, 440L, 2165L, 343L, 121L, 597L, 1163L, 1459L, 2093L, 1454L, 1761L, 1872L, 174L, 2040L, 552L,

Re: [R] Optimization problem

Hello, You cannot change the numerical accuracy, it's a built-in constant. To see it use ?.Machine .Machine$double.eps # smallest value different from zero Actually, .Machine$double.eps is the the smallest positive floating-point number x such that 1 + x != 1 You can try the following

[R] unsubscribe from mailing list help!

Dear all, I cant seem to unsubscribe from the mailing list even though it has removed my account, any tips? Best wishes, Lorna [[alternative HTML version deleted]] __ R-help@r-project.org mailing list

[R] %*%

What does these operators do: %*% Thanks -- Shane [[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

Re: [R] %*%

On 11-Apr-2013 10:25:17 Shane Carey wrote: What does these operators do: %*% Thanks -- Shane Enter the command ?%*% and you will see: Matrix Multiplication Description: Multiplies two matrices, if they are conformable. If one argument is a vector, it will be promoted to

Re: [R] %*%

On 11-04-2013, at 12:25, Shane Carey careys...@gmail.com wrote: What does these operators do: %*% Please read a manual. And at least try some searching in R: ?%*% Introduction to R (http://cran.r-project.org/doc/manuals/r-release/R-intro.html) section 5.7.1 Berend Thanks -- Shane

Re: [R] %*%

Thank you both On Thu, Apr 11, 2013 at 11:46 AM, Berend Hasselman b...@xs4all.nl wrote: On 11-04-2013, at 12:25, Shane Carey careys...@gmail.com wrote: What does these operators do: %*% Please read a manual. And at least try some searching in R: ?%*% Introduction to R (

Re: [R] %*%

On Apr 11, 2013, at 12:25 , Shane Carey wrote: What does these operators do: %*% Try help(%*%) -- Peter Dalgaard, Professor Center for Statistics, Copenhagen Business School Solbjerg Plads 3, 2000 Frederiksberg, Denmark Phone: (+45)38153501 Email: pd@cbs.dk Priv: pda...@gmail.com

Re: [R] unsubscribe from mailing list help!

On Apr 11, 2013, at 12:03 , McAusland, Lorna wrote: Dear all, I cant seem to unsubscribe from the mailing list even though it has removed my account, any tips? Patience? The mailing list software can't un-send mail so a handful of posts may still be in transit when you unsubscribe.

[R] Problem with caret

Hello all. I am using method bag() to develop a bagging with my dataset. When I do bag( dataP, dataP$Score, B=10) I obtain this error: error en bag.default(Score ~ ., data = dataP, B = 10) : entrada en evaluacion: recursivo por defecto o problemas anteriores? ( in english, evaluating

[R] (no subject)

Dear all, Is there a quick and easy way of converting utf characters to the \u form (necessary e.g. for packages)? I mean something working like this: utf2u(õäöü) [1] \u00f5\u00e4\u00f6\u00fc It is easy to program but perhaps someone already has implemented this. (I couldn't find

Re: [R] Cannot find ldfortran (R on Cygwin)

On 13-04-11 2:57 AM, Lalitha Viswanathan wrote: Hi I am new to Cygwin and Linux. I installed R under Cygwin as part of the setup I chose All during installation, for all packages. So I have the FULL installlation of cygwin up and running, including gfortran. *Under Cygwin, how do I check and

[R] Creating a GUI based workflow

Hello R-users, I am interested to know if there is a tool that allows us to build SAS Enterprise Miner like workflows using the R scripts. For example, let's say I have the R scripts/functions for the following tasks: 1. Data Cleaning 2. Model Building 3. Model Validation 4. Visualization I

Re: [R] Error when using fitdist function in R

Hi I think you get some answer from prof.Ripley. However when I did hist(k) plot(density(k)) it seems to me, that k is mixture of 2 normal distributions. So you probably cannot fit any single distribution to it. Regards Petr From: Paul Bernal [mailto:paulberna...@gmail.com] Sent: Wednesday,

Re: [R] (no subject)

On 13-04-11 7:56 AM, Kenn Konstabel wrote: Dear all, Is there a quick and easy way of converting utf characters to the \u form (necessary e.g. for packages)? I mean something working like this: utf2u(õäöü) [1] \u00f5\u00e4\u00f6\u00fc It is easy to program but perhaps someone

Re: [R] Creating a GUI based workflow

On Thu, Apr 11, 2013 at 8:09 AM, Ravishankar Rajagopalan viora...@gmail.com wrote: Hello R-users, I am interested to know if there is a tool that allows us to build SAS Enterprise Miner like workflows using the R scripts. For example, let's say I have the R scripts/functions for the

Re: [R] Table with n(%) formatting: A better way to do this?

Hello Dr. Murdoch, Thanks for your reply. Your tables package is interesting. Spent some time experimenting with it yesterday and it looks like something I can use. The n(%) format seems to be the standard for what I'm doing. So I can't directly use the output your package produces. It does

Re: [R] Table with n(%) formatting: A better way to do this?

On 11/04/2013 8:40 AM, Paul Miller wrote: Hello Dr. Murdoch, Thanks for your reply. Your tables package is interesting. Spent some time experimenting with it yesterday and it looks like something I can use. The n(%) format seems to be the standard for what I'm doing. So I can't directly use

Re: [R] Creating a GUI based workflow

You might look into Kepler: https://kepler-project.org/ Sarah On Thu, Apr 11, 2013 at 8:09 AM, Ravishankar Rajagopalan viora...@gmail.com wrote: Hello R-users, I am interested to know if there is a tool that allows us to build SAS Enterprise Miner like workflows using the R scripts. For

Re: [R] %*%

Hello, Just try ?%*% And the hel page title will give you the answer. Hope this helps, Rui Barradas Em 11-04-2013 11:25, Shane Carey escreveu: What does these operators do: %*% Thanks __ R-help@r-project.org mailing list

[R] Make barplot with error bars, anova out of a table

Helo everybody, I'm new to R and have some issues with my data in R. My raw data look like that: ID Day size 1 1 7 1 1 7.2 1 1 7.1 2 1 7.3 2 1 7.4 2 1 7.2 3 1 7 3 1 7.1 3 1 7.5 4 1 7.3 4 1 7.2 4 1 7.6 1 2 7 1 2 7.2 1 2 7.1 2 2 7.1 2 2 7.4 2 2 7.2 3 2 7.5 3 2 7.1 3 2 7.5 4 2 7.2 4 2 7.2 4 2 7.3

Re: [R] (no subject)

HI, May be this helps: library(Unicode) utf2uxxx- function(x){  a- tolower(as.u_char(utf8ToInt(x)))  a- gsub([+],,a)  a- paste(paste(\\,a,sep=),collapse=)  a-cat('',a,'',sep=)} utf2uxxx(õäöü) #\u00f5\u00e4\u00f6\u00fc sessionInfo() R version 3.0.0 (2013-04-03) Platform: x86_64-unknown-linux-gnu

Re: [R] (no subject)

Thanks a lot both! I finally used it (on windows) as library(Unicode) utf2uxxx- function(x){ x-iconv(x, to=utf-8) # doesn't look nice but seems necessary on Windows a- tolower(as.u_char(utf8ToInt(x))) a- gsub([+],,a) a- paste(paste(\\,a,sep=),collapse=) a-cat('',a,'',sep=) } It converts all

[R] Notes and warnings in building a package ('ms-dos style file name', 'no visible binding' and using .C(bincount, ..., PACKAGE = base))

Dear all, While buildign a package (R 2.15.2 on windows XP) I had three warnings/notes -- I would be grateful for any hints on what to do. 1. The first one is probably not really important if I'm not missing something: cygwin warning: MS-DOS style path detected: C:/x/y/z/package_1.0.tar.gz

[R] Cluster analysis

I am doing cluster analysis of my SNPs data. I have 2 questions: 1. I draw the cluster in hclust using the following codes.change direction to vertical. data - read.table(as.matrix(file.choose()), header=T, row.names = 1, sep=\t) plot(hclust(as.dist(data),method=complete)) it is horizontal,

[R] plotting the cumulative hazard function of a process whose intensity function is specified

I would like to simulate a stochastic process by specifying an intensity function, say lambda(t)=0.05 + 3  exp(−t) I have the following code, but it has a number of problems. One issue I have is I do not know where the time variable is. I would appreciate your help with this code.

Re: [R] Reshaping a table

-Original Message- Here's the relevant portion of the as.table Help file: But ... is an argument to table(), not to as.table; this part of the help file is not referring to as.table. The first argument to as.table is x, which is an 'arbitrary R object'. as.table is a generic; it

Re: [R] Reshaping a table - PS

as.table(dat1) #Error in as.table.default(dat1) : cannot coerce to a table You _can_ coerce to a matrix first: as.table( as.matrix(dat1) ) S Ellison *** This email and any attachments are confidential. Any use...{{dropped:8}}

[R] Connecting SSAS Cube From R

Hi All, Does anyone have a connection string example for connecting SqlServer Analysis Services cube (MS SQL Server 2005 or 2008). WITH R . For Connecting Sqlserver, i tried using RODBC/RJDBC Connection..but Still hard luck for me..Could you please provide me..which package i need to download

Re: [R] Creating a GUI based workflow

Or red-r.org.. Tom Sent from my iPhone On Apr 11, 2013, at 9:14 AM, Sarah Goslee sarah.gos...@gmail.com wrote: You might look into Kepler: https://kepler-project.org/ Sarah On Thu, Apr 11, 2013 at 8:09 AM, Ravishankar Rajagopalan viora...@gmail.com wrote: Hello R-users, I am

Re: [R] remove duplicates in data frame

-Original Message- How do I find and remove the two duplicate rows? Try ?unique() See also ?duplicated and consider the usage d[ !duplicated(d), ] where d is your data frame. But read the 'details' part of the help page closely... S Ellison

Re: [R] Reshaping a table

Yes. My error was previously pointed out to me and there was further offlist discussion. Maybe I should have kept it onlist, but I was reticent about displaying my stupidity. :( Mea Culpa. Cheers, Bert On Thu, Apr 11, 2013 at 9:11 AM, S Ellison s.elli...@lgcgroup.com wrote: -Original

Re: [R] Unbalance to balance panel

Is there any package to balance the unbalance panel data? If you mean you have an unbalanced experiment/observation data set, what balancing strategy would you be willing to accept? (for example, are you looking to delete data to achieve balance, or use imputation to invent new data to fill

[R] Dotchart per groups

Hi all, I would like to ask you for help. I did a dotplot - using dotchart function. There are two localites (loc) with values 75 or 56 in my data ZZ. The f column has 4 levels: P1, S1, S8, R6. The dataframe is ordered by N value, pchloc value is assign to use pch in plot. head(ZZ) loc f

Re: [R] Remove data 3 standard deviatons from the mean using R?

-Original Message- I have a very long list of data-points (+2300) and i know from my histogram that there are outliers which are affecting my mean. If extreme values are known to be unreliable*, take a look at robust statisical methods (or even the median) instead of the mean.

[R] Understanding contents of packages

Colleagues During the transition to R3.0.0 (OS X), one of the packages that I used -- SASxport -- did not work until I used the type=source option in install.packages. This led to an adventure: 1. I downloaded the package source:SASxport_1.2.4.tar.gz 2. in the

[R] group data

Hey, I have a dataset and I want to identify the records by groups for further use in ggplot. Here is a sample data: ID Value AL1 1 AL2 2 CA1 3 CA4 4 I want to identify all the records that in the same state (AL1 AND A2), group them as AL, and do the same for CA1 and CA4. How can I have

Re: [R] Understanding contents of packages

Hi Dennis, Out of curiosity, what is the purpose of this exercise? Best, Ista On Thu, Apr 11, 2013 at 2:34 PM, Fisher Dennis fis...@plessthan.com wrote: Colleagues During the transition to R3.0.0 (OS X), one of the packages that I used -- SASxport -- did not work until I used the

Re: [R] Understanding contents of packages

On Thu, Apr 11, 2013 at 2:42 PM, Fisher Dennis fis...@plessthan.com wrote: to learn how to use .C() Have you read the manual? http://cran.r-project.org/doc/manuals/R-exts.html#System-and-foreign-language-interfaces Best, Ista Dennis Fisher MD P (The P Less Than Company) Phone:

Re: [R] group data

Hello, Try the following. dat - read.table(text = ID Value AL1 1 AL2 2 CA1 3 CA4 4 , header = TRUE, stringsAsFactors = FALSE) dat$State - substr(dat$ID, 1, 2) Note that this dependes on having State being defined by the first two characters of ID. Hope this helps, Rui Barradas

Re: [R] group data

Thanks! But What if I have a very large data with 1000+rows, anyway to identify all AL1 and AL2 under ID and mark them as AL under new column State? Thanks. On Thu, Apr 11, 2013 at 11:54 AM, Rui Barradas ruipbarra...@sapo.pt wrote: Hello, Try the following. dat - read.table(text = ID

Re: [R] group data

Hi, dat1- read.table(text= ID  Value AL1  1 AL2  2 CA1  3 CA4  4 ,sep=,header=TRUE,stringsAsFactors=FALSE) dat2- dat1 dat1$State-gsub(\\d+,,dat1$ID) dat1 #   ID Value State #1 AL1 1    AL #2 AL2 2    AL #3 CA1 3    CA #4 CA4 4    CA #or library(stringr)

Re: [R] Understanding contents of packages

On 11-04-2013, at 20:34, Fisher Dennis fis...@plessthan.com wrote: Colleagues During the transition to R3.0.0 (OS X), one of the packages that I used -- SASxport -- did not work until I used the type=source option in install.packages. See below. There is a binary version available.

Re: [R] Understanding contents of packages

On Apr 11, 2013, at 1:34 PM, Fisher Dennis fis...@plessthan.com wrote: Colleagues During the transition to R3.0.0 (OS X), one of the packages that I used -- SASxport -- did not work until I used the type=source option in install.packages. This led to an adventure: 1. I

[R] Adding time series to time graphs

Hello: I have done this before but cannot figure out how to do it again. I would like to graph campaign evolution of news stories on certain topics. The campaign time period is as follows: campaign-seq.Date(from=as.Date('2011-09-06'), to=as.Date('2011-10-5'), by=1) I have a table of newspaper

Re: [R] group data

I think I just misinterpret~Thanks for your help~ On Thu, Apr 11, 2013 at 11:54 AM, Rui Barradas ruipbarra...@sapo.pt wrote: Hello, Try the following. dat - read.table(text = ID Value AL1 1 AL2 2 CA1 3 CA4 4 , header = TRUE, stringsAsFactors = FALSE) dat$State -

Re: [R] Adding time series to time graphs

-Original Message- From: r-help-boun...@r-project.org [mailto:r-help-bounces@r- project.org] On Behalf Of Simon Kiss Sent: Thursday, April 11, 2013 12:35 PM To: r-help@r-project.org Subject: [R] Adding time series to time graphs Hello: I have done this before but cannot figure out

[R] Batch read netcdf files and average one variable

Hi all, I'm a new R user. I now have daily netcdf data for year 1979 such as these: sm19790101.1.nc sm19790102.1.nc . . . sm19791231.1.nc I need to average a variable called sm to monthly resolution. I can now do this: glob2rx(sm197901*.1.nc)

[R] odfWeave: Some questions about potential formatting options

Hello All, Learning to use the odfWeave package. I really like the package. It has good documentation, makes some very nice looking tables, and seems to have lots of options for customizing output. There are a few things I'd like to do that don't seem to be covered in the documentation

Re: [R] Make barplot with error bars, anova out of a table

Hi On Thu, Apr 11, 2013 at 6:43 AM, Simza alex_steckba...@gmx.at wrote: Helo everybody, I'm new to R and have some issues with my data in R. My raw data look like that: ID Day size 1 1 7 1 1 7.2 1 1 7.1 2 1 7.3 2 1 7.4 2 1 7.2 3 1 7 3 1 7.1 3 1 7.5 4 1 7.3 4 1 7.2 4 1 7.6 1 2 7 1 2 7.2 1 2

Re: [R] Graphic window dimensions

Dear Brian, I have solved the problem using dev.size, as follows: I obtain dev.size()[2] and I take it as my height value. Then I obtain my width value as my height value divided by 1.41 (in order to reach the same aspect ratio as din A4). Thanks!. Eva --- El lun, 8/4/13, Prof Brian Ripley

[R] Read the data from a text file and reshape the data

I have a data set for different time intervals. The data has three comment lines before data for each time interval. For each time interval there are 500 data points. I want to change the dataset such that I have the following format: t1t2t3

Re: [R] Make barplot with error bars, anova out of a table

On Apr 11, 2013, at 4:43 AM, Simza wrote: Helo everybody, I'm new to R and have some issues with my data in R. My raw data look like that: ID Day size 1 1 7 1 1 7.2 1 1 7.1 2 1 7.3 2 1 7.4 2 1 7.2 3 1 7 3 1 7.1 3 1 7.5 4 1 7.3 4 1 7.2 4 1 7.6 1 2 7 1 2 7.2 1 2 7.1 2 2 7.1 2 2 7.4 2 2

[R] Ordination Plotting: Warning: Species scores not available

Hi, I am working with a species-by-trait .csv file (columns=traits, rows=species) and get the following warning message when trying to plot results of both metaMDS and pcoa:  Warning message: In ordiplot(x, choices = choices, type = type, display = display,  :   Species scores not available I

Re: [R] Read the data from a text file and reshape the data

Is this what you are looking for: input - readLines(C:\\Users\\Owner\\Downloads\\WAT_DEP.DAT) start - grep(N:SNAPSHOT, input) # find start of the data # add index of what would have been the last block start - c(start, tail(start, 1) + 53L) # now read in the data using 'text' parameter of

Re: [R] Read the data from a text file and reshape the data

Hi Jim, Yes, this is exactly what I was looking for. The code works perfect. Thanks a lot for your time and effort. Best, Janesh Devkota From: jim holtman [mailto:jholt...@gmail.com] Sent: Thursday, April 11, 2013 7:32 PM To: Janesh Devkota Cc: R mailing list Subject: Re: [R] Read

Re: [R] Solving an integral in R gives the error “The integral is probably divergent”

I don't get an error message (after I correct the missing line break after the comment b- sapply(a, Cfun, upper=1) b [1] 1.583458e-54 7.768026e-50 2.317562e-45 4.206260e-41 4.645737e-37 3.123801e-33 1.279358e-29 3.193257e-26 4.860876e-23 [10] 4.516582e-20 2.564400e-17

Re: [R] Read the data from a text file and reshape the data

Hi, May be this helps:  lines1- readLines(WAT_DEP.DAT.part) indx- which(grepl([*],lines1)) indx2-indx[seq(from=indx[2],length(indx),by=2)]+1 lines2-str_trim(lines1[indx2],side=left) dat1-read.table(text=lines2,sep=,header=FALSE) library(stringr) lst1-

[R] A strange behaviour in the graphical function curve

I thought the curve function was a very flexible way to draw functions. So I could plot funtions like the following: # I created a function to produce functions, for instance: fp - function(m,b) function(x) sin(x) + m*x + b # So I can produce a function like this ff - fp(-0.08, 0.2)

Re: [R] Problem with caret

Le 11/04/13 13:21, Nicolás Sánchez a écrit : Hello all. I am using method bag() to develop a bagging with my dataset. When I do bag( dataP, dataP$Score, B=10) I obtain this error: error en bag.default(Score ~ ., data = dataP, B = 10) : entrada en evaluacion: recursivo por defecto o

Re: [R] A strange behaviour in the graphical function curve

On 12-04-2013, at 05:15, Julio Sergio julioser...@gmail.com wrote: I thought the curve function was a very flexible way to draw functions. So I could plot funtions like the following: # I created a function to produce functions, for instance: fp - function(m,b) function(x) sin(x) +