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

[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 the other

[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

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

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 <jdnew...@dcn.davis.ca.us> wrote:

[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

Re: [R] Opposite color in R

see ww.cvrl.org under New CIE XYZ functions transformed from the CIE (2006) LMS functions best, Ken -- Kenneth 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

[R] help with ff matrix indexing and value assignment

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-project.org mailing list -- To UNSUBSCRIBE and more

[R] Are melt and dcast with ffdf possible?

(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 . Is it possible to use this method with ffdf? I appreciate any help Ken

[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

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

Andrew.Hoskins at 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

Re: [R] outer() problem

Marc Girondot marc_grt at 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

Re: [R] Generate Binary Matrix

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 PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide

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

Chen, George George.Chen at 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

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

Roland Deutsch roland.deutsch at 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

Re: [R] Select fixed number of elements

Alaios alaios at 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

Re: [R] colour code areas of a plot

Martin Batholdy batholdy at 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

Bos, Roger roger.bos at 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

Re: [R] speeding up

Ken Knoblauch ken.knoblauch at inserm.fr writes: Bos, Roger roger.bos at 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

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

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,

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

Martin Batholdy batholdy at 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

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

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 posting guide http://www.R-project.org/posting-guide.html and provide commented

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

] [,3] [1,]111 Quoting Bryan Hanson han...@depauw.edu: 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, 0.329, 0.3583) # D65 chromaticity

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 han...@depauw.edu: 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

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?

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épine 69500 Bron

Re: [R] measuring distances between colours?

. best, Ken Any suggestions would be appreciated. John John Fox Sen. William McMaster Prof. of Social Statistics Department of Sociology McMaster University Hamilton, Ontario, Canada http://socserv.mcmaster.ca/jfox/ -- Kenneth

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

Re: [R] plotting CIE chromaticity diagram?

ishi soichi soichi777 at 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,

Re: [R] plotting CIE chromaticity diagram?

ishi soichi soichi777 at 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 sermon And following up on my previous mail (diatribe), after having

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 han...@depauw.edu: I am, unfortunately, well-aware

[R] R intermittently crashes across cluster

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 the LORD, plans to prosper you

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

, starting where the variable xint is defined. best, Ken -- Kenneth 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 (0)6 84 10 64 10 http

[R] Wrong bin count number with 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 clue to get the correct count number per bin ? Regards, Ken

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

(an atmospheric pressure), rather than the name or index of the level. How to do that? TIA, Tom Roche Tom_Roche at pobox.com maybe, see ?strip.custom in lattice -- Ken Knoblauch Inserm U846 Stem-cell and Brain Research Institute Department of Integrative Neurosciences 18 avenue du Doyen Lépine 69500

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

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 gunter.ber...@gene.com: Ken: I would have thought so, too. However: x - 1:10; y - runif(1:10); f - factor(rep(1:2,5

Re: [R] history and readline, Mac OSX

://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 (0)6 84 10 64 10 http://www.sbri.fr

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

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.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE

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

Jessica da Silva jessica.m.dasilva at 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

[R] Memory errors using lmer

: 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 Daily Newshttp://www.daytondailynews.com/ | 1611 S. Main St. Dayton, Ohio 45409 Ph: 937-000-000 Cell: 937-225-2393 Email

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

for 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 Neurosciences 18 avenue du Doyen

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

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 Integrative

Re: [R] interpolation to montly data

monthly.price = ddply(dat, .(month), mean.var, var = price) monthly.price = cbind(monthly.price, month.freq = as.vector(table(df$month))) names(monthly.price) = c(month, average.price, month.freq) HTH, Ken __ R-help@r-project.org mailing list https

Re: [R] Multiple cbind according to filename

) unique.first.letter - unique(first.letter) l.files - list() for(i in 1:length(unique.first.letter)){ l.files[[i]] = as.data.frame(list.files[first.letter == unique.first.letter[i]]) } HTH, Ken __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman

Re: [R] Multiple cbind according to filename

) unique.first.letter - unique(first.letter) l.files - list() for(i in 1:length(unique.first.letter)){ l.files[[i]] = as.data.frame(list.files[first.letter == unique.first.letter[i]]) } HTH, Ken __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman

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 b...@xs4all.nl 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 t- file(C:\\Users\\Amin

Re: [R] Seeking help with LOGIT model

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 __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help

Re: [R] Seeking help with LOGIT model

0.0323 * `X 5`2.134544 2.570756 0.830 0.4064 Good luck. Ken Quoting Christofer Bogaso bogaso.christo...@gmail.com: 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

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

) 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 Institute Department of Integrative Neurosciences 18 avenue du

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

ken knoblauch ken.knoblauch at inserm.fr writes: Michael Bach phaebz at 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 vector. I looked into the apply

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

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, BP -- Ken Knoblauch Inserm U846 Stem-cell and Brain Research Institute

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

(a,b,c,d), v2 = c(a,b,e), v3 = c(a,f,g))) Many thanks -- 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

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

Stephanie Cooke cooke.stephanie at 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

Re: [R] Create subset using variable

is not created # analyses on subset take place j = j + 1 } Hope that helps, Ken __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html

Re: [R] 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 jwiley.ps...@gmail.com wrote: Hi Jonas, Look at Hadley Wickham's devtools package. It is designed with this sort of thing. That said

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

Begin forwarded message: From: Ken Hutchison vicvoncas...@gmail.com Date: January 15, 2012 8:54:49 PM EST To: Ben Bolker bbol...@gmail.com 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

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 jdnew...@dcn.davis.ca.us 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

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

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 marc_schwa...@me.com wrote

Re: [R] clear plot linear mixed model

Christof Kluß ckluss at email.uni-kiel.de writes: Am 02-01-2012 10:54, schrieb ken knoblauch: Christof Klußcklussat email.uni-kiel.de writes: lme- lme(conc ~ name/time - 1, random=conc~time|nr,method=ML,data=measurements) see plot.augPred in the nlme package thx, but how to set primary

Re: [R] clear plot linear mixed model

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 69500 Bron France tel: +33 (0)4 72 91 34 77 fax: +33 (0)4 72 91

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 riccardo.giacome...@gmail.com wrote: Hi, I need to maximize a quadratic function under constraints in R

Re: [R] R on Android

for food; else be prepared to name your phone Wilson. Ken Hutchison On Jan 1, 2012, at 11:28 AM, David Winsemius dwinsem...@comcast.net wrote: On Jan 1, 2012, at 11:07 AM, David Winsemius wrote: On Jan 1, 2012, at 10:05 AM, Federico J. Villatoro wrote: Hello,I am new

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 lig...@statistik.tu-dortmund.de 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 2 values. I am

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 rbun...@ufl.edu wrote: Hi, I am attempting to impute my data for missing values using

Re: [R] window manager interface commands for linux

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 rrast...@gmail.com wrote: How can i replicate this in Linux: source(file.choose()) I've tried source(tkgetOpenFile()) but with no luck

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

, if created, is going to be very hectic. HTH Ken Hutchison On Nov 28, 2554 BE, at 5:55 AM, cm mbnchakravar...@gmail.com 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

Re: [R] Question about randomForest

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 mattjamesfran...@gmail.com wrote: Thanks for the help. Let me explain in more detail how I think

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

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 deleted

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 rv

Re: [R] palettes for the color-blind

://statmath.wu.ac.at/~zeileis/papers/ Zeileis+Hornik+Murrell-2009.pdf HTH, Ken -- 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 (0)6 84 10 64 10

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 saschav...@gmail.com wrote: Hello I am working on my first attempt to impute missing data of a data set with systematically incomplete answers

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 flavi...@gmail.com wrote: Hi people, I´m trying to development

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 tomasm...@hotmail.com wrote: My apologies for my vague comment. My data comprises 400.000 x 21 (17 explanatory

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 jholt...@gmail.com wrote: It runs fine off a flash drive. On Thu, Oct 20, 2011 at 1:29 PM, Bogaso Christofer bogaso.christo...@gmail.com wrote: Dear all, the company I work for has Matlab installed

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 bps0...@auburn.edu wrote: This is just scientific notation, so 8.15e-01 is the same as: 8.15*10^-1 [1] 0.815 niki wrote: Dear all, i have done

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

, Ken Hutchison On Wed, Oct 19, 2011 at 7:22 PM, J Toll jct...@gmail.com 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 monthly simple returns of the CRSP value

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 j...@kin.umass.edu wrote: Greetings and gratitude, I have 19

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

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 jason.robe...@duke.eduwrote: I would like to build a forest of regression trees to see how well some covariates

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. chainsawti...@gmail.com wrote: This is one solution

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

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 francy.casal...@gmail.com wrote: Hi, I am having trouble understanding how

Re: [R] Plotting a polygon with xyplot

- with(Data_poly, which(z 240)) D_poly - rbind(Data_poly[fh, ], Data_poly[-rev(fh), ]) D_poly - rbind(D_poly, Data_poly[1, ]) plot(z ~ y, D_poly, type = n) with(D_poly, polygon(y, z, col = lightblue)) -- Ken Knoblauch Inserm U846 Stem-cell and Brain Research Institute Department of Integrative Neurosciences

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 tomdhar...@gmail.com

Re: [R] Quelplot

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 Research Development, LLC 965

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) ted.hard...@wlandres.net wrote: Hi

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: I just

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 On Sep

Re: [R] Weibull point process

) 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 t.h.er...@bio.uio.no Dear list, I'm looking for a function

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 anandro...@gmail.com 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 presence

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

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 Hutchison On Sep

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

[R] Fwd: help with by command

-- Forwarded message -- From: Ken Hutchison vicvoncas...@gmail.com Date: Wed, Aug 24, 2011 at 6:06 PM Subject: Re: [R] help with by command To: amalka ama...@gmail.com ?tapply or more specifically ?ave Hope this helps, Ken On Wed, Aug 24, 2011 at 2:51 PM, amalka ama

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

-- Forwarded message -- From: Ken Hutchison vicvoncas...@gmail.com Date: Wed, Aug 24, 2011 at 6:27 PM Subject: Re: [R] Importing data from MS EXCEL (.xls) to R To: Dan Abner dan.abne...@gmail.com save as csv. ?read.csv Ken On Wed, Aug 24, 2011 at 6:20 PM, Dan Abner

[R] require(dataset) for example.

, 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 commented, minimal

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 pauljoh...@gmail.com wrote: A client came into our consulting center with some data that had been damaged by somebody who opened it in MS Excel

Re: [R] Getting bootstrap statistic to work

, n.data.to.compute.from,replace=T) ##Insert code here to compute statistic stat.holder[count]=computed.statistic } It may be overtly simplistic, but you may find it helpful. Thanks, Ken On Fri, Aug 12, 2011 at 10:10 AM, Alex Olssen alex.ols...@gmail.com wrote: Hi R-help, I am

  1   2   3   >