[R] stepwise on generlized linear model

Hi- Can anybody give me a quick answer that whether step() can be directly applied to a glm model including both continuous variables and categorical variables (by factor())? Maybe my data is not suitable, I applied step() on my glm object and did not get any reduction (i.e., returned full mode

[R] stack overflow using step()

Hi, I'm using step() for forward regression with 680 main effects and the correspomding 230860 interactions. However, the calculation is stopped with the message: Error: protect(): stack overflow How can I specify the maximum size of the pointer protection stack when running R not from command lin

Re: [R] measurement error model - "simple" linear regression

Dear Troels, you might try McArdle, Limnol. Oceanogr., 48(3), 2003, 1363-1366 for some reading, and the pls package (newly updated) for tools. (I apologize if this information replicates previous postings, but Googling R-help McArdle 2005 draws a blank so I think I'm on safe ground!) Good luck!

[R] Error in lazyLoadDBfetch

R was working fine for a while. However, now after I get some data from the web and try to "write", "write.table" or "save.image", I get this error: Error in lazyLoadDBfetch(key, datafile, compressed, envhook) : file open failed I may have downloaded additional packages in between the time it

Re: [R] SEM with dichotomous indicators

Dear Renata, You can do this by combining the sem() and boot.sem() functions in the sem package with hetcor() in the polycor package; boot.sem() computes standard errors for parameter estimates by bootstrapping. See ?boot.sem for an example of a confirmatory factor-analysis model with ordinal ind

Re: [R] grid.edit problem

Hi Gabor Grothendieck wrote: > Thanks, again. I see its basically an origin problem so that, for > example, this would have worked too: > > grid.edit("L", x = grid.get("L")$x + unit(10, "native") - unit(0, "native")) Right. And your solution might be easier for people to read and write too :)

Re: [R] Need help on ARIMA (time series analysis)

Regarding the loop, the code below: > ws <- i > we <- i+w-1 > temp <- th16k[ws:we] can be written as: temp <- th16k[seq(i, length = w)] or you can get rid of the loop entirely using embed or using running in package gtools or rapply in the zoo package. > park wrote: > > H

Re: [R] acf.plot() question

Have you received a reply to this? I haven't seen one. If not, you might have received a quicker reply if you had noted that the example came from ch. 14 of Venables and Ripley (2002) Modern applied Statistics with S, 4th ed. (Springer, pp. 390-391, to be precise). You ask

[R] measurement error model - "simple" linear regression

Dear friends, I found the thread on this subject this summer but wonder whether it has been taken any further. I have an important medical problem where X is computed from a three independent and complicated measurements (exchangeable sodium and potassium and total body water - i.e. X = (Nae+Ke

Re: [R] Need help on ARIMA (time series analysis)

Have you received a reply to this post? I have not seen one. If you would still like some suggestions from this group, please provide more detail on your question, as requested in the posting guide (www.R-project.org/posting-guide.html), including which version of R under which oper

Re: [R] grid.edit problem

Thanks, again. I see its basically an origin problem so that, for example, this would have worked too: grid.edit("L", x = grid.get("L")$x + unit(10, "native") - unit(0, "native")) By the way, could unit.c, unit.rep and unit.length be named c.unit, rep.unit and length.unit (since then they would

Re: [R] boxplot: how to display x-labels vertically

On Sun, 9 Oct 2005, Steven T. wrote: > I need to draw a boxplot for a dataset consisting about 60 groups, the > problems is by default the x labels are drawn horizontally, so it is > impossible to display the names of all 60 groups. However, if the labels are > drawn vertically, it should be possi

Re: [R] finding roots of a polynomial

?polyroot, and watch that x^4 appears to be missing. On Fri, 7 Oct 2005, Vumani Dlamini wrote: > Dear listers, > I am trying to find the root of a polynomial of the form, > a+bx+cx^2+dx^3+ex^5+...=0 > where a,b,c,d,e are known. > I have looked at "quadprog" without any success. -- Brian

[R] [R-pkgs] pls version 1.1-0

Version 1.1-0 of the pls package is now available on CRAN. The pls package implements partial least squares regression (PLSR) and principal component regression (PCR). Features of the package include - Several plsr algorithms: orthogonal scores, kernel pls and simpls - Flexible cross-validation

[R] looking for a consultant/trainer with expertise in R

We are looking for a consultant/trainer with expertise in R. (Ideally, the person will also be familiar with “Microsoft SQL”). Training can be by phone, email, instant messenger and/or in person. (We are located in San Francisco.) Please respond to [EMAIL PROTECTED] if interested.

Re: [R] Hmisc latex function

Charles Dupont <[EMAIL PROTECTED]> writes: > For a temporary fix source the attached file after loading the Hmisc > library. I will fix this in source for the next version. Or, slightly shorter but somewhat dirtier: evalq(dQuote <- shQuote, environment(latex)) -- O__ Peter Dalgaard

[R] STATIS

Hi, i'm having trouble using Statis in ADE4 package. I want to study a matrix with 18x414...could any body help me?!?!' Urgent Thanks Miguel Ribeiro Portugal [[alternative HTML version deleted]] __ R-help@stat.math.ethz.ch mailing list https:/

[R] SEM with dichotomous indicators

Hello, I'd like to know if there is a way to fit a Structural equation model with dichotomous indicators (ex: problem with a phone solved/ or not) having effects on a ordinal variable. How I do that using R? Do you have an example with the code in R that you can send to me? Thanks a lot!

Re: [R] plot.augPred sorted and labelled according second factor

Hi, Christoph: Have you received a satisfactory reply to this? If no, you might wish to consider the following; with a few more hours, I could probably produce what you seem to be asking. If you try something further, please ping the authors and maintainers for augPred and xyplot (

[R] boxplot: how to display x-labels vertically

Hi, I need to draw a boxplot for a dataset consisting about 60 groups, the problems is by default the x labels are drawn horizontally, so it is impossible to display the names of all 60 groups. However, if the labels are drawn vertically, it should be possible to display all names. How can this b

[R] finding roots of a polynomial

Dear listers, I am trying to find the root of a polynomial of the form, a+bx+cx^2+dx^3+ex^5+...=0 where a,b,c,d,e are known. I have looked at "quadprog" without any success. Regards, Vumani __ R-help@stat.math.ethz.ch mailing list https://stat.eth

Re: [R] grid.edit problem

Hi Gabor Grothendieck wrote: > I am having a problem in editing a grob. It works ok if I try to > shift the grob using npc coordinates but if I do the same thing > using native coordinates the grob disappears. What is wrong? > > > library(grid) > grid.newpage() > > # create viewport > pushVi

Re: [R] repeated measures with random effects

Have you tried 'RSiteSearch("repeated measures with random effects")'? I just got 74 hits, some of which might interest you. Also, have you reviewed Pinheiro and Bates (2000) Mixed-Effects Models in S and S-Plus (Springer)? This book, the lme function in library(nlme) and

[R] TRAMO-SEATS confusion?

Hi, For what its worth I had to hack some horrendous old FORTRAN which I could not come close to understanding. The main issue was to allow dynmically allocated arrays. What I did was to run it through f2c to produce C++ code. Then I verified that it still worked. After that I replaced the pointer

Re: [R] TRAMO-SEATS confusion?

LICENSING AND OWNERSHIP CRAN wants the source code for any contributed package, both for security reasons (protection against viruses) and for consistency wiht the GNU license, which requires distributors of software using GNU software to extend the GNU restrictions to their derivativ

Re: [R] where is Defn.h and implementation

On Sun, 16 Oct 2005, jiesheng zhang wrote: > I am in the process of writing my first R extesion package. Have you read our posting guide? Please ask C programming questions on the apppropriate list -- R-devel. Have you read `Writing R Extensions'? It does cover this topic with copious warnin

Re: [R] Sorting a data frame by one of the variables

Leaf, using your example data as 'dat' below -- dat<-read.table("clipboard",header=T) dat XYZ 1 22.0 24.0 4.3 2 2.3 3.4 5.3 3 57.2 23.4 34.0 #to order the data frame by say X (for column 1)-- dat1<-dat[order(dat[,1]),] dat1 XYZ 2 2.3 3.4 5.3 1 22.0 24.0 4.3

Re: [R] asking the user for data

Use readline for a text dialog or use tkentry from the tcltk package as shown in the Edit Box example at: http://bioinf.wehi.edu.au/~wettenhall/RTclTkExamples/ On Windows other possibilities exist as well: winDialogString and in package svDialog there is guiDlg . Also check out qnorm . On 10/

Re: [R] asking the user for data

library(sm) > z <- ask(message="Please enter the z value for") Please enter the z value for: 1.96 > z [1] 1.96 I found this using RSiteSearch("ask", restrict="functions") . Iain Gallagher wrote: > Hello everyone. > > How do I get R to ask users for data to be entered? Specifically I want >

Re: [R] asking the user for data

Dear Iain, There's an ask() function in the sm package that does what you want, but you'll have to compose your message properly: ask(paste("Please enter the z value for", x)) Alternatively, eval(parse(text=readline(paste("Please enter the z value for ", x, ": ", sep="" will do what you want,

[R] asking the user for data

Hello everyone. How do I get R to ask users for data to be entered? Specifically I want to ask for a z score to be entered (the user would look this up in a table) and then use the entered data to compute a Dunn's post-hoc test (post kruskal.test). I've tried the "ask" function but it's not re

[R] where is Defn.h and implementation

Hi, all I am in the process of writing my first R extesion package. One of my function uses StrToInternal function(I copy it from other extension package). I noticed the StrToInternal was defined in "Defn.h". I have two questions: 1)should "Defn.h" be under $R_HOME/include. I could not find it

Re: [R] write R extenesion issue

Hi, duncan I remove the built and packaged line from DESCRIPTION and install from source. It worked right now Thanks -jason Duncan Murdoch wrote: > jiesheng zhang wrote: >> I found my attachment was discarded. >> My R package structure lis listed here >> - >> btRRTest

Re: [R] How to get to the varable in a list

Hi Lisa, Try the following -- there are probably better solutions, but these seem to work. Hank Stevens Miami University Oxford, OH, USA a <- vector("list",3) names(a) <- 1:3 a[[1]] <- data.frame(x1=3,x2=1) a[[2]] <- data.frame(x1=c(3,3),x2=c(2,2)) a[[3]] <- data.frame(x1=c(3,3), x2=c(3,3)) a u

Re: [R] enter a survey design in survey2.9

svydesign needs to be told what data frame to use, via the data= argument. I get a slightly different error message from you: Error in eval(expr, envir, enclos) : object "subdiv" not found but this may be because I'm using survey version 3.3.1, not 2.9. There are several new features (listed at h

Re: [R] Your order (#gzmgss-R-help) Is somebody's machine infected to zombie this trash??

Brian Lunergan wrote: > Dennis Ewing wrote: > >> Ciialis Softt Tabss actts up to 36 houurs! Simplyy disollve halff a >> tablet underr youur tonngue 15 minuutes beforee s.ex. >> >> Alll orderrs willl be processed and dispatched withiin 24hrs. >> 100 pills - 1.98

Re: [R] solve() versus ginv()

Werner Bier <[EMAIL PROTECTED]> writes: > Dear All, > > While inverting a matrix the following error appears on my console: > > Error in solve.default(my_matrix) : Lapack routine dgesv: system is exactly > singular > > With this respect, I have been replacing the solve() function with ginv(