Re: [R] Time complexity of functions in R

Hello Neha, You can try to measure those instructions time-complexiy by yourself. First, generate a benchmark dataset with increasing object size, i.e., set A. Have a look at how to use 'system.time' https://stat.ethz.ch/R-manual/R-devel/library/base/html/system.time.html Best, Mehmet On 24

Re: Alternative for numpy in Spark Mlib

You can use Breeze, which is part of spark distribution: https://github.com/scalanlp/breeze/wiki/Breeze-Linear-Algebra Check out the modules under import breeze._ On 23 May 2018 at 07:04, umargeek wrote: > Hi Folks, > > I am planning to rewrite one of my python

Re: [R] Comparing figures?

I suggest perceptual diff. You could write a wrapper around it. http://pdiff.sourceforge.net On Mon, 7 May 2018 16:49 Ramiro Barrantes, wrote: > Hello, > > I am working on tests to compare figures. I have been using ImageMagick, > which creates a figure

Re: [R-pkg-devel] Fwd: Collaboration Request: h2o R Package Function Cheatsheet

Dear Juan, A good start. A suggestion for versioning, instead of versioning on file names, maybe you can use git tags for release numbers. Github will create a release bundle with you release tag. RStudio has nice templates for cheatsheets too [1], I think you use their template and possibly

Re: [Rd] Why R should never move to git

On 31 January 2018 at 16:18, Barry Rowlingson wrote: >> > > Let the record also state that *gitlab* is an open source project and can be > downloaded and self-hosted, like gogs, but unlike github. Good to know. Nice one: https://github.com/gitlabhq/gitlabhq Best,

Re: [Rd] Best practices in developing package: From a single file

Dear Dr. Pfaff, Thank you for this, creating a package out of single file was my oriingal question, but not only creating and also maintaining it that way so R package is an artifact of the development process rather than "manually maintained" structure. I will have look at your sources. Best,

Re: [Rd] Why R should never move to git

ly event of GitHub / GitLab > closing business, that is _not_ such a big to any active project that > is hosted there. > > Gabor > > On Tue, Jan 30, 2018 at 11:07 PM, Suzen, Mehmet <mehmet.su...@gmail.com> > wrote: >> This might be off topic, but if R-core developmen

Re: [Rd] Why R should never move to git

This might be off topic, but if R-core development ever moves to git, I think it would make sense to have its own git service hosted by a university, rather than using github or gitlab. It is possible via https://gogs.io/ project. Just for the record. Best, -m

Re: [Rd] Best practices in developing package: From a single file

On 30 January 2018 at 21:31, Cook, Malcolm wrote: > > I think you want to see the approach to generating a skeleton from a single > .R file presented in: > > Simple and sustainable R packaging using inlinedocs > http://inlinedocs.r-forge.r-project.org/ > > I have not

Re: [Rd] Best practices in developing package: From a single file

entations, so richly democratised. Many regards, Mehmet On 30 January 2018 at 17:00, Suzen, Mehmet <mehmet.su...@gmail.com> wrote: > Dear R developers, > > I am wondering what are the best practices for developing an R > package. I am aware of Hadley Wickham's best practice > do

[Rd] Best practices in developing package: From a single file

Dear R developers, I am wondering what are the best practices for developing an R package. I am aware of Hadley Wickham's best practice documentation/book (http://r-pkgs.had.co.nz/). I recall a couple of years ago there were some tools for generating a package out of a single file, such as using

Re: [Rd] R CMD check warning about compiler warning flags

On 26 December 2017 at 00:00, Juan Telleria wrote: > Maybe I'm new, and forgive my ignorance, but maybe in the future (~ X years > from now) the R Project could be managed entirely from github, by doing I strongly disagree. Are you aware that github is a commercial company,

Re: [Rd] binary form of is() contradicts its unary form

On 30 November 2017 at 16:30, Iñaki Úcar wrote: > If you really believe that references should be needed to know what to > expect from a function call, then we work with different definitions A behaviour of a function call might be quite complex depending on the arguments

Re: [Rd] binary form of is() contradicts its unary form

On 30 Nov 2017 14:32, "Iñaki Úcar" wrote: >> >> Am I supposed to read every reference on a man page just to know what >> to expect from a function? >> > > If the reference is from John Chamber, you are supposed to read it. As a joke, it's funny. Not a joke. John Chambers

Re: [Rd] binary form of is() contradicts its unary form

On 30 November 2017 at 14:04, Iñaki Úcar wrote: > > Am I supposed to read every reference on a man page just to know what > to expect from a function? > If the reference is from John Chamber, you are supposed to read it. It is not always possible for maintainers to document

Re: [Rd] binary form of is() contradicts its unary form

On 30 November 2017 at 11:37, Iñaki Úcar <i.uca...@gmail.com> wrote: > 2017-11-30 3:14 GMT+01:00 Suzen, Mehmet <mehmet.su...@gmail.com>: >> My understanding is that there is no inconsistency. `is` does what it >> claims, from the documentation: >> >> ‘is’: Wi

Re: [Rd] binary form of is() contradicts its unary form

On 29 November 2017 at 21:45, Hervé Pagès wrote: > You're missing the point of my original post. Which is that > there is a serious inconsistency between the unary and binary > forms of is(). Maybe the binary form is right in case of My understanding is that there is no

Re: [Rd] binary form of is() contradicts its unary form

ot entirely sure. Best, -Mehmet On 29 November 2017 at 20:46, Hervé Pagès <hpa...@fredhutch.org> wrote: > Hi Mehmet, > > On 11/29/2017 11:22 AM, Suzen, Mehmet wrote: >> >> Hi Herve, >> >> I think you are confusing subclasses and classes. There is n

Re: [Rd] binary form of is() contradicts its unary form

Hi Herve, I think you are confusing subclasses and classes. There is no contradiction. `is` documentation is very clear: `With one argument, returns all the super-classes of this object's class.` Note that object class is always `data.frame` here, check: > class(data.frame()) [1] "data.frame"

Re: [R] run r script in r-fiddle

2017 at 15:09, Suzen, Mehmet <msu...@gmail.com> wrote: > On 31 October 2017 at 12:42, Martin Maechler <maech...@stat.math.ethz.ch> > wrote: >> Notably as I think it's been provided by a company that no >> longer exists under that name, and even if that'd be wrong, R-Fiddle

Re: [R] run r script in r-fiddle

On 31 October 2017 at 12:42, Martin Maechler wrote: > Notably as I think it's been provided by a company that no > longer exists under that name, and even if that'd be wrong, R-Fiddle > does not seem free software (apart from the R parts, I hope !). For the record,

Re: [R] run r script in r-fiddle

Note that, looks like r-fiddle runs R 3.1.2. __ R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide

Re: [R] run r script in r-fiddle

hler <maech...@stat.math.ethz.ch> wrote: >>>>>> Suzen, Mehmet <msu...@gmail.com> >>>>>> on Mon, 30 Oct 2017 11:16:30 +0100 writes: > > > Hi Frank, You could upload your R source file to a public > > URL, for example to git

Re: [R] run r script in r-fiddle

Hi Frank, You could upload your R source file to a public URL, for example to github and read via RCurl, as source do not support https as far as I know. Here is a working example. library('RCurl') tmatrix <- getURL("https://raw.githubusercontent.com/msuzen/isingLenzMC/master/R/isingUtils.R;)

Re: [R] Regarding Principal Component Analysis result Interpretation

Usually, PCA is used for a large number of features. FactoMineR [1] package provides a couple of examples, check for temperature example. But you may want to consult to basic PCA material as well, I suggest a book from Chris Bishop [2]. [1]

Re: RDD order preservation through transformations

Hi Johan, DataFrames are building on top of RDDs, not sure if the ordering issues are different there. Maybe you could create minimally large enough simulated data and example series of transformations as an example to experiment on. Best, -m Mehmet Süzen, MSc, PhD | PRIVILEGED

Re: RDD order preservation through transformations

On 14 September 2017 at 10:42, wrote: > val noTs = myData.map(dropTimestamp) > > val scaled = scaler.transform(noTs) > > val projected = (new RowMatrix(scaled)).multiply(principalComponents).rows > > val clusters = myModel.predict(projected) > > val result =

Re: RDD order preservation through transformations

of partitions in mapPartition? On 13 Sep 2017 19:54, "Ankit Maloo" <ankitmaloo1...@gmail.com> wrote: > > Rdd are fault tolerant as it can be recomputed using DAG without storing the > intermediate RDDs. > > On 13-Sep-2017 11:16 PM, "Suzen, Mehmet" <

Re: [Rd] establishing a Code of Conduct for R

On 13 September 2017 at 13:22, Brian G. Peterson wrote: > I am not an official representative of the R team, so this is only my > opinion. > Thank you. > It seems to me that you are trying to create a solution to a problem > which does not exist. I am not trying to create

Re: RDD order preservation through transformations

y a map operation can change sequence across a > partition as partition is local and computation happens one record at a > time. > > On 13-Sep-2017 9:54 PM, "Suzen, Mehmet" <su...@acm.org> wrote: > > I think the order has no meaning in RDDs see this post, specia

Re: RDD order preservation through transformations

I think the order has no meaning in RDDs see this post, specially zip methods: https://stackoverflow.com/questions/29268210/mind-blown-rdd-zip-method - To unsubscribe e-mail: user-unsubscr...@spark.apache.org

Re: [R] (no subject)

Hello David, As error message says you have a version dependency not satisfied. "error: Need GSL version >= 1.12". If you are using Ubuntu for example you could do; sudo apt-get install libgsl2 Or you can compile by yourself, I am sure there are people in LRZ can help you on this:) Best, Mehmet

Re: [R] comparition of occurrence of multiple variables between two dataframes

; > CV43 > > CV44 > > CV51 > > CV52 > > IN11 > > IN12 > > IN13 > > 4728 > > 0 > > 0 > > 0 > > 1 > > 0 > > 0 > > 0 > > 3 > > 0 > > 0 > > 0 > > 0 > > 0 > > 0 > > 0 &

Re: [R] comparition of occurrence of multiple variables between two dataframes

Do you have a simplified example with a code? It is not clear to me what do you mean by tree but if you refer to tree data structure, maybe you could change the data structure to tree (https://cran.r-project.org/web/packages/data.tree/vignettes/data.tree.html) and try to write comparison of two

Re: [Rd] Suggestion: Create On-Disk Dataframes

It is not needed. There is a large community of developer using SparkR. https://spark.apache.org/docs/latest/sparkr.html It does exactly what you want. On 3 September 2017 at 20:38, Juan Telleria wrote: > Dear R Developers, > > I would like to suggest the creation of a new

Re: [Rd] readLines() segfaults on large file & question on how to work around

Jennifer, Why do you try Sparkr? https://spark.apache.org/docs/1.6.1/api/R/read.json.html On 2 September 2017 at 23:15, Jennifer Lyon wrote: > Thank you for your suggestion. Unfortunately, while R doesn't segfault > calling readr::read_file() on the test file I

Re: [R] Block comment?

AFAIK block comment is not possible it needs to be implemented in R interpreter and defined in the parser.'If' solution is not elegant. On 2 September 2017 at 14:09, Uwe Ligges wrote: > > > On 02.09.2017 11:40, Christian wrote: >> >> I consider it quite worth

Re: Training A ML Model on a Huge Dataframe

...@gmail.com=22> > > On Wed, Aug 23, 2017 at 2:59 PM, Suzen, Mehmet <su...@acm.org> wrote: > >> It depends on what model you would like to train but models requiring >> optimisation could use SGD with mini batches. See: >> https://spark.apache.org/docs/latest/

Re: Training A ML Model on a Huge Dataframe

It depends on what model you would like to train but models requiring optimisation could use SGD with mini batches. See: https://spark.apache.org/docs/latest/mllib-optimization.html#stochastic-gradient-descent-sgd On 23 August 2017 at 14:27, Sea aj wrote: > Hi, > > I am

Re: [R] Directional Forecast

I suggest, you read: Forecasting: principles and practice from Hyndman-Athana­sopou­los https://www.otexts.org/fpp __ R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the

Re: [R] Nested cross validation with lapply

Hi Jesús, Do you have a code you tried without lapply? Why don't you post that here too? There are a couple of packages supporting nested CV; TANDEM, blkbox you may want to check their code. Also, `cvTools` package may help you to write one. On 7 August 2017 at 15:21, Jesús Para Fernández

Re: How can i remove the need for calling cache

On 3 August 2017 at 03:00, Vadim Semenov wrote: > `saveAsObjectFile` doesn't save the DAG, it acts as a typical action, so it > just saves data to some destination. Yes, that's what I thought, so the statement "..otherwise saving it on a file will require

Re: How can i remove the need for calling cache

On 3 August 2017 at 01:05, jeff saremi wrote: > Vadim: > > This is from the Mastering Spark book: > > "It is strongly recommended that a checkpointed RDD is persisted in memory, > otherwise saving it on a file will require recomputation." Is this really true? I had the

Re: [R] How export data set (available in the package) from R?

I also suggest you Hadley's optimized package for interoperating xls files with R: https://github.com/tidyverse/readxl https://cran.r-project.org/web/packages/readxl/index.html __ R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see

Re: [Rd] matrices with names

Not always, see what happens with lapply: > x<-matrix(12,1,1) > names(x)<-"one" > y<-matrix(1,1,1) > names(y)<-"one" > dput(lapply(x,`+`,e2=y)) structure(list(one = structure(13, .Dim = c(1L, 1L))), .Names = "one") >dput(lapply(x,`+`,e2=1)) structure(list(one = 13), .Names = "one") Prof. Ripley

Re: A tool to generate simulation data

I suggest RandomRDDs API. It provides nice tools. If you write wrappers around that might be good. https://spark.apache.org/docs/latest/api/scala/index.html#org.apache.spark.mllib.random.RandomRDDs$ - To unsubscribe e-mail:

Re: [R] R package for scorecard development

I suggest you to have a look at this R document: https://cran.r-project.org/doc/contrib/Sharma-CreditScoring.pdf On 28 June 2017 at 13:26, Nikhil Abhyankar wrote: > Hello all, > > Is there any R package that can develop a scorecard model for a binary > target

Re: [R] Nash equilibrium and other game theory tools implemented in networks using igraph or similar

Hello Chris, I was implying you are capable enough to implement it, while you have already identify a research paper. If there is no package out there, uploading to CRAN would help future user too. I am more than happy to help if you want to implement from scratch. Best, Mehmet On 27 June 2017

Re: [R] Nash equilibrium and other game theory tools implemented in networks using igraph or similar

Why don't you implement and uplad the package to CRAN? On 27 Jun 2017 17:45, "Chris Buddenhagen" wrote: Does anyone know of some code, and examples that implement game theory/Nash equilibrium hypothesis testing using existing packages like igraph/statnet or similar?

Re: Do we anything for Deep Learning in Spark?

There is a BigDL project: https://github.com/intel-analytics/BigDL On 20 June 2017 at 16:17, Jules Damji wrote: > And we will having a webinar on July 27 going into some more details. Stay > tuned. > > Cheers > Jules > > Sent from my iPhone > Pardon the dumb thumb typos :)

Re: [R] Latin Hypercube Sampling when parameters are defined according to specific probability distributions

No it is an R programming questions. Nelly specifically asked you: "how can I use your code to apply my model to each of the 50 rows of the data frame “tabLHS”?" __ R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see

Re: [R] (Somewhat?) Off topic: Containerization software

This is a nice summary addressing the same with R: https://arxiv.org/pdf/1410.0846.pdf On 30 May 2017 at 17:43, Bert Gunter wrote: > Folks: > > This is **off topic**, but I thought it might be informative to this > community. Consequently: please **no on list public

Re: [Rd] as.POSIXct behaviour

On 9 March 2017 at 01:29, Arunkumar Srinivasan wrote: > The time info is lost on the first index as well. And it happens *silently*. Yes, because it assumes homogeneous format on the entire vector. You may want to do two passes with different formats or run a

partition size inherited from parent: auto coalesce

Hello List, I was wondering what is the design principle that partition size of an RDD is inherited from the parent. See one simple example below [*]. 'ngauss_rdd2' has significantly less data, intuitively in such cases, shouldn't spark invoke coalesce automatically for performance? What would

partition size inherited from parent: auto coalesce

Hello List, I was wondering what is the design principle that partition size of an RDD is inherited from the parent. See one simple example below [*]. 'ngauss_rdd2' has significantly less data, intuitively in such cases, shouldn't spark invoke coalesce automatically for performance? What would

Re: [R] Social Network Simulation

Dear Professor Haenlein, Have you solved this issue yet? I found this eally interesting problem I was wondering if it is possible to wrapper "objective function" around igraph's 'sample_pa' and 'sample_smallworld'. If you have an example data set, I can have a look at this. Viele Gruesse aus

Re: [Rd] Inconsistency when naming a vector

There is no inconsistency. Documentation of `names` says ...value should be a character vector of up to the same length as x... In the first definition your character vector is not the same length as length of x, so you enforce NA by not defining value[2] x - 1:2 value-c(a) value[2] [1] NA where

Re: [gmx-users] Using thermostats to create temperature gradient in system

Hi Agnivo, Temperature gradient means non-equilibrium MD (NEMD) See notes from Prof. Martini: https://nanohub.org/resources/7582/download/Martini_L10_NonequilibruimMD.pdf What observable would you like to measure? Lets say you want to measure observable A. One procedure I can think of: 1.

Re: [gmx-users] Using thermostats to create temperature gradient in system

Hi Agnivo; Yes, the idea of freezing the solvent or the protein in many times is to sample the non-equilibrium thermal process. It will remain in the target temperature on average, over many fixed configurations obtained by different NVT runs (equilibrium runs). But you may need to run this many

Re: [R] lm() funtion

try lm.ridge from MASS package. __ R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented,

Re: [R] Multinomial Fitting Distrbution

mixtools package has mixture of Gaussian fitting, maybe that might help? __ R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide

Re: [R] Random Forest in Caret

Can you post your memory profile and codes? __ R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide

Re: [R] Cost-effectiveness Analysis using R

Yes. __ R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained,

Re: [R] any way to write sas7bdat with R

I didn't try this but there is an experimental package from Dr. Shotwell. http://cran.r-project.org/web/packages/sas7bdat/index.html if it can read, maybe you can modify to write as well? __ R-help@r-project.org mailing list -- To UNSUBSCRIBE and more,

Re: [R] Cost-effectiveness Analysis in R

Do you have specific example that you have tried to implement in R? Can you post your codes too? There are high quality package BCEA and BayesTree, that could be helpful; http://cran.r-project.org/web/packages/BCEA/index.html http://cran.r-project.org/web/packages/BayesTree/index.html

Re: [Numpy-discussion] IDE's for numpy development?

Spyder supports C. Thanks for correcting this. I wasn't aware of it. How was your experience with it? Best, -m ___ NumPy-Discussion mailing list NumPy-Discussion@scipy.org http://mail.scipy.org/mailman/listinfo/numpy-discussion

Re: [Numpy-discussion] IDE's for numpy development?

Hi Chuck, Spider is good. If you are coming from Matlab world. http://spyder-ide.blogspot.co.uk/ I don't think it supports C. But Maybe you are after Eclipse. Best, -m ___ NumPy-Discussion mailing list NumPy-Discussion@scipy.org

Re: [R] glmnet: converting coefficients back to original scale

This is interesting, can you post your lm.ridge solution as well? I suspect in glmnet, you need to use model.matrix with intercept, that could be the reason. -m __ R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see

Re: [R] weighted network centrality measures by network size

Hi Jenny, Have you tried igraph before? See, http://igraph.org/r/doc/ There are couple of centrality measures there. Best, -m On 6 August 2014 02:50, Jenny Jiang jiangyun...@y7mail.com wrote: Dear R-help, My name is Jenny Jiang and I am a Finance Honours research student from the

Re: [R] Using R to analyze multiple MRI studies

Did you inspect the CRAN view for Medical imaging? http://cran.r-project.org/web/views/MedicalImaging.html On 3 July 2014 17:09, moleps islon mole...@gmail.com wrote: I need to analyze multiple T1 contrast enhanced MRI studies from different patients. They are all in DICOM format. I see that

Re: [R-SIG-Finance] Kalman Filter Implementation in R

I suggest you to read the paper by Fernando Tusell from University of Basque Country, Kalman Filtering in R, JSS Vol. 39, Issue 2, Mar 2011 On 16 June 2014 11:21, Manuj Goel mg...@st-andrews.ac.uk wrote: Hello everyone, I am an applied statistics post-graduate student and am doing my

Re: [R] Logistic Regression

You might want to read this vignette: http://cran.r-project.org/web/packages/HSAUR/vignettes/Ch_logistic_regression_glm.pdf On 14 June 2014 19:53, javad bayat j.bayat...@gmail.com wrote: Dear all, I have to use Zelig package for doing logistic regression. How can I use Zelig package for

Re: [R] Defining default method for S3, S4 and R5 classes

There is a nice tutorial on this: http://adv-r.had.co.nz/OO-essentials.html For an in depth guide, have a look at the book from John Chambers, Software for data analysis programming with R. On 13 June 2014 12:20, Luca Cerone luca.cer...@gmail.com wrote: Dear all, I am writing a script

Re: [R] copula fitting

yes you can. On 7 June 2014 16:04, mudit gupta muditf...@gmail.com wrote: Hi guys, can i fit a copula to two marginal distributions with different sample size? like one has 2340 observations and other has 1912. thanks Mudit [[alternative HTML version deleted]]

Re: [R] copula fitting

Have you checked out 'copula' package? On 11 June 2014 00:36, Suzen, Mehmet msu...@gmail.com wrote: yes you can. On 7 June 2014 16:04, mudit gupta muditf...@gmail.com wrote: Hi guys, can i fit a copula to two marginal distributions with different sample size? like one has 2340

Re: [R] How to Optimize two functions together in R

Use defaul values initially, to see if you got reasonable results. See here for the details of nsga2 http://dx.doi.org/10.1109/4235.996017 On 19 May 2014 16:42, Mingxuan Han han...@purdue.edu wrote: I try to use NSGAII function in the mco but I am kind of confusing about the numbers of input

Re: [R] How to Optimize two functions together in R

This deals with the multi-objective optimisation. Try MCO and emoa packages. http://stats.stackexchange.com/questions/77580/optimization-of-multiple-objective-functions-with-constraints On 15 May 2014 17:47, Mingxuan Han han...@purdue.edu wrote: I am trying to minimize two functions with same

Re: [R] Calculating transition probabilties

This looks like this is your homework about Markov chains. not an R question actually. But have a look at the markovchain package from CRAN: http://cran.r-project.org/web/packages/markovchain/vignettes/an_introduction_to_markovchain_package.pdf On 13 May 2014 16:49, Baba Bukar bbu...@nda.edu.ng

Re: [R] problem in r-code

Wrong list. This is an R list not Bugs. You may want to consult Bugs materials: http://www2.mrc-bsu.cam.ac.uk/bugs/weblinks/webresource.shtml On 8 May 2014 11:36, thanoon younis thanoon.youni...@gmail.com wrote: dear all members is there anyone explain to me the code below and how can i

Re: [R] How can I make this nested loop faster?

Your code is not re-producable. Can you provide a working example using a standard dataset from R? But, you could first try to use compiler package, see ?enableJIT. Another option would be to use doMC/foreach packages if you can run your assignment in the nested loop in parallel, see %dopar%. On

Re: [R] uniform number

WTF? Is that a R package from you? On 5 May 2014 09:27, Rolf Turner r.tur...@auckland.ac.nz wrote: On 05/05/14 17:05, Ragia Ibrahim wrote: Dear group, How to generate uniform probability choosing p to be 2% and 5%, in separate trials for 100 times. No idea WTF you are talking about.

Re: [R] uniform number

That paper you cite is about Social networks. You may want to use igraph or sna packages On 5 May 2014 10:54, Ragia Ibrahim ragi...@hotmail.com wrote: thanks for replying in the following paper http://www.cs.cornell.edu/home/kleinber/kdd03-inf.pdf page 6 third paragraph the author

Re: [gmx-users] Silica Nanoparticles

You may try to use special potential designed for amorphous silica: A numerical investigation of the liquid–vapor coexistence curve of silica Yves Guissani and Bertrand Guillot J. Chem. Phys. 104, 7633 (1996) On 22 April 2014 17:12, Kazem Sepehrinia ksepehri...@gmail.com wrote: Hi Dear All,

Re: [R] inverse normal distribution function

You may want to read about generalized linear modelling and link functions for forming appropriate categorical variable/link function. See documentations in R: ?glm, ?family and ?inverse.gaussian. Also look at the original paper of Nelder, John; Wedderburn, Robert , it is available freely with

Re: [R] inverse normal distribution function

Not sure how would you do that but there is a package SEM on CRAN for structural equation models. On 20 April 2014 01:10, thanoon younis thanoon.youni...@gmail.com wrote: thank you so much Suzen i want to use bayesian analysis in structural equation models with ordered categorical data and i

Re: [R] operating website through R

This looks not so elegant, while normally data provider must have a nice accessing API, anyway, for example you can do this: myAdd -

Re: [R] operating website through R

You just need to pass the parameters on Giovanni_cgi.pl with action=ASCII+Output On 11 April 2014 17:19, eliza botto eliza_bo...@hotmail.com wrote: Dear Users of R, I wanted to operate certain slots of this website (http://disc2.nascom.nasa.gov/Giovanni/tovas/TRMM_V7.3B42_daily.2.shtml)

Re: [R] Sampling according to type

If I understood correctly, you need weighted sampling. Try 'prob' argument from 'sample'. For your example: n - 10 ntype - rbinom(n, 1, 0.5) myProbs - rep(1/10, 10) # equally likely myProbs[ which(ntype == 0)] - 0.75/7 # Divide so the sum will be 1.0 myProbs[ which(ntype == 1)] - 0.25/3

Re: [R] Sampling according to type

myProbs[ which(ntype == 0)] - 0.75/7 # Divide so the sum will be 1.0 myProbs[ which(ntype == 1)] - 0.25/3 Here of course you need to divide by number of 0s and 1s, 7 and 3 were was just an example. __ R-help@r-project.org mailing list

Re: [R] Shortest connected path in a matrix

You may want to check bioconductor packages doing graph algorithms. Maybe this one: http://www.bioconductor.org/packages/release/bioc/manuals/RBGL/man/RBGL.pdf See for example ?dijkstra.sp On 5 March 2014 18:44, McCloskey, Bryan bmcclos...@usgs.gov wrote: Here is some example data (hopefully

Re: [gmx-users] Computing melting point

Lindemann criterion might be easier. See For example, * Materials science: Melting from within Nature 413, 582-583 (11 October 2001) | doi:10.1038/35098169 Robert W. Cahn On 14 January 2014 15:35, Golshan Hejazi golshan.hej...@yahoo.com wrote: Hello everyone! I would like to compute the

Re: [R] Season's Greetings (and great news ... )!

I wouldn't blame R for floating-point arithmetic and our personal feeling of what 'zero' should be. options(digits=20) pi [1] 3.141592653589793116 sqrt(pi)^2 [1] 3.1415926535897926719 (pi - sqrt(pi)^2) 1e-15 [1] TRUE There was a similar post before, for example see:

Re: [R] Functional Programming patterns

Have you checked the r.lambda package of Brian Lee Yung Rowe ? http://cran.r-project.org/web/packages/lambda.r/index.html On 20 November 2013 10:02, mohan.radhakrish...@polarisft.com wrote: Hi, ' Not specific to 'R'. I search for patterns and found http://patternsinfp.wordpress.com/ which

Re: [R] Sending a matrix in an email

On 18 November 2013 05:37, Ira Fuchs irafu...@gmail.com wrote: I have a matrix which has colnames and I would like to send this matrix using sendmailR. How can I convert this simple matrix My 1 cent; In case of large objects or full session, suitable for attachment; RData might be more

Re: [R] Fitting arbitrary curve to 1D data with error bars

question has nothing to do with drawing error bars on a plot. What I want is to do a curve fit to a data with error bars. Best, e. On 14 Nov 2013, at 04:21, Suzen, Mehmet wrote: If you are after adding error bars in a scatter plot; one example is given below : #some example data

Re: [R] Fitting arbitrary curve to 1D data with error bars

If you are after adding error bars in a scatter plot; one example is given below : #some example data set.seed(42) df - data.frame(x = rep(1:10,each=5), y = rnorm(50)) #calculate mean, min and max for each x-value library(plyr) df2 - ddply(df,.(x),function(df)

Re: [R] computation of hessian matrix

On 1 November 2013 11:06, IZHAK shabsogh ishaqb...@yahoo.com wrote: below is a code to compute hessian matrix , which i need to generate 29 number of different matrices for example first You may consider using Numerical Derivatives package for that instead, see:

Re: [R] Revo R for Arima Implementation

On 28 October 2013 14:26, Anindita Chattopadhyay anindit...@mu-sigma.com wrote: We need to understand how we can implement this in Revo R. Most of the people here contribute to community of R not Revo R. I think it is unfair of you to request from this list to solve your Revo R issue.

Re: [R] Incorporate Julia into R

On 17 October 2013 15:38, Timo Schmid timo_sch...@hotmail.com wrote: I have some code in R with a lot of matrix multiplication and inverting. R can be very slow for larger matrices like 5000x5000. I have seen the new programming language Julia (www.julialang.org) which is quite fast in doing

Re: [R] RStudio with Revolution-R

On 15 October 2013 01:27, Maxim Linchits mlinch...@gmail.com wrote: Hello, Is it possible to use Revolution-R's multithreading capability with RStudio as the IDE? Apparently, RevoR is available for Ubuntu, Wrong list! But for reference:

  1   2   >