Re: [R] r-help volcano plot

2011-09-05 Thread Hasan Diwan
On 6 September 2011 08:01, KnifeBoot wrote: >  Can't installe packag maDB or limma. Which R version, and what platform are you using? -- Sent from my mobile device Envoyait de mon telephone mobil __ R-help@r-project.org mailing list https://stat.ethz.

[R] r-help volcano plot

2011-09-05 Thread KnifeBoot
Can't installe packag maDB or limma. Error is shown as Using R version 2.13.1, biocinstall version 2.8.4. Installing Bioconductor version 2.8 packages: Is there any other way to draw volcano plot ? Thanks [[alternative HTML version deleted]] __

Re: [R] capturing a figure to PDF or Image

2011-09-05 Thread Eran Eidinger
The text seems to be white, reformatted: On Tue, Sep 6, 2011 at 8:53 AM, Eran Eidinger wrote: > Thank you Sarah and David, > > I only used a simple plot, and remembered to dev.off(). > attached is the session: > > > plot(c(1,2,3),c(3,2,4))> jpeg("test.jpg")> dev.off()RStudioGD > 2 > ses

Re: [R] capturing a figure to PDF or Image

2011-09-05 Thread Eran Eidinger
Thank you Sarah and David, I only used a simple plot, and remembered to dev.off(). attached is the session: > plot(c(1,2,3),c(3,2,4))> jpeg("test.jpg")> dev.off()RStudioGD 2 > sessionInfo()R version 2.13.1 (2011-07-08) Platform: x86_64-pc-mingw32/x64 (64-bit) locale: [1] LC_COLLATE=Hebre

Re: [R] plm package, R squared, dummies in panel data

2011-09-05 Thread Daniel Malter
Hi, I answered this question before in this post: http://r.789695.n4.nabble.com/Regressions-with-fixed-effect-in-R-td2173314.html , specifically in my message from May 11, 2010; 4:30pm. However, I believe the newer version of plm shows an R-squared, which should be the within R-squared. Why the pr

Re: [R] function censReg in panel data setting

2011-09-05 Thread Arne Henningsen
Dear Igors On 5 September 2011 23:58, Igors wrote: > I have a problem estimating Random Effects model using censReg function. > > small part of code: > > UpC <- censReg(Power ~ Windspeed, left = -Inf, right = > 2000,data=PData_In,method="BHHH",nGHQ = 4) > > Error in maxNRCompute(fn = logLikAttr,

Re: [R] In optim() function, second parameter in par() missing

2011-09-05 Thread colstat
Thanks, David. You suggestion worked very well. The par() in optim() only takes one argument, so I combine it into a vector. Now, I will run it with my actual code and see what happens. Colin -- View this message in context: http://r.789695.n4.nabble.com/In-optim-function-parameter-in-par-miss

[R] htmlParse hangs or crashes

2011-09-05 Thread Simon Kiss
Dear colleagues, each time I use htmlParse, R crashes or hangs. The url I'd like to parse is included below as is the results of a series of basic commands that describe what I'm experiencing. The results of sessionInfo() are attached at the bottom of the message. The thing is, htmlTreeParse a

[R] Receive "unable to load shared object RNetCDF.o" during R INSTALL of RNetCDF

2011-09-05 Thread David Brown
On a Red Hat Linux cluster I am seeing the following after multiple other packages were successfully installed. The error seems to suggest that RNetCDF.o was not copied to the appropriate lib folder. The admin user performing the install has the required privileges to perform the install. Words of

Re: [R] p values greater than 1 from lme4

2011-09-05 Thread RTSlider
Hi Ben, I've attached the data and Tinn-R file I've been using to run it, thank you for the fast response. Please let me know if you see something odd, I have no idea where to begin looking -Rob On Mon, Sep 5, 2011 at 5:30 PM, bbolker [via R] wrote: > RTSlider gmail.com> writes: > >> >> Hello, >

[R] multivariate bandwidth: regband

2011-09-05 Thread Soberon Velez, Alexandra Pilar
Hello to everybody, Somebody know how I can apply the instruction (regband) of the package (locfit) to a multivariate covariate. In the instructions states that the formula has to have only one predictor. How can I use that when I have an unknown function with two predictors, m(X,Z)? Is there

[R] function censReg in panel data setting

2011-09-05 Thread Igors
Hello all, I have a problem estimating Random Effects model using censReg function. small part of code: UpC <- censReg(Power ~ Windspeed, left = -Inf, right = 2000,data=PData_In,method="BHHH",nGHQ = 4) Error in maxNRCompute(fn = logLikAttr, fnOrig = fn, gradOrig = grad, hessOrig = hess, : N

Re: [R] Bayesian functions for mle2 object

2011-09-05 Thread Ben Bolker
On 11-09-05 02:13 PM, Billy wrote: > Dear Dr. Ben Bolker and 'JLucke', > > Thanks for your comments, but I'm still facing some problems. > For example, using the gls() function, I receive an error message and > I'm not sure I'm writing the arguments in the right way. Well, (if you still need he

Re: [R] Getting the values out of histogram (lattice)

2011-09-05 Thread Rolf Turner
On 05/09/11 21:26, Deepayan Sarkar wrote: 1. The `official' way to get panel arguments is trellis.panelArgs(); e.g., p<- histogram(~rnorm(100) | gl(2, 50), type = "density") str(trellis.panelArgs(p, 2)) List of 5 $ x : num [1:50] 0.277 1.144 1.13 -0.912 -0.892 ... $ breaks

Re: [R] p values greater than 1 from lme4

2011-09-05 Thread Ben Bolker
RTSlider gmail.com> writes: > > Hello, > I'm running linear regressions using the following script where I have > separated out species using the "IDtotsInLn" identifier > > x<-read.csv('tbl02TOTSInLn_ENV.csv', header=T) > x > attach (x)

Re: [R] Need more information about VGLM

2011-09-05 Thread Steve Lianoglou
Hi, On Mon, Sep 5, 2011 at 6:46 AM, privat NDOUTOUME wrote: > Hi, > > I'm working on multiple logistic regression. I used the function vglm > (Package VGAM) in R. Now, i'd like to implent this function (vglm) in C++. > Could someone help me or send me the algorithm. Download and extract the sou

Re: [R] what is wrong with my quicksort?

2011-09-05 Thread Daniel Nordlund
> -Original Message- > From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org] > On Behalf Of peter dalgaard > Sent: Sunday, September 04, 2011 11:09 AM > To: warc > Cc: r-help@r-project.org > Subject: Re: [R] what is wrong with my quicksort? > > > On Sep 4, 2011, at 13:1

Re: [R] How to create R executable?

2011-09-05 Thread S Ellison
Do you mean 'how do you execute the script?'? If so, the executable is called R. You might mean R -f .R , assuming R is on your path. S > -Original Message- > From: r-help-boun...@r-project.org > [mailto:r-help-boun...@r-project.org] On Behalf Of Pariksheet > Sent: 05 September 2011

[R] package gmp installation problem

2011-09-05 Thread Giacomo Bezzi
Hello everybody, Trying to install the package gmp I get the following errors and fail to install: Error in dyn.load(file, DLLpath = DLLpath, ...) : unable to load shared object '/usr/local/lib/R/site-library/gmp/libs/gmp.so': libgmp.so.10: cannot open shared object file: No such file or dire

Re: [R] In optim() function, second parameter in par() missing

2011-09-05 Thread David Winsemius
On Sep 5, 2011, at 11:41 AM, David Winsemius wrote: On Sep 5, 2011, at 10:30 AM, colstat wrote: Hi, First time using the optim(), can someone please tell me what I am doing wrong? The error looks like this Error in .Internal(pnorm(q, mean, sd, lower.tail, log.p)) : 'sd' is missing Yo

Re: [R] How to create R executable?

2011-09-05 Thread Uwe Ligges
On 05.09.2011 12:04, Pariksheet wrote: Hi , I have created .R file which connects to Teradata database and then does some manipulation and produces the output graph. How to create executable for .R file? You cannot. Or do you mean you want to create an R package? In any case, R code is alw

Re: [R] capturing a figure to PDF or Image

2011-09-05 Thread David Winsemius
On Sep 5, 2011, at 11:14 AM, Eran Eidinger wrote: Sorry guys, just removed, the "at" prameters. The warnings are gone now, but I still get an empty image, no matter which function i use (pdf, jpeg, bmp...) Any clue? The most common error is failing to call dev.off() after the plotting is

[R] SAS code in R

2011-09-05 Thread Bonnett, Laura
Dear all, I was wondering if anyone can help? I am an R user but recently I have resorted to SAS to calculate the probability of the event (and the associated confidence interval) for the Cox model with combinations of risk factors. For example, suppose I have a Cox model with two binary vari

[R] p values greater than 1 from lme4

2011-09-05 Thread RTSlider
Hello, I'm running linear regressions using the following script where I have separated out species using the "IDtotsInLn" identifier x<-read.csv('tbl02TOTSInLn_ENV.csv', header=T) x attach (x) library(lme4) rInLn<-lm

Re: [R] In optim() function, second parameter in par() missing

2011-09-05 Thread David Winsemius
On Sep 5, 2011, at 10:30 AM, colstat wrote: Hi, First time using the optim(), can someone please tell me what I am doing wrong? The error looks like this Error in .Internal(pnorm(q, mean, sd, lower.tail, log.p)) : 'sd' is missing You should be using a textbook. Which one are you consul

Re: [R] Hessian Matrix Issue

2011-09-05 Thread Ben Bolker
Uwe Ligges statistik.tu-dortmund.de> writes: > > I have not really looked into the details of the lengthy and almost > unreadable code below. In any case, there are good reasons why numerics > software typically uses Fisher scoring / IWLS in order to fit GLMs. > > And if your matrix is th

Re: [R] sorry, [WinBUGS] question

2011-09-05 Thread Uwe Ligges
On 05.09.2011 16:08, Kehl Dániel wrote: Dear Community, I know this is not the place to ask WinBUGS questions, but I did not get any answers on other lists. I am rather new to the BUGS language and to bayesian modeling, excuse me for probably simple questions. I have to conduct a bayesian meta

Re: [R] help with glmm.admb

2011-09-05 Thread Ben Bolker
Elizabeth C Eadie unm.edu> writes: > > > R glmmADMB question > I am trying to use glmm.admb (the latest alpha version > from the R forge website 0.6.4) to model my count data > that is overdispersed using a negative binomial family but > keep getting the following error message: > > Error

Re: [R] capturing a figure to PDF or Image

2011-09-05 Thread Sarah Goslee
Hi Eran, To be able to help, we need at a minimum a reproducible example and the output of sessionInfo(). Sarah On Mon, Sep 5, 2011 at 11:14 AM, Eran Eidinger wrote: > Sorry guys, just removed, the "at" prameters. The warnings are gone now, but > I still get an empty image, no matter which func

Re: [R] capturing a figure to PDF or Image

2011-09-05 Thread Eran Eidinger
Sorry guys, just removed, the "at" prameters. The warnings are gone now, but I still get an empty image, no matter which function i use (pdf, jpeg, bmp...) Any clue? Thanks, Eran. On Mon, Sep 5, 2011 at 6:09 PM, Eran Eidinger wrote: > Hello, > > I've been using jpeg(), bmp() and pdf() to captur

[R] capturing a figure to PDF or Image

2011-09-05 Thread Eran Eidinger
Hello, I've been using jpeg(), bmp() and pdf() to capture plots. I've used the parameter "at" in a plot, to define the tickmarks. It works fine on screen, but when I try to print it to a file, it gives a warning: "at" is not a graphical parameter and prints an empty figure. Can you help? Than

[R] KernSmooth: dpill

2011-09-05 Thread Soberon Velez, Alexandra Pilar
Hello, I trying to find a way to calculate a bivariate bandwidth for a bivariate local polynomial regression and I have a question. Somebody know if I can use the instruction "dpill" to do it? When I read the instructions I see that "x" must be a vector of variables, but if I state "x" is

Re: [R] Bayesian functions for mle2 object

2011-09-05 Thread Billy
Dear Dr. Ben Bolker and 'JLucke', Thanks for your comments, but I'm still facing some problems. For example, using the gls() function, I receive an error message and I'm not sure I'm writing the arguments in the right way. Instead, I thought about my original models and realized that I was modelli

[R] Need more information about VGLM

2011-09-05 Thread privat NDOUTOUME
Hi,  I'm working on multiple logistic regression. I used the function vglm (Package VGAM) in R. Now, i'd like to implent this function (vglm) in C++.  Could someone help me or send me the algorithm.  Thanks in advance,  Quentin  [[alternative HTML version deleted]] ___

[R] Dealing with NA's in a data matrix

2011-09-05 Thread Adam Hesse
hello... I am trying to get this code to work, but as I get to the predict command, it displays an error due to the length of the data sets from the removal of the NA's. Here is the data, and the code that I am using so far, if you run it, you'll see the error pop up please help me to get aro

[R] In optim() function, second parameter in par() missing

2011-09-05 Thread colstat
Hi, First time using the optim(), can someone please tell me what I am doing wrong? The error looks like this Error in .Internal(pnorm(q, mean, sd, lower.tail, log.p)) : 'sd' is missing An example of the error dat = c(20, 19, 9, 8, 7, 4, 3, 2) dat_mu=mean(dat) dat_s=sd(dat) max.func = funct

[R] ggplot2-grid/viewport and PNG

2011-09-05 Thread ashz
Dear All, The following code save my graphs as pdf: pdf("j:/mix.pdf", width = 18, height = 16) grid.newpage() pushViewport(viewport(layout = grid.layout(3,1))) vplayout <- function(x, y) viewport(layout.pos.row = x, layout.pos.col = y) print(Aplot, vp = vplayout(1, 1)) print(Bplot, vp

[R] How to create R executable?

2011-09-05 Thread Pariksheet
Hi , I have created .R file which connects to Teradata database and then does some manipulation and produces the output graph. How to create executable for .R file? Thanks Pariksheet -- View this message in context: http://r.789695.n4.nabble.com/How-to-create-R-executable-tp3790883p3790883.htm

[R] Implent the function vglm in C++

2011-09-05 Thread Quentin
Hi, I'm working on multiple logistic regression. I used the function vglm (Package VGAM) in R. Now, i'd like to implent this function (vglm) in C++. Could someone help me or send me the algorithm. Thanks in advance, Quentin -- View this message in context: http://r.789695.n4.nabble.com/Impl

Re: [R] Alternatives to integrate?

2011-09-05 Thread . .
Sorry about the indentation, it seems I'm doing some mistake when sending the e-mail. I've tried to indent here like the Google's style, but it becomes a garbage in the e-mail. On Mon, Sep 5, 2011 at 11:27 AM, . . wrote: > Hi, continuing the improvements... > > I've prepared a new code: > > ddae

Re: [R] Alternatives to integrate?

2011-09-05 Thread . .
Hi, continuing the improvements... I've prepared a new code: ddae <- function(individuals, frac, sad, samp="pois", trunc=0, ...) { dots <- list(...) Compound <- function(individuals, frac, n.species, sad, samp, dots) { print(c("Size:", length(individuals), "Compound individuals:", individ

[R] Power analysis in hierarchical models

2011-09-05 Thread Tom Wilding
Dear All I am attempting some power analyses, based on simulated data. My experimental set up is thus: Bleach: main effect, three levels (control, med, high), Fixed. Temp: main effect, two levels (cold, hot), Fixed. Main effect interactions, six levels (fixed) For each main-effect combination I ha

[R] sorry, [WinBUGS] question

2011-09-05 Thread Kehl Dániel
Dear Community, I know this is not the place to ask WinBUGS questions, but I did not get any answers on other lists. I am rather new to the BUGS language and to bayesian modeling, excuse me for probably simple questions. I have to conduct a bayesian meta-analysis of some data. We have collecte

Re: [R] help with installing tar.gz package

2011-09-05 Thread Kay Cecil Cichini
i finally found the XML.zip: http://www.stats.ox.ac.uk/pub/RWin/bin/windows/contrib/2.14/ and thanks to you pointing me at the URL of the RCurl.zip: http://www.stats.ox.ac.uk/pub/RWin/bin/windows/contrib/2.13/RCurl_1.6-10.1.zip i was able to to install the wanted RGooglesDoc package successfully

Re: [R] generating multiple dataset and applying function and output multiple output dataset......

2011-09-05 Thread Sarah Goslee
Hi, On Sun, Sep 4, 2011 at 9:25 AM, John Clark wrote: > Dear R experts: > > Here is my problem, just hard for me... > > I want to generate multiple datasets, then apply a function to these > datasets and output corresponding output in single or multiple dataset > (whatever possible)... > > # my e

Re: [R] help with installing tar.gz package

2011-09-05 Thread Uwe Ligges
On 05.09.2011 14:52, Kay Cecil Cichini wrote: internet connection exists - dont' know why, but i can't set to CRAN mirror - maybe because i run R from an usb-stick, firewall or whatever. Probably not correctly configured proxy setting? thus, install.packages("RCurl_0.91-0.tar.gz") won't wor

Re: [R] glht (multcomp): NA's for confidence intervals using univariate_calpha (fwd)

2011-09-05 Thread Torsten Hothorn
fixed @ R-forge. New version should appear on CRAN soon. Thanks for the report! Torsten -- Forwarded message -- Date: Sat, 3 Sep 2011 23:56:35 +0200 From: Ulrich Halekoh To: "r-help@r-project.org" Subject: [R] glht (multcomp): NA's for confidence intervals using univaria

Re: [R] help with installing tar.gz package

2011-09-05 Thread Kay Cecil Cichini
internet connection exists - dont' know why, but i can't set to CRAN mirror - maybe because i run R from an usb-stick, firewall or whatever. thus, install.packages("RCurl_0.91-0.tar.gz") won't work. a zip-file does not exist, only the tar.gz file. so i tried to install from the source that i

Re: [R] help with installing tar.gz package

2011-09-05 Thread Uwe Ligges
On 05.09.2011 13:52, Kay Cecil Cichini wrote: hi, i'd like to install the package "RGoogleDocs ". i downloaded to path "E:/R/R-2.13.0/library/RCurl_0.91-0.tar.gz" i run R from an usb-stick and can't get the install.packages() prompt to run correctly - can anyone help with this? 1. Why not

[R] help with installing tar.gz package

2011-09-05 Thread Kay Cecil Cichini
hi, i'd like to install the package "RGoogleDocs ". i downloaded to path "E:/R/R-2.13.0/library/RCurl_0.91-0.tar.gz" i run R from an usb-stick and can't get the install.packages() prompt to run correctly - can anyone help with this? thanks, kay sessionInfo() R version 2.13.0 (2011-04-13

[R] plm package, R squared, dummies in panel data

2011-09-05 Thread Cecilia Carmo
Hi R-helpers, I have two questions I hope you could help me with them: In the plm package how can I calculate the R2 within, R2 between and R2 overall? Is there any special reason to not display these values? When using first differences do I need to have some special care with dummies

[R] plot regression line, log-transformation

2011-09-05 Thread Johannes Radinger
Hello UseRs, I've somehow general questions. I've got a dataset which shows signs of heteroscedasticity and non-normality in errors if I do a normal linear regression of the form Y~X. So to things came into my mind, either transforming the variables (log or log10) or using robust regression. So

Re: [R] about raw type

2011-09-05 Thread Jim Holtman
it is useful if you need to process your data 'byte by byte'. I have used it to read in a large file and then translate one byte value to another. you can learn by creating some small script and processing example data with them. do you have a specific application tha you want to use it on?

[R] glm

2011-09-05 Thread Samuel Okoye
Dear all, I am using glm with quasibinomial. What does the following error message mean: Error in eval(expr, envir, enclos) : y values must be 0 <= y <= 1 Does it mean that the predictor variable should only have zero and one or it is also possible to have continuous values between zero and one

Re: [R] calling R from PHP error

2011-09-05 Thread Uwe Ligges
On 05.09.2011 09:40, Katerina Karayianni wrote: Yes, you were right about the make install command for R, now it is ok. The new issue is that when the R script is trying to generate a plot file, the following error occurs Error in function (file = ifelse(onefile, "Rplots.pdf", "Rplot%03d.pdf"

Re: [R] Getting the values out of histogram (lattice)

2011-09-05 Thread Deepayan Sarkar
On Thu, Sep 1, 2011 at 10:29 AM, Rolf Turner wrote: > > 'Scuse me, but I don't see anything in your example relating to what the OP > asked for. She wanted to get at the ``actual data defining the histogram'', > which > I interpret as meaning the bar heights (the percentages, density values, or >

[R] Stemming functions only work on the last word of plain text documents

2011-09-05 Thread Christian Timmermann
Hello, I want to use the SnowballStemmer on a collection of plain text documents. However, when I apply it to my corpus using the tm_map function it only stems the last word of each document (The problem is the for wordStem and stemDocument does not work at all).  An example: > path <- c("c:

Re: [R] calling R from PHP error

2011-09-05 Thread Katerina Karayianni
Yes, you were right about the make install command for R, now it is ok. The new issue is that when the R script is trying to generate a plot file, the following error occurs Error in function (file = ifelse(onefile, "Rplots.pdf", "Rplot%03d.pdf"), : cannot open file 'Rplots.pdf' Calls: plot .

Re: [R] Hysteresis modeling and simulation

2011-09-05 Thread Yvonnick Noel
Hi Bill, I once modelled a hysteresis phenomenon (on binary data) with a simple logistic model. I am not sure I understand how this pattern appears in your data, but in my previous analyses, it appeared as an order effect: The response increased in probability later with increasing than with

Re: [R] Gradients in optimx

2011-09-05 Thread Rubén Roa
Well, I guess this doesn't make necessary for me to prepare a report with the CatDyn package. However, I am available to test a new optimx R-forge version with my package. Cheers Rubén -- Dr. Ruben H. Roa-Ureta Senior Researcher, AZTI Tecnalia, Marine Research Division, Txatxarramendi Ugartea z