Re: [R] Filling empty List in a FOR LOOP

Ok. I appreciate your help. i don't understand what you guys mean by reproducible example, if you explained better, I will be more than glad to post it. Regards, -- View this message in context: http://r.789695.n4.nabble.com/Filling-empty-List-in-a-FOR-LOOP-tp4522694p4524532.html Sent from the

Re: [R] Filling empty List in a FOR LOOP

Hello, michaelyb wrote Ok. I appreciate your help. i don't understand what you guys mean by reproducible example, if you explained better, I will be more than glad to post it. Regards, Try this: R A - list(a=NULL, b=NULL) R A R dput(A) Now, the output of 'dput' is a reproducible

Re: [R] simulate correlated binary, categorical and continuous variable

Hello David Duffy-2, I see that you just proved using rmvnorm and then dichotomize/categorize them should work. Thanks but please take a look at this link; http://biostat.mc.vanderbilt.edu/twiki/bin/view/Main/CatContinuous and this article; Analysis by Categorizing or Dichotomizing Continuous

[R] Selecting files with a particular string in filename

I am trying to find an efficient way to select certain text files from a folder with thousands of text files. I used file.list function with a specific pattern to obtain 70 files for each of the i's in 1:100. Next, I need to select some files from this list of 70 files. The files to be selected do

[R] How to export plot outputs as gif

Hi, I really need to export plot image as gif. I know that other image forms are available with functions like png(), jpeg() etc. However, The reason that I only need the gif form is gimage function in gWidgetstcltk library only works with gif or pnm. I want to load some

Re: [R] Degrees of Freedom for lme.

thanks. well, basically. i ran this analysis once before, with just general linear model, but nothing stuck after multiple comparisons. the question is: do reading scores predict volume change over time (in canonical reading regions). so i tried again, using linear mixed effects, adding

[R] get data from Perl arrays

Hi, a very simple doubt, but I do not know how to manage this. I want to plot a histogram for all data in 'datos.txt'. a) by using R: datos-scan(datos.txt) pdf(xh.pdf) hist(datos) dev.off() b) How could I invoke R inside Perl to do the same?? #!/usr/bin/perl open(DAT,datos.txt); while (DAT) {

[R] NaN - trouble fixing NaN

Hi R-listers, I am using the package plyr. I am just trying to get the hatching success mean of each nesting event and have typed in the following and received the below results: tapply(HSuccess, Aeventexhumed, mean) AB C 0.2156265 0.1288559

[R] Non-linear least squares

Greetings, I am having some troubles with the nls() function in R V 2.14.2. I am doing some modelling where I want to predict the mass of leaf litter on the forest floor (X) as a function of time since fire (t). Fortunately, I have a differential equation that I can fit to the data which is

Re: [R] simulate correlated binary, categorical and continuous variable

Hello Greg, Sorry for the confusion. Lets say, I have a population. I have 6 variables. They are correlated to each other. I can get you pearson correlation, tetrachoric or polychoric correlation coefficients. 2 of them continuous, 2 binary, 2 categorical. Lets assume following conditions; Co1

Re: [R] How to adjust the stack size of R

On Mon, 2 Apr 2012, Xiaobo Gu wrote: You would need do compile your own copy of R and increase the stack size, Can we do this at runtime? That depends on your unstated OS. Uwe gave you an answer for Windows. On many other OSes, the stack size is set up by the OS when the process is

[R] linear-by-linear association model in R?

Dear all, can somebody give me some pointer how I can fit a linear-by-linear association model (i.e. loglinear model for the ordinal variables) in R? A brief description can be found here 'https://onlinecourses.science.psu.edu/stat504/node/141'. Thanks for your help

Re: [R] How to adjust the stack size of R

Hi Uwe, you said 'You need at least 76252 obs and that means the design matrix needs 46 Gbyte! ' Can you please explain me how you come up with these 2 numbers? Is the number '76252 = 76251 + 1 (1 for the intercept in the model)' somehow related with the Saturated model? Thanks, On Sun, Apr

[R] How to adjust fontsize of direct.label annotations in ggplot2?

Hello! I'm using direct.label() to annotate a ggplot, but I can't get it to change fontsize of annotations even tho I'm running the example from ?directlabels::dlcompare. Like this: library(lattice) library(ggplot2) oldopt - lattice.options(panel.error=NULL) dts -

[R] how can i get AIC value for AR model

attached file:- http://r.789695.n4.nabble.com/file/n4525253/1A1I_A.txt 1A1I_A.txt arr-ar(t$V1) str(arr) arr$aic 0 1 2 3 4 5 6 4.273791 4.397484 2.258405 3.008703 1.884487 0.00 1.476055 7 8 9

[R] [R-pkgs] New package joineR

Dear All, The 'joineR' package for the joint analysis of repeated measurements and time-to-event outcomes is now available on CRAN. The package contains utilities for creating and manipulating 'jointdata' objects, graphical summaries, a variogram function for estimating correlation structure,

Re: [R] get data from Perl arrays

On 4/2/2012 8:20 AM, José Ramón Blas Pastor wrote: Hi, a very simple doubt, but I do not know how to manage this. I want to plot a histogram for all data in 'datos.txt'. b) How could I invoke R inside Perl to do the same?? #!/usr/bin/perl open(DAT,datos.txt); while (DAT) { chomp;

Re: [R] How to adjust the stack size of R

You would need do compile your own copy of R and increase the stack size, Can we do this at runtime? That depends on your unstated OS. Uwe gave you an answer for Windows. On many other OSes, the stack size is set up by the OS when the process is started, and your OS documentation will tell

[R] date-time series plot

Hi all, i have a date-time series plot like this : Datum-Ura;padavine;P-kum-15dni;BLPV-1N;T-BLPV-1N;BLPV-4N;T-BLPV-4N;BLPV-5N;T-BLPV-5N;pretok;N-min;N-max 1.1.2011 18:00;0.0;0.0;174.62;14.60;174.56;14.10;174.64;14.20;299.26;174.18;181.78 2.1.2011

[R] R2wd: Insert MS Word field

Hello Everyone, the basic idea is to generate an initial (MS Word) report by R (plots, tables and outline structure), pass this document around so everyone can add / change content (except plots and tables) and later to be able to update the R plot / table logic and update the already modified

Re: [R] Selecting files with a particular string in filename

Doesn't grep() solve your problem? cheers, Rolf Turner On 02/04/12 10:51, Kumar Mainali wrote: I am trying to find an efficient way to select certain text files from a folder with thousands of text files. I used file.list function with a specific pattern to obtain 70 files for

[R] Make package out of own function

Hello, I already posted that on stackoverflow[1], but although it's crossposting, I think this question can probably easier to be answered by other R-users on this list, which maintain packages etc. I would like to make a package out of a function. The function is working in a script, but when

[R] partial likelihood device

Hello, I would like to get the partial likelihood device. fit - coxph(Surv(stop,event)~rx+size+number,data=bladder) Does fit$loglik give the partial likelihood device? Many thanks, Dunia [[alternative HTML version deleted]] __

[R] how to model check in cozigam

Hi to all, I'd like to check my cozigam model with a continuous response. If I only had a GAM I could do gam.check(b), but with cozigam I'm not getting a way to check for residuals, fitted values, qq plot.. Even if I have to do it separately for zeros and non-zeros part, how can I assure that the

[R] table: output: all variables in rows

Dear people, I would like to create a table out of a data.frame. How can I determine, which variables are put in the rows and which in the columns? I would like to get all the variables in the ROWS: I am including a simple example: D-data.frame(age=c(8,9,10),county=c(B,W,W)) the output should

[R] How to put latex mdframed side by side (separated by space) ?

Hi, I have one query regarding latex. I am working on latex and i need to put two mdframes side (separated by space) on A4 paper size. How can i implement it? Thanks -- View this message in context:

[R] how to read netcdf file in R

Dear R Users, I am using R 2.14.1 on windows How to read netcdf files in R ? Which packeges do we need to install for this, and what commands are used for reading netcdf files. Thanks, Best Regards, Yogesh Tiwari -- Yogesh K. Tiwari (Dr.rer.nat), Scientist, Centre for Climate Change

[R] Bootstrapped Tobit regression - get standard error 0...

I am trying to work out a bootstrapped Tobit regression model. I get the coefficients all right, but they all have standard error zero. And I am unable to figure out why. I know the coefficients are correct because that's what I get when do a Tobit (without bootstrapping). Here's my code: #

Re: [R] Bootstrapped Tobit regression - get standard error 0...

On Mon, 2 Apr 2012, Ravi Kulkarni wrote: I am trying to work out a bootstrapped Tobit regression model. I get the coefficients all right, but they all have standard error zero. And I am unable to figure out why. I know the coefficients are correct because that's what I get when do a Tobit

Re: [R] linear-by-linear association model in R?

On 4/2/2012 3:34 AM, Christofer Bogaso wrote: Dear all, can somebody give me some pointer how I can fit a linear-by-linear association model (i.e. loglinear model for the ordinal variables) in R? A brief description can be found here 'https://onlinecourses.science.psu.edu/stat504/node/141'.

Re: [R] how to read netcdf file in R

Yogesh, The package I most commonly use for reading and working with netcdfs is called ncdf. It has a collection of commands for working with netcdf files, but the open command is open.ncdf in this package Adrienne -- Adrienne Wootten Graduate Research Assistant State Climate Office of North

Re: [R] ff usage for glm

Thomas, I tried biglm and it does not work see http://r.789695.n4.nabble.com/unable-to-get-bigglm-working-ATTN-Thomas-Lumley-td2276524.html#a2278381 . There are other posts from people who cannot get biglm working and others who get strange results. Please, advise if you can help. I have row

Re: [R] NaN - trouble fixing NaN

Without seeing data, it's only a guess, but you should pass na.rm = TRUE to mean() so it will pass over NaNs in the data : tapply(HSuccess, Aeventexhumed, mean, na.rm=TRUE) Cheers, Michael On Sun, Apr 1, 2012 at 7:27 PM, Jhope jeanwaij...@gmail.com wrote: Hi R-listers, I am using the package

Re: [R] How to export plot outputs as gif

caTools::write.gif() Michael On Mon, Apr 2, 2012 at 12:41 AM, mrzung mrzun...@gmail.com wrote: Hi, I really need to export plot image as gif. I know that other image forms are available with functions like png(), jpeg() etc. However, The reason that I only need the gif form is            

Re: [R] ff usage for glm

Did you try the example described on the ff man page? On Monday, April 2, 2012, Bond, Stephen wrote: Thomas, I tried biglm and it does not work see http://r.789695.n4.nabble.com/unable-to-get-bigglm-working-ATTN-Thomas-Lumley-td2276524.html#a2278381 . There are other posts from people

Re: [R] How do I get a rough quick utility plot of a time series?

Thanks much for the reply. I did the ?ts and it seems to be for equal interval time series. Our time series are hardly ever equal interval. The time is in the first column used for the abscissa. Probably need the easiest simplest tier plot. When I say plot, it wants limits and stuff. I guess it

Re: [R] linear-by-linear association model in R?

dear Christofer, Try the following d-expand.grid(a=1:3,b=1:4) d$freq-rpois(12,5) o-glm(freq~factor(a)+factor(b)+I(a*b), family=poisson, data=d) vito Il 02/04/2012 9.34, Christofer Bogaso ha scritto: Dear all, can somebody give me some pointer how I can fit a linear-by-linear association

Re: [R] How do I get a rough quick utility plot of a time series?

Hi, On Mon, Apr 2, 2012 at 9:49 AM, Hurr hill0...@umn.edu wrote: Thanks much for the reply. I did the ?ts and it seems to be for equal interval time series. Our time series are hardly ever equal interval. The time is in the first column used for the abscissa. Probably need the easiest

Re: [R] How do I get a rough quick utility plot of a time series?

Convert your data to a zoo object and plot() that. ?zoo There's a whole bunch of documentation for plot.zoo available. Michael On Mon, Apr 2, 2012 at 9:49 AM, Hurr hill0...@umn.edu wrote: Thanks much for the reply. I did the ?ts and it seems to be for equal interval time series. Our time

[R] Unwanted page break in Rd2pdf

Hi, I want to create a pdf of my Man-pages from my package. Therefore I run in the terminal Rd2pdf on the package and a pdf of all the pages is created. After the titlepage there is the general package page, which includes Description and Details etc. Unfortunately after the Subtitle Details

Re: [R] A function like sum but with functions other than '+'

Michael Sumner mdsumner at gmail.com writes: Try ?Reduce Thanks Michael! --Sergio. __ 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

Re: [R] Unwanted page break in Rd2pdf

On 02/04/2012 9:57 AM, Johannes Radinger wrote: Hi, I want to create a pdf of my Man-pages from my package. Therefore I run in the terminal Rd2pdf on the package and a pdf of all the pages is created. After the titlepage there is the general package page, which includes Description and Details

[R] comparison of stochastic matrices

Dear List Members, i am looking for a statistical method or test wich helps me to verify the equality of two stochastic matrices (the sums in the rows = 1). Could you help me? Thanks! regards, galla -- View this message in context:

[R] Problem with fa.poly in psych package

I am attempting to use the psych package to do EFA on a set of data that mixes dichotomous and 5-value ordinal variables. I tried using fa.poly with disappointing results. When I use a data subset that consists solely of dichotomous or solely of 5-value variables, fa.poly works. However,

[R] gamm: tensor product and interaction

Hi list, I'm working with gamm models of this sort, using Simon Wood's mgcv library: gm- gamm(Z~te(x,y),data=DATA,random=list(Group=~1)) gm1-gamm(Z~te(x,y,by=Factor)+Factor,data=DATA,random=list(Group=~1)) with a dataset of about 7 rows and 110 levels for Group in order to test whether tensor

[R] Limitation on number of variables and data points in clustering

Hi All, Is there a limit on number of variables and the number of data points for clustering in R? Thanks for your help in advance! -- View this message in context: http://r.789695.n4.nabble.com/Limitation-on-number-of-variables-and-data-points-in-clustering-tp4525874p4525874.html Sent from the

Re: [R] how can i get AIC value for AR model

If you look at the documentation for ar, you will find the following description of aic: aic The differences in AIC between each model and the best-fitting model So this means that your model order is 5 (look at the pacf). If you want the actual AIC of your model: (I call your

Re: [R] How do I get a rough quick utility plot of a time series?

I tried ?zoo, but it's not available. Do I do something to make it available? Thanks Sarah, it's incomprehensible to me too. That's why I showed you the data and asked. -- View this message in context:

Re: [R] How do I get a rough quick utility plot of a time series?

Yes, zoo is a package for irregular time series but it needs to be installed (only once) and loaded (again every new session) install.packages(zoo) # Will install library(zoo) # Do this every session Michael On Mon, Apr 2, 2012 at 10:31 AM, Hurr hill0...@umn.edu wrote: I tried ?zoo, but it's

Re: [R] table: output: all variables in rows

On Apr 2, 2012, at 7:11 AM, Marion Wenty wrote: Dear people, I would like to create a table out of a data.frame. How can I determine, which variables are put in the rows and which in the columns? I would like to get all the variables in the ROWS: I am including a simple example:

Re: [R] Limitation on number of variables and data points in clustering

You should probably read the posting guide, but the answer is Yes. -- David L Carlson Associate Professor of Anthropology Texas AM University College Station, TX 77843-4352 -Original Message- From: r-help-boun...@r-project.org

Re: [R] How do I get a rough quick utility plot of a time series?

Thanks Michael, sounds right, will try later today when I have time. -- View this message in context: http://r.789695.n4.nabble.com/How-do-I-get-a-rough-quick-utility-plot-of-a-time-series-tp4522709p4526310.html Sent from the R help mailing list archive at Nabble.com.

[R] fatal error prevents multi-line functions or loop to run

Dear all, I have written some function that can go up to many lines. Two of the lines give warning errors like: Error in print()[c(1, : error in evaluating the argument 'x' in selecting a method for function 'print': Error in ... regions with no neighbours found now obviously this

[R] \Sexpr{}

Did anyone know whether the problem it doesn't evaluate \Sexpr{} in Sweave has been fixed? I checked r help page and realized it was asked a couple of years ago. But I didn't find solutions. Thanks for your help in advance! Hua [[alternative HTML version deleted]]

[R] gnm and gnlr3

Hi, I am quite new to R and would like to do nonlinear regressions with Poisson distributed data. I would like to estimate paramters of an equation of this type: FR = [c*NO * exp(a+b*NO)] / [(c+NO)*(1+exp(a+b*NO))] a,b and c are parameters, NO are input values I found both the gnm and gnlr3

Re: [R] Non-linear least squares

On 2012-04-01 17:31, n.surawski wrote: Greetings, I am having some troubles with the nls() function in R V 2.14.2. I am doing some modelling where I want to predict the mass of leaf litter on the forest floor (X) as a function of time since fire (t). Fortunately, I have a differential

Re: [R] NaN - trouble fixing NaN

Also, mean(x) will generate NaN's if both Inf and -Inf are in x. E.g., group - c(A,A,A, B,B,B, C,C, D,D) when - c( 1,2,3, -Inf,4,Inf, NA,5,NaN,6) tapply(when, group, mean) A B C D 2 NaN NA NaN tapply(when, group, mean, na.rm=TRUE) A B C D

Re: [R] rgl package broke with R 2.14.2

On Mar 28, 2012, at 10:15 AM, Duncan Murdoch wrote: On 12-03-27 6:31 PM, Grimes Mark wrote: Dear People I can't figure out how to fix this problem: rgl won't run under R 2.14.2 (it was working for me before under 2.14.0). The error message is: rgl is currently changing fairly rapidly.

Re: [R] how to read netcdf file in R

On Mon, Apr 2, 2012 at 1:54 AM, Yogesh Tiwari yogesh@googlemail.com wrote: How  to read netcdf files in R ? Which packeges do we need to install for this, and what commands are used for reading netcdf files. Typical code would look something like this: library('ncdf'') filename -

[R] Error: (subscript) logical subscript too long

Hello, I am trying to perform a logistic regression using counts. For example: cedegren - read.table(http://www.cloudstat.org/index.php?do=/attachment/download/id_95 /, header=T) attach(cedegren) ced.del - cbind(sDel, sNoDel) ced.logr - glm(ced.del ~ cat + follows + factor(class),

[R] summaryBy: transformed variable on RHS of formula?

Hi Folks, I'm trying to cut my data inside the summaryBy function. Perhaps formulas don't work that way? I'd like to avoid adding another column if possible, but if I have to, I have to. Any ideas? Thanks, Allie require(doBy) df = dataframe(a - rnorm(100), b -rnorm(100))

Re: [R] \Sexpr{}

Hi Hua, try options(SweaveSyntax=SweaveSyntaxNoweb) before sweaving your file (see A.14 in the FAQ section of the manual). In most cases R2HTML interferes with Sweave, causing your problem hth. Am 02.04.2012 15:12, schrieb Liang, Hua: Did anyone know whether the problem it doesn't evaluate

Re: [R] summaryBy: transformed variable on RHS of formula?

Sorry, no - you will have to do something like the latter suggestion. Regards Søren -Original Message- From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org] On Behalf Of Alexander Shenkin Sent: 2. april 2012 17:59 To: r-help@r-project.org Subject: [R] summaryBy:

Re: [R] summaryBy: transformed variable on RHS of formula?

On Mon, Apr 2, 2012 at 11:58 AM, Alexander Shenkin ashen...@ufl.edu wrote: Hi Folks, I'm trying to cut my data inside the summaryBy function.  Perhaps formulas don't work that way?  I'd like to avoid adding another column if possible, but if I have to, I have to.  Any ideas? Thanks, Allie

Re: [R] How to put latex mdframed side by side (separated by space) ?

Hi, as long as you don't want to stetch the frames across a page break, would minipges do? \documentclass[12pt]{article} \usepackage{lipsum} \begin{document} \begin{minipage}[t]{0.45\textwidth} \lipsum[1] \end{minipage} \hfill \begin{minipage}[t]{0.45\textwidth} \lipsum[2]

Re: [R] simulate correlated binary, categorical and continuous variable

On Sun, Apr 01, 2012 at 06:00:43PM -0700, Burak Aydin wrote: Hello Greg, Sorry for the confusion. Lets say, I have a population. I have 6 variables. They are correlated to each other. I can get you pearson correlation, tetrachoric or polychoric correlation coefficients. 2 of them

Re: [R] Nested brew call yields Error in .brew.cat(26, 28) : unused argument(s) (26, 28)

Many thanks for this. I have a follow-up question. The output that I have from the nested brew call includes output like this: NANANANANANANANAN ... then a graph or a table ... then more NANANANANANANANANNANANANANANANANANNANANANANANANANANNANANANANANANANAN ... etc. It only occurs in the

Re: [R] How to adjust the stack size of R

On 02.04.2012 10:04, Christofer Bogaso wrote: Hi Uwe, you said 'You need at least 76252 obs and that means the design matrix needs 46 Gbyte! ' Can you please explain me how you come up with these 2 numbers? Is the number '76252 = 76251 + 1 (1 for the intercept in the model)' somehow related

Re: [R] Error: (subscript) logical subscript too long

On 2012-04-02 08:52, IOANNA wrote: Hello, I am trying to perform a logistic regression using counts. For example: cedegren- read.table(http://www.cloudstat.org/index.php?do=/attachment/download/id_95 /, header=T) attach(cedegren) ced.del- cbind(sDel, sNoDel) ced.logr- glm(ced.del ~ cat +

[R] Prediction on link- and responselevel for glm

Dear R-ers, We are doing prediction of risk for aortic aneurysm in screened subjects based on different riskfactors. When we do the prediction on the linklevel and transform the prediction and the limits of the 95%- confidence intervall (with exponentiation and then transformation from

Re: [R] R process taking over memory

Oops, yes, I am using CentOS Linux 6.0, R 2.14.1 and nlme 3.1-103 I looked at the problem more carefully. For some datasets applied to nlme, nlme gets stuck in one of the iterations and the memory usage just grows and grows. nlme works by alternating between solving two conditional

Re: [R] Cox model approximaions (was comparing SAS and R survival....)

I have a question about Cox's partial likelihood approximations in coxph function of survival package (and in SAS as well) in the presence of tied events generated by grouping continuous event times into intervals. I am processing estimations for recurrent events with time-dependent covariates in

Re: [R] unable to move temporary installation

OK - so I followed the following steps, which I think rule out those causes 1) I uninstalled all remaining versions of R, and then deleted all the directories in c:\progra~1\R 2) I restarted the computer 3) I installed 2.14.2, and attempted to install the Rcmdr package. Same error message for

Re: [R] Error, Variable is missing

Thanks Michael , its works perfectly now. -- View this message in context: http://r.789695.n4.nabble.com/Error-Variable-is-missing-tp4515418p4526683.html Sent from the R help mailing list archive at Nabble.com. __ R-help@r-project.org mailing list

[R] lattice: multiple plots on the same page across pages

Dear All, I have two sets of plots: plot.1 is a set of qqplots that contains 3 pages and 1 panel per page; plot.2 is a set of xyplot that also contains 3 pages and 1 panel per page. These two sets of plots are conditioning on the same variable so I want to generate a display that shows both the

[R] Error in gamma(delta + (complex(0, 0, 1) * (x - mu))/alpha) : unimplemented complex function

I am trying to obtain the grafic of a pdf . but this error keeps showing . Here is the code MXN.fd = function(x,alpha,beta,mu,delta) { A = (2*cos(beta/2))^(2*delta) B = 2*alpha*pi*gamma(2*delta) C = (beta*(x-mu))/alpha D = abs(gamma(delta + (complex(0,0,1)*(x-mu))/alpha)^2) M = A/B*exp(C)*D M

[R] Requirement for Java Developer [REQ:104605]

[1]Click here to unsubscribe if you no longer wish to receive our emails Dear Partner, Here is our Direct client requirement which can be filled immediately. Kindly respond to this requirement with your consultant resume, contact and current location info to speed

[R] nls() error

Hello, I am running a simple nls model (which a friend ran OK) but I get the following error: Error in nls(y ~ R * (1 - (x/K)^2), data = nls.dat, start = list(R = 0.3, : object 'R_nls_iter' not found Does anyone know what the 'R_nls_iter' error is? The data are: x=1:8 ;

Re: [R] nls() error

On 02-04-2012, at 20:48, Jeff Breiwick wrote: Hello, I am running a simple nls model (which a friend ran OK) but I get the following error: Error in nls(y ~ R * (1 - (x/K)^2), data = nls.dat, start = list(R = 0.3, : object 'R_nls_iter' not found Does anyone know what the

Re: [R] \Sexpr{}

just FYI: the knitr package supports HTML as well and will not be interfered by R2HTML. See http://yihui.name/knitr/demo/minimal/ Regards, Yihui -- Yihui Xie xieyi...@gmail.com Phone: 515-294-2465 Web: http://yihui.name Department of Statistics, Iowa State University 2215 Snedecor Hall, Ames, IA

[R] Parallel writes in R

R-helpers: I'm curious what support R has for parallel writes to a binary file? If I want to use snow to have each node write different rows of a flat binary file (possibly out of sequence), are there any tricks/issues I should be aware of? --j -- Jonathan A. Greenberg, PhD Assistant Professor

Re: [R] nls() error

On 02-04-2012, at 21:20, Jeff Breiwick wrote: Thanks for that. Sorry, I did forget to give some basics: You should also reply to the list so that other can follow the thread. I am also sending this the list. Running R-1.15.0 (x64) on Windows 7. I assume you mean R2.15.0 The data I

[R] How to hide console window

Hi, all I'm wondering the code that hides or minimizes console window. What I'm making is gui for some function and that is saved as gui.RData that is started with gui window when I open the gui.RData. In this sitiuation I don't want to see console window but just gui window. How can I solve my

[R] sampling rows from a list

Hi: I'm sure this seems like a rudimentary question, but I am not well versed with R syntax for lists. I have a ragged array from which I've removed records (entire rows) with missing data. The functions I used to remove the missing cases resulted in the generation of an R list class object,

[R] STL decomposition of time series with multiple seasonalities

Hi all, I have a time series that contains double seasonal components (48 and 336) and I would like to decompose the series into the following time series components (trend, seasonal component 1, seasonal component 2 and irregular component). As far as I know, the STL procedure for

Re: [R] sampling rows from a list

?? Something like: lapply(mydata, function(x){ nr - nrow(x) x[sample(seq_len(nr),nr,rep=TRUE),] }) maybe. The idea is to use the sampled rows as your row index. -- Bert On Mon, Apr 2, 2012 at 11:24 AM, Bcampbell99 briand.campb...@ec.gc.ca wrote: Hi: I'm sure this seems like a rudimentary

Re: [R] sampling rows from a list

## recreating your data mydata-list(matrix(1:9, nrow=3, byrow=T), matrix(10:15, nrow=2, byrow=T), matrix(16:30, nrow=5, byrow=T)) ## get the shortest matrix in your list n - min(unlist(lapply(mydata, nrow))) ## subset the list into random samples of length n

Re: [R] nls() error

Thanks to Berend for his help with this. The problems seems to have been due to an environment that wasn't clean. I also gave x y that were untransformed. But I am able to get results now. Jeff __ R-help@r-project.org mailing list

Re: [R] Error in gamma(delta + (complex(0, 0, 1) * (x - mu))/alpha) : unimplemented complex function

Guaramy Guaramy at hotmail.com writes: [SNIP] i think the problem is the gamma function, does anyone know how to compute gamma with imaginary numbers? Try lngamma_complex() in the gsl package? Ben Bolker __ R-help@r-project.org mailing

Re: [R] Is it possible to de-select with sqlQuery from the RODBC library?

I've never heard of a an SQL de-select, but if there is such a thing it shouldn't be too hard to find via some web searches. In the meantime, I would probably just do a select * to get all the fields from the database tables, and then drop the unwanted ones afterwards in R. I think this will give

[R] Default parameter values in R functions?

Hi all, I have a newbie question: If I have a function with the following documentation: ca.jo(x, type = c(eigen, trace), ecdet = c(none, const, trend), K = 2, spec=c(longrun, transitory), season = NULL, dumvar = NULL) Let's take type as an example... if I omit this parameter when calling the

Re: [R] Default parameter values in R functions?

The first one -- this is implemented by match.arg(). Michael On Mon, Apr 2, 2012 at 5:30 PM, Michael comtech@gmail.com wrote: Hi all, I have a newbie question: If I have a function with the following documentation: ca.jo(x, type = c(eigen, trace), ecdet = c(none, const, trend), K =

Re: [R] Default parameter values in R functions?

Yes, it's a bit subtle. Argument matching is usually done through match.arg(), so see ?match.arg for an explanation.. (the first chosen). -- Bert On Mon, Apr 2, 2012 at 2:30 PM, Michael comtech@gmail.com wrote: Hi all, I have a newbie question: If I have a function with the following

Re: [R] fatal error prevents multi-line functions or loop to run

Sounds like you tried everything I know. However, I, and anyone else on this list whose mind-reading powers are undeveloped, really can't be sure because you didn't show us what you tried. Perhaps you should (re-)read the posting guidelines and try showing us reproducible sample code that

Re: [R] Default parameter values in R functions?

Okay! Thank you both for your help! So it chooses the 1st one by default... Thank you! On Mon, Apr 2, 2012 at 4:43 PM, Bert Gunter gunter.ber...@gene.com wrote: Yes, it's a bit subtle. Argument matching is usually done through match.arg(), so see ?match.arg for an explanation.. (the

[R] How to layout the output nicely into webpage?

How to layout the output nicely into webpage? Hi all, I have data which is about 6 x 1 divided into trunks and in the following for loop I am doing data analysis trunk by trunk. The goal is to layout the data and analysis summaries side-by-side in a clear manner. I guess the best way to

Re: [R] Default parameter values in R functions?

Just to point out another reference, which is the R Language Definition, where in section 4.3.2: http://cran.r-project.org/doc/manuals/R-lang.html#Argument-matching entitled Argument Matching, there are some pointers to the use of match.arg(), leading you to the help page pointed to by Bert

Re: [R] How to layout the output nicely into webpage?

Also, is there a way to add a plot onto each trunk, so it will become 3 parts: Left: data (xtrunk: 40 rows) Mid: plot of this data Right: lm summary outputs (around 10 rows) Any thoughts? Thanks a lot! On Mon, Apr 2, 2012 at 5:30 PM, Michael comtech@gmail.com wrote: How to layout the

[R] Calculating NOEL using R and logistic regression - Toxicology

Hello, I used the glm function in R to fit a dose-response relationship and then have been using dose.p to calculate the LC50, however I would like to calculate the NOEL (no observed effect level), ie the lowest dose above which responses start occurring. Does anyone know how to do this?

Re: [R] Calculating NOEL using R and logistic regression - Toxicology

dear Danielle, The NOEL is a threshold value or breakpoint in the range of dose. Have a look to the package segmented to estimate a GLM with unknown breakpoints. The code (untested) should be something like library(segmented) o-glm(y~1, family=binomial) os-segmented(o, ~dose, psi=starting_psi)

  1   2   >