Re: [R] classes and methods

Hi Edna -- Quoting Edna Bell : I have put together a test class and methods setClass("foo",representation(x="numeric")) [1] "foo" This creates an *S4* class which HAS A 'slot' x that is of type numeric. setMethod("plot","foo",function(x,y,...)boxplot(x,...)) [1] "plot" x <- rnorm(100)

[R] Howto Supress Extra Blank Page in gridBase

Dear all, I have a simple plot using "gridBase" like this. The problem occurs whenever I execute this code there is always a blank page created before the actual plot. How can we disable that blank page? I am using: R version 2.7.2 (2008-08-25) and gridBase version: 0.4-3 __ BEGIN__ library(g

[R] classes and methods

I have put together a test class and methods > setClass("foo",representation(x="numeric")) [1] "foo" > setMethod("plot","foo",function(x,y,...)boxplot(x,...)) [1] "plot" > x <- rnorm(100) > class(x) <- "foo" > plot(x) > y <- new("foo",x=rnorm(50)) > class(y) [1] "foo" My question is: should I

Re: [R] functions within a list

On Mar 19, 2009, at 9:01 PM, Altaweel, Mark R. wrote: Hi, I am trying to perform various functions on a list with a number of elements. For example. I would like to take the mean of different variable values in the entire list. As an example, say I have a list with 1000 elements and variab

Re: [R] how to make aggregation in R ?

On Thu, Mar 19, 2009 at 8:40 PM, jim holtman wrote: > Try this technique.  I use it with large data objects since it is > sometime faster, and uses less memory, by using indices: > > x <- read.table(textConnection("  v1 v2 n1 n2 > 1   a a1  1 21 > 2   a a1  2 22 > 3   a a1  3 23 > 4   a a2  4 24 >

[R] CCA - manual selection

Hello, I am trying to obtain f-values for response (independent) variables from a CCA performed in vegan package, to see which ones of them have significative influence in my dependent variables (like the manual selection in canoco), but I can't find any function (or package) that do such a thing

Re: [R] print of objects in R

xtable in the library xtable does a good job if you are using latex to write up your results. For example: > xtable(matrix(rnorm(20),5,5)) % latex table generated in R 2.8.0 by xtable 1.5-4 package % Fri Mar 20 13:48:53 2009 \begin{table}[ht] \begin{center} \begin{tabular}{rr} \hline & 1 &

Re: [R] which.na

> ?is.na > x <- c(NA, 3, 4, 5, NA) > which(is.na(x)) [1] 1 5 Charles Annis, P.E. charles.an...@statisticalengineering.com phone: 561-352-9699 eFax: 614-455-3265 http://www.StatisticalEngineering.com -Original Message- From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project

[R] print of objects in R

Hello, I have been watching my output as I create functions and do other things in r. One thing I don't like is the [1,] type notation at the beginning of a line. I have been able to change that to a number such as 1 2 etc. using as.data.frame(object). How can I stop the printing of a lin

[R] which.na

Hi R- users I was wondering if there is any function equivalent to which.na used in S+? Thanks much in advance! Regards, Santosh [[alternative HTML version deleted]] __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r

Re: [R] how to make aggregation in R ?

Try this technique. I use it with large data objects since it is sometime faster, and uses less memory, by using indices: x <- read.table(textConnection(" v1 v2 n1 n2 1 a a1 1 21 2 a a1 2 22 3 a a1 3 23 4 a a2 4 24 5 a a3 5 25 6 b b1 6 26 7 b b1 7 27 8 b b2 8 28 9 b b2

Re: [R] Extreme AIC or BIC values in glm(), logistic regression

Maggie Wang wrote: Hi, Dieter, Gad, and all, Thank you very much for your reply! So here is my data, you can copy it into a file names "sample.txt" Hi Maggie, With this data (allowing for more iterations) I get: lr <- glm(fo, family=binomial(link=logit), data=matrix, control=glm.contr

Re: [R] how to make aggregation in R ?

Here are two solutions: > aggregate(testDF[c("n1", "n2")], testDF[c("v1", "v2")], sum) v1 v2 n1 n2 1 a a1 6 66 2 a a2 4 24 3 a a3 5 25 4 b b1 13 53 5 b b2 27 87 6 c c1 11 31 7 c c2 39 99 8 c c3 15 35 9 d d1 16 36 10 d d2 17 37 11 d d3 18 38 12 d d4 39 79 > library(sqldf

[R] randomForest

Hi! I am dealing with random forest using R. Is there a way to sample a fixed no.of rows from a dataset for use with different trees in random Forest. To be more clear, my data set contains 1500 rows, and I am growing 500 trees in Random Forest Is it possible to sample only 500 rows of data from

Re: [R] how to make aggregation in R ?

Hi, I think I found the solution. Using doBy library, I got: testDF.result2 <- summaryBy(n1+n2 ~ v1+v2, data = testDF, FUN=sum) > testDF.result2 v1 v2 n1.sum n2.sum 1 a a1 6 66 2 a a2 4 24 3 a a3 5 25 4 b b1 13 53 5 b b2 27 87 6 c c1

Re: [R] converting to strings

On 20/03/2009, at 12:25 PM, Peter Palenchar wrote: I have a list: x<-c("abcg", "bcah", "defi") I want the third and fourth thing in each item of the list ("bc", "ca", "ef"). ?substr ## Attention:\ This e-mail me

[R] how to make aggregation in R ?

Hi, I am trying to aggregate the sum of my test data.frame as follow: testDF <- data.frame(v1 = c("a", "a", "a", "a", "a", "b", "b", "b", "b", "b", "c", "c", "c", "c", "c", "d", "d", "d", "d", "d"), v2 = c("a1", "a1", "a1", "a2", "a3", "b1", "b1", "b2", "b2", "b2", "c1", "c2"

Re: [R] run r scripts (tinn) via command

Thank you! Works perfectly well! jholtman wrote: > > You just create the "superior" script as a set of 'source' statements > reading in the files that you want to execute in order. > > On Thu, Mar 19, 2009 at 4:13 AM, thoeb wrote: >> >> The scripts are written with the TINN editor (***.r) and

[R] converting to strings

I have a list: x<-c("abcg", "bcah", "defi") I want the third and fourth thing in each item of the list ("bc", "ca", "ef"). Peter __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting gu

[R] Mean-replacing NAs in a 3d array

Hi all I have a 3d array containing missing values. > Xa , , 1 [,1] [,2] [1,]13 [2,] NA4 , , 2 [,1] [,2] [1,]57 [2,] NA NA , , 3 [,1] [,2] [1,]9 11 [2,] 10 12 I want to replace the missing values with the mean, but the mean of each 'page' i

[R] Package HDF5

The package works fine, but it seems not to provide access to fields, but only to the whole data in a certain file (function hdf5load(file, load = TRUE, verbosity = 0, tidy = FALSE)) Since hdf5 organizes the data and metadata in a hierarchical structure, we must explore in our problem (file>7GB),

[R] functions within a list

Hi, I am trying to perform various functions on a list with a number of elements. For example. I would like to take the mean of different variable values in the entire list. As an example, say I have a list with 1000 elements and variables called Id and Tick. What I would like to do is take the

Re: [R] Difference between gam() and loess().

Rolf Turner wrote: It seems that in general gam(y~lo(x)) # gam() from the gam package. and loess(y~x) give slightly different results (in respect of the predicted/fitted values). Most noticeable at the endpoints of the range of x. Can anyone enlighten me about the reason for this d

Re: [R] How do I add a variable to a text file?

On 19/03/2009 4:55 PM, Guillaume Filteau wrote: Hello all, I have a 2.0 GB dataset that I can't load into R, due to memory issues. The dataset itself is in a tab-delimited .txt file with 25 variables. I have a variable I'd like to add to the dataset. How do I do this? Best, Read a block, ad

[R] GLMM: non-linear fixed effect

Dear List, The formula expression in a glmm call is typically linear, which I believe is a requirement in order to preserve the assumption of odds proportionality. Is glmm capable of fitting a non-linear function of a covariate? In medical research a popular nonlinear function is the Emax model

Re: [R] simple "for loop" program for merging datasets?

You say you want to merge the two datasets on Country but the last step in you unsuccessful effort tries to match the results of the nested loop (necessarily with values of Measure2 from B) to B which has 3 rows. You never told us what the final results of this effort should be. but surely

Re: [R] Aggregating multiple columns

Try this: coef2 <- function(i) cov(x[i,"order"], x[i, "y"]) / var(x[i, "order"]) tapply(1:nrow(x), x$subject, coef2) On Thu, Mar 19, 2009 at 5:41 PM, Adam D. I. Kramer wrote: > Dear colleagues, > >        Consider the following data frame: > > x <- data.frame(y=rnorm(100),order=rep(1:10,10),subj

[R] Phase Spectrum

Hello, I have trouble to interpret the phase spectrum correctly. Here is the code (with different data) I used to get the phase spectrum: mfdeaths.spc <- spec.pgram(ts.union(mdeaths, fdeaths), spans = c(3,3)) plot(mfdeaths.spc, plot.type = "phase") I would like to know whether positive or nega

[R] Aggregating multiple columns

Dear colleagues, Consider the following data frame: x <- data.frame(y=rnorm(100),order=rep(1:10,10),subject=rep(1:10,each=10)) ...it is my goal to aggregate x to compute a linear effect of order for each subject. So, ideally, result would be a vector containing a single number f

[R] simple "for loop" program for merging datasets?

Dear R community, I would like to merge two datasets based on the categorical predictor “country”. Dataset A: Country Measure1 Afganistan1 Afganistan1 Russia 5 Poland 3 Poland 2 Dataset B: Country Measure 2 Russia

Re: [R] How to set R_PROFILE conditional on batch or interactive mode

Yes, thank you. That does the trick. Mark Lyman -Original Message- From: Rolf Turner [mailto:r.tur...@auckland.ac.nz] Sent: Thursday, March 19, 2009 2:42 PM To: Lyman, Mark Cc: r-help@r-project.org Subject: Re: [R] How to set R_PROFILE conditional on batch or interactive mode On 20/03

[R] Difference in client vs. server graphics defaults

Hello, I am having trouble with the difference between default graphic settings on my client machine and the instance of R on our company's server. I created a script locally that output graphs, but when I run it on the server the output graphs have titles running past the margins, legends imprope

[R] Difference between gam() and loess().

It seems that in general gam(y~lo(x)) # gam() from the gam package. and loess(y~x) give slightly different results (in respect of the predicted/fitted values). Most noticeable at the endpoints of the range of x. Can anyone enlighten me about the reason for this difference?

[R] How do I add a variable to a text file?

Hello all, I have a 2.0 GB dataset that I can't load into R, due to memory issues. The dataset itself is in a tab-delimited .txt file with 25 variables. I have a variable I'd like to add to the dataset. How do I do this? Best, Guillaume __ R-help@r-

Re: [R] How to set R_PROFILE conditional on batch or interactive mode

On 20/03/2009, at 9:17 AM, Lyman, Mark wrote: Thanks for the suggestion Rolf. Unfortunately, this uses the BATCH profile when I use for example R CMD INSTALL. I would like the BATCH profile to be used only when R CMD BATCH is used. Sorry, about the vagueness of my question. I have found a way

Re: [R] How to set R_PROFILE conditional on batch or interactive mode

Thanks for the suggestion Rolf. Unfortunately, this uses the BATCH profile when I use for example R CMD INSTALL. I would like the BATCH profile to be used only when R CMD BATCH is used. Sorry, about the vagueness of my question. I have found a way that approximates what I want to do. If I create a

[R] subtract values

Dear R-help I am using R version 2.6.2. I am trying to subtract specific values from a larger data frame. I feel this should be straightforward, but I am struggling. I have a dataframe "Bk" as follows: DateTime cumPrecip 01/01/2008 00:00 348 01/01/2008 01:00 348

Re: [R] function question

Not surprising at all -- expected! is.function(vector) ##TRUE -- Bert -Original Message- From: Wacek Kusnierczyk [mailto:waclaw.marcin.kusnierc...@idi.ntnu.no] Sent: Thursday, March 19, 2009 1:03 PM To: Bert Gunter Cc: Edna Bell; R help Subject: Re: [R] function question Bert Gunter

Re: [R] function question

Bert Gunter wrote: > Edna: > > Please make at least a minimal effort to answer such questions before > posting. > > is.vector(function(x)x) ## FALSE > as.vector(function(x)x) ## try it > > or even > > is.vector(plot) > > ... or even the surprizing is.vector(vector) vQ ___

Re: [R] nth root

x^(1/n) # ? exp(log(x)/n) # ? Either of those could be rolled into a function with the desired name: On Mar 19, 2009, at 11:21 AM, Martin Biuw wrote: Hi, Is there a function in R to calculate the nth root, similar to the MATLAB function NTHROOT()? Thanks, Martin Biuw David Winsemius, MD He

Re: [R] nth root

Matlab's Nthroot calculates the real nth root. For positive a, you can use a^(1/b); for negative a, b must be odd for the result to be real, and you can use -abs(a)^(1/b). So you could write: nthroot <- function(a,b) ifelse( b %% 2 == 1 | a >= 0, sign(a)*abs(a)^(1/b), NaN) This is a so-called V

Re: [R] zip.file.extract

Duncan Murdoch wrote: > On 3/19/2009 2:30 PM, Li, Xuejun wrote: >> Hello everyone, >> >> >> >> Can anyone give me some hint/suggestion for unzipping a file in R? The >> help file for "zip.file.extract" function seems to be difficult to >> understand. >> >> >> For example, the name of my zip f

Re: [R] subtract values

I think the subject heading should be "looping one day at a time". I am guessing that one of the zoo functions has already solved that problem but it's not a package that I have much (or any) depth. If you create a date vector and pass it through as.POSIXct(), you will get a series of "mi

[R] Matching failure in merge()

Hi! I've done a merging between 2 dataframes using merge(): delme <- merge(miDUNS50peqB,Bnomscodmunicipis,by.x="POBLACION",by.y="NOMMUNI",all.x=T,sort=F) After realizing some problems in the resulting dataset, I've found that the problem was that, in some cases, there was no match between the b

Re: [R] zip.file.extract

On 3/19/2009 2:30 PM, Li, Xuejun wrote: Hello everyone, Can anyone give me some hint/suggestion for unzipping a file in R? The help file for "zip.file.extract" function seems to be difficult to understand. For example, the name of my zip file is "refGene.txt.gz", which is located in

[R] How to rank risks and opportunities;Tornado chart plot

I would like to use Tornado charts/diagrams to rank risks and opportunities in a variable (for example, Cost) that depends on a few other variables which have risks using R and R-packages. I did some RSiteSearchs but did not find Tornado chart plots and how to rank risks and opportunities. Any idea

Re: [R] two plots side-by-side with persp3d(rgl)

On 3/19/2009 2:27 PM, Camarda, Carlo Giovanni wrote: Dear R-users, I would like to place two 3D plots side-by-side in a rgl-setting. It would nice to have something like "par(mfrow=c(1,2))" for basic plots, or an array framework for wireframe(lattice) (see example below). I only manage

Re: [R] How to set R_PROFILE conditional on batch or interactive mode

On 20/03/2009, at 5:46 AM, Lyman, Mark wrote: I am running R 2.8.1 on SUSE Linux. I would like to be able to set R_PROFILE to two different possibilities. One when R is called by R CMD BATCH and another any other time. Is this possible? I think so. In your .Rprofile have code like unto:

Re: [R] function question

Patrick Burns wrote: Duncan Murdoch wrote: On 3/19/2009 12:49 PM, Edna Bell wrote: Dear R Gurus: I read somewhere that functions are considered vectors. Is this true, please? Your question is a little ambiguous (you probably did read that, and probably someone does consider them to be vect

[R] Re : Randomly splitting a data frame in half

selected<-sample(nrow(data),19811) randomsample1 <- data[selected,] randomsample2 <- data[-selected,] # for non select But I think is good to have a variable that indicate selected and not selected cases in same data frame. You can try this selected<-rep(0,39622) selected[sample(1:39622,39622/2

[R] zip.file.extract

Hello everyone, Can anyone give me some hint/suggestion for unzipping a file in R? The help file for "zip.file.extract" function seems to be difficult to understand. For example, the name of my zip file is "refGene.txt.gz", which is located in "C:\Documents and Settings\Desktop\test1".

Re: [R] R editor that will work with Rcmdr

Thank you both, Daniel and Alain. I tried eclipse as well and it seems as if the combination with eclipse and StatEt was exactly what I was looking for. Regards, Jannis Am 08.03.2009 22:05, schrieb Tal Galili: Daniel, thank you very much for this posting!I felt similar "inspiration" as you ca

[R] nth root

Hi, Is there a function in R to calculate the nth root, similar to the MATLAB function NTHROOT()? Thanks, Martin Biuw [[alternative HTML version deleted]] __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE d

[R] two plots side-by-side with persp3d(rgl)

Dear R-users, I would like to place two 3D plots side-by-side in a rgl-setting. It would nice to have something like "par(mfrow=c(1,2))" for basic plots, or an array framework for wireframe(lattice) (see example below). I only managed to overlap two persp3d plots. My final idea would be

Re: [R] Randomly splitting a data frame in half

Thanks that makes perfect sense! Chris On 3/19/09 1:17 PM, Dimitris Rizopoulos wrote: well, you need to keep track of the rows you sampled, e.g., dat <- data.frame(x = rnorm(20), y = rnorm(20), w = rnorm(20)) ii <- seq_len(nrow(dat)) ind1 <- sample(ii, 10) ind2 <- ii[!ii %in% ind1] dat[ind1, ]

Re: [R] Randomly splitting a data frame in half

well, you need to keep track of the rows you sampled, e.g., dat <- data.frame(x = rnorm(20), y = rnorm(20), w = rnorm(20)) ii <- seq_len(nrow(dat)) ind1 <- sample(ii, 10) ind2 <- ii[!ii %in% ind1] dat[ind1, ] dat[ind2, ] I hope it helps. Best, Dimitris Christopher David Desjardins wrote: I

[R] Randomly splitting a data frame in half

I have a data frame in long format and I would like to randomly divide this data frame in half. The data frame consists of 39622 rows and I initially tried ... randomsample1 <- data[sample(nrow(data),19811), ] Where allows me to randomly select half of the rows and assign them to randomsample

Re: [R] find the index of the smallest or biggest number in a vetor or data.frame

Try this: # for the minimum a$price[which.min(a$n)] # for the maximum a$price[which.max(a$n)] On Thu, Mar 19, 2009 at 10:15 AM, tedzzx wrote: > > Dear R experts, > > How to find out the index of minimum or maxmum number in a vetor or > data.frame? > > For example, > > a= > > n price > 1

[R] sunflowerplot error

A sunflowerplot crossing two categorial variables with NAs fails: ### sample: start ### set.seed(20) a <- c(letters[1:4]) z <- c(letters[23:26]) fa <- factor(sample(rep.int(a, 1000), 100, replace=T), levels=a, ordered=T) fz <- factor(sample(rep.int(z, 1000), 100, replace=T), levels=z, ordered

Re: [R] Help with Scviews

Steve Sidney wrote: Many thanks Phillipe I was a little confused with web site and must have missed the fact that Komodo Edit has replaced SciViews. I will try and follow the instructions on your email and if I run into any more problems I will communicate with you. Some additional questio

Re: [R] function question

In mathematics that is the case. That is vector v = (v_1, v_2, ..., v_n) can be considered to be the function which maps k to v_k but there is no real relation to this in terms of R classes and objects. On Thu, Mar 19, 2009 at 12:49 PM, Edna Bell wrote: > Dear R Gurus: > > I read somewhere that

Re: [R] function question

Duncan Murdoch wrote: On 3/19/2009 12:49 PM, Edna Bell wrote: Dear R Gurus: I read somewhere that functions are considered vectors. Is this true, please? Your question is a little ambiguous (you probably did read that, and probably someone does consider them to be vectors), but I think the

[R] Re : function question

Functions and vectors are both objects in R, but a functions are not vectors try this foo<-function() print("I'm the foo function") is.vector(foo) as.vector(foo)  Justin BEM BP 1917 Yaoundé Tél (237) 99597295 (237) 22040246 De : Edna Bell À : r-h...@stat.mat

Re: [R] function question

Edna: Please make at least a minimal effort to answer such questions before posting. is.vector(function(x)x) ## FALSE as.vector(function(x)x) ## try it or even is.vector(plot) -- Bert Gunter Genentech -Original Message- From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-pr

Re: [R] using Search engine from help.start()

On 3/19/2009 12:57 PM, Erin Hodgess wrote: Dear R People: I ran into a weird situation...when using the help.start() on a Windows XP, (Firefox browser), the help files would not open. Has anyone else ran into that please? If so, how did you solve it, please? That's a Firefox bug. If you upg

Re: [R] function question

On 3/19/2009 12:49 PM, Edna Bell wrote: Dear R Gurus: I read somewhere that functions are considered vectors. Is this true, please? Your question is a little ambiguous (you probably did read that, and probably someone does consider them to be vectors), but I think the right answer is no: w

Re: [R] Import R-output into Java

#I used sink ?sink #Thomas Maxl18 schrieb: Hello, I want to import R-output via Rserve to Java, especially for the function ctree from the package party. Rserve is working properly. Yet, I only get the predictions with the Java code try{ RConnection c = new RConnection();

[R] using Search engine from help.start()

Dear R People: I ran into a weird situation...when using the help.start() on a Windows XP, (Firefox browser), the help files would not open. Has anyone else ran into that please? If so, how did you solve it, please? Thanks, Erin -- Erin Hodgess Associate Professor Department of Computer and

[R] How to set R_PROFILE conditional on batch or interactive mode

I am running R 2.8.1 on SUSE Linux. I would like to be able to set R_PROFILE to two different possibilities. One when R is called by R CMD BATCH and another any other time. Is this possible? Mark Lyman, Statistician Engineering Systems & Integration, ATK _

[R] function question

Dear R Gurus: I read somewhere that functions are considered vectors. Is this true, please? thanks Edna Bell __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/po

Re: [R] Run fortran from R

On Thu, 19 Mar 2009 16:30:18 +0100 Uwe Ligges wrote: > > > boem paukeslag wrote: > > Hi there, > > > > I've got a fortran simplex-routine that I would like to run from R. > > Is this possible? > > See the Writing R Extensions manual. > > Uwe Ligges > Hi Boem, Yes, you can. This is one of

Re: [R] changing the description of the x-axi

johnhj wrote: Hii, Is it possible to change the description ot the x-axis ? I need some barplots with high x-axis values. After plotting the graph, the x-axis values are shown in this form: 0e+00 1e+06 2e+06 3e+06 4e+06 5e+06 6e+06 Can I change this description to numerary values, like 53

Re: [R] Adding labels to heatmaps from image()

draja...@bidmc.harvard.edu wrote: Hi, I have been trying to add labels to the rows of a heatmap produced using image() function. It is simply not working. Here is what I did. A2Rplot.hclust(hcc,k=length(num),col.up="black",col.down=num,lty.up=2,lty.down=1,lwd.up=1,lwd.down=2,show.labels=FALS

Re: [R] permutations in R

Yes, it would have been better if I had mentioned in my original reply that it was an omegahat package. I found it by using help.search() (or ??) where that package was installed on the computer where I did the search. Looking at the help for setRepositories and packageStatus show one way to s

Re: [R] coefficient graph

Please read the posting guide. 1. We do not have the data and hence cannot reproduce. 2. Is this the most recent version of R and the packages that you are using? Uwe Ligges Sachi Ito wrote: Dear R list members, I'd like to make a graph of coefficients of the intercept, variable 1, and vari

Re: [R] How to keep attributes when dropping factor levels?

Dear Dieter, did you try it? I kew this post and below I extended my example. The result seems to be identical. Maybe that other attributes of the _DF_ don't get lost, but the label of the factor does still disappear. ff <- factor(substring("statistics", 1:10, 1:10), levels=letters) attribute

Re: [R] Run fortran from R

boem paukeslag wrote: Hi there, I've got a fortran simplex-routine that I would like to run from R. Is this possible? See the Writing R Extensions manual. Uwe Ligges _ [[elided Hotmail spam]] [[alternative HTML ver

Re: [R] Import R-output into Java

Maxl18 wrote: Hello, I want to import R-output via Rserve to Java, especially for the function ctree from the package party. Rserve is working properly. Yet, I only get the predictions with the Java code try{ RConnection c = new RConnection(); ... c.voidEval("modell <- c

[R] Import R-output into Java

Hello, I want to import R-output via Rserve to Java, especially for the function ctree from the package party. Rserve is working properly. Yet, I only get the predictions with the Java code try{ RConnection c = new RConnection(); ... c.voidEval("modell <- ctree(...)");

[R] Possible bug in charting: shaded region? (quantmod)

Dear Jeff and all, I think there is bug in plotting shaded regions: the region is one bar off on the left side. Can somebody please confirm this? By the way, I like this quantmod package a lot! Thanks for the great work! Regards, --Chang Lin Sent from my Verizon Wireless BlackBerry ___

[R] Run fortran from R

Hi there, I've got a fortran simplex-routine that I would like to run from R. Is this possible? _ [[elided Hotmail spam]] [[alternative HTML version deleted]] __ R-help@r-project

Re: [R] xlsReadWrite library

See: http://wiki.r-project.org/rwiki/doku.php?id=tips:data-io:ms_windows&s=excel On Thu, Mar 19, 2009 at 10:45 AM, Pascal Candolfi wrote: > Any idea why this library was removed and where could I find it for Windows > (only Unix in the Archive) ? > Much thanks in advance. > Pascal > >        [[al

Re: [R] find the index of the smallest or biggest number in a vetor or data.frame

?which On Thu, Mar 19, 2009 at 2:15 PM, tedzzx wrote: > > Dear R experts, > > How to find out the index of minimum or maxmum number in a vetor or > data.frame? > > For example, > > a= > > n   price > 1     50 > -2    100 > 0    200 > -1   300 > .. > > > I want to find out the row which the n

[R] xlsReadWrite library

Any idea why this library was removed and where could I find it for Windows (only Unix in the Archive) ? Much thanks in advance. Pascal [[alternative HTML version deleted]] __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listin

[R] find the index of the smallest or biggest number in a vetor or data.frame

Dear R experts, How to find out the index of minimum or maxmum number in a vetor or data.frame? For example, a= n price 1 50 -2100 0200 -1 300 .. I want to find out the row which the n is the smallest or largestest and extract the price. Thanks Ted -- View this messag

Re: [R] R with MKL

Both gcc and MKL are moving targets. So if you are using a more recent version of gcc than the one used to compile MKL then you are potentially in trouble. In this case, you still have the option of using Intel threads (libmkl_intel_thread) instead of the GNU thread library. In this case you must

Re: [R] How to keep attributes when dropping factor levels?

Heinz Tuechler gmx.at> writes: > > to drop unused factor levels two ways are outlined in R-help. In both > cases a label attribute is lost. Brian Ripley has posted this http://markmail.org/message/pl2odydwzv64v3u3 Dieter __ R-help@r-project.org ma

[R] Using unicode character for 'middle dot' in options(OutDec)

When drawing a graph, I'd like the unicode character 'middle dot' (or something else similar to \cdot in latex) to be used when writing numbers. Something like the following works for me: x <- 1:10 y <- runif(length(x)) par(las=1, bty='n') plot(x,y, ylim=c(0,1), yaxt='n') p <- pretty(y) axis(2,

Re: [R] Help with Scviews

Many thanks Phillipe I was a little confused with web site and must have missed the fact that Komodo Edit has replaced SciViews. I will try and follow the instructions on your email and if I run into any more problems I will communicate with you. Some additional questions. 1) Is it your int

[R] Generalized Extreme Value Distribution (LMOM package) and Frechet Distribution

Dear R helpers I have some data and through some other software, it is understood that I can fit the Frechet Distribution to it. However, I need to fit the distribution using R code only. I have searched many R packages and one R helper has suggested some sites too, but unfortunately paramete

Re: [R] yet another package building question please

Edna Bell wrote: Hello again. I'm trying to use package.skeleton to build my package. However, my package will contain a Fortran subroutine. Can you use package.skeleton with that subroutine, please or do you need to add it manually? The latter. Uwe Ligges Thanks again, Sincerely, Edn

Re: [R] run r scripts (tinn) via command

You just create the "superior" script as a set of 'source' statements reading in the files that you want to execute in order. On Thu, Mar 19, 2009 at 4:13 AM, thoeb wrote: > > The scripts are written with the TINN editor (***.r) and are organised within > a project. I want to create another scrip

Re: [R] problem with Dates

Try this (and read R News 4/1): Lines <- "6/4/1992 12:00:00 AM 2/13/1992 12:00:00 AM 6/19/1992 12:00:00 AM 2/11/1992 12:00:00 AM 6/22/1992 12:00:00 AM 10/3/1991 12:00:00 AM 5/12/1992 12:00:00 AM 7/11/1991 12:00:00 AM 2/6/1992 12:00:00 AM 10/4/1991 12:00:00 AM 1/31/1992 12:00:00 AM 6/9/1992 12:00:0

[R] How to keep attributes when dropping factor levels?

Dear All, to drop unused factor levels two ways are outlined in R-help. In both cases a label attribute is lost. The same happens, when using car:::recode. Is there a simple way to avoid losing attributes? Thanks, Heinz ## example ff <- factor(substring("statistics", 1:10, 1:10), levels=lett

Re: [R] Minimum cutsets

In the 'igraph' package, there are two (private) functions that can do this. I.e. install and load igraph, and call igraph:::find.all.min.cutsets() or igraph:::kCutsets() They might be quite slow if your graph is big. Check the source code for some comments. Best, Gabor On Thu, Mar 19, 2009 at

[R] correlation with ties & multiple comparisons problem

Hello! I am currently working on a big data set. There are different subjects that have been measured in about 40 different variables (in this case a measure of brain metabolism). Now I want to look at possbile significant correlations between those variables. As some of the variables are not

Re: [R] Updated R on Debian testing machine...

Hi Petar, I use Debian Lenny, these are the lines I added to my sources.list: # Add R-repository deb http://cran.at.r-project.org/bin/linux/debian etch-cran/ deb-src http://cran.at.r-project.org/bin/linux/debian etch-cran/ cheers, Paul Petar Milin wrote: Thanks! However, that is my point: it

Re: [R] the quote problem with readLines()

Hi Jim, Finally, I got the result I want, thanks a lot! Best, Dongyan jholtman wrote: > > Check out this reference: > > http://tolstoy.newcastle.edu.au/R/e2/help/07/02/9709.html > > > > On Wed, Mar 18, 2009 at 11:16 AM, Dongyan Song > wrote: >> >> Hi Jim, >> >> Thank you very much! I wil

[R] Minimum cutsets

The minCut function in RBGL package returns only a value or the minimum cut. I would be really greatful if any knows of any R function or package available for finding all minimal cut sets ( i.e., components whose failure will results in a network failure) between any given pair of vertices.

Re: [R] Simple Plot with Grid's Viewport

On 19 Mar 2009, at 07:22, Dieter Menne wrote: Gundala Viswanath gmail.com> writes: I have the following code that try to plot simple sinus curve into 2x2 grid in 1 page. But this code of mine create 4 plots in 1 page each. What's wrong with my approach? ... library(lattice) library(grid)

  1   2   >