Re: [R] Undocumented behaviour of diag when replacing the diagonal of a matrix?

i, i)] <- value } x I fiddled a bit around with expressions, but to no avail. Could you (or anybody else) maybe still shed some more light on this? Thx a lot in advance! Best -- Gerrit --------- Dr. Gerrit Eichner

Re: [R] Undocumented behaviour of diag when replacing the diagonal of a matrix?

-- Gerrit - Dr. Gerrit Eichner Mathematical Institute, Room 215 gerrit.eich...@math.uni-giessen.de Justus-Liebig-University Giessen Tel: +49-(0)641-99-32104 Arndtstr. 2, 35392 Giessen, Germany https://www.uni-giessen.de/ma

[R] Undocumented behaviour of diag when replacing the diagonal of a matrix?

ng the code of `diag<-` line by line, it throws errors. So why does it work? diag(A[-1, ]) <- 2; A diag(A[-5, -1]) <- 3; A diag(A[-5, -(1:2)]) <- 4; A Any ideas? TIA and best regards -- Gerrit -----

Re: [R] Puzzled by results from base::rank()

- Dr. Gerrit Eichner Mathematical Institute, Room 215 gerrit.eich...@math.uni-giessen.de Justus-Liebig-University Giessen Tel: +49-(0)641-99-32104 Arndtstr. 2, 35392 Giessen, Germany http://www.uni-giessen.de/eichner

Re: [R] Get data from a list of data frames

on4 = Station4, Station5 = Station5) select.stat <- paste0("Station", c(1, 4, 5)) select.sens <- "thermometer" sapply(Total[select.stat], function(x, sens) x$code[x$sensor == sens], sens = select.sens) Hth -- Gerrit -

Re: [R] How to create density ellipses with R

Hi Paul, take a look at base R's function contour (and the Examples section of its help page), and follow maybe the hints under "See also". Hth -- Gerrit --------- Dr. Gerrit Eichner Mathematical

Re: [R] Replacing certain rows with values from a different column

Hi, James, if I understand you correctly, maybe, with(firstdf, ifelse(Province %in% seconddf$Country, Province, Country) ) does what you want? Hth -- Gerrit - Dr. Gerrit Eichner

Re: [R] Plotting Coxph model with an interaction.

xph(Surv(futime, fustat) ~ age+rx+age*rx, data = ovarian) plot(survfit(coxfit, newdata=data.frame(age=60,rx=2))) Or do I missunderstand your question? Hth -- Gerrit --------- Dr. Gerrit Eichner Mathematical Institut

Re: [R] Group by and duplicate a value/dplyr

Hello, Elahe, you were, of course, supposed to insert my suggested code-snippet into you code and test it therein ... Regards -- Gerrit - Dr. Gerrit Eichner Mathematical Institute, Room 212 gerrit.eich

Re: [R] Group by and duplicate a value/dplyr

Homework? Try maybe mutate(MinValue = min(Value[Value != 0]) ) or mutate(MinValue = sort(unique(Value))[2]) Hth -- Gerrit - Dr. Gerrit Eichner Mathematical Institute, Room 212 gerrit.eich...@math.uni

Re: [R] cox.zph

Yes. :-) Best regards -- Gerrit - Dr. Gerrit Eichner Mathematical Institute, Room 212 gerrit.eich...@math.uni-giessen.de Justus-Liebig-University Giessen Tel: +49-(0)641-99-32104 Arndtstr. 2

Re: [R] Calculating column differences

Dear Jeff, read diff's help page, and you'll find out what is wrong with your expectation. What do think diff(df$Score) should give for the first element in df$Score?? Hth -- Gerrit - Dr. Gerr

Re: [R] Color in stripchart

ur problem using stripchart, but you may want to look at the examples of function beeswarm in package beeswarm. Hth -- Gerrit - Dr. Gerrit Eichner Mathematical Institute, Room 212 gerrit.eich...@math.uni-

Re: [R] Fwd: Concatenation?

f course, then, e.g., do.call(paste0, data.frame(X)) would work. Best -- Gerrit --------- Dr. Gerrit Eichner Mathematical Institute, Room 212 gerrit.eich...@math.uni-giessen.de Justus-Liebig-University Giessen Tel:

[R] Fwd: Concatenation?

[I forgot to keep it on the list.] Weitergeleitete Nachricht Betreff: Re: [R] Concatenation? Datum: Wed, 17 Feb 2021 22:14:13 +0100 Von: Gerrit Eichner Organisation: JLU Gießen An: Parkhurst, David Hi David, checkout, e. g., base-R's paste0(site, depth) or, if stor

Re: [R] [effects] Wrong xlevels in effects plot for mixed effects model when multiline = TRUE

r support and the brilliant effects package in general! :-) Best regards -- Gerrit ----- Dr. Gerrit Eichner Mathematical Institute, Room 212 gerrit.eich...@math.uni-giessen.de Justus-Liebig-University Giessen Tel: +49-(0)641-99-32104

[R] [effects] Wrong xlevels in effects plot for mixed effects model when multiline = TRUE

.4 nloptr_1.2.2.2 [11] Matrix_1.2-18boot_1.3-25 splines_4.0.2statmod_1.4.34 lme4_1.1-23 [16] tools_4.0.2 survival_3.2-3 yaml_2.2.1 compiler_4.0.2 colorspace_1.4-1 [21] mitools_2.4 insight_0.9.5nnet_7.3-14 -----------

Re: [R] how to generate this kind of graph

Hi, John, look for the package beeswarm. Regards -- Gerrit - Dr. Gerrit Eichner Mathematical Institute, Room 212 gerrit.eich...@math.uni-giessen.de Justus-Liebig-University Giessen Tel: +49-(0)641-99

Re: [R] nearest lower and higher integers to a multiple of another integer

Hi Stefano, maybe floor(x / n) * n and ceiling(x / n) * n does what you want? Best regards -- Gerrit - Dr. Gerrit Eichner Mathematical Institute, Room 212 gerrit.eich...@math.uni-giessen.de Justus-Liebig

Re: [R] how to find number of unique rows for combination of r columns

guess. Regards -- Gerrit - Dr. Gerrit Eichner Mathematical Institute, Room 212 gerrit.eich...@math.uni-giessen.de Justus-Liebig-University Giessen Tel: +49-(0)641-99-32104 Arndtstr. 2, 35392 Giessen

Re: [R] how to find number of unique rows for combination of r columns

Hi, Ana, doesn't udt <- unique(dt[c("chr", "pos", "gene_id")]) nrow(udt) get close to what you want? Hth -- Gerrit --------- Dr. Gerrit Eichner Mathematical Institute, Room

Re: [R] [effects] allEffects does not accept integer value for xlevels

ring to xlevels[[name]] <- levs where xlevels is an integer. (But that doesn't matter anymore because my conjecture pointed to the wrong function anyway.) Best regards -- Gerrit --------- Dr. Gerrit Eichner

[R] [effects] allEffects does not accept integer value for xlevels

rspace_1.4-1 mitools_2.4 nnet_7.3-12 --------- Dr. Gerrit Eichner Mathematical Institute, Room 212 gerrit.eich...@math.uni-giessen.de Justus-Liebig-University Giessen Tel: +49-(0)641-99-32104 Arndtstr. 2, 35392 Giessen

Re: [R] Partition a vector into select groups with fixed length

Hi, Christofer, try something along len <- 5 split(Vec, rep(seq(ceiling(length(Vec)/len)), each = len)) Hth -- Gerrit - Dr. Gerrit Eichner Mathematical Institute, Room 212 gerrit.eich...@math.

Re: [R] Integer Sample with Mean Dependent on Size

s -- Gerrit ----- Dr. Gerrit Eichner Mathematical Institute, Room 212 gerrit.eich...@math.uni-giessen.de Justus-Liebig-University Giessen Tel: +49-(0)641-99-32104 Arndtstr. 2, 35392 Giessen, Germany http://www.uni-giess

Re: [R] Fwd: How to retain the NA as a value for splitting a dataframe

Hello, Jun, try split(df, f = factor(df$C, exclude = NULL)) For more info see ?factor, of course. Regards -- Gerrit - Dr. Gerrit Eichner Mathematical Institute, Room 212 gerrit.eich...@math.uni

Re: [R] Diff'ing 2 strings

Don't you mean ?Rdiff ? Hth -- Gerrit - Dr. Gerrit Eichner Mathematical Institute, Room 212 gerrit.eich...@math.uni-giessen.de Justus-Liebig-University Giessen Tel: +49-(0)641-99-32104 Arn

Re: [R] summary function does not work with Westfall correction in multcomp with 27 comparisons

t, or wait/hope for a more knowledgeable list member to "jump in". Hth -- Gerrit --------- Dr. Gerrit Eichner Mathematical Institute, Room 212 gerrit.eich...@math.uni-giessen.de Justus-Liebig-University G

Re: [R] verInd= and HorInd= arguments to pairs() function

Am 08.06.2018 um 12:02 schrieb Martin Maechler: Martin Maechler on Fri, 8 Jun 2018 11:13:24 +0200 writes: [..] >> Thank you, Chris, for the report and >> Gerrit for your proposed fix !! >> >> It looks good to me, but I will test some more (also with >> 'r

Re: [R] verInd= and HorInd= arguments to pairs() function

Hi, Chris, had the same problem (and first thought it was my fault), but there seems to be a typo in the code of pairs.default. Below is a workaround. Look for two comments (starting with #) in the code to see what I have changed to make it work at least the way I'd expect it in one of your e

Re: [R] Manipulation of data.frame into an array

Why not use as.matrix(Imp) in this case? Regards -- Gerrit Am 24.05.2018 um 17:04 schrieb Bert Gunter: This is one of those instances where a less superficial knowledge of R's technical details comes in really handy. What you need to do is convert the data frame to a single (numeric) vector

Re: [R] Wired result when I convert from Character to Numeric

Hi, Christofer, try print(as.numeric(x), digits = 10) # or another number of digits to increase the printed(!) precision. Hth -- Gerrit Am 23.05.2018 um 12:42 schrieb Christofer Bogaso: Hi, Below is my simple result in R x = "1282553.821000" as.numeric(x) [1] 1282554 Any idea whe

Re: [R] Further questions

. Third, just to be sure: This list is not for home work. Two more remarks see inline below. Hth -- Gerrit - Dr. Gerrit Eichner Mathematical Institute, Room 212 gerrit.eich...@math.uni-giessen.de Justus

Re: [R] Problem with regression line

Hi, Anne, assign Age and Bloodpressure in the correct order to the axes in your call to plot as in: plot(y = Age, x = BloodPressure) abline(SimpleLinearReg1) Hth -- Gerrit - Dr. Gerrit Eichner

Re: [R] effects & lme4: error since original dataframenotfoundWASeffects: error when original data frame is missing

-checking the intended name and if neccessary changing it (somehow ...) automatically. (Have to think about that ...) Best regards -- Gerrit - Dr. Gerrit Eichner Mathematical Institute, Room 212 gerrit.eich

Re: [R] effects & lme4: error since original data frame notfoundWASeffects: error when original data frame is missing

m <- lmer(Reaction ~ Days + (Days | Subject), data = X) Effect("Days", fm) } Hth -- Gerrit --------- Dr. Gerrit Eichner Mathematical Institute, Room 212 gerrit.eich...@math.uni-giessen.de Justus-Lieb

Re: [R] effects & lme4: error since original data frame not found WASeffects: error when original data frame is missing

3.4.2 compiler_3.4.2 [17] colorspace_1.3-2 mgcv_1.8-22nnet_7.3-12 quantreg_5.33 - Dr. Gerrit Eichner Mathematical Institute, Room 212 gerrit.eich...@math.uni-giessen.de Justus-Liebig-University Giessen Tel: +49-(0)641-99-

[R] effects: error when original data frame is missing

n a list like fm.plus.origdata <- list(fm, myair = myair)) to be able to restore the original environemt (or at least the needed opriginal data frame) of the time when fm was created? Thx for any hint! Regards -- Gerrit ------------

Re: [R] barplot_add=TRUE

ot the cause of your "problem". I suggest you check the structure of m: str(m) I suspect you'll see that the assumption of 52 trapcatch values is somehow wrong. Hth -- Gerrit --------- Dr. Gerrit Eichner

Re: [R] barplot_add=TRUE

Hi, Sibylle, since you write '"mathematically" add', does barplot(rbind(m$trapcatch, w$trapcatch)) do what you want (modulo layout details)? Hth -- Gerrit --------- Dr. Gerrit Eichner Ma

Re: [R] Help avoiding setting column type two times

t) save(result, file = "folkeafstemning2009.Rdata") Maybe two loops simplify this a little bit for you (not tested): for(v in c("StedType", )) result[[v]] <- factor(result[[v]]) for(v in c("StemmerAntal", )) result[[v]] <- as.integer(result[[v]]) Hth -- Gerrit ---

Re: [R] adding predictor to linear model without changingexistingcoefficients

Hello, Urs, you may have seen Wolfgang Viechtbauer's answer already which offers an R-technical solution, but this may leave the mathematical grounds of linear models. See inline below for my concern and a hint. Am 17.05.2017 um 09:12 schrieb Urs Kleinholdermann: Dear list members, I want to a

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

(untested). Hth -- Gerrit --------- Dr. Gerrit Eichner Mathematical Institute, Room 212 gerrit.eich...@math.uni-giessen.de Justus-Liebig-University Giessen Tel: +49-(0)641-99-32104 Arndtstr. 2, 35392 Giessen, Germany Fax: +49-(0)641-99-32109htt

Re: [R] plot empty when drawing from custom function

Hi, Luigi, you are probably missing a call to print() around the call to the latticeExtra plotting function useOuterStrips() you use inside your function printer(). Hth -- Gerrit - Dr. Gerrit Eichner

Re: [R] Arguments imply differing number of rows: 1, 0

_ 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, reproducible code. -- Hth -- Gerrit -----

Re: [R] Why is options(digits = 1) displays two digits?

Hi, Tal, in print.default it says: digits: a non-null value for digits specifies the __minimum__ number of __significant__ digits to be printed in values. Maybe this clarifies your observation. Hth -- Gerrit - Dr. Gerrit

Re: [R] How to stop as.integer removing dimenions

... or: net1 <- array(0L, dim=c(5,5)) Note the difference between 0 and 0L. Gerrit - Dr. Gerrit Eichner Mathematical Institute, Room 212 gerrit.eich...@math.uni-giessen.de Justus-Liebig-Univers

Re: [R] How to stop as.integer removing dimenions

Hi, Thomas, maybe mode(net1) <- "integer" does what you want? Hth -- Gerrit ----- Dr. Gerrit Eichner Mathematical Institute, Room 212 gerrit.eich...@math.uni-giessen.de Justus-Liebig-Universit

Re: [R] Power of t test for unequal variances?

Hello, Mauricio, maybe pwr.t2n.test() in package pwr or/and n.ttest() in package samplesize do what you need/want. Hth -- Gerrit - Dr. Gerrit Eichner Mathematical Institute, Room 212 gerrit.eich

[R] update of formula with conditional rhs introduces unwantedparentheses

how to extract various components of a formula and how to "re-build" formulae from such single components, but I would like to avoid that. Thx for any hint -- Gerrit - Dr. Gerrit Eichner Mathematical

Re: [R] make a right subset!

Hello, Elahe, look at ?match and check if df2$Count[ match( Matched, df2$Serial)] does what you want/need. Hth -- Gerrit - Dr. Gerrit Eichner Mathematical Institute, Room 212 gerrit.eich...@math.uni

Re: [R] Where to get quantiles (1st parameter q) for pnorm?(Normalapproximation to binomial distribution)

Hi, Matti, use ?pnorm and read about qnorm. Hth -- Gerrit On Wed, 20 Jan 2016, mviljamaa wrote: I'm doing Normal approximation to binomial distribution. My variables are generated by rbinom. Here: http://msemac.redwoods.edu/~darnold/math15/spring2013/R/Activities/ApproxBinomWithNorm.

Re: [R] Scaling rows of a large Matrix::sparseMatrix()

Hello, Dirk, maybe I'm missing something, but to avoid your for-loop-approach doesn't M <- M/Matrix::rowSums(M) do what you want? Hth -- Gerrit ----- Dr. Gerrit Eichner Mathematical Instit

Re: [R] plot several lines programmaticaly

Hi, Petr, from your example it's not doubtlessly clear to me sure how the sequence of ifs should really be continued, but couldn't a nested loop help? Something like: for (i in 1:npks) { plot(res[[i]]$x, res[[i]]$y, pch=20) for (j in 1:i) { lines(res[[i]]$x, res[[i]]$fitpk[j,], col=j, lw

Re: [R] vector of constant values

Hi, Thanoon! I want to simulate a vector of constants values with dimention = 200x2 and all values of this vector are 1. You mean "I want to _construct_ a _matrix_ of dimension 200 x 2 with all entries equal to 1", don't you? Homework? Take a look at ?matrix Hth -- Gerrit Any help

Re: [R] rbind - names in dataframe. Beginner Question

Hello, Sarah, take a look at the online help page of the function cut() (and apply the function to the Hunger-column). Hth -- Gerrit On Thu, 29 Oct 2015, Dagmar wrote: Dear Sarah, (dear All), Thank you for trying to help! You are right, that I want to add a column to my dataframe with a

Re: [R] extract all dataframes from list

Hello, Maicel, if you only want to extract info from those data frames, maybe it is enough to attach the list to the search path using attach (and afterwards detach()); see ?attach. Otherwise the online help page of assign() might be a starting point for you. Hth -- Gerrit Hi Or basica

Re: [R] Fisher's Test 5x4 table

d*** p-values [of] larger than 2 by 2 tables." (Similarly, control and hybrid are not needed either here.) Regards -- Gerrit --------- Dr. Gerrit Eichner Mathematical Institute, Room 212 gerrit.eich...@math.u

Re: [R] Fisher's Test 5x4 table

initab result for the Chi squared test as proof (I know that getting very low p values are highly unlikely but sometimes it happens). Seeing is believing i suppose! Regards, Paul On Fri, Aug 28, 2015 at 8:56 AM, Gerrit Eichner < gerrit.eich...@math.uni-giessen.de> wr

Re: [R] Fisher's Test 5x4 table

Dear Paul, quoting the email-footer: "PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code." So, what exactly did you try and what was the actual problem/error message? Besides that, have you noted th

Re: [R] replacing then summing by values from another dataframe

Hello, Xianming, I have changed your (particular) data structure: use matrices because you have only numeric scores and effects, use NA instead of -1 as missing value (as usual), don't use columns for ids or row/column names (except for the easy of reading the data structures), increase your s

Re: [R] Recursive looping of a list in R

Hi Evans, not many people (incl. me) are going to guess the building law for your recursive structure from the -- in fact at first sight not so clear -- picture, but I have some comments inline below. Hth -- Gerrit I am trying to creat a list from a loop such that once you loop the value

Re: [R] numeric of length 1

Hi, Ragia, use names( pr_sub$vector[ 1]) Regards -- Gerrit On Tue, 23 Jun 2015, Ragia Ibrahim wrote: Dear group, I have the following object pr_sub$vector[1] 14 0.07782809 class( pr_sub$vector[1]) [1] "numeric" > length( pr_sub$vector[1]) [1] 1 how can I separate pr_sub$

Re: [R] [FORGED] Re: Dunnett Test in 'multicomp' package

Hello, everyone, aside from Rolf's hint (and Richard's warning!) you could also consider relevel(): viagraData$dose <- relevel( viagraData$dose, ref = "placebo") Hth -- Gerrit On Fri, 5 Jun 2015, Rolf Turner wrote: On 05/06/15 11:08, Jim Lemon wrote: Hi James, You can change the order

Re: [R] generate a list as follows: 1, 1, 1, 1, 2, 2, 2, 2, 3, 3, 3, 3, . . . . ., n, n, n, n

Hi, John, doesn't n <- lapply( 1:n, rep, each = 4) do what you need? Hth -- Gerrit On Sat, 18 Apr 2015, John Sorkin wrote: Windows 7 64-bit R 3.1.3 RStudio 0.98.1103 I am trying to generate a list of length 4n which consists of the integers 1 to n repeated in groups of four, i.e.

Re: [R] check a list that a sublist exists

Hello, Alaios, try is.null( Mylist[[i]]$parameters) It returns TRUE if your list doesn't have a component named "parameters" or if that component contains NULL. HtH -- Gerrit On Tue, 3 Mar 2015, Alaios via R-help wrote: Hi all,I have a list that has the following fields. $`80` [1] "Er

Re: [R] from expand.grid to a list (lapply)

Hie Alaiso, try to simply apply lapply() to your data frame. It should work w/o any further operation since data frames are lists (check with is.list()). Hth -- Gerrit On Wed, 25 Feb 2015, Alaios via R-help wrote: Hi all,I would like to use an expand.grid functionality that would give me

Re: [R] i need help for var.test()

Hello, Sait, take, for example, a look at bartlett.test(), but keep in mind that Bartlett's test is quite sensitive to deviations from normality. Levene's test (e.g., in leveneTest() in package car) is said to be more robust. Hth -- Gerrit On Thu, 8 Jan 2015, sait k wrote: Dear Sir or M

Re: [R] factor levels > numeric values

Hello, David, take a look at the beginning of the "Warning" section of ?factor. Hth -- Gerrit Hi everybody, I have another question (to which I could not find an answer in my r-books. I am sure, it's not a great issue, but I simply lack of a good idea how to solve this: One of my variable

Re: [R] qq-plot in R

Hello, ron, have you looked at the help page of qqplot and consequently tried, e.g., the following? xy <- qqplot( rt(300, df = 5), rt(300, df = 5), xlim = c(-10, 10), ylim = c(0, 10)) str( xy) Hth -- Gerrit On Tue, 11 Nov 2014, Ron Michael wrote: Hi, I am some questions

[R] [survMisc]: error message in examples of comp()

ind regards -- Gerrit --------- Dr. Gerrit Eichner Mathematical Institute, Room 212 gerrit.eich...@math.uni-giessen.de Justus-Liebig-University Giessen Tel: +49-(0)641-99-32104 Arndtstr. 2, 35392 Giessen, Germany Fax: +49-(0)6

Re: [R] Compare data in two rows and replace objects in data frame

"-1" = "Bht")[ paste( x, collapse = "")]) presumably does what you want (and can certainly be improved). Hth -- Gerrit - Dr. Gerrit Eichner Mathematical Institute, Room 212 gerrit.eich...@math.uni-giessen.de

Re: [R] Filling in missing values in a column based on previousandfollowingvalues

column, not all: Something like yourdataframe$ID <- na.locf( yourdataframe$ID) should replace the ID-column with the modified version you want. Regards -- Gerrit Kind regards, Florian Am 04.08.2014 um 11:31 schrieb Gerrit Eichner : Hello, Florian, function na.locf() from package

Re: [R] Filling in missing values in a column based on previousandfollowing values

Hello, Florian, function na.locf() from package zoo mightdo what you want. Hth -- Gerrit - Dr. Gerrit Eichner Mathematical Institute, Room 212 gerrit.eich...@math.uni-giessen.de Justus-Liebig-University

[R] format of output of residual standard errors of manova()

x64 (64-bit) locale: [1] LC_COLLATE=German_Germany.1252 LC_CTYPE=German_Germany.1252 [3] LC_MONETARY=German_Germany.1252 LC_NUMERIC=C [5] LC_TIME=German_Germany.1252 attached base packages: [1] stats graphics grDevices utils datasets methods base other attached packages: [1] f

Re: [R] Sweave: cat() in a chunk with option results=texdoesn'tproducelinebreaksatthe end of a character string anymore

lse "locally" in case of need, i.e., as the option for the respective code chunks (which also have results=tex, and from which I want "only" linebreaks in the TeX-output). Actually, embarrassingly simple. Thanks again -- Gerrit On Wed, 5 Mar 2014, Duncan Murdoch wrote:

Re: [R] Sweave: cat() in a chunk with option results=texdoesn'tproducelinebreaks at the end of a character string anymore

arks ... Any further ideas? Regards -- Gerrit On Wed, 5 Mar 2014, Duncan Murdoch wrote: On 05/03/2014 7:32 AM, Gerrit Eichner wrote: Hello, everyone, I am struggling with an Sweave-problem that didn't occur sofar (and I have no clue what I might have changed in my system; see below). Th

[R] Sweave: cat() in a chunk with option results=tex doesn't producelinebreaks at the end of a character string anymore

.1252 attached base packages: [1] graphics grDevices utils datasets stats methods base other attached packages: [1] fortunes_1.5-2 loaded via a namespace (and not attached): [1] grid_3.0.2 lattice_0.20-27 tools_3.0.2 -

Re: [R] calculate an value in dependence of another column

Hi, Mat, e. g., try something like ORDER$VALUE * c( NY = 3, BER = 2, FRA = 1, WAS = 4)[ ORDER$DESTINY] Hth -- Gerrit On Fri, 17 Jan 2014, Mat wrote: Hello together, i have a little problem, to create a new column, in a data.frame. I know i can calculate one column with as a example 2 lik

Re: [R] How to verify char variables contain at least one value

e: Hi Gerrit, Thank you for the suggestion. Unfortunatedly the line of code you suggest - i.e. apply( subset( d2, subset = V13:V239), 1, function( x) any( x != "")) returns the same error. Any other suggestion? Danke schôn und auf wiederhoren Luca 2014/1/2 Gerrit Eichner

Re: [R] How to verify char variables contain at least one value

Hello, Luca, also a happy new year! It's not quite clear to me what you want to do, but note first that the ":"-operator is a short-cut for seq() with by = 1 (look at ?seq), and that it usually (!) does not work on columns of data frames. Exception: when used for the argument subset of functi

Re: [R] filter a data.frame

Hello, Mat, see below. hello together, i want to filter a data.frame. My problem is, that i want to filter 2 numbers. My data.frame look like this one. No. text 1 abc 2 def 3 ee 4 ff 5 gg I want now to filter No. 2 and 3, so my solution should be look like this on

Re: [R] how can i write the function into a file c:/mytest.R withcatfunction?

Hello, Ë®¾²Á÷Éî! mytest<-function(x,f){ sum(x*f)/sum(f) } cat(mytest,file="c:/mytest.R") Error in cat(list(...), file, sep, fill, labels, append) : argument 1 (type 'closure') cannot be handled by 'cat' how can i write the mytest function into a file c:/mytest.R with cat

Re: [R] Splitting a vector

Hello, Long Vo, take a look at the help page of split or directly at str(Y) They tell you that Y is a list, and list components are indexed using "[[": mean(Y[[4]]) should do what you want. Regards -- Gerrit This does what I needed. However, as the output is a list object, is there

Re: [R] If-statement in for-loop

Hello, "gncl dzgn", your problem has a flavor of homework which is usually not delt with on this list. However, a few comments: 0. The description of your problem is rather vague, in particular, the meaning of "input" in the description of your "conditions" is unclear! (By the way, your main

Re: [R] lattice: superposed boxplots with same colors forrectanglesandumbrellas and filled boxes

mycolors, alpha = myalpha)) ) Rich On Mon, Dec 9, 2013 at 5:17 PM, Gerrit Eichner wrote: Dear R-list, I've been trying to produce a sort of an interaction plot wherein colored stripplots and boxplots are superposed. In particul

[R] lattice: superposed boxplots with same colors for rectanglesand umbrellas and filled boxes

ox.[dot, rectangle, umbrella], which does not give what I need. I may have overlooked the obvious and would be quite grateful if somebody could give me a hint where to look further. Or is a workaround necessary (and available)? Thanks a lot in advance! Best Regards -- Gerrit -

Re: [R] interaction plot with SE bar

Hi, Kristi, I have adapted the interaction plot example of function panel.average() of the lattice package and modified the code of panel.average() to a new function panel.loc_and_scale(). See below. (There might be a much simpler solution, though.) Remark on the side: I can't resist to poin

Re: [R] remap values from one vector to the other

Alex, it is a bit unclear what you mean by "remap" etc., but maybe y0 <- 880e6; y1 <- 1020e6 x0 <- 1; x1 <- ncol(x) y0 + (x0-1):x1 * (y1 - y0)/(x1 - (x0-1)) gives what you want. Hth -- Gerrit On Thu, 31 Oct 2013, Alaios wrote: Hi everyone, I am plotting some legend and I am usi

Re: [R] omitting integer(0) rows from data frame

Hi, Jack, well, I disagree: What do you expect to grab out of a bucket (= data frame) if you do not at all grab into it (indexing with an _empty_ index, i.e. with nothing)? And changing the sign of nothing is still nothing ... Hth -- Gerrit On Wed, 30 Oct 2013, Jack Tanner wrote: I'm not

Re: [R] Select fixed number of elements

Hello, Alaois, if x is your vector maybe n <- length( x) positions <- trunc( quantile( seq( n), prob = 0:5/5)) x[ positions] comes close to what you want. Hth -- Gerrit Hi all, I have in my code some vectors that are not of equal size. I would like to be able for each of these vectors sel

Re: [R] Nonparametric k-way ANOVA

Hello, Vincent, you may want to take a look at "Nonparametric methods in factorial designs" by Edgar Brunner and Madan L. Puri in Statistical Papers 42, 1-52 (2001). There is the R-package nparcomp for one-way layouts, but the paper goes further (and mentions another software) and is maybe a

Re: [R] Where is element 30?

Hi, Alaios, check out ?which and in particular its "Examples" and "See Also" section. Hth -- Gerrit On Tue, 22 Oct 2013, Alaios wrote: Hi I have a vector like that readCsvFile$V1  [1]  30  31  32  33  34  35  36  37  38  39 310 311 312 313 314 315 316 317 318 [20] 319 320 321  20  21  22 

Re: [R] bug(?) in str() with strict.width = "cut" when appliedtodataframe with numeric component AND factor or character component withlongerlevels/strings

252 LC_NUMERIC=C [5] LC_TIME=German_Germany.1252 attached base packages: [1] stats graphics grDevices utils datasets methods base other attached packages: [1] fortunes_1.5-0 loaded via a namespace (and not attached): [1] tools_3.0.2 On Wed, 16 Oct 2013, Gerrit Eichner wrote: Thanks, D

Re: [R] bug(?) in str() with strict.width = "cut" when appliedtodataframe with numeric component AND factor or character component withlongerlevels/strings

just forgot to take that possibility into consideration. :( Regards -- Gerrit On Tue, 15 Oct 2013, Duncan Murdoch wrote: On 15/10/2013 7:53 AM, Gerrit Eichner wrote: Dear list subscribers, here is a small artificial example to demonstrate the problem that I encountered when looking at the

Re: [R] extract regression coefficients

ic to give you any hint. So, "PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code." Regards -- Gerrit On 15 October 2013 14:59, Gerrit Eichner wrote: Hello, Catalin, check out

Re: [R] extract regression coefficients

Hello, Catalin, check out ?coef Regards -- Gerrit On Tue, 15 Oct 2013, catalin roibu wrote: Hello all! I have a problem with R. I want to extract regression coefficients from summary and use it for compute the theoretical values. How can I do that in R? thank you! best regards, -- ---

[R] bug(?) in str() with strict.width = "cut" when applied to dataframe with numeric component AND factor or character component with longerlevels/strings

4.4 grid_3.0.2 lattice_0.20-23 MASS_7.3-29 [5] nnet_7.3-7 rpart_4.1-3 stats4_3.0.2tools_3.0.2 --------- Dr. Gerrit Eichner Mathematical Institute, Room 212 gerrit.eich...@math.uni-giessen.de Justus-Liebig-U

Re: [R] How to write a wrapper function which can honor defaultvalues, when the target function has ones?

Hello, Adam, I'm rather uncertain about your goal (and consequently even more so about how to reach it), but anyway, maybe the function match.call() with its argument expand.dots is of some help for you. From its help page: "match.call is most commonly used in two circumstances: To recor

Re: [R] substituting string variable in expression with its value

Hi, Jannis, maybe plot( 1, 1, main = bquote( paste( .(a), " [", .(b), "]"))) comes close to what you want, but I think you may even have to use the following to get a varying exponent really printed elevated: a <- "speed" b <- "m * s" cc <- -2 plot( 1, 1, main = bquote( paste( .(a), " [",

  1   2   >