[R] Re : BOOTSTRAP_CROSS VALIDATION

2009-03-16 Thread justin bem
See the boot package  Justin BEM BP 1917 Yaoundé Tél (237) 99597295 (237) 22040246 De : per243 jose.perezsua...@csiro.au À : r-help@r-project.org Envoyé le : Lundi, 16 Mars 2009, 0h12mn 03s Objet : [R] BOOTSTRAP_CROSS VALIDATION I need a script that works

Re: [R] Help_Unable to run glmer

2009-03-16 Thread Dieter Menne
David Winsemius dwinsemius at comcast.net writes: Maybe there was a run on the British bit and the bit banks are a wee short? Runs are by definition in the order of TeraBits. Bits are peanuts and freely available. D __ R-help@r-project.org mailing

Re: [R] [OT] two question about color space.

2009-03-16 Thread Thomas Lumley
The XYZ to RGB conversion is already in convertColor(), and also in the colorspace package. -thomas On Sun, 15 Mar 2009, baptiste auguie wrote: I've put together a rough R port of that C code [*] in a package on r-forge:

Re: [R] primitives again

2009-03-16 Thread Wacek Kusnierczyk
Berwin A Turlach wrote: On Sun, 15 Mar 2009 14:23:40 +0100 Wacek Kusnierczyk waclaw.marcin.kusnierc...@idi.ntnu.no wrote: Edna Bell wrote: How do I find the functions which are primitives, please? you can scan the whole search path for functions that are primitives:

Re: [R] Package download problem {was Help_Unable to run glmer}

2009-03-16 Thread Martin Maechler
DW == David Winsemius dwinsem...@comcast.net on Sun, 15 Mar 2009 07:32:09 -0400 writes: DW It's a 404 message. Try a repository that does have DW it. The CMU repository has a copy but the one you picked DW does not seem to have one at the moment. Maybe there was DW a run on

Re: [R] How to combine xtable and minipage with Sweave ?

2009-03-16 Thread Ptit_Bleu
Dear Duncan, Thanks for the tips, but I haven't found a way to modify your solution so that I can include a table and not a graph (includegraphics works for graphs, I need an 'includetable', I think). I found scripts like yours explaining how to include graphs into table but this is not my case:

Re: [R] How to combine xtable and minipage with Sweave ?

2009-03-16 Thread Ptit_Bleu
Dear Duncan, Thanks for the tips, but I haven't found a way to modify your solution so that I can include a table and not a graph (includegraphics works for graphs, I need an 'includetable', I think). I found scripts like yours explaining how to include graphs into table but this is not my case:

Re: [R] [OT] two question about color space.

2009-03-16 Thread baptiste auguie
Thanks for pointing this out, i did check the available packages before dwelling on this but i was focussed on the first part of the problem (spectrum to xyz). The available convertColor() function is especially welcome that it was the most inefficient part of the code I produced yesterday.

Re: [R] [OT] two question about color space.

2009-03-16 Thread Jinsong Zhao
Thanks for pointing this out, i did check the available packages before dwelling on this but i was focussed on the first part of the problem (spectrum to xyz). The available convertColor() function is especially welcome that it was the most inefficient part of the code I produced yesterday.

[R] Creating variables with different levels depending on conditions

2009-03-16 Thread Miriam Gade
Dear R-help-list, I started using R coming from the former use of SPSS, mainly writing my own syntax. I do classical experimental cognitive psychology and am interested in the cognitive control of task sequences (i.e. classification task). Especially I am interested in the performance of the

Re: [R] How to combine xtable and minipage with Sweave ?

2009-03-16 Thread Duncan Murdoch
Ptit_Bleu wrote: Dear Duncan, Thanks for the tips, but I haven't found a way to modify your solution so that I can include a table and not a graph (includegraphics works for graphs, I need an 'includetable', I think). Sorry, I misread your request. (I saw the letters t, a, b, l, e and

Re: [R] Creating variables with different levels depending on conditions

2009-03-16 Thread Jim Lemon
Miriam Gade wrote: Dear R-help-list, I started using R coming from the former use of SPSS, mainly writing my own syntax. I do classical experimental cognitive psychology and am interested in the cognitive control of task sequences (i.e. classification task). Especially I am interested in the

[R] Updates to plotrix and clinsig

2009-03-16 Thread Jim Lemon
Hi everybody, An update to plotrix (2.5-3) is now on CRAN. Apart from the usual polishing, this includes an interesting new function (hierobarp) that was suggested on this very list a few weeks ago. This function will plot nested breakdowns of numeric values by one or more factors. It seems

[R] mixture models

2009-03-16 Thread Rita Gaio
Dear R users, I want to apply latent class models (mixture models) with covariates in the case of continuous and categorical variables. My questions are the follow: a) the only package in R that can handle mixture models with data of mixed type is moc (which unfortunately was left orphaned

Re: [R] gz netCDF files

2009-03-16 Thread jbaiao
I am experiencing the same problem with ncdf bzipped (bz2) files. After many tries I still can't find a way to open them... Does any one know how to do it? I would really appreciate your help. Fernando Lima Dept. Biol. Sciences Univ. South Carolina Magdalena Lucini wrote: Hello, I am

Re: [R] hazard function in a Cox model

2009-03-16 Thread Terry Therneau
-- begin included message --- I am hoping for some advice regarding obtaining the values for the hazard function in a cox regression that I have undertaken. I have a model in the following form, analysed with the package survival (v. 2.34-1) and a log-log plot obtained using Design (v. 2.1-2).

Re: [R] mixture models

2009-03-16 Thread Ingmar Visser
Dear Rita, The depmixS4 package handles such models. It is still in beta testing so I would appreciate it if you share your experiences with it. Also let me know if you have further questions. Best, Ingmar On 16 Mar 2009, at 10:19, Rita Gaio wrote: Dear R users, I want to apply latent

[R] XYplot simple question

2009-03-16 Thread AllenL
Hello R friends, Simple question today: I am desiring to do an xyplot with the below code, which graphs time series across different experimental Plots- xyplot(AbvBioAnnProd~Year|Plot)### Plots each monoculture biomass vs time xyplot(AbvBioAnnProd~Year|Plot,panel=function(x,y){

Re: [R] XYplot simple question

2009-03-16 Thread Sundar Dorai-Raj
Convert Plot to factor: xyplot(AbvBioAnnProd ~ Year | Plot, type = c(b, r), pch = 16) Also note that using the type argument with multiple values prevents the necessity of a custom panel function. HTH, --sundar On Mon, Mar 16, 2009 at 5:54 AM, AllenL allen.laroc...@gmail.com wrote: Hello R

Re: [R] XYplot simple question

2009-03-16 Thread David Winsemius
You should look again at the first example on the xyplot help page, run it, and pay attention to how the strip function works. All should be revealed. I suspect you want something like: xyplot(AbvBioAnnProd~Year|Plot,panel=function(x,y){ panel.xyplot(x,y,type=b,pch=16, strip =

[R] changing the description of the x-axi

2009-03-16 Thread johnhj
Hii, Is it possible to change the description ot the x-axis ? I need some barplots with high x-axis values. After plotting the graph, the x-axis values are shown in this form: 0e+00 1e+06 2e+06 3e+06 4e+06 5e+06 6e+06 Can I change this description to numerary values, like 5354352 ?

Re: [R] XYplot simple question

2009-03-16 Thread Sundar Dorai-Raj
Sorry, I should have xyplot(AbvBioAnnProd ~ Year | factor(Plot), type = c(b, r), pch = 16) On Mon, Mar 16, 2009 at 6:17 AM, Sundar Dorai-Raj sdorai...@gmail.com wrote: Convert Plot to factor: xyplot(AbvBioAnnProd ~ Year | Plot, type = c(b, r), pch = 16) Also note that using the type

[R] JGR and system()

2009-03-16 Thread Liviu Andronic
Dear all, JGR seems unable to report the results of system() commands. In a command-line R session: system(pwd) /home/liviu system(acpitool) AC adapter : on-line Thermal zone 1 : ok, 67 C system(asdf) sh: asdf: command not found In JGR: system(pwd) system(ls -F) system(acpitool)

Re: [R] Package download problem {was Help_Unable to run glmer}

2009-03-16 Thread David Winsemius
On Mar 16, 2009, at 4:26 AM, Martin Maechler wrote: DW == David Winsemius dwinsem...@comcast.net on Sun, 15 Mar 2009 07:32:09 -0400 writes: DW It's a 404 message. Try a repository that does have DW it. The CMU repository has a copy but the one you picked DW does not seem to have

Re: [R] Fw: Fitting GUMBEL Distribution - CDF function and P P Plot

2009-03-16 Thread J. R. M. Hosking
Congratulations, you have discovered a bug in the implementation of function cdfgum. A corrected version of package lmom will be submitted today and should be available on CRAN within a day or two. Meanwhile, as a workaround define cdfgum in your own workspace, thus: cdfgum -

[R] Uniroot and Newton-Raphson Anomaly

2009-03-16 Thread Doran, Harold
I have the following function for which I need to find the root of a: f - function(R,a,c,q) sum((1 - (1-R)^a)^(1/a)) - c * q To give context for the problem, this is a psychometric issue where R is a vector denoting the percentage of students scoring correct on test item i in class j, c is the

[R] FW: Select a random subset of rows out of matrix

2009-03-16 Thread Haenlein.Michael
Dear all, I have a large dataset (N=100,000 with 89 variables per subject). This dataset is stored in a 100.000 x 89 matrix where each row describes one individual and each column one variable. What is the easiest way of selecting a subset of let's say 1.000 individuals out of that whole

[R] svDialogs

2009-03-16 Thread Antje
Hi there, I was using Open/Save-dialogs from the package svDialogs (SciViews). But now the package has dissapeared? How do I have to set up my R-installation to further use these dialogs??? (beside copying my old packages to the new installation). Ciao, Antje

Re: [R] FW: Select a random subset of rows out of matrix

2009-03-16 Thread David Winsemius
On Mar 16, 2009, at 10:34 AM, Haenlein.Michael wrote: Dear all, I have a large dataset (N=100,000 with 89 variables per subject). This dataset is stored in a 100.000 x 89 matrix where each row describes one individual and each column one variable. What is the easiest way of selecting a

Re: [R] FW: Select a random subset of rows out of matrix

2009-03-16 Thread Paul Hiemstra
Hi, You could try something like this: nrow = 1000 ncol = 10 nosample = 50 m = matrix(runif(1), nrow = nrow, ncol = ncol) # Randomly select 50 row numbers and extract the data m_subset = m[sample(1:nrow, size = nosample),] cheers, Paul Haenlein.Michael wrote: Dear all, I have a large

[R] How to read large file

2009-03-16 Thread Dongyan Song
Dear all, I want to read a 1.2 G file with 136,736,218 numbers with scan function, and then want make a histogram with those numbers, but it showed Error: cannot allocate vector of size 2.0 Gb. How can I import the file then? The version I used: sessionInfo() R version 2.8.1 (2008-12-22)

Re: [R] How to read large file

2009-03-16 Thread Gabor Grothendieck
Example 6e of http://code.google.com/p/sqldf/#Example_6._File_Input shows how to use the sqldf package to read a sample of rows from a file without reading the file itself into R. On Mon, Mar 16, 2009 at 11:11 AM, Dongyan Song yzhsk...@hotmail.com wrote: Dear all, I want to read a 1.2 G file

Re: [R] Uniroot and Newton-Raphson Anomaly

2009-03-16 Thread Ravi Varadhan
Harold, Here is an approach that might work well for all the students (it doesn't use derivatives): require(BB) # I use the dfsane() root-finder from the B package f - function(a,R,c1,q1) sum((1 - (1-R)^a)^(1/a)) - c1 * q1 # I have re-named the variables

[R] Dear R help,

2009-03-16 Thread Elfatih Abdel-Rahman
Dear R help, I really appreciate if someone could tell me how I do leave-one-out cross validation with randomForest for regression, in fact I could do 10 fold with cv-k. Kind regards Elfatih - Elfatih M. Abdel-Rahman PhD student School

[R] [solved] How to combine xtable and minipage with Sweave ?

2009-03-16 Thread Ptit_Bleu
Thanks to you all (and a bit to some websites dedicated to latex), I finally managed to have a table and a graph side-by-side !!! I'm not sure it will interested people, but here is the script working for me (before the graph was below the table only because I have an empty line between the 2

Re: [R] read.xls question

2009-03-16 Thread Hans-Peter Suter
2009/3/15 David Scott d.sc...@auckland.ac.nz: On Sat, 14 Mar 2009, Mathew, Abraham T wrote: [snip] Just checking on CRAN xlsReadWrite is not currently available. There is an archived version available however. Thanks for mentioning. The reason is that CRAN is strict(er) now with binary

[R] Scatterplot

2009-03-16 Thread Shreyasee
Hello, I am having a dataset in which the values range from 0 - 70. When I tried to plot a scatterplot using following command plot(x,y) I got the plot where the x-axis and y-axis were displayed as 0e+00, 1e+05, 2e+05, ... I saw the dataset, but there is no erroneous value. Can anybody

[R] readline() is reading the entire line and not just the user input

2009-03-16 Thread Farrel Buchinsky
I just discovered readline to provide user input. It is working as one would expect in Rgui but not in Rterm. In Rterm it reads the entire line of code not just the user input. How do I fix and what is the issue that I should read about? Here is my two line script: ANSWER - readline(prompt=Are

Re: [R] readline() is reading the entire line and not just the user input

2009-03-16 Thread jim holtman
Seems to work fine with RTERM under my Windows: ANSWER - readline(prompt=Are you a satisfied R user? ) Are you a satisfied R user? yes ANSWER [1] yes sessionInfo() R version 2.8.1 (2008-12-22) i386-pc-mingw32 On Mon, Mar 16, 2009 at 12:44 PM, Farrel Buchinsky fjb...@gmail.com wrote: I just

[R] dpaste, pastebin etc. Syntax highlighting for R?

2009-03-16 Thread Jose Quesada
-BEGIN PGP SIGNED MESSAGE- Hash: SHA1 Hi All, There are many sites to get bite-sized chuncks of code. The nice thing is that they use syntax highlighting. Examples: http://dpaste.com/ http://pastie.org/ http://pastebin.com/ Twitter with a link to a syntax highlighting site would do

[R] R-code in html help pages: syntax highlighting

2009-03-16 Thread Jose Quesada
-BEGIN PGP SIGNED MESSAGE- Hash: SHA1 Hi, A follow up to my previous post. It'd be good to have syntax highlighting in the html help pages. There is highlight (http://www.andre-simon.de/) which seems to be open source. Since code is delimited in the help file format, it should not be

Re: [R] R-code in html help pages: syntax highlighting

2009-03-16 Thread Romain Francois
Hi Jose, highlight is fairly easy to use from command line. I have something like this when I use R from a console that supports color characters: highlight - function( fun ){ fun - match.fun( fun ) tf - tempfile() dump( fun, tf ) system( sprintf( highlight -A -S R %s, tf ) ) unlink( tf

[R] Best way to do temporal joins in R?

2009-03-16 Thread Jonathan Greenberg
I've been playing with zoo a bit, and it seems ok except it doesn't support non-unique time stamps when performing joins. I have two databases which contain a dataframe of a Date object (with the time, not just MM/DD/YY), e.g.: DB 1: UniqueID,Date1,Data 1,Data 2 DB 2: Date2, Station, Data 3

[R] errors when install RSQLite

2009-03-16 Thread Dongyan Song
Dear all, I am trying to install RSQLite package since I want to install sqldf, and I used install.packages(RSQLite) first, which gave Error message as below: make: *** [RS-DBI.o] Error 1 chmod: cannot access `/usr/lib/R/library/RSQLite/libs/*': No such file or directory ERROR: compilation

Re: [R] R-code in html help pages: syntax highlighting

2009-03-16 Thread Duncan Murdoch
On 3/16/2009 1:49 PM, Jose Quesada wrote: -BEGIN PGP SIGNED MESSAGE- Hash: SHA1 Hi, A follow up to my previous post. It'd be good to have syntax highlighting in the html help pages. This needs two things: making the highlight program available to R, and then adding it to the Rd

Re: [R] R-code in html help pages: syntax highlighting

2009-03-16 Thread Romain Francois
Duncan Murdoch wrote: On 3/16/2009 1:49 PM, Jose Quesada wrote: -BEGIN PGP SIGNED MESSAGE- Hash: SHA1 Hi, A follow up to my previous post. It'd be good to have syntax highlighting in the html help pages. This needs two things: making the highlight program available to R, and

Re: [R] R-code in html help pages: syntax highlighting

2009-03-16 Thread Duncan Murdoch
On 3/16/2009 2:50 PM, Romain Francois wrote: Duncan Murdoch wrote: On 3/16/2009 1:49 PM, Jose Quesada wrote: -BEGIN PGP SIGNED MESSAGE- Hash: SHA1 Hi, A follow up to my previous post. It'd be good to have syntax highlighting in the html help pages. This needs two things: making

Re: [R] Best way to do temporal joins in R?

2009-03-16 Thread Jonathan Greenberg
Sorry for the immediate follow-up, but Phil Spector correctly reminded me this is a lot easier for the community I provide some sample data, so I'm attaching 3 small CSVs to this email: species_data_Rexample.csv contains the field data (which species was ID'd and what time it was ID'd),

[R] Implementation of ZIGP in R

2009-03-16 Thread kohi001
Please any one who can help, i have count data which have many zero and variance larger than mean. I need to applyZero inflated generalized poisson regression. I have installed the ZIGP package but I failed to implement. My file look loke count Species treatment 1 duiker

Re: [R] permutations in R

2009-03-16 Thread Greg Snow
Look at the permutations function in the Combinations package. Using the fun argument may accomplish what you want. If not, there are references on the help page to other code that may work for you. Hope this helps, -- Gregory (Greg) L. Snow Ph.D. Statistical Data Center Intermountain

[R] - help - LARS LASSO based on glm

2009-03-16 Thread eitan lavi
Hello I'm trying to use LARS on a glm: log.1 = glm(Train.Data$TL~(.),Train.Data,family = binomial) LARS.1 = lars(model.matrix(log.1),Train.Data$TL,type = lar,max.steps = ncol(model.matrix(log.1))-1) self.fit = predict.lars(LARS.1,model.matrix(log.1),type=fit,mode=step)$fit

[R] plot lines at 0,0 axis

2009-03-16 Thread Mohan Singh
Hii I am trying to plot lines at (0,0) axis where my scatterplot goes between positive n negative values for x and y axis i can plot point using points(0,0) but if i want to draw lines along it, can't seem to get it right Mohan __

Re: [R] Testing for Inequality à la select case

2009-03-16 Thread Stavros Macrakis
On Sun, Mar 15, 2009 at 11:46 PM, diegol diego...@gmail.com wrote: ...Steve, ... Actually Stavros (ΣΤΑΥΡΟΣ), not Stephen/Steve (ΣΤΕΦΑΝΟΣ). Both Greek, but different names. I still don't understand the analogy. I agree that in this case the R approach is vectorized. However, your function

Re: [R] matrix multiplication, tensor product, block-diagonal and fast computation

2009-03-16 Thread Camarda, Carlo Giovanni
Dear Chuck, thanks a lot for your reply. Nevertheless, and though not so elegant, the for-loop solution seems faster (see below). Best, Carlo Giovanni Camarda m - 50 n - 40 nl - 1000 A - array(1:(m*m*n), dim=c(m,m,n)) M - matrix(1:(m*n), nrow=m, ncol=n) # for-loop solution

Re: [R] plot lines at 0,0 axis

2009-03-16 Thread David Winsemius
x- rnorm(100) y - rnorm(100) plot(x,y) abline(h=0) abline(v=0) box(col=grey) On Mar 16, 2009, at 3:55 PM, Mohan Singh wrote: Hii I am trying to plot lines at (0,0) axis where my scatterplot goes between positive n negative values for x and y axis i can plot point using points(0,0)

[R] : put grids in a plot at specific x points

2009-03-16 Thread Laura Rodriguez Murillo
Hi dear all! I wonder if anybody can help me with this: I have a plot: allchr - read.table(allchrtog.txt, header=F) attach(allchr) names(allchr) [1] V1 V2 V3 V4 plot(V1,V3, type=n) lines(V1,V3) and I want to add grids to the plot but just at specific x points. So for example, the first

Re: [R] - help - LARS LASSO based on glm

2009-03-16 Thread Max Kuhn
how can I fix it so that the LARS predicts values between 0 and 1 ? I don't know, but you can use glmnet to accomplish the same thing. -- Max __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the

Re: [R] : put grids in a plot at specific x points

2009-03-16 Thread David Winsemius
grid refers you to abline whenmore fine tuning is required, i.e. when the situation is not a regular grid. ?abline abline(v=286.16, lty=dotted) If you had a vector, vec, then this should do the trick: abline(v=vec, lty=dotted) -- David Winsemius On Mar 16, 2009, at 4:25 PM, Laura

[R] sort data frame columns according to a list

2009-03-16 Thread Jie Zhang
HI,   I have searched the R-help achive, trying to find a way to sort columns in a data frame according to a specific list, but not be able to find a general way to do it.   My data frame contains columns, with each one named from a sample ID, and each sample was measured in a specific time,

Re: [R] R-code in html help pages: syntax highlighting

2009-03-16 Thread hadley wickham
It would be pretty easy to use the output from the R parser (which is never wrong, is it?), and dump some markup out of it. For example the showTree function in codetools dumps an R expression as Lisp, this is not too far from generating html, or any other markup. As this sounds like fun,

Re: [R] Scatterplot

2009-03-16 Thread David Winsemius
options(scipen=10) On Mar 16, 2009, at 12:34 PM, Shreyasee wrote: Hello, I am having a dataset in which the values range from 0 - 70. When I tried to plot a scatterplot using following command plot(x,y) I got the plot where the x-axis and y-axis were displayed as 0e+00, 1e+05,

Re: [R] : put grids in a plot at specific x points

2009-03-16 Thread Laura Rodriguez Murillo
Thank you! that worked nicelly! Laura 2009/3/16 David Winsemius dwinsem...@comcast.net: grid refers you to abline whenmore fine tuning is required, i.e. when the situation is not a regular grid. ?abline abline(v=286.16, lty=dotted) If you had a vector, vec,  then this should do the

[R] Cointegration Vectors

2009-03-16 Thread EDUARDO LONGO, BANCO SAFRA S.A.
Hi, I am trying to test the cointegration among 5 time series, grouped in pairs. I would like to save in a table the cointegration vectors for the 10 tests. I used the urca package, but I dont know how to extract the data only for the cointegration vector. Thanks in advance for help !

Re: [R] sort data frame columns according to a list

2009-03-16 Thread David Winsemius
So your second dataframe, say DF2, has sample IDs in one of its columns, say colnames, (or are they in in the names?) ( Examples and output of str() are requested to avoid exactly such ambiguity.) If it is the former, then; DF1[ , DF2colnames ] # could do it. -- David Winsemius On Mar

Re: [R] sort data frame columns according to a list

2009-03-16 Thread David Winsemius
On Mar 16, 2009, at 5:05 PM, David Winsemius wrote: So your second dataframe, say DF2, has sample IDs in one of its columns, say colnames, (or are they in in the names?) ( Examples and output of str() are requested to avoid exactly such ambiguity.) If it is the former, then; DF1[ ,

Re: [R] Save the elements of an atomic vector to a text fil

2009-03-16 Thread johnhj
I can save a atomic vector with the write.table function so far... But how can I save 2 vectors to a text file ? I tried to do it by merging 2 vector to 1 vector like this here: summary_vector-c(vector1, vector2) and tried to write the contents of summary_vector to a textfile with:

Re: [R] How to read large file

2009-03-16 Thread Darren Norris
Try putting memory limit linux or mem.limits in the nabble help. Or type ?mem.limits in the R gui. Also have a look at: http://www.nabble.com/RAM%2C-swap%2C-Error%3A-cannot-allocate-vector-of-size%2C-Linux%3A-tt11680213.html#a11680213 Think you also need to check / understand how R deals with

Re: [R] R-code in html help pages: syntax highlighting

2009-03-16 Thread Romain Francois
hadley wickham wrote: It would be pretty easy to use the output from the R parser (which is never wrong, is it?), and dump some markup out of it. For example the showTree function in codetools dumps an R expression as Lisp, this is not too far from generating html, or any other markup. As this

Re: [R] Save the elements of an atomic vector to a text fil

2009-03-16 Thread David Winsemius
http://cran.r-project.org/doc/manuals/R-data.pdf -- David Winsemius On Mar 16, 2009, at 4:35 PM, johnhj wrote: I can save a atomic vector with the write.table function so far... But how can I save 2 vectors to a text file ? I tried to do it by merging 2 vector to 1 vector like this

Re: [R] plot lines at 0,0 axis

2009-03-16 Thread p_connolly
Quoting Mohan Singh mohan.si...@ucd.ie: Hii I am trying to plot lines at (0,0) axis where my scatterplot goes between positive n negative values for x and y axis i can plot point using points(0,0) but if i want to draw lines along it, can't seem to get it right ?abline maybe abline(h

[R] Logistic regression and mixed effects / hierarchial structure

2009-03-16 Thread Kjetil Halvorsen Løland
Hi, I’ve tried to find a general approach to my problem without any success, although this might very well be due to my inexperience with R help resources (and statistics in general). My general problem is a straightforward 2 by 2 table (“Belonging to the upper quartile” vs “not-belonging

[R] Get user system name

2009-03-16 Thread Etienne Bellemare Racine
I would like to get the name of the user form the system. Is it possible ? Something like system.user() returning something like [1] etber12 Thanks, Etienne [[alternative HTML version deleted]] __ R-help@r-project.org mailing list

Re: [R] Get user system name

2009-03-16 Thread Sundar Dorai-Raj
Assuming USER is defined on your system then Sys.getenv(USER) ought to work. --sundar On Mon, Mar 16, 2009 at 3:04 PM, Etienne Bellemare Racine etienn...@gmail.com wrote: I would like to get the name of the user form the system. Is it possible ? Something like  system.user() returning

Re: [R] Get user system name

2009-03-16 Thread Etienne Bellemare Racine
It seems like it is not defined. I get an empty string. Is there a workaround or another solution ? -- Etienne Sundar Dorai-Raj a écrit : Assuming USER is defined on your system then Sys.getenv(USER) ought to work. --sundar On Mon, Mar 16, 2009 at 3:04 PM, Etienne Bellemare Racine

Re: [R] Get user system name

2009-03-16 Thread Sundar Dorai-Raj
If that's the case, it's best to post the output from version to this thread. Type version from the R command prompt and paste the results here. On Mon, Mar 16, 2009 at 3:14 PM, Etienne Bellemare Racine etienn...@gmail.com wrote: It seems like it is not defined. I get an empty string. Is there a

Re: [R] Get user system name

2009-03-16 Thread Etienne Bellemare Racine
This one is the one working for me. Thanks ! -- Etienne Gábor Csárdi a écrit : Etienne, Sys.info()[user] user csardi Best, Gabor On Mon, Mar 16, 2009 at 11:04 PM, Etienne Bellemare Racine etienn...@gmail.com wrote: I would like to get the name of the user form the

Re: [R] Scatterplot

2009-03-16 Thread Shreyasee
Thanks a lot. It works... On Tue, Mar 17, 2009 at 4:53 AM, David Winsemius dwinsem...@comcast.netwrote: options(scipen=10) On Mar 16, 2009, at 12:34 PM, Shreyasee wrote: Hello, I am having a dataset in which the values range from 0 - 70. When I tried to plot a scatterplot using

[R] R with MKL

2009-03-16 Thread Cristián
Hi, I have seen a lot of problems from people trying to compile R with MKL. So I am writing my experience in case it helps and to ask one question. I installed R-2.8.1.patched in Ubuntu 9.04 (gcc 4.3.3) using MKL 10.1.1.019. I configured correctly (following MKL userguide) with : sudo

Re: [R] permutations in R

2009-03-16 Thread Stavros Macrakis
Greg, Thanks for helping this user. I assume you mean the permn function in the combinat package? For a new user (including me), it is not obvious how to get from the permutations function in the Combinations package to that. I see there is also a function gtools::permutations. The gtools

Re: [R] sort data frame columns according to a list

2009-03-16 Thread Jun Shen
Jie, df[match(df.t$ID,df$ID),] should work. df.t$ID is your target sequence. df$ID is your current sequence. Jun On Mon, Mar 16, 2009 at 3:41 PM, Jie Zhang jieu...@yahoo.com wrote: HI, I have searched the R-help achive, trying to find a way to sort columns in a data frame according to a

Re: [R] permutations in R

2009-03-16 Thread David Winsemius
On Mar 16, 2009, at 6:43 PM, Stavros Macrakis wrote: Greg, Thanks for helping this user. I assume you mean the permn function in the combinat package? For a new user (including me), it is not obvious how to get from the permutations function in the Combinations package to that. I see there

Re: [R] Best way to do temporal joins in R?

2009-03-16 Thread Gabor Grothendieck
There was nothing attached. On Mon, Mar 16, 2009 at 3:11 PM, Jonathan Greenberg greenb...@ucdavis.edu wrote: Sorry for the immediate follow-up, but Phil Spector correctly reminded me this is a lot easier for the community I provide some sample data, so I'm attaching 3 small CSVs to this email:

Re: [R] Get user system name

2009-03-16 Thread Duncan Murdoch
On 16/03/2009 6:09 PM, Sundar Dorai-Raj wrote: Assuming USER is defined on your system then Sys.getenv(USER) ought to work. In Windows you need to look at Sys.getenv(USERNAME) instead. Duncan Murdoch __ R-help@r-project.org mailing list

Re: [R] Best way to do temporal joins in R?

2009-03-16 Thread Jonathan Greenberg
Weird -- the email was sent through my gmail account, looks like the .csvs got intercepted somewhere along the way. At any rate, I placed them on a website: http://cstars.ucdavis.edu/~jongreen/temp/temporal_join_R/ --j Gabor Grothendieck wrote: There was nothing attached. On Mon, Mar 16,

Re: [R] R-code in html help pages: syntax highlighting

2009-03-16 Thread Duncan Murdoch
On 16/03/2009 5:06 PM, Romain Francois wrote: hadley wickham wrote: It would be pretty easy to use the output from the R parser (which is never wrong, is it?), and dump some markup out of it. For example the showTree function in codetools dumps an R expression as Lisp, this is not too far from

[R] Splitting variables

2009-03-16 Thread Rickard Grönholm
Dear R Users,I have a data set with several variables. One of these variables is the number of farms per county and I'd like to split this variable into two domains: (a) counties with fewer than 600 farms and (b) counties with 600 or more farms I've tried solving this problem for hours, but I'm

Re: [R] Splitting variables

2009-03-16 Thread jim holtman
y - x[x600] z - x[x=600] 2009/3/16 Rickard Grönholm gronholm.rich...@gmail.com: Dear R Users,I have a data set with several variables. One of these variables is the number of farms per county and I'd like to split this variable into two domains: (a) counties with fewer than 600 farms and (b)

[R] Non-Linear Optimization - Query

2009-03-16 Thread Lars Bishop
Dear All, I couple of weeks ago, I’ve asked for a package recommendation for nonlinear optimization. In my problem I have a fairly complicated non-linear objective function subject to one non-linear equality constrain. I’ve been suggested to use the *Rdonlp2* package, but I did not get any

[R] - help - predicting with glmnet/lars for dataframes with different nrow then the train set

2009-03-16 Thread eitan lavi
Hello I'm having trouble using lars and glmnet functions to predict on a new data set with different nrow then the original : for instance: = log.1 = glm(temp.data$TL~(.),temp.data,family = binomial,x=TRUE,y=TRUE) nrow(test.data) != nrow(temp.data # == TRUE Val.frame =

[R] Calculate across columns

2009-03-16 Thread timbean
I'd like to perform some calculations across a series of columns. Is there an easy way to do so? For example, I'd like to calculate the maximum value from columns 6-42. Although I can use, say: pmax(data[,6],data[,7],data[,8], etc.), I'd prefer to be able to do something like pmax(data[,6:42]).

[R] [S] ASA Southern California Chapter Applied Statistics Workshop

2009-03-16 Thread Madeline Bauer
The workshop organizing committee of the Southern California Chapter of the American Statistical Association announces the 28th Annual Workshop in Applied Statistics. Professor Colin Cameron from the department of Economics at UC Davis will give a one-day workshop titled Advances in Count Data

[R] the effect of blocking on the size of confidence intervals - analysis using lme and lmer

2009-03-16 Thread Lars Kunert
This is a follow-up mail of the effect of blocking on the size of confidence intervals - analysis using aov. In both mails I pursue the idea of using blocking factors in order to reduce the width of confidence intervals. My dataset comprises, a quantitative response variable, namely: response,

[R] the effect of blocking on the size of confidence intervals - analysis using aov

2009-03-16 Thread Lars Kunert
Hi, I am a method developer in drug discovery. I have developed a similarity searching method for small (drug-like) molecules. I want to compare the performance of my method to the performance of other methods. Similarity searching methods are commonly assessed by their ability to (re)discover a

Re: [R] Calculate across columns

2009-03-16 Thread stephen sefick
?apply look at the MARGIN argument On Mon, Mar 16, 2009 at 8:56 PM, timbean timb...@gmail.com wrote: I'd like to perform some calculations across a series of columns. Is there an easy way to do so? For example, I'd like to calculate the maximum value from columns 6-42. Although I can use,

[R] R package to automatically produce combination plot?

2009-03-16 Thread Jason Rupert
By any chance is there an R package that automatically produces the plot shown at the following link: http://n2.nabble.com/Can-R-produce-this-plot--td2489288.html That is an R package to produce on plot that has the following: (a) a vertically oriented histogram, (b) associated barplot, and

[R] Bonacich Power (memory Error) (Igraph)

2009-03-16 Thread Sur Nathan
Hello, I am trying to run Bonpower for my graph.I have posted the R objects along with this file. There is a object called G(graph) that I created for running bonpower.ummary information shows the number of Vertices and edges. When I try to run Bonacich Power I get out of memory error.I did

Re: [R] Non-Linear Optimization - Query

2009-03-16 Thread Paul Smith
On Tue, Mar 17, 2009 at 12:09 AM, Lars Bishop lars...@gmail.com wrote: I couple of weeks ago, I’ve asked for a package recommendation for nonlinear optimization. In my problem I have a fairly complicated non-linear objective function subject to one non-linear equality constrain. I’ve been

[R] set intercept in glm to non-0

2009-03-16 Thread Peter Palenchar
I know how to make the intercept 0 in glm, but I actually want to force my intercept to be 1. Is there a way I can do this? Peter __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide

[R] sweep?

2009-03-16 Thread rkevinburton
I am having a hard time understanding just what 'sweep' does. The documentation states: Return an array obtained from an input array by sweeping out a summary statistic. So what does it mean weeping out a summary statistic? Thank you. Kevin __

[R] bigglm() results different from glm()

2009-03-16 Thread Francisco J. Zagmutt
Dear all, I am using the bigglm package to fit a few GLM's to a large dataset (3 million rows, 6 columns). While trying to fit a Poisson GLM I noticed that the coefficient estimates were very different from what I obtained when estimating the model on a smaller dataset using glm(), I wrote a

  1   2   >