Re: [R] MLE Estimation of Gamma Distribution Parameters for data with 'zeros'

There's at least one package that can do zero-inflated gamma regression (Rfast2::zigamma). I'm not sure it's ML, though. On Thu, Jan 19, 2023 at 10:17 AM Jeff Newmiller wrote: > Beware of adding a constant... the magnitude of the constant used can have > an outsized impact on the answer obtain

[R] startup loading issue

happens upon starting R. I’ve also created a new .Rprofile that successfully executes but does not prevent the load problem unless the .Rprofile includes a quit() command. Any suggestions? Ken [[alternative HTML version deleted]] __ R-help@r

Re: [R] generate random numeric

That's all right. Thanks. On Sat, Oct 30, 2021 at 12:29 AM Marc Schwartz wrote: > Ken Peng wrote on 10/29/21 2:39 AM: > > I saw runif(1) can generate a random num, is this the true random? > > > >> runif(1) > > [1] 0.8945383 > > > > What's

[R] generate random numeric

I saw runif(1) can generate a random num, is this the true random? > runif(1) [1] 0.8945383 What's the other better method? Thank you. [[alternative HTML version deleted]] __ R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see h

Re: [R] Error with RMySQL

It seems like mysql server doesn't open the right authentication for your access, such as user or host permissions. Thanks. On Thu, Oct 28, 2021 at 1:42 PM Ashim Kapoor wrote: > Dear R - users, > > I have 2 databases on a MySQL server. I am able to access the old one > but not the freshly creat

Re: [R] Help needed in double bar plot made using ggplot2

Thanks Dr. Burradas too. i also had the same question. regards August 20, 2021 6:02 AM, "bharat rawlley via R-help" wrote: > Thank you, Dr. Burradas! > That resolved my query > Have a great rest of your day > On Thursday, 19 August, 2021, 04:47:42 pm GMT-4, Rui Barradas > wrote: > __

Re: [R] R matrix multiplication slowdown

I've done more searching and found a problem with the data in one of the matrices that was corrupting the calculation. Data now fixed and problem is solved. Apologies if anyone wasted time on this. Ken On Fri, Jan 27, 2017 at 8:29 AM, Jeff Newmiller wrote: > You are asked by the

[R] R matrix multiplication slowdown

o B, the same command takes 2.4 seconds; at 1050 rows in B, the command is up to almost 4 seconds. Just trying to understand why the big slowdown is occurring, especially since the matrices I actually want to multiply have 1000 and 5000 rows, not 100 and 1017. Thanks, Ken (Macbook Pro runn

[R] problem with Rcpp and boost threadpool

I'm having a problem calling a local library through Rcpp with R Studio Server. It's a bit perplexing, since I have no issues when I call it from R at the command line. I've written an analytics library which uses boost's threadpool functionality for running multiple threads. I've stripped eve

Re: [R] Opposite color in R

the results are unlikely to match a given individual's vision. On top of that, decisions made when this norm was specified are such that it deviates from human vision for short wavelengths so that you would be better off using a corrected version like that proposed by Judd in the 1950's or

[R] help with ff matrix indexing and value assignment

g cbind = still getting only zeros. Tried a few other methods I just cant seem to wrap my head around how to do this. I had a look at ffindexset, but that looks like assigning whole rows at a time. I appreciate any help !! Thanks Ken __ R-help@r-p

[R] Are melt and dcast with ffdf possible?

5 515 6 616 c.melt.ff <- melt(comb.ff, id.vars=c("Var1","Var2")) > c.melt.ff value NA NA 1 1 1 1 2 2 1 2 3 3 1 3 4 4 1 4 5 5 1 5 6 6 1 6 7 7 1 7 .

[R] data.table cant find function melt?

Hi all, trying to use the melt function in data.table and I'm getting an error Anyone seen this before or know how to fix it? Thanks str(Distdata) Classes ‘data.table’ and 'data.frame': 828451 obs. of 3 variables: $ Poly1 : int 50088 50088 50088 50088 50088 50088 50088 50088 50088 500

Re: [R] Help with GLM starting values in user defined link function

csiro.au> writes: > I'm trying to fit a binomial GLM with user defined link function (negative exponential), however I seem to > be unable to find the correct starting values to initialise such a model. I've tried taking starting > values from a logistic and log models fit to the same data an

Re: [R] outer() problem

Marc Girondot yahoo.fr> writes: > > outer(0:1, 0:1, FUN=function(x, y) {1}) > Erreur dans outer(0:1, 0:1, FUN = function(x, y) { : >dims [produit 4] ne correspond pas à la longueur de l'objet [1] Because whatever the dimensions of your 2 input vectors, this function simply returns the value

Re: [R] Generate Binary Matrix

t; > table(rowSums(dimMat)) Wht don't you sample from the distribution of row sums for each row and then distribute that many 1's randomly among the columns. Ken __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLE

Re: [R] How do I perform conditional annotation of lattice panel plots?

Chen, George roswellpark.org> writes: > My apologies for asking this question that may have been asked before. I am trying to plot activity > dependent on time conditioned by the subject. Code for sample data below. > So I have something like this > xyplot(Activity~Time|Subject). > This works fi

Re: [R] Adding the complementary log-link to binomial() and make.link()

Roland Deutsch tuwien.ac.at> writes: > in my research I frequently work with binomial response models, which > are of course part of the generalized linear models. While I do use > common link functions such as the logit, probit and cloglog, I often > have the need of invoking the lesser-kno

Re: [R] Select fixed number of elements

Alaios yahoo.com> writes: > I have in my code some vectors that are not of equal size. I would like to be able for each of these vectors > select 6 elements that are (almost) equally spaced. So the first one would be at (or close) to the beginning > the last one at (or close) to the end and the

Re: [R] colour code areas of a plot

Martin Batholdy googlemail.com> writes: > I would like to colour different areas of a plot. > But I don't know how to do this efficiently. > > here an example: > (t = time) > > t <- 1:100 > bg_colors <- c(rep('green',20), rep('yellow',10), rep('green',20), rep('red',5), > rep('yellow',45)) >

Re: [R] speeding up

Ken Knoblauch inserm.fr> writes: > > Bos, Roger rothschild.com> writes: > > I am using a sum of squared differences in the > objective function of an optimization problem I am > doing and I > > have managed to speed it up using the > outer function ve

Re: [R] speeding up

Bos, Roger rothschild.com> writes: > I am using a sum of squared differences in the objective function of an optimization problem I am doing and I > have managed to speed it up using the outer function versus the nested for loops, but my suspicion is that > the calculation could be done even q

Re: [R] Looking for package to solve for exponent using newton's method

Thanks! That's just what I needed. __ 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, reproducibl

[R] Looking for package to solve for exponent using newton's method

Hi, I'm looking for an R function/package that will let me solve problems of the type: 13 = 2^x + 3^x. The answer to this example is x = 2, but I'm looking for solutions when x isn't so easily determined. Looking around, it seems that there is no algebraic solution for x, unless I'm mistaken. Do

Re: [R] import function without overwriting function with the same name

Martin Batholdy googlemail.com> writes: > I have to import multiple R-files. > Each file consists of several functions with the same function name across the R-files. > When I import all files one by one (with source()) I overwrite the function definition of the previous file > until only the ve

Re: [R] grDevices::convertColor XYZ space is it really xyY?

If they sum to 1 then they are one and the same. Look at how chromaticity coordinates are defined in terms of the tristimulus values. Quoting Bryan Hanson : Thank you Ken. 90% of my problem was missing the factor of 100. I was just inputing xyY as a test, I wasn't sure whether the

Re: [R] grDevices::convertColor XYZ space is it really xyY?

uot;XYZ", to = "sRGB") [,1] [,2] [,3] [1,]1 11 Quoting Bryan Hanson : Ken, I followed your suggestion and perhaps I don't understand what to expect from convertColor or maybe I'm not using it correctly. Consider the following tests: D65 <- c(0.3127,

Re: [R] grDevices::convertColor XYZ space is it really xyY?

Z are tristimulus values and xyY are chromaticity coordinated and the luminance which is the Y tristimulus value for the CIE 1931 standard observer. Ken __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the

Re: [R] measuring distances between colours?

I'd have to look it up and I'm not home at the moment. Can see later on. I would have thought that it would be normalized to have a jnd equal to 1 but I'm not sure. Ken Sent from my iPhone ___ Ken Knoblauch Inserm U846 Stem-Cell and Brain Research Institute 18 av du Doyen Lé

Re: [R] measuring distances between colours?

Hi John, Out of curiosity and if it is not much trouble, I would be curious if Luv worked any better than Lab. I think that Luv is supposed to be preferred for monitors and Lab for surfaces but they are generally pretty similar. Best, Ken Sent from my iPhone ___ Ken Knoblauch Inserm U846

Re: [R] measuring distances between colours?

so should be better at finding the nearest color, I would think. Untried though. Don't forget that these are based on a standard color calibration, sRGB, I think in the case of R and your display may or may not be close to that, so expect some error from that as well. best, Ken >

[R] Using Rscript in version 3.0.0

Greetings, I have just installed version 3.0.0. I am trying to use code that I have used numerous times in previous versions of R. My code executes correctly until I try to call makePSOCKcluster. I issue the following command and get the following error: > cluster <- makePSOCKcluster(nodes, p

Re: [R] plotting CIE chromaticity diagram?

an interpolation algorithm, to reproduce what is on the wikipedia page. Of course, I'm kicking myself for encouraging this be done given the opportunity it provides for misguidance and misrepresentation. Ken Quoting Bryan Hanson : I am, unfortunately, well-aware of the limitations tha

Re: [R] plotting CIE chromaticity diagram?

ishi soichi gmail.com> writes: > Has anyone plotted or is it possible to plot > > CIE *xy* chromaticity diagram > > http://en.wikipedia.org/wiki/File:CIE1931xy_blank.svg > > I need this plot in color. > > ishida And following up on my previous mail (diatribe), after having contemplated th

Re: [R] plotting CIE chromaticity diagram?

ishi soichi gmail.com> writes: > Has anyone plotted or is it possible to plot > > CIE *xy* chromaticity diagram > > http://en.wikipedia.org/wiki/File:CIE1931xy_blank.svg > > I need this plot in color. > > ishida > I think that plotting the spectral locus and the line of purples is trivial,

[R] R intermittently crashes across cluster

as to what might be going on? I have run this code successfully many times when I do not use the loop. I have a lot of data to process and recreating the cluster every time that I want to run my function is a waste of time. Thanx, Ken "For I know the plans I have for you," declares

Re: [R] how to suppress the intercept in an lm()-like formula method?

; method fails when called > > normally, but works if I explicitly use -1 in the formula. I could hack > > the result of model.matrix(), > > but maybe there's an easier way? Have a look at the polr function in MASS where this same problem is handled, I think, around lines 1

[R] Wrong bin count number with hist() ?

ght gray",main="wrong - hist()") truehist(c(1:15,B,50),breaks=c(0,15,A,50),main="correct - truehist()") # I was not able to find an explanation online or in the R documentation.  # The option "include.lowest" doesn't help in such cases. # Any

Re: [R] [lattice] how to label panels with variable value (not name)?

list. I think by selecting the one with which.panel, it may have been looking for the second one in the vector, but each time the vector was of length 1. Just idle speculation though... Ken Quoting Bert Gunter : Ken: I would have thought so, too. However: x <- 1:10; y <- runif(1:10)

Re: [R] [lattice] how to label panels with variable value (not name)?

lattice with the _value_ of the level (an > atmospheric pressure), rather than the name or index of the level. > How to do that? > > TIA, Tom Roche pobox.com> maybe, see ?strip.custom in lattice -- Ken Knoblauch Inserm U846 Stem-cell and Brain Research Institute Department of Integ

Re: [R] history and readline, Mac OSX

vice? Try looking at http://R.research.att.com/libs/ > > TIA -- Christian > -- Ken Knoblauch Inserm U846 Stem-cell and Brain Research Institute Department of Integrative Neurosciences 18 avenue du Doyen Lépine 69500 Bron France tel: +33 (0)4 72 91 34 77 fax: +33 (0)4 72 91 34 61 portable: +33

Re: [R] Calculating the mean in one column with empty cells

blank cells with NA before importing. If you try to import an data frame with empty cells, you usually get an error using read.table(). But since you seem to have already got you data into R, that may not be the problem. HTH, Ken __ R-help@r-project.o

Re: [R] blank plot----how do I make symbols appear

Jessica da Silva gmail.com> writes: > I am trying to create a scatterplot, coding each point to one of 5 > populations. I was successful when I did this for one set of data, yet > when I try plotting other data a blank plot appears (although the axes are > labelled and I can fit the regression

[R] Memory errors using lmer

help(memory.size) 3: In model.matrix.default(mt, mf, contrasts) : Reached total allocation of 2187Mb: see help(memory.size) 4: In model.matrix.default(mt, mf, contrasts) : Reached total allocation of 2187Mb: see help(memory.size) Thanks in advance, Ken Ken McCall | Database reporter Dayton D

Re: [R] self-starter functions for y = a + b * c^x

this. It was designed to fit gamma functions to the luminance vs frame buffer values measured on CRT screens. But the functional form is similar. > thx > Christof > > best, Ken -- Kenneth Knoblauch Inserm U846 Stem-cell and Brain Research Institute Department of Integrative Neuro

[R] [R-pkgs] new package MPDiR version 0.1-11

es several data sets from published psychophysical experiments using detection and rating scale measures for estimating signal detection parameters, psychometric functions, classification images, etc. Best, Ken -- Ken Knoblauch Inserm U846 Stem-cell and Brain Research Institute Department of I

Re: [R] interpolation to montly data

uot;, "29/03/09","26/04/09","24/05/09", "28/06/09", "26/10/08","23/11/08","21/12/08","18/01/09", "15/02/09","16/03/09",&qu

Re: [R] Multiple cbind according to filename

le name just makes sure that the strings in filenames is in the same order as the files are in list.files: a1 <- data.frame("col1" = seq(1,10, 1)) a2 <- data.frame("col2" = seq(11,20, 1)) b1 <- data.frame("col3" = seq(21,30, 1)) b2 <- data.f

Re: [R] Multiple cbind according to filename

le name just makes sure that the strings in filenames is in the same order as the files are in list.files: a1 <- data.frame("col1" = seq(1,10, 1)) a2 <- data.frame("col2" = seq(11,20, 1)) b1 <- data.frame("col3" = seq(21,30, 1)) b2 <- data.f

Re: [R] Replacing cretin value in a file

Perhaps you are in the first circle of the R Inferno? Ken Hutchison On May 13, 2012, at 9:31 AM, Berend Hasselman wrote: > > On 13-05-2012, at 15:08, Jonsson wrote: > >> Dear All, >> >> I am trying to replace a value of 528.8933 to - in my file >&

Re: [R] Seeking help with LOGIT model

638766 2.141 0.0323 * `X 5`2.134544 2.570756 0.830 0.4064 Good luck. Ken Quoting Christofer Bogaso : Thanks Ken for your reply. No doubt your english is quite tough!! I understand something is not normal with the 5th explanatory variable (se:2872.17069!) However could not unde

Re: [R] Seeking help with LOGIT model

kes with my above implementation? I s it just because, I have too > less number of '0' in my response Variable? > Look at the output of summary, especially the standard errors. You seem to be getting complete separation on X5 and X4 doesn,'t look so hot either. Ken ___

Re: [R] Apply function to every 'nth' element of a vector

ken knoblauch inserm.fr> writes: > > Michael Bach gmail.com> writes: > > how do I e.g. square each second element of a > vector with an even > > number of elements? Or more generally to > apply a function to every > > 'nth' element of a vec

Re: [R] Apply function to every 'nth' element of a vector

example: > v <- c(1, 2, 3, 4) > mysquare <- function (x) { return (x*x) } > w <- applyfun(v, mysquare, 2) > then w should be c(1, 4, 3, 16) > Michael Bach Hi Michael, v^(2 - seq_along(v) %% 2) [1] 1 4 3 16 Ken -- Ken Knoblauch Inserm U846 Stem-cell and Brain Research

Re: [R] User defined link function with extra parameters

t; Is there any way to tell glm() to add this > parameter in the estimation or do I have to write my own estimator with > optim()? If the parameter cannot be made into a coefficient of the linear predictor, then I'm afraid that you will have to roll your own. > Thanks, > &g

Re: [R] How to get intersection of multiple vectors?

How can I do then? > > What I know is only for 2 vectors via "intersect" function, but don't know how to deal with multiple vectors. > Reduce(intersect, list(v1 = c("a","b","c","d"), v2 = c("a","b&quo

Re: [R] Placing a Shaded Box on a Plot

Stephanie Cooke gmail.com> writes: > > Hello, > > I would like to place shaded boxes on different areas of a > phylogenetic tree plot. Since I can not determine how to find axes on > the phylogenetic tree plot I am not able to place the box over certain > areas. Below is example code for the sh

Re: [R] Create subset using variable

# if null values are given as NA df.data.sub = df.data.sub[df.data.sub >= 0] # if null values are < 0 #subset dataset is not created # analyses on subset take place j = j + 1 } Hope that helps, Ken __ R-help@r-project.or

Re: [R] howto test a package without installation

Did not fully read the without installing it part. Mea Culpa, Ken On Jan 15, 2012, at 8:56 PM, Jeff Newmiller wrote: > I don't believe you can. However, you need not install it into a system-wide > library directory... your personal library (e.g. > /home/jonas/R/x86_6

[R] Fwd: Trouble installing packages on R2.14.1

Begin forwarded message: > From: Ken Hutchison > Date: January 15, 2012 8:54:49 PM EST > To: Ben Bolker > Subject: Re: [R] Trouble installing packages on R2.14.1 > > Check browser proxy settings and run R.exe with the proper flags to use them > from cmd. > Ho

Re: [R] howto test a package without installation

install.packages('pathtotargz', repos=FALSE) I believe that will get you going. Hope that's helpful, Ken On Jan 15, 2012, at 8:33 PM, Joshua Wiley wrote: > Hi Jonas, > > Look at Hadley Wickham's devtools package. It is designed with this > sort of thing.

Re: [R] Problem Installing R to SuSE 10 via RPM

convert the deb package to rpm locally to install after downloading the deb, although building from source shouldn't be that challenging unless you have ridiculous Fortran compiler problems like I do with Linux. Hope that helps, Ken On Jan 13, 2012, at 1:13 PM, Marc Schwartz

Re: [R] clear plot linear mixed model

Christof Kluß email.uni-kiel.de> writes: > Am 02-01-2012 10:54, schrieb ken knoblauch: > > Christof Kluß email.uni-kiel.de> writes: > >> lme<- lme(conc ~ name/time - 1, > >> random=conc~time|nr,method="ML",data=measurements) > > see pl

Re: [R] quadratic programming-maximization instead of minization

I don't have experience with this in R and I'm not sure I understand the question that well but maybe something like nearPD()? Ken Hutchison On Jan 2, 2012, at 6:36 AM, riccardo24 wrote: > Hi, I need to maximize a quadratic function under constraints in R. > For minimization

Re: [R] clear plot linear mixed model

me) and > seperate colors for the measurements (nr). > > How would you do that? > > thx > Christof > see plot.augPred in the nlme package -- Ken Knoblauch Inserm U846 Stem-cell and Brain Research Institute Department of Integrative Neurosciences 18 avenue du Doyen Lépine

Re: [R] R on Android

you, I'd advise looking for food; else be prepared to name your phone Wilson. Ken Hutchison On Jan 1, 2012, at 11:28 AM, David Winsemius wrote: > > On Jan 1, 2012, at 11:07 AM, David Winsemius wrote: > >> >> On Jan 1, 2012, at 10:05 AM, Federico J. Villatoro w

Re: [R] boot.ci: [Error: cannot allocate vector of size 1.5 Gb]

the maxreps that are right for your machine. Ken Hutchison 2011/12/20 Uwe Ligges > > > On 20.12.2011 06:47, Vikram Bahure wrote: > >> Dear R users, >> >> I am getting following error while using boot.ci. I have int.inc function >> with 2000

Re: [R] Help with the Mice Function

Try reducing the maximum iterations. Probably won't make your call instantaneous, but might be worth the slack gained over a day or so. Ken Hutchison On Dec 9, 2554 BE, at 1:59 PM, "Richard J. Buning" wrote: > Hi, > > I am attempting to impute my data for missi

Re: [R] Principal componet plot from lower triangular matrix file

created, is going to be very hectic. HTH Ken Hutchison On Nov 28, 2554 BE, at 5:55 AM, cm wrote: > Hi, > > I have a comma separated file with element names in first column like shown > below : > > Name_1,0 > Name_2,0.8878,0 > Name_3,0.6777,0.7643,0 > Na

Re: [R] window manager interface commands for linux

our R terminal settings so that organization and dir() will get you what you want. Hope that helps, Ken Hutchison On Nov 28, 2554 BE, at 6:07 AM, Ana wrote: > How can i replicate this in Linux: > source(file.choose()) > > > I've tried source(tkgetOpe

Re: [R] Question about randomForest

uld get more votes on the error. Please correct me if I am wrong, Hopefully a specialist will come along and clear this up, Ken Hutchison On Nov 27, 2554 BE, at 3:21 AM, Matthew Francis wrote: > Thanks for the help. Let me explain in more detail how I think that > ran

Re: [R] append to PDF file

PDF files contain information at the end of them and so you cannot append without altering the file (universally true for PDF). Perhaps pdf() your plots and use external tools to convert the PDFs to .ps then re-merge. Might not be the best way, but an effective one. Ken Hutchison

[R] NA handling in tree package

r and possibly more robust ways to handle NAs within tree objects? Note: I only would like information on tree( ) objects,random forest imputation et cetera are not wanted in this application. Thank you in advance, Ken [[alternative HTML version de

Re: [R] installing java on ubuntu 11.10 installation

Which is better is a matter of opinion: Try: sudo apt-get install sun-java6-jre sun-java6-plugin sun-java6-fonts if it doesnt work, update/add repos. Also, miss google knows this; lastly: not really a geRmaine topic for the R-list. HTH Ken Hutchison On Wed, Nov 9, 2011 at 4:40 PM, ravi

Re: [R] palettes for the color-blind

ibm.com/dx/proceedings/ pravda/truevis.htm who was (is) quite concerned with this issue, as well, as the excellent article by Zeileis, Hornik and Murrell http://statmath.wu.ac.at/~zeileis/papers/ Zeileis+Hornik+Murrell-2009.pdf HTH, Ken -- Ken Knoblauch Inserm U846 Stem-cell and Brain Research

Re: [R] Imputing Missing Data: A Good Starting Point?

Hope this helps: http://rss.acs.unt.edu/Rdoc/library/randomForest/html/rfImpute.html Ken Hutchison On Nov 1, 2554 BE, at 5:29 PM, Sascha Vieweg wrote: > Hello > > I am working on my first attempt to impute missing data of a data set with > systematically incomplete answers (school

Re: [R] glm-poisson fitting 400.000 records

Your memory shouldn't be capped there, try ?memory.size and ?memory.limit. Background less things. Good luck, Ken Hutchison On Oct 21, 2554 BE, at 11:57 AM, D_Tomas wrote: > My apologies for my vague comment. > > My data comprises 400.000 x 21 (17 explanatory variables,

Re: [R] Arima Models - Error and jump error

Perhaps: require(forecast) ?auto.arima # Or look into package fitAR. The first performs seasonal optimization so it is likely better for your application. Ken Hutchison On Oct 21, 2554 BE, at 1:59 PM, Flávio Fagundes wrote: > Hi people, > > I´m trying to development a simple r

Re: [R] Running R with browser without installing anything

Try Googling "R Portable" Ken Hutchison On Oct 20, 2554 BE, at 2:13 PM, jim holtman wrote: > It runs fine off a flash drive. > > On Thu, Oct 20, 2011 at 1:29 PM, Bogaso Christofer > wrote: >> Dear all, the company I work for has Matlab installed for >> sta

Re: [R] ar() - AIC and BIC

gaRds, Ken Hutchison On Wed, Oct 19, 2011 at 7:22 PM, J Toll wrote: > Hi, > > I'm slowly working through Tsay's "Analysis of Financial Time Series" > 3rd ed. I'm trying to replicate Table 2.1 on p.47, which gives PACF, > AIC, and BIC for the mont

Re: [R] p value in R - beginners question

data as you saw it: than that's not so exciting after all. Hope that was helpful, Ken Hutchison On Tue, Oct 18, 2011 at 2:35 PM, B77S wrote: > This is just scientific notation, so > 8.15e-01 is the same as: > > 8.15*10^-1 > [1] 0.815 > > > > > > niki wrote:

Re: [R] Party package: varimp(..., conditional=TRUE) error: term 1 would require 9e+12 columns

rtance for each variable when the others are held out (inferential only) Weak I know, but I hope it helps! Ken Hutchison On Fri, Oct 14, 2011 at 12:06 PM, Jason Roberts wrote: > I would like to build a forest of regression trees to see how well some > covariat

Re: [R] aov(variable~group*(speed*person)) but How to get 95% confidence intervals ?

Hey, If I understand correctly, library(gplots) plotmeans(). You might also try TukeyHSD() to see if that gets you where you are trying to go. Good luck! Ken Hutchison On Fri, Oct 14, 2011 at 12:11 PM, Jebb Remelius wrote: > Greetings and gratitude, > > I have 19 person

Re: [R] Perform 20 x one-way anova in 1 go

Hey, Doesn't this give you a ridiculous Type 1 error? Maybe randomly select one result and trust it. Try bestglm or stepwise regression maybe. Hope that's helpful, Ken Hutchison On Oct 14, 2554 BE, at 12:39 AM, "C.H." wrote: > This is on

Re: [R] Permutation or Bootstrap to obtain p-value for one sample

will let you know) you should be able to make inference from that using parametric methods (once) which will fit the truth a bit better than a t.test. Hope that's helpful, Ken Hutchison On Sat, Oct 8, 2011 at 10:04 AM, francy wrote: > Hi, > > I am having troubl

Re: [R] fgrep with caret (^) meta-character in system() call

man awk? I've used awk for similar tasks (if I am reading the post correctly.) Google-Fu should turn up some useful examples. Also awk should be on your linux installation in some form or another. Regards, Ken Hutchison On Oct 4, 2554 BE, at 10:52 PM, "Tom D. Harray" wro

Re: [R] Plotting a polygon with xyplot

ot;http://r.789695.n4.nabble.com/file/n3870788/111004_Lode_Outlines.csv";, header = TRUE,sep = ",",) par(mfrow = c(1, 2), pty = "s") plot(z ~ y, Data_poly, type = "l") fh <- with(Data_poly, which(z > 240)) D_poly <- rbind(Data_poly[fh, ], Data_poly[-rev(fh), ]) D_

Re: [R] Quelplot

t may be a few weeks though until I can get started on any of this though, but I will keep you informed. Regards, Ken Kenneth M. Goldberg, Ph.D. Associate Director, Nonclinical Statistics Johnson & Johnson Pharmaceutical Researc

Re: [R] Limitations of audio processing in R

Also with Linux you can add more swap memory(which I'm pretty sure R spills into if it hasn't reached it's internal limits on 32 bit installations). Windows pagefile is kind of obnoxious. Ken Hutchison On Sep 21, 2554 BE, at 5:05 PM, (Ted Harding) wrote: > Hi Ulisses

Re: [R] open source editor for r for beginners

Rstudio and Rcmdr are very popular and for good reason, find a good book while the latter is installing though. IF you are using linux, rkward is fantastic to woRk with. Ken Hutchison On Tue, Sep 20, 2011 at 11:21 AM, Steve Lianoglou < mailinglist.honey...@gmail.com> wrote: &g

Re: [R] Where to put tryCatch or similar in a very big for loop

What type of singularity exactly, if you're working with counts is it a special case? If using a Monte Carlo generation scheme, there are various workarounds such as while(sum(vec)!=0) {sample} for example. More info on the error circumstances would help. Good luck! Ken Hutchison O

Re: [R] Weibull point process

ogLambda) Point.Process.Counts=rpois(T,Lambda) return(Point.Process.Counts) } I haven't actually tried this code (may contain clerical errors) but I hope it gets you on the right track. Good luck, Ken Hutchison 2011/9/14 Torbjørn Ergon > Dear list, > > I'm looking for a fun

Re: [R] how to remove NA/NaN/Inf in a matrix??

Try: library(Hmisc) ?na.delete Ken Hutchison On Sep 11, 2554 BE, at 5:38 AM, anand m t wrote: > Hi all.. > I'm very new R, and i'm analyzing microarray data using Bioconductor.. > Recently i was given microarray data to analyze. The problem is whenever i > run MAS5 pr

Re: [R] Generating data when mean and 95% CI are known

ble to anything because there is probably more to the story distributionally. If you can assume, say a normal distribution you are ?rnorm and ?quantile away from Monte Carlo-ing a good part of the story yourself for conclusions. Best of luck, and sorry for the bad R jokes. Ken Hutchiso

Re: [R] Test for Random Walk and Makov Process

For random walk, there are entropy based tests (Robinson 1991), or you could empirically test the hypothesis by generating random normal data with the same mean and standard deviation and looking at the distribution of your quantiles. You could make generic statements also about whether or not t

[R] Fwd: Importing data from MS EXCEL (.xls) to R XXXX

-- Forwarded message -- From: Ken Hutchison Date: Wed, Aug 24, 2011 at 6:27 PM Subject: Re: [R] Importing data from MS EXCEL (.xls) to R To: Dan Abner save as csv. ?read.csv Ken On Wed, Aug 24, 2011 at 6:20 PM, Dan Abner wrote: > Hello everyone, > > Wh

[R] Fwd: help with "by" command

-- Forwarded message -- From: Ken Hutchison Date: Wed, Aug 24, 2011 at 6:06 PM Subject: Re: [R] help with "by" command To: amalka ?tapply or more specifically ?ave Hope this helps, Ken On Wed, Aug 24, 2011 at 2:51 PM, amalka wrote: > Hello, > > I

[R] require(dataset) for example.

Thank you, Ken [[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 and provide com

Re: [R] seeking advice about rounding error and %%

How about something like: If(round(x)!=x){zap} not exactly working code but might help Ken On Aug 13, 2554 BE, at 3:42 PM, Paul Johnson wrote: > A client came into our consulting center with some data that had been > damaged by somebody who opened it in MS Excel. The columns were >

Re: [R] Automating an R function call

seconds (and don't need to run another process in said console) Ken On Fri, Aug 12, 2011 at 3:41 PM, peter dalgaard wrote: > > On Aug 12, 2011, at 21:26 , Duncan Murdoch wrote: > > > On 12/08/2011 2:03 PM, RobertJK wrote: > >> Any way to run an R functi

  1   2   3   >