Re: [R] Sensitivity and Specificity

MKclass::perfMeasures(predict_testing, truth = testing$case, namePos = 1) should also work and computes 80 performance measures. Best Matthias Am 25.10.22 um 06:42 schrieb Jin Li: Hi Greg, This can be done by: spm::pred.acc(testing$case, predict_testing) It will return both sensitivity

Re: [R] What are the pros and cons of the log.p parameter in (p|q)norm and similar?

Response to 1You need the log version e.g. in maximum likelihood, otherwise the product of the densities and probabilities can become very small. Ursprüngliche Nachricht Von: r-help-requ...@r-project.org Datum: 04.08.21 12:01 (GMT+01:00) An: r-help@r-project.org Betreff:

Re: [R] density with weights missing values

. Best wishes, Matthias PS. Sorry for the HTML email. I’ve given up trying to fix such behavior. Von: Martin Maechler Gesendet: Dienstag, 13. Juli 2021 09:09 An: Matthias Gondan Cc: r-help@r-project.org Betreff: Re: [R] density with weights missing values >>>>>

Re: [R] density with weights missing values

You're right, of course. Extrapolating your argument a bit, the whole practice of na.rm is questionable, since there's always a reason for missingness (that is not in x and rarely elsewhere in the data)Best wishes Matthias  Ursprüngliche Nachricht Von: Jeff Newmiller Datum

Re: [R] density with weights missing values

?  Ursprüngliche Nachricht Von: Bert Gunter Datum: 12.07.21 16:25 (GMT+01:00) An: Matthias Gondan Cc: r-help@r-project.org Betreff: Re: [R] density with weights missing values The behavior is as documented AFAICS.na.rmlogical; if TRUE, missing values are removed from x. If FALSE anymissing values

Re: [R] density with weights missing values

a.rm=TRUE) [1] 2 Von: Richard O'Keefe Gesendet: Montag, 12. Juli 2021 13:18 An: Matthias Gondan Betreff: Re: [R] density with weights missing values Does your copy of R say that the weights must add up to 1? ?density doesn't say that in mine. But it does check. On Mon, 12 Jul 2021 at 22:42,

[R] density with weights missing values

the weights, the check for sum(weights) == 1 might trigger false positive warnings since the weights might not add up to 1 anymore Best wishes, Matthias [[alternative HTML version deleted]] __ R-help@r-project.org mailing list -- To UNSUBSCRIBE

[R] Embedded R: Test if initialized

, or should I choose another solution? Having said this, it may be a good idea to expose a function Rf_R_initialized that performs such a test. Thank you for your consideration. Best regards, Matthias *The use case is an R library that connects to swi-prolog and allows the „embedded“ swi

Re: [R] confidence intervals for the difference between group means

you could try: library(MKinfer) meanDiffCI(a, b, boot = TRUE) Best Matthias Am 04.08.20 um 16:08 schrieb varin sacha via R-help: Dear R-experts, Using the bootES package I can easily calculate the bootstrap confidence intervals of the means like in the toy example here below. Now, I am

Re: [R] ncol() vs. length() on data.frames

should have added: dim(x)[2L] -> length(x) Am 31.03.20 um 16:21 schrieb Prof. Dr. Matthias Kohl: Dear Ivan, if I enter ncol in the console, I get function (x) dim(x)[2L] indicating that function dim is called. Function dim has a method for data.frame; see methods("dim"). T

Re: [R] ncol() vs. length() on data.frames

) Hence, it calls length on the provided data.frame. In addition, some "magic" with .row_names_info is performed, where base:::.row_names_info function (x, type = 1L) .Internal(shortRowNames(x, type)) Best Matthias Am 31.03.20 um 16:10 schrieb Ivan Calandra: Thanks Ivan for the answer

Re: [R] How to rum Multiple ANOVA and Multiple T-test between the same groups?

Have a look at Bioconductor package genefilter, especially functions colttests and colFtests. Best Matthias Am 10.02.19 um 10:35 schrieb AbouEl-Makarim Aboueissa: Dear All: good morning *Re:* How to rum Multiple ANOVA and Multiple T-test between the same groups. Your help

Re: [R] Strange lazy evaluation of default arguments

So, I’ll stick with R. Still 25 years or so until retirement, but I’ll survive, even without crossreferenced default arguments. Best wishes, Matthias Von: S Ellison Gesendet: Dienstag, 5. September 2017 16:17 An: Matthias Gondan; r-help@r-project.org Betreff: RE: [R] Strange lazy evaluation of de

Re: [R] Strange lazy evaluation of default arguments

mu=0.53, sigma2=4.3^2) # instead of l=u And maybe also „in-place“ changes of values… Best regards, Matthias Von: William Dunlap Gesendet: Samstag, 2. September 2017 19:41 An: Rui Barradas Cc: Matthias Gondan; r-help@r-project.org Betreff: Re: [R] Strange lazy evaluation of default argume

[R] Strange lazy evaluation of default arguments

later place, and L is divided twice by 4.3. Is this behavior intended? It seems strange that the result depends on a debugging message. Best wishes, Matthias > sessionInfo() R version 3.4.1 (2017-06-30) Platform: x86_64-w64-mingw32/x64 (64-bit) Running under: Windows >= 8 x64 (build 9200)

[R] IF LOOOP

why... Thanks for the help r <- as.vector(lw) count=0 for (i in r) { if(i == 0.990956) { break } print(i) } Best Matthias [[alternative HTML version deleted]] __ R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see http

Re: [R] Paired sample t-test with mi.t.test

, not the variables themselves. Best Matthias Am 06.04.2017 um 18:32 schrieb Joel Gagnon: Dear all, It is my first time posting on this list so forgive me for any rookie mistakes I could make. I want to conduct t-tests on a dataset that has been imputed using the mice package: imput_pps <- mi

Re: [R] Quantiles with ordered categories

could be 1 if the data is from ordered categories. - Or both. It is probably a little thing to fix, but I lack the skills to do this myself. Best wishes, Matthias Von: Bert Gunter Gesendet: Dienstag, 14. März 2017 21:34 An: matthias-gon...@gmx.de Cc: r-help@r-project.org Betreff: Re: [R] Quantiles

[R] Quantiles with ordered categories

? Is there an elegant workaround? Thank you. Best wishes, Matthias [[alternative HTML version deleted]] __ 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

[R] summing up a matrix

Hello together, is it possible, to summing up a matrix? I have the following matrix at the moment: [,1] [,2] [,3][,4] [,5] [,6] 2016-1120200100 50 100 30

Re: [R] fPortfolio dont work in R 3.3.1

Andre, you need to make sure you got a C/C++ compiler as well as a Fortran compiler to compile a package from source that makes use of these language. Many R packages use one of those languages under the hood to speed things up. gcc / gfortran are common and free choices for such compilers.

[R] filter a data.frame in dependence of a column value

Hello togehter, i have short question, maybe anyone can help me. I have a data.frame like this one: NO ORDER 1 1530 for Mr. Muller (10.0 -> 11.2) 2 1799 for Mr Giulani 3 1888 for Mr. Marius (11.2 -> 12) I need a solution, which only contains the values in

[R] sustring in dependence of values

Hello togehter, maybe anyone can help me with a little problem. I have the following column in an data.frame (called TEST) VERSION 1abc (9.2 -> 10.0) 2def (10.2 -> 11.0) 3ghi (7.4 -> 8.4) 4jkl (10.2 -> 10.4) 5mno (8.1 -> 9.0) I now need the code for the following

[R] (no subject)

ependently). Thank you for the help!! If you need to create the plots you can simply change the "runmean_10yr_Nino_3_Index_17th_century" to any other vector. Best Matthias __ R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see htt

[R] hello everyone and happy new year

Matthias aa <- OHC_d[40:70] bb <- OHC_d[50:80] cc <- OHC_d[173:203] dd <- OHC_d[205:235] ee <- OHC_d[273:303] ff <- OHC_d[292:322] gg <- OHC_d[302:332] hh <- OHC_d[370:400] ii <- OHC_d[381:411] s3 <- data.frame(aa,bb,cc,dd,ee,ff,gg,hh,ii) y2 <- ggplot(s3,aes(x=time2,y=

Re: [R] How to install pkg when "not found" ?

-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code. -- Prof. Dr. Matthias Kohl www.stamats.de __ R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see https://stat.ethz.ch/mailman/listinfo/r

[R] Compute Regressions with R-Function

ing heat content. The error I got was the following: Error: unexpected symbol in "Trend <- function(x,y,z) { lm x" Thanks a lot for the help! Best Matthias [[alternative HTML version deleted]] __ R-help@r-project.

[R] ..lsmeans and coxme..

Dear list-eners, I run into the following problem when I want to get contrasts from a coxme model using the lsmeans package: A call to lsmeans on the coxme model throws the following error: Error in if (adjustSigma object$method == ML) V = V * object$dims$N/(object$dims$N - : missing value

Re: [R] calculate value in dependence of target value

-Ursprüngliche Nachricht- Von: David L Carlson [mailto:dcarl...@tamu.edu] Gesendet: Montag, 9. März 2015 16:08 An: Matthias Weber; r-help@r-project.org Betreff: RE: calculate value in dependence of target value It is very hard to figure out what you are trying to do. 1. All of the VALUEs

[R] calculate value in dependence of target value

Hello together, i have a litte problem. Maybe anyone can help me. I have to calculate a new column in dependence of a target value. As a example: My target value is 100.000 At the moment I have a data.frame with the following values. IDVALUE 1 111 2 125 3 13

[R] move date-values from one line to several lines

Hello together, i have a data.frame with date-values. What I want is a data.frame with a several lines for each date. My current data.frame looks like this one: ID FROM TOREASON 1 2015-02-27 2015-02-28Holiday 1 2015-03-15 2015-03-20Illness 2

[R] merge 2 data.frames in dependence of 2 values

Hello togehter, i have a little problem. Maybe anyone can help me. I have 2 data.frames, which look like as follows: First: NAMEMONTH BONUS 1 Andy 2014-10 100 2 Pete 2014-10200 3 Marc2014-10300 4 Andy2014-11400

[R] split a field in dependence of the semicolon

Hello togehter, i have a little problem, maybe anyone can help me. I have a data.frame, which look like this one: ID Members 1 1 ; abc; def; ghi 2 2 ; abc; 3 3 ;def; How can I create another column for each value between 2 semicolons? The

Re: [R] dotplot with library lattice

a lot. Best regards. Mat -Ursprüngliche Nachricht- Von: Jim Lemon [mailto:j...@bitwrit.com.au] Gesendet: Montag, 27. Oktober 2014 09:47 An: Matthias Weber Cc: r-help@r-project.org Betreff: Re: AW: [R] dotplot with library lattice On Mon, 27 Oct 2014 08:53:51 AM Matthias Weber wrote: Hi

Re: [R] dotplot with library lattice

- Von: Jim Lemon [mailto:j...@bitwrit.com.au] Gesendet: Samstag, 25. Oktober 2014 03:27 An: Matthias Weber Cc: r-help@r-project.org Betreff: Re: AW: [R] dotplot with library lattice On Fri, 24 Oct 2014 12:49:39 PM Matthias Weber wrote: I want always to see the Process of the green button

Re: [R] dotplot with library lattice

11:29 An: r-help@r-project.org Cc: Matthias Weber Betreff: Re: [R] dotplot with library lattice On Thu, 23 Oct 2014 05:57:27 PM Matthias Weber wrote: Hello together, i have a short question. Maybe anyone can help me to create a barplot in R with the package lattice. I have the following data

[R] dotplot with library lattice

Hello together, i have a short question. Maybe anyone can help me to create a barplot in R with the package lattice. I have the following data as ouput values and the following code: Data (d): KOST BudgetIST 1060 -2.18 0 1080

[R] sort a data.frame in a different way

Hello together, i have a little problem. Maybe anyone can help me. I have a data. frame which look like this one: 1000 1001 10021003 15 6 1211 24 3 81 What i need is a data.frame, which looks like this one. The Column

[R] strange behavior of the compare operator

I had a strange behavior of a function written a few days ago. I pointed the problem down to the following minimal example. can anyone explain why the following comparisons don't reply the samecorrect answer? Thanks for your reply! Matthias R version 3.1.1 (2014-07-10) -- Sock it to Me

[R] ..nlme: start values necessary even though using self-start function?..

Dear list. I am stuck with an subscript out of bounds error when I play with an nlme-example from the yellow Pinheiro/Bates book. It is about using nlme() using the formula interface for non-linear mixed models. The following code snippet is inspired by the book (in section 8.2), using the

Re: [R] Could someone recommend a package for time series?

. This is a pretty nice applied course that uses R for illustration: https://stat.ethz.ch/education/semesters/ss2012/atsa apart from this course, I think Prof. Rob Hyndman (who also has a blog) has a lot of useful stuff to say when it comes time series analysis with R best Matthias Bannert KOF

[R] Propensity Score Matching with Restrictions (Matching package)

): restrict the matching in way that only allows for matches within the same year, so that only propensity scores within the same year are considered for a respective match. Note that I don't what to estimate the years separately. Is there a way to do so? best regards matthias --- Matthias

Re: [R] C.D.F

/ (will appear soon) resp. a previous version at http://arxiv.org/abs/1006.0764 hth Matthias Am 11.08.2014 um 10:17 schrieb pari hesabi: Hello everybody, Can anybody help me to write a program for the CDF of sum of two independent gamma random variables ( covolution of two gamma distributions

Re: [R] legendre quadrature

you could use package distrEx: library(distrEx) GLIntegrate(function(x) x^2, lower = -1, upper = 1, order = 50) hth Matthias On 01.05.2014 09:43, pari hesabi wrote: Hello everybody I need to approximate the amount of integral by using legendre quadrature. I have written a program which

Re: [R] joint distribution

plot(x4) For more details on the computations see http://arxiv.org/abs/1006.0764 hth, Matthias On 22.04.2014 13:11, Ms khulood aljehani wrote: Hello i have two independent variablesx1 from normal (0,1)x2 from bernoulli (o.75) i need the density estimation of(b1*x1) + (b2*x2) where b1 and b2

Re: [R] Mistakes in date conversion for future date/time (POSIXct)

but not in Linux. So this seems to be a problem related to Windows. Kind regards Matthias   -Ursprüngliche Nachricht- Von: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org] Im Auftrag von David McPearson Gesendet: Sonntag, 6. April 2014 12:19 An: r-help@r-project.org Betreff: Re: [R

[R] Mistakes in date conversion for future date/time (POSIXct)

this with R 3.0.3, it showed the same problem. Thank you for your help! Kind regards, Matthias [[alternative HTML version deleted]] __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide

[R] computationally singular

help me with that. All the best, Matthias // [[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 http://www.R-project.org/posting-guide.html

Re: [R] MPICH2 Rmpi and doSNOW

Hi I have managed to install MPICH2 and Rmpi on my Windows 7 machine. I can also run the following code library(Rmpi) mpi.spawn.Rslaves() 4 slaves are spawned successfully. 0 failed. master (rank 0, comm 1) of size 5 is running on: MyMaster slave1 (rank 1, comm 1) of size 5

[R] MPICH2 Rmpi and doSNOW

Hi I have managed to install MPICH2 and Rmpi on my Windows 7 machine. I can also run the following code library(Rmpi) mpi.spawn.Rslaves() 4 slaves are spawned successfully. 0 failed. master (rank 0, comm 1) of size 5 is running on: MyMaster slave1 (rank 1, comm 1) of size 5 is running

Re: [R] Creating new instances from original ones

see function SMOTE in package DMwR hth Matthias On 31.03.2013 10:46, Nicolás Sánchez wrote: I have a question about data mining. I have a dataset of 70 instances with 14 features that belong to 4 classes. As the number of each class is not enough to obtain a good accuracy using some classifiers

Re: [R] create bar chart with different totals in a bar

1 0 2 1 02 Thanks for your help Von: John Kane [mailto:jrkrid...@inbox.com] Gesendet: Freitag, 8. März 2013 16:29 An: Matthias Weber Betreff: RE: AW: [R] create bar chart with different totals in a bar The image did not come

Re: [R] integrate function

use pmin instead of min. hth Matthias On 12.02.2013 16:41, dan wang wrote: Hi All, Can any one help to explain why min and max function couldn't work in the integrate function directly. For example, if issue following into R: integrand - function(x) {min(1-x, x^2)} integrate(integrand, lower

Re: [R] problem with any function

at the command prompt to see what type of unintended things happen. Actually, 1:4 + 1:2 is running silently. Best wishes, Matthias Am 17.11.2012 19:56, schrieb Rui Barradas: Hello, First of all, that's not an error message, it's a warning. As for your condition, it's too complicated. You

Re: [R] survreg gompertz

On 15/11/12 21:22, David Winsemius wrote: On Nov 15, 2012, at 5:38 AM, Matthias Ziehm wrote: Hi all, Sorry if this has been answered already, but I couldn't find it in the archives or general internet. A Markmail/Rhelp search on: gompertz survreg ...brings this link to a reply by Terry

[R] survreg gompertz

distribution. Many thanks! Matthias __ 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/posting-guide.html and provide commented, minimal, self-contained, reproducible

Re: [R] Two-way Random Effects with unbalanced data

.. Matthias. On Mon, Oct 29, 2012 at 10:42 PM, asafwe as...@wharton.upenn.edu wrote: Hi there, I am looking to fit a two-way random effects model to an *unblalanced* layout, y_ijk = mu + a_i + b_j + eps_ijk, i=1,...,R, j=1,...,C, k=1,...,K_ij. I am

[R] ylim with only one value specified

, Inf)[as a replacement for NULL/ autoselection] and ylim=c(Inf, -Inf)[autoselection, reversed y axis] should be handled correctly. I would find such a feature useful. Do you think it would interfere with other functions? Thank you for your consideration. Best wishes, Matthias Gondan

Re: [R] transform RTs

you might to do something like library(SuppDists) t = runif(100, 100, 500) # original RT t_IG = qinvGauss(ecdf(t)(t)-0.5/length(t), 1, 16) plot(density(t_IG)) but what is the purpose of it? Usually reaction times are thought to follow a certain kind of distribution (e.g. an inverse Gaussian

[R] Font size in geom_dl (using ggplot2)

Hey everyone, I am an R-newby... so sorry for bothering you with simple-to-solve questions;) I have the following issue: trying to add labels to my scatterplots (with geom_dl in ggplot2). Everything works fine, but after checking every resource I do not find a way to change the font size of my

Re: [R] What makes R different from other programming languages?

My vote: 1. Symbolic function arguments: fn = function(a, b) { a/b } fn(b=10, a=2) 2. Names for elements of a vector and matrices v = c(a=1, b=2) v['a'] = v['a'] * 2 same for matrices 3. about 10,000 user-contributed packages on CRAN 4. weird things like a = numeric(10) a[1:10] = 1:2

Re: [R] Convolve 2 uniform distributed variables

take a look at the distr package library(distr) U - Unif() U2 - U + U # or U2 - convpow(U, 2) plot(U2) # or curve(d(U2)(x), from = 0, to = 2) Best Matthias On 24.05.2012 08:29, c...@mail.tu-berlin.de wrote: Hi, I want to convolve 2 uniform distributed [0,1] variables, so that I get

Re: [R] ggplot2 - geom_bar

On 23.04.2012 20:55, Brian Diggs wrote: On 4/23/2012 9:24 AM, Matthias Rieber wrote: Hello, I've some problem with the ggplot2. Here's a small example: [...] Is it wrong to use geom_bar with that kind of data? I could avoid this issue when I cast the data.frame, but I like to avoid

Re: [R] Compare String Similarity

you should also look at Bioconductor Package Biostrings hth Matthias Am 19.04.2012 18:01, schrieb R. Michael Weylandt: Though if you do decide to use Levenstein, it's implemented here in R: http://finzi.psych.upenn.edu/R/library/RecordLinkage/html/strcmp.html I'm pretty sure this is in C code

Re: [R] Sweave UFT8 problem

try: Sweave(sim_pi.Rnw, encoding = utf8) Best, Matthias On 15.04.2012 11:41, Philippe Grosjean wrote: Hello, Have you tried to put that command in a comment: %\usepackage[utf8]{inputenc} I haven't tested it in this particular case, but it works in some other situations. Best, Philippe

[R] line width in legend of interaction.plot

leg.lwd Best wishes, Matthias Gondan __ 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/posting-guide.html and provide commented, minimal, self-contained

[R] Reference for dataset colon (package survival)

Dear R team, dear Prof. Therneau, library(survival) data(colon) ?colon gives me only a very rudimentary source (only a name). Is there a possibility to get a reference to the clinical trial these data are taken from? Many thanks in advance. With best wishes, Matthias Gondan

Re: [R] calculate quantiles of a custom function

Dear Gerhard, you could also use package distr; e.g. library(distr) ## use generating function AbscontDistribution D - AbscontDistribution(d = function(x) dbeta(x, 2, 6) + dbeta(x,6,2), low = 0, up = 1, withStand = TRUE) ## quantiles q(D)(seq(0,1,0.1)) Best Matthias On 03.01.2012 19:33

[R] Fwd: Re: Poisson GLM using non-integer response/predictors?

, your example seems incorrect, because it is a predictor, not the dependent variable tnoise_sqrt ~ lengthfeeding_log Best wishes, Matthias Am 30.12.2011 16:29, schrieb Lucy Dablin: Great lists, I always find them useful, thank you to everyone who contributes to them. My question

Re: [R] Fwd: Re: Poisson GLM using non-integer response/predictors?

Hi Ben, Thanks for clarifying this, I used a misleading word, model the observation time sounds as if observation time were the dependent variable - which it is not, of course, instead, in the scenario described, the parrot counts are modeled. Best wishes, Matthias Am 30.12.2011 20:50

[R] arrow egdes and lty

, Matthias -- __ 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/posting-guide.html and provide commented, minimal, self-contained, reproducible code.

[R] image problem in 2.13.1

and save as metafile causes a segmentation fault. Everything worked fine on 2.13.0 Best wishes, Matthias -- __ 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/posting

Re: [R] Defining functions inside loops

Dear Eduardo, try: f - function(x){} s - 0.2 body(f) - substitute({x^2 + s}, list(s = s)) Best, Matthias On 15.02.2011 16:53, Eduardo de Oliveira Horta wrote: Thanks... but I guess I didn't make myself clear. What I was trying to do was precisely to store inside the function the number

Re: [R] Package distr and define your own distribution

Dear Gabriel, it is possible. You can define a new class or just an object of an already existing class. Did you look at: library(distr) vignette(newDistributions) as well as ?AbscontDistribution ?DiscreteDistribution Please let me know if you have further questions! Best Matthias

[R] JGR and font antialiasing

installation? Has someone else noticed that issue? Regards, Matthias __ 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/posting-guide.html and provide commented

[R] random sequences for rnorm and runif

0.8830174 Best wishes, Matthias Gondan -- GMX DSL Doppel-Flat ab 19,99 Euro/mtl.! Jetzt mit gratis Handy-Flat! http://portal.gmx.net/de/go/dsl __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting

Re: [R] [OT] R on Atlas library

/libptf77blas.a /usr/local/lib/libatlas.a -lpthread make make install Using your benchmark example, all CPUs seem to be used. It is now working at such a high speed that I could not entirely be sure, top refreshes only every 5 s. Cheers Matthias Am 09.08.2010 11:51, schrieb Allan Engelhardt: I don't

Re: [R] optimization subject to constraints

try this (package Rsolnp) library(Rsolnp) g- function(x) { return(x[1]^2+x[2]^2) } # constraint f- function(x) { return(x[1]+x[2]) } # objective function x0 = c(1, 1) solnp(x0, fun=f, eqfun=g, eqB=c(1)) Am 10.08.2010 14:59, schrieb Gildas Mazo: Thanks, but I still cannot get to

Re: [R] optimization subject to constraints

try command solnp in package Rsolnp Am 09.08.2010 18:56, schrieb Dwayne Blind: Hi ! Why not constrOptim ? Dwayne 2010/8/9 Gildas Mazogildas.m...@curie.fr Dear R users, I'm looking for tools to perform optimization subject to constraints, both linear and non-linear. I don't mind which

[R] [OT] R on Atlas library

. Searching around the internet was not very encourageing. Some people wrote that it is not so simple to have Atlas fully exploit a multicore computer. I hope this is not too unspecific. Best wishes, Matthias -- Dr. rer. nat. Matthias Gondan Institut für Psychologie Universität Regensburg D-93050

Re: [R] how to generate a random data from a empirical distribition

/distr.pdf Best, Matthias On 27.07.2010 10:37, Dennis Murphy wrote: Hi: On Mon, Jul 26, 2010 at 11:36 AM, xin weixin...@stat.psu.edu wrote: hi, this is more a statistical question than a R question. but I do want to know how to implement this in R. I have 10,000 data points. Is there any way

Re: [R] Robust Asymptotic Statistics (RobASt)

(as in the case of package RobLox) are unfortunately not yet implemented. So, if you are willing to accept these drawbacks I will be glad to give you a hand to get our estimators to work for your data. Best, Matthias On 06.06.2010 06:47, Alex Weslowski wrote: Hi all, Other than: http://www.stamats.de

[R] Problem with format(,%G-%V) - Segfault

- format(test.data, format=%G-%V) Error: segfault from C stack overflow this happens with a self compiled R version and the debian packages. The same operation works with R 2.9.x and R 2.10.x I've attached the backtrace. kind regards, Matthias Core was generated by `/opt/R/lib64/R/bin/exec/R

Re: [R] Nonparametric generalization of ANOVA

Brunner et al. available on his website http://www.ams.med.uni-goettingen.de/de/sof/ld/index.html But they seem not to be working with current R versions. Best regards, Matthias Gondan -- Sicherer, schneller und einfacher. Die aktuellen Internet-Browser - jetzt kostenlos herunterladen! http

[R] Accessing named elements of a vector

, vx$a is not functional. Would it break existing compatibility if the $ would be allowed to access elements of a vector, as well? Best regards, Matthias __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read

[R] [R-pkgs] new package RFLPtools

sequence alignment. To get a first impression try: vignette(RFLPtools) as well as library(RFLPtools) example(RFLPplot) example(RFLPrefplot) Best regards, Fabienne Alexandra Matthias -- Dr. Matthias Kohl www.stamats.de ___ R-packages mailing list r-packa

[R] [R-pkgs] new version of RobASt-family of packages

in ExamplesEstimation.R in folder scripts Best Peter Matthias -- Dr. Matthias Kohl www.stamats.de ___ R-packages mailing list r-packa...@r-project.org https://stat.ethz.ch/mailman/listinfo/r-packages __ R

[R] [R-pkgs] new version of distr-family of packages

an example with framed code in vignette Best Peter Matthias -- Dr. Matthias Kohl www.stamats.de ___ R-packages mailing list r-packa...@r-project.org https://stat.ethz.ch/mailman/listinfo/r-packages __ R-help@r

[R] plot filled.contour over continent map

of Europe, without the contourplot. When commenting the map statements I can see the contourplot. So I am doing something wrong, but I really have no idea what? Anybody could help me out here? Thanks in advance, Matthias - Department of Earth

[R] COURSE: Introduction to Metabolomics and biomarker research

Introduction to Metabolomics and biomarker research. The course will take place in Innsbruck, 16th to 17th November 2009, and will be held in German. For more details see: http://www.gdch.de/vas/fortbildung/kurse/fortbildung2009.htm#1175 -- Dr. Matthias Kohl www.stamats.de

Re: [R] programming to calculate variance

I think it should be var(y[i-3:i-1,]) instead of var(x[i-3:i-1,]) otherwise the values of the vector are overwritten Best wishes, Matthias marlene marchena schrieb: Dear R-user Suppose I have the following data y=c(2,1,5,8,11,3,1,7,50,21,33,7,60) x=data.frame(y) for(i in 4:nrow(x

[R] Example data for analysis of a highly pseudoreplicate mixed-effects experiment

give me some advice on this problem ? Or should I just stick to reducing the data from each condition to a single parameter, as explained in my first email below? Thank you again! Matthias Gralle wrote: Hello everybody, I have been trying for some weeks to state the correct design of my

Re: [R] Example data for analysis of a highly pseudoreplicate mixed-effects experiment

, how do I correctly account for pseudoreplication, avoiding the artificially high number of df ? Thank you, Matthias Matthias Gralle wrote: Hello everybody, it may be better to have sample data. I have provided data with less levels of gene and day and only ca. 400 data points per condition

[R] Analysis of a highly pseudoreplicate mixed-effects experiment

) on Ubuntu 8.04 on a linux 2.6.24-24-generic kernel on different Intel systems. I am using the lme4 that came with R 2.8.0. Thank you very much for your time! -- Matthias Gralle, PhD Dept. Evolutionary Genetics Max Planck Institute for Evolutionary Anthropology Deutscher Platz 6 04103 Leipzig

[R] Test for stochastic dominance, non-inferiority test for distributions

wishes, Matthias PS. This one would be ok, as well: F(t) G(t), for all t null hypothesis: F(t) = G(t), for some t. -- GRATIS für alle GMX-Mitglieder: Die maxdome Movie-FLAT! __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman

Re: [R] help with median for each row

if your matrix has many rows you might want to consider rowMedians from Bioconductor package Biobase. hth, Matthias Greg Hirson schrieb: Edward, See ?apply x = matrix(rnorm(100), ncol = 10) apply(x, 1, median) Hope this helps, Greg Edward Chen wrote: Hi, I tried looking through google

Re: [R] how to change the quantile method in bwplot

for teaching purposes I wrote a corresponding function; cf. qbxp.stats (as well as qboxplot ...) in package MKmisc. hth, Matthias Deepayan Sarkar schrieb: On Tue, Jul 21, 2009 at 7:47 AM, Jun Shenjun.shen...@gmail.com wrote: Uwe, Thank you for your reply. I am still not very clear about

Re: [R] trouble using optim for maximalisation of 2-parameter function

), llik.nor, x=x, control = list(fnscale = -1)) hth, Matthias Anjali Sing schrieb: Hello, I am having trouble using optim. I want to maximalise a function to its parameters [kind of like: univariate maximum likelihood estimation, but i wrote the likelihood function myself because of data issues

Re: [R] implementing Maximum Likelihood with distrMod when only the PDF is known

and it would of course be better to have a valid definition for r, d, p and q. Best, Matthias guillaume.martin schrieb: Dear Mathias, That's pretty amazing, thanks a lot ! I'll have to look all this through because I don't easily understand why each part has to be set up, in particular

Re: [R] implementing Maximum Likelihood with distrMod when only the PDF is known

Dear Guillaume, thanks for your interest in the distrMod package. Regarding your question I took up your example and put a file under: http://www.stamats.de/distrModExample.R Hope that helps ... Don't hesitate to contact me if you have further questions! Best, Matthias guillaume.martin

  1   2   >