[R] include C functions from nmath in my own C functions

Hi: I followed the README in src/nmath/standalone/ to make the use the command make shared to make the libRmath.so file. I also add the directories containg libRmath.so to LD_LIBRARY_PATH by using command export LD_LIBRARY_PATH=$LD_LIBRARY_PATH:$/home/zhliu/Backup/R-2.0.1/src/nmath/standalon

Re: [R] connexion problem getHdata (HMisc)

Sorry about the missing info, I'm just trying to do too many things at the same time! I installed R on my Windows2000 computer (no choice there) . version platform i386-pc-mingw32 arch i386 os mingw32 system i386, mingw32 status major2 minor1.1 year 2005 month

Re: [R] background colors in image()

Le 11.08.2005 00:59, array chip a écrit : Thanks for the suggestion! It works in a way that the entire graph window is in the background color, is there a way to only have the plotting area (i.e. the area within the axis box in the background color, but leave the area outside the axes to be

Re: [R] connexion problem getHdata (HMisc)

Hi Anne, does update.packages() work, or does it just hang? cheers! Sean On 11/08/05, [EMAIL PROTECTED] [EMAIL PROTECTED] wrote: Sorry about the missing info, I'm just trying to do too many things at the same time! I installed R on my Windows2000 computer (no choice there) . version

Re: [R] include C functions from nmath in my own C functions

First, please read the posting guide: it says C programming questions should go to R-devel so it seems you have not done so. Second, you have not told us your OS nor what you did to run your code, so we have to guess. But sin etc are not 'R functions' but C mathematical functions, and the

Re: [R] connexion problem getHdata (HMisc)

Hi Sean, Your intuition is spots on: update.package also hang: Message d'avis : connexion à 'cran.ch.r-project.org' impossible sur le port 80 Avis : unable to access index for repository http://cran.ch.r-project.org/bin/windows/contrib/2.1 (sorry for the french...) I did not try it before as

[R] how to disable RODBC log file

I am using RODBC with Postgresql. I have disabled the log file in ~/.odbc.ini, but nevertheless I still get a logfile in /tmp, which can grow to gigabytes, so I would like to disable it. It appears that the logfile is generated by R: % lsof | grep /tmp/psqlodbc_tpapp7482.log R 7482

Re: [R] connexion problem getHdata (HMisc)

Hi Sean, Thanks for the tip, I will try to set up the command line in the target field as recommended in the faq Anne Sean O'Riordain [EMAIL PROTECTED] 11.08.2005 08:36 Veuillez répondre à [EMAIL PROTECTED] A [EMAIL PROTECTED] [EMAIL PROTECTED] cc Objet Re: [R] connexion problem

[R] conjoint optim

Hi, have anybody a starting point or experince how i could optimize the results from a conjoint-ananlysis. I have the idea to optimize this with 2 function , one for the parth-worth utilities which have to maximize and a secondconcurrent for the cost's of the parth-worth utilities

Re: [R] off-topic ODBC Q (was how to disable RODBC log file)

RODBC does not have a log file, and this is nothing to do with RODBC or R. What you are talking about here is an *ODBC* log file, and for that you need to consult your (unstated) device manager documentation for your (unstated) OS. The name of the log file (psqlodbc_*) should have given you a

[R] getting xpinch and ypinch

Hello! I would like to know if it is possible (and how) to get xpinch and ypinch values that R gets (from the system) for drawing plots. If I understand correctly; with this values I could save the same picture as I see on computer 1 on computer 2 as a let say wmf file. Thank you in advance,

Re: [R] rpart plot question

Petr Pikal wrote: Dear all I am quite confused by rpart plotting. Here is example. set.seed(1) y - (c(rnorm(10), rnorm(10)+2, rnorm(10)+5)) x - c(rep(c(1,2,5), c(10,10,10)) fit - rpart(x~y)## NB should be y~x plot(fit) text(fit) Text on first split says x 3.5 and on the second

[R] sub set selection

hi all is there a package that undertakes subset selection but BASED ON AIC or any other information criteria. i've seen the subselect and the leaps package but i have not played around with them yet. thanx__ R-help@stat.math.ethz.ch mailing list

[R] tcltk programming guide

Dear List, I'm looking for some documentation about the R tcltk package The one I found in the help doesn't look exaustive, I need information on the use of the single tk widget, maybe with some examples thank you, simone gabbriellini __

Re: [R] tcltk programming guide

Go to http://bioinf.wehi.edu.au/~wettenhall/RTclTkExamples/. There are good examples. X-Original-To: [EMAIL PROTECTED] Delivered-To: [EMAIL PROTECTED] X-Virus-Scanned: by amavisd-new at stat.math.ethz.ch Mime-Version: 1.0 (Apple Message framework v733) To: R-help@stat.math.ethz.ch From: Simone

[R] Error in autoloader

Hi, After installing the latest versions of lme4, Matrix, VR on a Debian Box, I have run into a problem. When I use library(lme4) on R Version 2.1.0 (2005-04-18) I get the following error message: library(lme4) Loading required package: Matrix Error in autoloader(name = confint, package =

Re: [R] sub set selection

in what model, glm or gam? I believe you can use aic in both. On 8/11/05, Clark Allan [EMAIL PROTECTED] wrote: hi all is there a package that undertakes subset selection but BASED ON AIC or any other information criteria. i've seen the subselect and the leaps package but i have not played

[R] translation Hmisc Design

** This email and any files transmitted with it are confidenti...{{dropped}} __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the

Re: [R] repeated - R package - Compilation Error

On Wed, 10 Aug 2005 [EMAIL PROTECTED] wrote: Thompson's Manual to Accompany Agresti's book refers to a package named repeated. It's not on CRAN from what I can see. I have seen rpm's for it. Where is the best place to download this package? See the FAQ, Q5.1.5. As to the `best' place, it

Re: [R] repeated - R package - Compilation Error

On Wed, 10 Aug 2005 [EMAIL PROTECTED] wrote: Thompson's Manual to Accompany Agresti's book refers to a package named repeated. It's not on CRAN from what I can see. I have seen rpm's for it. Where is the best place to download this package? See the FAQ, Q5.1.5. As to the `best' place, it

Re: [R] sub set selection

just normal linear model: multiple regression in particular Wensui Liu wrote: in what model, glm or gam? I believe you can use aic in both. On 8/11/05, Clark Allan [EMAIL PROTECTED] wrote: hi all is there a package that undertakes subset selection but BASED ON AIC or any other

Re: [R] repeated - R package - Compilation Error

When I go into binnest.f, here is what I see: C Calculate First Derivative Do kk = 1, t1 Gradient(kk) = Gradient(kk) + (D1_beta(kk) / Li) End do Do kk = 1, t2 Gradient(t1+kk) = Gradient(t1+kk) + (D1_Sig1(kk) / Li) End do

Re: [R] repeated - R package - Now rmutil Compilation Errors

In my quest to install the repeated package, I have also to install the rmutil package. (BTW, I'm running Fedora Core 4 and R 2.1.1.) But I now get several compilation errors for rmutil and I'm afraid to try to fix them: [EMAIL PROTECTED] R]# R CMD INSTALL rmutil * Installing *source* package

Re: [R] repeated - R package - Compilation Error

On Thu, 11 Aug 2005 [EMAIL PROTECTED] wrote: When I go into binnest.f, here is what I see: C Calculate First Derivative Do kk = 1, t1 Gradient(kk) = Gradient(kk) + (D1_beta(kk) / Li) End do Do kk = 1, t2 Gradient(t1+kk) =

Re: [R] sub set selection

check stepAIC in MASS package. On 8/11/05, Clark Allan [EMAIL PROTECTED] wrote: hi all is there a package that undertakes subset selection but BASED ON AIC or any other information criteria. i've seen the subselect and the leaps package but i have not played around with them yet.

Re: [R] p-values

Spencer, Here is an example from rayner and best 2001 and the script sent by Felipe. This can be done as follows using the function durbin.grupos() in the attached file ###Ice cream example from Rayner and Best 2001 . Chapter 7 judge - rep(c(1:7),rep(3,7)) variety -

Re: [R] repeated - R package - Compilation Error

You know what spams are? It is highly recommended to try such basic things and to get the depencies right _before_ posting to mailing lists. Regards Michael [EMAIL PROTECTED] wrote: When I go into binnest.f, here is what I see: C Calculate First Derivative Do kk = 1, t1

Re: [R] repeated - R package - Now rmutil Compilation Errors

On Thu, 11 Aug 2005 [EMAIL PROTECTED] wrote: In my quest to install the repeated package, I have also to install the rmutil package. (BTW, I'm running Fedora Core 4 and R 2.1.1.) But I now get several compilation errors for rmutil and I'm afraid to try to fix them: [EMAIL PROTECTED] R]# R

Re: [R] repeated - R package - Now rmutil Compilation Errors

Please note that you are using gfortran rather than g77 - they are not the same, and very possibly differ in how they react to code. The 1e301 constants do overflow the maximum real value (of the order of 3e+38), but 1d301 is less than the maximum double value of the order of 1.7e+308. So

[R] Kouros Owzar is out

I will be out of the office starting 08/10/2005 and will not return until 08/15/2005. __ R-help@stat.math.ethz.ch 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] sub set selection

On Thu, 11 Aug 2005, Clark Allan wrote: hi all is there a package that undertakes subset selection but BASED ON AIC or any other information criteria. i've seen the subselect and the leaps package but i have not played around with them yet. The leaps package finds a best model of each

[R] How to insert a certain model in SVM regarding to fixed kernels

Dear David, Dear R Users , Suppose that we want to regress for example a certain autoregressive model using SVM. We have our data and also some fixed kernels in libSVM behinde e1071 in front. The question: Where can we insert our certain autoregressive model ? During creating data frame ? Or

Re: [R] tcltk programming guide

thank you, I knew that link, but I need something more document- oriented, more specific, i.e. if I want to know how to use tkadd, where should I look? or tkinsert and so on... thanx, simone Il giorno 11/ago/05, alle ore 14:18, John Zhang ha scritto: Go to

Re: [R] tcltk programming guide

Dear Simone, tcltk functions correspond closely to Tcl and Tk commands, so documentation for the latter, available at http://wiki.tcl.tk/3109, is helpful. I also found Welsch's Practical Programming in Tcl and Tk useful. I expect that you've already seen Peter Dalgaard's two R News articles on

Re: [R] Creating new columns inside a loop

Ok, here's an english pseudo coded version of what I'd like to do... 10 columns of (somedata) names(somedata): C1 C2 C3 C4 C5 Loop through each column FOR ColName = C1 through C5 Compute a new column, named ColNameA = some result for each row Compute another new Column named ColNameB = some

[R] question on creating a new logical variable

I have two variables. A - rep(c(1:9), 2) B - rep(c(2:10),2) -- Renuka Sane http://www.nyx.net/~rsane [[alternative HTML version deleted]] __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read

[R] question on creating a new logical variable

I think there was an incomplete mail that was accidently sent by me. I apologize for the inconvenience. Here is the full text. I have two variables. A - rep(c(1:9), 2) B - rep(c(2:10),2) I want to know the the value for A==1 and B==1 If I do A== 1 | B==1, I get FALSE TRUE 16 2 This is

[R] strips according to groups in xYplot or xyplot

Dear all, using the follwing code, i have been trying group the individual plots from trellis graphics according to genus. adding a groups=genus colours the lines on the individual plots accordingly. However, i need the strips to be coloured according to genus with the lines in each plot

Re: [R] Please don't blame your tools (was Error in autoloader)

The problem is **NOT** in MASS. confint is in stats, not MASS, and has been (in base then stats) since R 1.7.0. You could (and should) have checked that from the (O)NEWS files of R and MASS. Nor is there an error in autoloader, as your title stated. Some other piece of software on your

Re: [R] question on creating a new logical variable

table(ifelse(A==2B==3, TRUE, FALSE)) On 8/11/05, Renuka Sane [EMAIL PROTECTED] wrote: I think there was an incomplete mail that was accidently sent by me. I apologize for the inconvenience. Here is the full text. I have two variables. A - rep(c(1:9), 2) B - rep(c(2:10),2) I want to

Re: [R] question on creating a new logical variable

On Thu, 11 Aug 2005, Renuka Sane wrote: I have two variables. A - rep(c(1:9), 2) B - rep(c(2:10),2) I want to know the the value for A==1 and B==1 [SNIP] To solve the problem I therefore create a new variable C - c(A, B) and then do table(C==1) which gives me FALSE TRUE 34 2 Is

Re: [R] Creating new columns inside a loop

On Thu, 11 Aug 2005, David L. Van Brunt, Ph.D. wrote: Ok, here's an english pseudo coded version of what I'd like to do... 10 columns of (somedata) names(somedata): C1 C2 C3 C4 C5 Loop through each column FOR ColName = C1 through C5 Compute a new column, named ColNameA = some result for each

[R] Converting strings with internal delimiters into lists

Hi UserRs, I know that there has to be an easy way to do this in R (probably easy enough that once someone clues me in I'll smack myself on the forehead for not figuring it out myself), but my searches on my own have not yielded any hints. I have many fields in my dataset that participants

Re: [R] background colors in image()

Romain, Thanks for the code. It worked perfectly! --- Romain Francois [EMAIL PROTECTED] wrote: Le 11.08.2005 00:59, array chip a écrit : Thanks for the suggestion! It works in a way that the entire graph window is in the background color, is there a way to only have the plotting area

[R] easier way to print heatmap on multiple pages?

Hi All, I've worked on some code to take a heatmap with 1000 row entries, and split this up into 20 pages, each with 50 rows from the original heatmap. I want to preserve the row order such that all 20 pages, if put together, would comprise the original heatmap. Here's what I've done:

Re: [R] Converting strings with internal delimiters into lists

??? I guess I don't get it. Note that testtable- as.data.frame(cbind(c(Joe,Mary,Jane,Mary),c(Fred,Joe,Pete,Joe,Mary,Fr ed))) is probably not what you want since cbind expects vectors of equal length. By default, shorter vectors are recycled to the length of the longest one, which I doubt is

Re: [R] Converting strings with internal delimiters into lists

On Thu, 11 Aug 2005, shelby berkowitz wrote: Hi UserRs, I know that there has to be an easy way to do this in R (probably easy enough that once someone clues me in I'll smack myself on the forehead for not figuring it out myself), but my searches on my own have not yielded any hints. I

Re: [R] Converting strings with internal delimiters into lists

yes, that does it. thanks! --- Thomas Lumley [EMAIL PROTECTED] wrote: On Thu, 11 Aug 2005, shelby berkowitz wrote: Hi UserRs, I know that there has to be an easy way to do this in R (probably easy enough that once someone clues me in I'll smack myself on the forehead for not

[R] A coding question involving variable assignments

Dear List, I have the following code that does what I want: x - replicate(5,replicate(10,sum(rnorm(rpois(1,10) How might one change it such that the maximum value generated by rnorm(rpois(1,10)) can be retrieved for later use? __

[R] getting xpinch and ypinch

Hello! I would like to know if it is possible (and how) to get xpinch and ypinch values that R gets (from the system) for drawing plots. If I understand correctly; with this values I could save the same picture as I see on computer 1 on computer 2 as a let say wmf file. Thank you in advance,

[R] signal handling

Is ther a signal handling model in R? similar to Perl's %SIG hash. I want to do fast clean up in my R code before exit when a kill signal is issued. __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the

[R] clustering or homegenity approaches?

Hi, there: I have a question on the following dataset rbind(t2[which(t40.3),][1:3,], t2[1:3,]) # don't worry about what this line means [,1] [,2] [,3] [,4] [,5] [1,] 34.216166 96.928587 330.125990 330.183222 330.201215 [2,] 2.819183 8.134491 8.275841

[R] error tracing

Hi, I am running some Cox modeling on large number of variables (thousands) using apply. For some reasons, some of the variables have problems in Cox regression, so the the run was stopped automatically. Is there a way to keep the running for all the variables to complete and let us, when done, to

[R] scatter plot

I'd like to do a simple scatter plot but instead of using the variable values on the X axis I would like to plot the percentiles. I searched in the manual for percentiles but did not find what I was looking for. I've been using SAS for several years but I new to R. -- Dean Sonneborn

Re: [R] scatter plot

Here is a hack par(mfrow=c(1,2)) x - rnorm(1000) y - x + rnorm(1000) xp - (rank(x)-1)/(length(x)-1) plot( x, y ) plot( xp, y ) But do notice that by using percentiles, it spreads the 'x' values evenly. This may be important for points at the extremes. Regards, Adai On Thu,

[R] scatter plot

I'm doing this scatter plot and it's almost what I want. What I would like to change is to have the X axis (the lipid_adj_pcb_cent variable) be plotted in percentiles not the variable values. Here's the code I currently using: xyplot(zawgt_sk_cent ~ lipid_adj_pcb_cent, panel=function(x,y){

Re: [R] signal handling

Omar Lakkis wrote: Is ther a signal handling model in R? similar to Perl's %SIG hash. I want to do fast clean up in my R code before exit when a kill signal is issued. I don't know Perl's %SIG hash. However, there are several things you can do: the on.exit() function lets a function do

Re: [R] A coding question involving variable assignments

On 12/08/05 05:55, xpRt.wannabe wrote,: Dear List, I have the following code that does what I want: x - replicate(5,replicate(10,sum(rnorm(rpois(1,10) How might one change it such that the maximum value generated by rnorm(rpois(1,10)) can be retrieved for later use? set.seed(99) x

Re: [R] A coding question involving variable assignments

Jim and List, Thank you for the prompt reply. Perhaps I should have been more specific in the way I phrased the question. The code you gave would return the max value just one time. I was interested in getting as many max values generated by rnorm(rpois(1,10)) as specified by:

Re: [R] A coding question involving variable assignments

Ted and List, In your code that produced 'mx', you dropped sum() from my original code though. As a result, the 10 x 5 max's are of the same value. Unfortunately, that's not what I need. On 8/11/05, ecatchpole [EMAIL PROTECTED] wrote: On 12/08/05 05:55, xpRt.wannabe wrote,: Dear List,

Re: [R] A coding question involving variable assignments

Sorry, I don't follow. What's wrong with this? Ted. set.seed(99) x - replicate(5,replicate(10,sum(rnorm(rpois(1,10) set.seed(99) mx - replicate(5,replicate(10,max(rnorm(rpois(1,10) x [,1][,2][,3] [,4] [,5] [1,] -2.0071674 -7.0883335

[R] access assigned objects

Dear list, I am trying to write a function that will make a matrix of each row of a data frame (4 columns), calculate a fisher.test on each resulting matrix and assign a vector of the p-values. I have gotten through making the matrices, but cannot calculate the fisher.test. fm-function(x)

Re: [R] A coding question involving variable assignments

Ted and List, What I need is I need to know what max of rnorm(rpois(1,10)) is before R does sum(), replicate(10, ...) and replicate(5, ...). The fact that you have set.seed(99) twice, does that mean, say, entry [1,1] 0.4896243 in 'mx' is one of the z number of values generated by

Re: [R] A coding question involving variable assignments

On 12/08/05 13:27, xpRt.wannabe wrote,: Ted and List, What I need is I need to know what max of rnorm(rpois(1,10)) is before R does sum(), replicate(10, ...) and replicate(5, ...). The fact that you have set.seed(99) twice, does that mean, say, entry [1,1] 0.4896243 in 'mx' is one of the

[R] Help converting a function from S-Plus to R: family$weight

Hi all I am converting an S-Plus function into R. The S-Plus code uses some of the glm families, and family objects. The family objects in S-Plus and R have many different features, for example: In R: names(Gamma()) [1] family link linkfunlinkinvvariance [6]

Re: [R] A coding question involving variable assignments

Thank you! That was helpful. Another thing I learned was that I would need to do set.seed(99) not once but twice in this context. On 8/11/05, ecatchpole [EMAIL PROTECTED] wrote: On 12/08/05 13:27, xpRt.wannabe wrote,: Ted and List, What I need is I need to know what max of

[R] chisq warning

Hi I am running chisq as below and getting a warning. Can anyone tell me the significance or the warning? chisq.test(c(10 ,4 ,2 ,6 ,5 ,3 ,4 ,4 ,6 ,3 ,2 ,2 ,2 ,4 ,7 ,10 ,0 ,6 ,19 ,3 ,2 ,7 ,2 ,2 ,2 ,1 ,32 ,2 ,3 ,10 ,1 ,3 ,9 ,4 ,10 ,2 ,2 ,4 ,5 ,7 ,6 ,3 ,7 ,4 ,3 ,3 ,7 ,1 ,4 ,2 ,2 ,3 ,3 ,5 ,5 ,4