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 and

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: R-help

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 >>>>> Matthias Gonda

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 Newmill

Re: [R] density with weights missing values

thing?  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 v

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

[R] density with weights missing values

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

substitute, 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“

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 will be

Re: [R] Strange lazy evaluation of default arguments

happen? [it’s documented behavior, but still…] 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.

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

er 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

... 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

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

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 2016-12

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. D

[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 br

[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 resul

[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" ?

osting guide http://www.R-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.e

[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.org ma

[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 valu

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

t regards. Mat -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. Al

[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 1330

[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 Se

[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 resul

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

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

Re: [R] dotplot with library lattice

Jim Lemon [mailto:j...@bitwrit.com.au] Gesendet: Freitag, 24. Oktober 2014 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 m

[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 9

[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 na

[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 same"correct" answer? Thanks for your reply! Matthias R version 3.1.1 (2014-07-10) --

[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 Oran

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)

d like to do (in addition): 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 r

Re: [R] C.D.F

g/v59/i04/ (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 dist

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

t;- 3*x1 + 2*x2 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) whe

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

ot 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: R

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

quot;2038" "2039" "2040" "2041" "2042" "2043" "2044" "2045" "2046" "2047" "2048" "2049" [91] "2050" "2051" "2052" "2053" "2054&qu

[R] computationally singular

. I would be thankful if you could 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 htt

[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 runnin

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

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(i

Re: [R] problem with "any" function

+ 1:2 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 t

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

[R] survreg & gompertz

Same questions for the gompertz-makeham 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 commente

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

.. Matthias. On Mon, Oct 29, 2012 at 10:42 PM, asafwe 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

[R] ylim with only one value specified

lim=c(-Inf, 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, Matth

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 dis

[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 lab

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

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 thi

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

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 si

[R] line width in legend of interaction.plot

interaction.plot source: legend(xleg, yleg, legend = ylabs, col = col, pch = if (type %in% c("p", "b")) pch, lty = if (type %in% c("l", "b")) lty, bty = leg.bty, bg = leg.bg) <- here I woul

[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 Mat

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.

[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. > > >

[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

"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-p

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

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 Ma

[R] JGR and font antialiasing

em with my 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

[R] random sequences for rnorm and runif

yrnorm(4)) [1] 0.2875775 0.7883051 0.4089769 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/

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

l/lib/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

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 sol

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 Mazo Dear R users, I'm looking for tools to perform optimization subject to constraints, both linear and non-linear. I don't mind which algorithm may be used, m

[R] [OT] R on Atlas library

be helpful. 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 Regensbu

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 wei 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 to generate a empirical

Re: [R] Robust Asymptotic Statistics (RobASt)

clearly reduce computation time (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 Weslowsk

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

] 13888 14166 14188 14189 14189 ... > result <- 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 rega

Re: [R] Nonparametric generalization of ANOVA

are some R scripts from 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 koste

[R] Accessing named elements of a vector

ter of taste, but I like the $ very much. Unfortunately, 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.e

[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] [R-pkgs] new version of RobASt-family of packages

ce check time ... - some minor corrections 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.

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

vignette: * included an example with escape sequences in vignette * included an example with framed code in vignette Best Peter Matthias -- Dr. Matthias Kohl www.stamats.de ___ R-packages mailing list r-packa.

[R] plot filled.contour over continent map

range of lat/lon values. If I run the above-mentioned, I just get the 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

[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

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

4.6742 0.0309 In the original data set, denDF is >200 000. Once again, 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

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

it is the wrong formula. Can anybody 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

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

8-10-20) 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

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

some t. Best 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.eth

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 Shen wrote: Uwe, Thank you for your reply. I am still not very clear about the meanings of the

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

m optim(c(1,5), 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 mys

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

his is rather a quick and dirty solution 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 unders

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.m

  1   2   >