[R] fitting an interaction term

Hello! IŽm fitting a model with glm(family binomial). The best model counts 9 Variables and includes an interaction term that was generated by the product of to continuous variables (a*b). All variables are correlated under a value of 0.7 (Spearman rank order) While the estimates of both main

Re: [R] Substring and strsplit

you can also use substring(), e.g., substring(x3, 1:nchar(x3), 1:nchar(x3)) Best, Dimitris Dimitris Rizopoulos Ph.D. Student Biostatistical Centre School of Public Health Catholic University of Leuven Address: Kapucijnenvoer 35, Leuven, Belgium Tel: +32/(0)16/336899 Fax: +32/(0)16/337015

[R] Help on apply() function

Respected Sir/Madam, I have a problem with apply function. I have to two matrices of dimension of one column but n rows. I have to check whether one matrix is greater than other by going thru each row (ie) using if condition to check one matrix with another matrix. I like to use apply()

[R] Installation of SrcStatConnectorSrv on Windows

I am trying to install SrcStatConnectorSrv(2) and rcom from local zip files. I have successfully downloaded both files (from CRAN for rcom CRAN Other for SrcStatConnectorSrv) and installed rcom. However, I get the following error message when I try to install SrcStatConnectorSrv (the version

Re: [R] Substring and strsplit

If you are using 'only' English then str - dog strsplit(str,NULL)[[1]] works perfectly and it is fast. But if you also dealing with Unicode character have a look at http://wiki.r-project.org/rwiki/doku.php?id=tips:data- strings:decomposestring Cheers, Hans you can also use substring(),

[R] MCMClogit

Hi, I am using MCMCpack and the MCMClogit function to create logistic regression models in a medical (adverse event) study. My question is, is there a way where I can directly create the estimated probabilities of the adverse outcome, moreover the confidence interval for the estimated

[R] How to put title Vertically

Dear all R users, Suppose, Dear all R users, Suppose, pauto.cor = pacf(lh, plot=F) max.lag = max(pauto.cor$lag) min.lag = min(pauto.cor$lag) centre = (max.lag - min.lag)/2 pauto.cor = pauto.cor$acf pauto.cor = pauto.cor[-1] par(mar=c(3,0,1,1)) barplot(pauto.cor,

[R] Barplot

Dear all, I have a dataset. I want to make barplot from this data. Zero1 - V1 V2 V3 V4 V5 V6 V7 V8 V9 1 1 0 0 0 1 0 0 0 Positive 2 0 0 1 0 1 0 1 1 Negative 3 0 0 1 0 0 0 1 1 Positive 4 0 1 0 1 1 1 0 1 Negative 5 0 0 1 0 1 1 0 0 Positive 6 0

Re: [R] Installation of SrcStatConnectorSrv on Windows

That is not an R package. It looks like the sources of the server described on http://cran.r-project.org/contrib/extra/dcom/RSrv200.html which tells you about the approriate mailing list to ask about it. You probably want to use RSrv200.exe to install the server (and other tools and

Re: [R] Barplot

barplot(t(sapply(split(z1[,1:8], z1$V9),colSums)), beside=T) --- Jacques VESLOT CNRS UMR 8090 I.B.L (2ème étage) 1 rue du Professeur Calmette B.P. 245 59019 Lille Cedex Tel : 33 (0)3.20.87.10.44 Fax : 33 (0)3.20.87.10.31

[R] MCMClogit

Hi, I am using MCMCpack and the MCMClogit function to create logistic regression models in a medical (adverse event) study. My question is, is there a way where I can directly create the estimated probabilities of the adverse outcome, and the confidence interval for the estimated probabilities?

[R] density() with from, to or cut and comparrison of density()

Hi the function density() does normally integrate to one - I've checked it and it works and I also read the previous threads. But I realised that it does not integrate to one if I use from, to or cut. My scenario: I simulated densities of a plants originating from an sseed source at distance

Re: [R] Help on apply() function

On 8/30/06, Mark Lyman [EMAIL PROTECTED] wrote: I have a problem with apply function. I have to two matrices of dimension of one column but n rows. I have to check whether one matrix is greater than other by going thru each row (ie) using if condition to check one matrix with another

Re: [R] Barplot

Muhammad Subianto wrote: ... I have tried to learn ?xtabs ?table and ?ftable but I can't figure out. I need a barplot for all variables and the result maybe like | | | | | | | | | || | | |pos|neg| |pos|neg||pos|neg| |

Re: [R] Firefox extension fo R Site Search

May be it's not a bug, but I tried to search for the package rpanel and I was not find. At the r-project's site that package is available. How to explain it? Rob __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help

Re: [R] Firefox extension fo R Site Search

Robert Mcfadden a écrit : May be it's not a bug, but I tried to search for the package rpanel and I was not find. At the r-project's site that package is available. How to explain it? Rob Hi Rob, If it's not there : http://finzi.psych.upenn.edu/R/library/ it's not on the extension

Re: [R] How to put title Vertically

stat stat wrote: Dear all R users, Suppose, Dear all R users, Suppose, pauto.cor = pacf(lh, plot=F) max.lag = max(pauto.cor$lag) min.lag = min(pauto.cor$lag) centre = (max.lag - min.lag)/2 pauto.cor = pauto.cor$acf pauto.cor = pauto.cor[-1] par(mar=c(3,0,1,1))

Re: [R] density() with from, to or cut and comparrison of density()

Rainer M Krug wrote: Hi the function density() does normally integrate to one - I've checked it and it works and I also read the previous threads. But I realised that it does not integrate to one if I use from, to or cut. My scenario: I simulated densities of a plants originating from

Re: [R] Firefox extension fo R Site Search

Romain Francois wrote: Robert Mcfadden a écrit : May be it's not a bug, but I tried to search for the package rpanel and I was not find. At the r-project's site that package is available. How to explain it? Rob Hi Rob, If it's not there : http://finzi.psych.upenn.edu/R/library/

Re: [R] Firefox extension fo R Site Search

On 08/30/06 11:55, Romain Francois wrote: Robert Mcfadden a ?crit : May be it's not a bug, but I tried to search for the package rpanel and I was not find. At the r-project's site that package is available. How to explain it? Rob Hi Rob, If it's not there :

Re: [R] Barplot - thanks

Dear all, Many Thanks to Jacques VESLOT and Jim Lemon for their helps. Best, Muhammad Subianto #Jacques VESLOT barplot(t(sapply(split(z1[,1:8], z1$V9),colSums)), beside=T) #Jim Lemon barplot(sapply(z1[1:8],by,z1[9],sum),beside=TRUE) On this day 30/08/2006 11:43, Muhammad Subianto wrote:

[R] [R-pkgs] Version 1.2-0 of the Rcmdr package

I've submitted a new, and substantially enhanced, version (1.2-0) of the Rcmdr package to CRAN. Some highlights (from the CHANGES) file: o Added ability to import from Excel, Access or dBase files (contributed by Samir Messad, Renaud Lancelot and Matthieu Lesnoff). o Added ability to

[R] Datetime

Hi all, I'm getting confused by date handling. I wish to read a date from a file wich is a number of seconds since 1970 (POSIXct). Then i wish to convert this date to a human readable form (POSIXlt) By example : ctDate-1132963200 #Wed Aug 30 14:24:37 2006 is(ctDate) [1] numeric

Re: [R] Handling realisations in geoRglm

Chris Barker cmbarker at ucdavis.edu writes: However, after reading the packages' documentation and searching the mailing lists and other sources, it seems that the functions in geoRglm do not acknowledge the existence of multiple realisations per site. I see that the likfit() function

Re: [R] MCMClogit

Andras, At this point you need to write your own function to take the posterior density sample (stored in a coda mcmc object) and covariates of interest to get a Monte Carlo estimate of these probabilities. Best, ADM On Aug 30, 2006, at 5:10 AM, Andras Treszl wrote: Hi, I am using

[R] Create a vector from another vector

Dear list Suppose I have the following vector: x - c(3,4,2,5,6) Obviously, this sums to 20. Now, I want to have a second vector, call it x2, that sums to x where 5 = x = 20, but there are constraints. 1) The new vector must be same length as x 2) No element of the new vector can be 0 3)

Re: [R] Create a vector from another vector

Dear Harold package partitions does almost this: library(partitions) x - 1+restrictedparts(15,5) x[,1:10] [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] 16 15 14 13 12 11 10 [2,]1234567 [3,]1111111 [4,]111

[R] Optimization

Dear R-list, I'm trying to estimate the relative importance of 6 environmental variables in determining clam yield. To estimate clam yield a previous work used the function Yield = (A^a*B^b*C^c...)^1/(a+b+c+...) where A,B,C... are the values of the environmental variables and the weights a,b,c...

[R] converting decimal - hexadecimal

Hi, do you know, a method to convert an decimal value (integer) to the corresponding hexadecimal value ? thinks for help. Romain -- Lorrillière Romain UMR 8079 Laboratoire Ecologie, Systématique et Evolution Bât. 362 Université Paris-Sud 91405 Orsay cedex France tel : 01 69 15

Re: [R] Create a vector from another vector

On 8/30/2006 9:49 AM, Doran, Harold wrote: Dear list Suppose I have the following vector: x - c(3,4,2,5,6) Obviously, this sums to 20. Now, I want to have a second vector, call it x2, that sums to x where 5 = x = 20, but there are constraints. 1) The new vector must be same length

Re: [R] Create a vector from another vector

I think I've got it now. If I understand your question, try: x - do.call(cbind,lapply(5:7,function(i){restrictedparts(i, 5,include.zero=FALSE)})) acceptable - function(x){all(x=c(3,4,5,2,6))} x[,apply(x,2,acceptable)] [,1] [,2] [,3] [,4] [1,]1232 [2,]11

Re: [R] Barplot

Try this. First we reduce the data to a frequency matrix and then plot it using classic and then lattice graphics: zm - as.matrix(rowsum(z1[-9], z1[,9])) barplot(zm, beside = TRUE, col = grey.colors(2)) legend(topleft, legend = levels(z1[,9]), fill = grey.colors(2)) library(lattice)

Re: [R] converting decimal - hexadecimal

Try: sprintf(%x, 109) On 9/30/06, Romain Lorrilliere [EMAIL PROTECTED] wrote: Hi, do you know, a method to convert an decimal value (integer) to the corresponding hexadecimal value ? thinks for help. Romain -- Lorrillière Romain UMR 8079 Laboratoire Ecologie, Systématique et

Re: [R] converting decimal - hexadecimal

use sprintf sprintf(%x,123) [1] 7b On 9/30/06, Romain Lorrilliere [EMAIL PROTECTED] wrote: Hi, do you know, a method to convert an decimal value (integer) to the corresponding hexadecimal value ? thinks for help. Romain -- Lorrillière Romain UMR 8079 Laboratoire Ecologie,

Re: [R] converting decimal - hexadecimal

?sprintf ex : sprintf('%X',10) [1] A hih __ 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 and provide commented, minimal,

[R] function which gives the hessian matrix of the log-likelihood of a nonlinear mixed model?

Hi, Can anyone tell me which function in R gives the hessian matrix of the log-likelihood of a nonlinear mixed model? fdHess is for scarlar function only. Thanks in advance! Hongmei [[alternative HTML version deleted]] __

Re: [R] converting decimal - hexadecimal

how about learning to use help.search? (In may you already asked a similiar question...) help.search(hexadecimal) which would lead you to format.hexmode ?format.hexmode Romain Lorrilliere schrieb: Hi, do you know, a method to convert an decimal value (integer) to the corresponding

[R] working with summarized data

The data sets I am working with all have a weight variable--e.g., each row doesn't mean 1 observation. With that in mind, nearly all of the graphs and summary statistics are incorrect for my data, because they don't take into account the weight. For example median is incorrect, as the

Re: [R] Create a vector from another vector

Hi Duncan Here is a bit more detail, this is a bit tough to explain, sorry for not being clear. Ordering is not important because the vector I am creating is used as a sufficient statistic in an optimization routine to get some MLEs. So, any combination of the vector that sums to X is OK. But,

Re: [R] converting decimal - hexadecimal

An other way would be: a - 123 class(a) - hexmode a [1] 7b On 30 Aug 2006, at 16:26, mel wrote: ?sprintf ex : sprintf('%X',10) [1] A hih __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read

Re: [R] converting decimal - hexadecimal

Romain Lorrilliere [EMAIL PROTECTED] wrote in message news:[EMAIL PROTECTED] Hi, do you know, a method to convert an decimal value (integer) to the corresponding hexadecimal value ? Starting in R 2.1.0, sprintf can be used: x - c(0, 65535, 65536, 305419896, 2^31-1) y - sprintf(0x%X, x) y

Re: [R] working with summarized data

In each case, look around (help.search, RSiteSearch) to see if you can find a function that handles weights. For the case you mention, medians, it can be done via quantile regression: x - w - 1:5 library(quantreg) coef(rq(x ~ 1, weight = w)) On 8/30/06, Rick Bischoff

Re: [R] Create a vector from another vector

OK, With this extra detail, a third solution follows, which may be closer in spirit to your application. It may or may not be faster than the other two, depending on the exact parameters used: library(partitions) 1+blockparts(n=15,y=c(3,4,2,5,6)-1,include.fewer=T) (720 distinct solutions)

Re: [R] Create a vector from another vector

maybe something like this could be of help: max.score - c(3,4,3) # max score for each item all.pats - as.matrix(expand.grid(lapply(max.score, :, 1))) all.pats[rowSums(all.pats) == 5, ] Best, Dimitris Dimitris Rizopoulos Ph.D. Student Biostatistical Centre School of Public Health Catholic

Re: [R] working with summarized data

Unfortunately, it seems that most(all?) of R's graphics and summary statistic functions don't take a weight or frequency argument. (Fortunately the models do...) I have been been meaning to add this functionality to my graphics package ggplot (http://had.co.nz/ggplot), but unfortunately

Re: [R] working with summarized data

On Wed, 30 Aug 2006, Rick Bischoff wrote: Unfortunately, it seems that most(all?) of R's graphics and summary statistic functions don't take a weight or frequency argument. (Fortunately the models do...) I have been been meaning to add this functionality to my graphics package ggplot

[R] lmer applied to a wellknown (?) example

Dear all, During my pre-R era I tried (yes, tried) to understand mixed models by working through the 'rat example' in Sokal and Rohlfs Biometry (2000) 3ed p 288-292. The same example was later used by Crawley (2002) in his Statistical Computing p 363-373 and I have seen the same data being

[R] bootstrap for group and subgroup

Dear R-friends, I have a table data structured by group, subgroups, records and attributes. For each group and subgroup I have differente number of records (more than 200). I need bootstrap 100 records for each group/subgroup combinations and repeat it a big number of times. Could

Re: [R] Bootstraping for groups and subgroups and joing with other table

I have a table with following collumns: State, SamplePlot, Species and BodySize. I sampled bird species at 34 SamplePlots and 5 States (regions) monthly during two years. On each bird record I measured bodysize and identified the species. So I have many records of each species (about 150

Re: [R] Barplot

Dear all, To Gabor Grothendieck, (again) thanks you very much for your help. Now, I can play around with lattice package. Best, Muhammad Subianto #Gabor #reduce the data to a frequency matrix and #then plot it using classic and then lattice graphics: zm - as.matrix(rowsum(z1[-9], z1[,9]))

[R] Ranking and selection statistical procedure

Dear R helpers I would like to know if the Ranking and Selection statistical procedure has been implemented in R. I made a quick search in the R packages list but I could not find it. Thanks in advance Prasanna __ R-help@stat.math.ethz.ch mailing list

Re: [R] Producing R demos

One option is to use VNC along with vncrec to do the recording (see the website: http://www.sodan.org/~penny/vncrec/). I think there are some other recorders also available for vnc, so you might try a google search. -- Gregory (Greg) L. Snow Ph.D. Statistical Data Center Intermountain

Re: [R] density() with from, to or cut and comparrison of density()

You may want to look at the logspline package, it uses a different technique than density does, but it estimates densities and allows you to tell the routine that there is a minimum value and that the density does not extend beyond there. Hope this helps, -- Gregory (Greg) L. Snow Ph.D.

Re: [R] working with summarized data

There are functions to do weighted summary statistics in the Hmisc package (wtd.quantile, ...). For more complicated analyses (but not plots yet) the biglm package has a bigglm function that expects the data in chunks, you could write a function that expand parts of the dataset at a time. Hope

[R] Debugging with gdb

I tried to run gdb in linux with emacs But could not even run a simple example in the writing extensions tutorial. The execution history is as follows. Gdb worked fine for other debugging such as C++ codes. Thanks Han I started R at echo of emacs by typing (also tried other methods mentioned

Re: [R] lattice and several groups

hadley wickham a écrit : I would like to use the lattice library to show several groups on the same graph. Here's my example : ## the data f1 - factor(c(mod1,mod2,mod3),levels=c(mod1,mod2,mod3)) f1 - rep(f1,3) f2 - factor(rep(c(g1,g2,g3),each=3),levels=c(g1,g2,g3)) df -

[R] .Rprofile under Windoze.

I am (for my sins) having to do some work using R under Windoze. I wanted to set up a .Rprofile to control my set-up. The docs on .Rprofile say that it can/should be placed in ``the user's home directory''. ``An Introduction to R'' observes lucidly that this concept needs to be clarified under

Re: [R] lattice and several groups

Gabor Grothendieck a écrit : Try this: xyplot(val ~ x, data = df, type = p, col = as.numeric(df$f1), pch = as.numeric(df$f2)) key1 - list(border = TRUE, colums = 2, text = list(levels(df$f1)), points = list(pch = 1:nlevels(df$f1)) ) key2 - list(border = TRUE, colums = 2, text =

Re: [R] lattice and several groups

Gabor Grothendieck a écrit : Note that before entering this you need: library(lattice) library(grid) # to access the viewport function On 8/29/06, Gabor Grothendieck [EMAIL PROTECTED] wrote: Try this: xyplot(val ~ x, data = df, type = p, col = as.numeric(df$f1), pch =

Re: [R] .Rprofile under Windoze.

Under Windows mine is located here C:\Program Files\R\R-2.3.1\library\base\R The file name, however is not .Rprofile, but rather Rprofile Charles Annis, P.E. [EMAIL PROTECTED] phone: 561-352-9699 eFax: 614-455-3265 http://www.StatisticalEngineering.com -Original Message- From:

Re: [R] .Rprofile under Windoze.

On 8/30/06, Rolf Turner [EMAIL PROTECTED] wrote: I am (for my sins) having to do some work using R under Windoze. I wanted to set up a .Rprofile to control my set-up. The docs on .Rprofile say that it can/should be placed in ``the user's home directory''. ``An Introduction to R'' observes

Re: [R] lattice and several groups

In fact, my problem is to fit the data for every level of the f2 factor, showing the levels of the f1 factor and that for several surveys . Here's an example closer to my actual data : Then maybe you want: qplot(x, y, . ~ surveys, data=df, type=line, colour=f1, id=f2, size=f2) (which doesn't

Re: [R] .Rprofile under Windoze.

And... If you have/use shortcuts to R, you may also save an .Rprofile to whatever directory you name in the Start in: field of the shortcut. This allows one to have many profiles. jab -- John Bollinger, CFA, CMT www.BollingerBands.com If you advance far enough, you arrive at the

[R] how to read just a column

Hi, how can I read, using for example read.table() or scan(), just one column from a text file that has more columns without any header? Thanks, bye. __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read

Re: [R] .Rprofile under Windoze.

-Original Message- From: [EMAIL PROTECTED] [mailto:r-help- [EMAIL PROTECTED] On Behalf Of Rolf Turner Sent: Wednesday, August 30, 2006 2:50 PM To: r-help@stat.math.ethz.ch Subject: [R] .Rprofile under Windoze. I am (for my sins) having to do some work using R under Windoze. I

Re: [R] .Rprofile under Windoze.

On 8/30/2006 5:49 PM, Rolf Turner wrote: I am (for my sins) having to do some work using R under Windoze. I wanted to set up a .Rprofile to control my set-up. The docs on .Rprofile say that it can/should be placed in ``the user's home directory''. ``An Introduction to R'' observes lucidly

Re: [R] how to read just a column

you can read them all in and delete the ones you don't want. Or check out 'what' on 'scan' or colClasses on 'read.table' On 8/30/06, Carlo Trimarchi [EMAIL PROTECTED] wrote: Hi, how can I read, using for example read.table() or scan(), just one column from a text file that has more columns

Re: [R] lattice and several groups

To handle conditioning on survey we provide a panel function that subsets col and pch: # define test data - df # note that your val2 and val3 lines had a syntax # so we have commented them out and # replaced them as shown. n - 18 x1 - seq(1,n) val1 - -2*x1+50 # val2 - (-2*(x1-8)2)+100 val2 -

Re: [R] .Rprofile under Windoze.

Rolf The last invocation indicated that the name of the file was *really* ``.Rprofile.txt'' --- although I'd tried to save it as (simply) ``.Rprofile''. Is that the problem? If so, how can I persuade Windoze NOT to stick that damned .txt tag on the end? The easiest way is to use a smarter

[R] NaN when using dffits, stemming from lm.influence call

Hi all I'm getting a NaN returned on using dffits, as explained below. To me, there seems no obvious (or non-obvious reason for that matter) reason why a NaN appears. Before I start digging further, can anyone see why dffits might be failing? Is there a problem with the data? Consider:

Re: [R] lattice and several groups

Or maybe this is what you are looking for where pnl below was created by modifying source to the panel.plot.default in the zoo package (there might be a simpler way): pnl - function (x, y, subscripts, groups, col, pch, type, ...) { for (g in levels(groups)) { idx - g ==

Re: [R] lattice and several groups

In thinking about this a bit more we can use panel.superpose/panel.groups to shorten it: # define data -- df # note that your val2 and val3 lines had a syntax # so we have commented them out and # replaced them as shown. n - 18 x1 - seq(1,n) val1 - -2*x1+50 # val2 - (-2*(x1-8)2)+100 val2 -

[R] [R-pkgs] New package 'random' for non-deterministic random number

Dear useRs, A few days ago, the initial version 0.1.0 of a new package 'random' was uploaded to CRAN. The random packages provides convenient access to the non-deterministic random numbers provided by the random.org site created by Mads Haahr (http://www.random.org). While certain hardware

[R] Combine 'overlapping' dataframes, respecting row names

Hi, I've examined the archives and found quite a few questions on concatenating dataframes, but none that really addressed my issue, I'm afraid. I've also examined the cbind and rbind documentation but nonetheless here I am writing to r-help ;) This is what I have (the row names are dates

Re: [R] Combine 'overlapping' dataframes, respecting row names

If you are converting them to 'its' anyways then after the conversion to 'its' use the 'its' union command. On 8/31/06, James Howison [EMAIL PROTECTED] wrote: Hi, I've examined the archives and found quite a few questions on concatenating dataframes, but none that really addressed my issue,

Re: [R] Combine 'overlapping' dataframes, respecting row names

'merge' is the key here. You say you want to merge, but it seems did not try merge() (res - merge(cvsFrame, downloadsFrame, by=row.names, all=TRUE)) Row.names cvsactions downloads 1 2002-11-15 412 2 2002-12-15 9 8 3 2003-01-15 5NA 4

Re: [R] how to read just a column

Carlo Trimarchi a écrit : Hi, how can I read, using for example read.table() or scan(), just one column from a text file that has more columns without any header? Thanks, bye. afaik, you have to read all the table and then you select the column you want. eg read.table(blabla)[3] to get the