[R] Accidental misuse of 'R-packages' list

This should never have been on R-packages. I (as moderator) apologize the glitch in moderation. This is *not* CC'ed to that list, as it really really is only for *announcements* of new or vastly changed packages available from CRAN. Apropos: Please, package authors, keep those low volume: A good

Re: [R] issues with R Library on a Server

On Tue, May 18, 2010 at 4:03 AM, Daisy Englert Duursma daisy.duur...@gmail.com wrote: Hello, I am a bit over my head on this issue. My colleagues and I are running R off of our server. We all have admin rights and prior to yesterday we all had our own libraries. Our main system administrator

[R] GUI commands to call for a protein from protein data bank

What I am trying to do is use GUI function, traitr, and to call for a pdb file and save it and then display it. I want to call for it by taking it from the user and then displaying it on the screen. I am having problems with that. The line pdb - read.pdb(ProteinCode) where proteincode should be

[R] A problem in allocation of vector of size

Hi, r-users I happen to a problem in allocation of vector of size. When I run my R script, an error appears: Error: cannot allocate vector of size 450 Mb Could anyone happen to the same problem? Thank you for your help. Lee [[alternative HTML version deleted]]

Re: [R] Robust SE Heteroskedasticity-consistent estimation

Thanks all of you. -Message d'origine- De : Arne Henningsen [mailto:arne.henning...@googlemail.com] Envoyé : mardi 18 mai 2010 06:00 À : RATIARISON Eric Cc : r-help@r-project.org; Ott-Siim Toomet; Achim Zeileis Objet : Re: [R] Robust SE Heteroskedasticity-consistent estimation On 17

Re: [R] lme4 + R 2.11.0 + mac unavailable

David Winsemius dwinsem...@comcast.net on Mon, 17 May 2010 17:44:00 -0400 writes: On May 17, 2010, at 5:01 PM, Adam November wrote: I believe I'm working with the newest version of R (2.11.0) and I've tried a few of the most recent versions... No luck yet. Luck

Re: [R] Compiling R with --enable-R-shlib for rpy2 error

Thanks for the help! :) Regards, Padma Tan Genome Institute of Singapore 60 Biopolis Street, Genome #02-01 Singapore 138672 DID : 6478 8671 Fax : 6478 9058 email: ta...@gis.a-star.edu.sg

[R] how to select rows per subset in a data frame that are max. w.r.t. a column

Hi, I'd like to select one row in a data frame per subset which is maximal for a particular value. I'm pretty close to the solution in the sense that I can easily select the maximal values per subset using aggregate, but I can't really figure out how to select the rows in the original data frame

Re: [R] how to select rows per subset in a data frame that are max. w.r.t. a column

Hi, Maybe it's just me but I don't understand what you're trying to do. Isn't maxValues what you need? Providing a reproducible example with your data (using the function dput), the desired output, and the code you've tried would really help! Ivan Le 5/18/2010 10:58, Tim Van den Bulcke a

Re: [R] Create counter variable for subsets without a loop

take a look at the by, ave, aggregate and apply functions, perhaps one suits your needs Bart -- View this message in context: http://r.789695.n4.nabble.com/Create-counter-variable-for-subsets-without-a-loop-tp2220663p2220925.html Sent from the R help mailing list archive at Nabble.com.

Re: [R] Change order of columns in table?

On 05/17/2010 10:46 PM, someone wrote: I'm an R noob and have a (maybe) stupid question... I have a table where I have the weekdays and a number for each weekday of entries: Thats what the table looks like... Now I want to have an pie3D plot of this, but obviously the order of the weekdays are

[R] [RGL] Need help to modify current plot

Dear folks, I have created a plot on RGL device :   x = 1:6 y = seq(-12, 5, by=1) z = matrix(0, length(y), length(x)) z[13,3] = 1; z[13,4] = 1.011765 surface3d(x, y, t(z), col=rainbow(1000)) grid3d(c(x-, y-, z))   Now I want to draw 2 lines along x=3 x=4, over the surface (with different

Re: [R] [RGL] Need help to modify current plot

Megh Dal wrote: Dear folks, I have created a plot on RGL device : x = 1:6 y = seq(-12, 5, by=1) z = matrix(0, length(y), length(x)) z[13,3] = 1; z[13,4] = 1.011765 surface3d(x, y, t(z), col=rainbow(1000)) grid3d(c(x-, y-, z)) Now I want to draw 2 lines along x=3 x=4, over the surface (with

Re: [R] Graphing Titration Curves, how to graph in R

On 05/18/2010 06:54 AM, Amitoj S. Chopra wrote: I am working on a script with R to calculate pKa values. I want to be able then to graph pKa values to Energy, basically two variables which will give me titration curves. It is basically a graph with two variables and help with just graphing it

Re: [R] Graphing Titration Curves, how to graph in R

Hi Amitoj, Also look at the AquaEnv and titan packages. Jim __ R-help@r-project.org 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,

Re: [R] [RGL] Need help to modify current plot

On 18/05/2010 6:46 AM, Megh Dal wrote: Thanks Duncan for your reply. This could definitely be an answer of my query however I wanted something else. What I want is to draw 2 lines over the surface which pass through x=3,4 That's even more ambiguous than your original request, so I won't

Re: [R] [R-sig-DB] Package RPostgreSQL : Problem with dbWriteTable

On Wed, May 12, 2010 at 2:44 AM, PALMIER Patrick - CETE NP/INFRA/TRF patrick.palm...@developpement-durable.gouv.fr wrote: Hello, I have a probem with dbWriteTable method of package RPostrgreSQL The table is well added in the database but R doesn't succeed in inserting rows But If I send

Re: [R] how to select rows per subset in a data frame that are max. w.r.t. a column

Hi: Here are a couple of ways using the doBy and plyr packages: library(doBy) library(plyr) # doBy: subsetBy(~ wool + tension, subset = breaks == max(breaks), data = warpbreaks) breaks wool tension A|H 43A H A|L 70A L A|M 36A M B|H 28B

Re: [R] Create counter variable for subsets without a loop

Here are four solutions: data - cbind(state.region,as.data.frame(state.x77))[,1:2] # ave data2 - data[order(data$state.region, -data$Population), ] data2$rank - ave(data2$Population, data2$state.region, FUN = seq_len)) # by f - function(x) cbind(x[order(-x$Population), ], rank = 1:nrow(x))

[R] Query on linear mixed model

Hi R Forum I am a newbie to R and I have been amazed by what I can get my team to accomplish just by implementing Scripting routines of R in all my team's areas of interest.. Recently i have been trying to adopt R scripting routine for some analysis with longitudanal data.. I am presenting my R

Re: [R] Where has the stats-rosuda-devel mailing list gone?

Oliver, Apologies for the confusion, there was a server upgrade in the computer centre here which gave us some grief. The list should be fine now. Best regards Antony Antony Unwin Professor of Computer-Oriented Statistics and Data Analysis, Mathematics Institute, University of Augsburg,

[R] Fwd: [R-sig-DB] Package RPostgreSQL : Problem with dbWriteTable

-- Forwarded message -- From: PALMIER Patrick - CETE NP/INFRA/TRF patrick.palm...@developpement-durable.gouv.fr Date: Tue, May 18, 2010 at 8:05 AM Subject: Re: [R-sig-DB] Package RPostgreSQL : Problem with dbWriteTable To: Davis, Sean (NIH/NCI) [E] sdav...@mail.nih.gov Hello,

Re: [R] GUI commands to call for a protein from protein data bank

Amitoj S. Chopra amitojc at gmail.com writes: What I am trying to do is use GUI function, traitr, and to call for a pdb file and save it and then display it. I want to call for it by taking it from the user and then displaying it on the screen. I am having problems with that. The line pdb

[R] lattice::panel.levelplot.raster too picky with unequal spacing

Dear all, I got a couple of warnings using panel.levelplot.raster, In panel.levelplot.raster(..., interpolate = TRUE) : 'y' values are not equispaced; output will be wrong although I was quite sure my data were equally spaced (indeed, I created them with seq()). A closer look at the source

Re: [R] Query on linear mixed model

Hi Vijayan, You are really asking for this list to serve as your statistical consultant, which is not its purpose. If you have a specific problem (and if you know how to ask for help -- see the posting guide) this list is a tremendous resource. But it is not a replacement for a statistician.

Re: [R] best polynomial approximation

I guess you may be looking for the Remez algorithm. AFAIK there is no implementation in one of the R packages. You can find FORTRAN code in the Collected Algorithms of the ACM (no. 604) which probably could be called from R. There appears to exist a discrete, equi-distant(?) version as function

Re: [R] timing a function

Thank you Alexander 2010/5/17 Alexander Shenkin ashen...@ufl.edu You could also put the call to system.time inside the function itself: f = function(x) { system.time({ ... #function's code ret_val = ... }); flush.console(); return ret_val; } i s'pose you'd miss

Re: [R] [RGL] Need help to modify current plot

Thanks Duncan for your reply. This could definitely be an answer of my query however I wanted something else. What I want is to draw 2 lines over the surface which pass through x=3,4   Any better idea?   Thanks, --- On Tue, 5/18/10, Duncan Murdoch murdoch.dun...@gmail.com wrote: From: Duncan

Re: [R] lattice::panel.levelplot.raster too picky with unequal spacing

On Tue, May 18, 2010 at 6:32 PM, baptiste auguie baptiste.aug...@googlemail.com wrote: Dear all, I got a couple of warnings using panel.levelplot.raster, In panel.levelplot.raster(..., interpolate = TRUE) :  'y' values are not equispaced; output will be wrong although I was quite sure my

[R] Using the zero-inflated binomial in experimental designs

I'm trying to use the inflated binomial distribution of zeros (since 75% of the values are zeros) in a randomized block experiment with four quantitative treatments (0, 0.5, 1, 1.5), but I'm finding it difficult, since the examples available in VGAM packages like for example, leave us unsure of

[R] doubt with auto.arima

hello, i have a doubt with this function, i need get the returns values because i have to save in a variables,something like this: Invernadero-ts(x2) test-auto.arima(x2) x2.pred.ar31-predict(arima(x2,order=c(p,d,q)),n.ahead=10)$pred can i get p,d and q from auto.arima?this is my question.

[R] looking for .. dec if vector if element x

Hi to all, I am just looking for more efficient ways ;-) is there a better way instead a loop to decrease x if greater y test - c(1,3,5,7,9) decrease if greater 1 to test2 - c(1,2,4,6,8) Kind regards Knut __ R-help@r-project.org mailing list

Re: [R] looking for .. dec if vector if element x

On Tuesday 18 May 2010, Knut Krueger wrote: Hi to all, I am just looking for more efficient ways ;-) is there a better way instead a loop to decrease x if greater y test - c(1,3,5,7,9) decrease if greater 1 to test2 - c(1,2,4,6,8) Does this help? test - c(1, 3, 5, 7, 9) test[test

Re: [R] looking for .. dec if vector if element x

On 18/05/2010 7:34 AM, Knut Krueger wrote: Hi to all, I am just looking for more efficient ways ;-) is there a better way instead a loop to decrease x if greater y test - c(1,3,5,7,9) decrease if greater 1 to test2 - c(1,2,4,6,8) test2 - ifelse( test 1, test-1, test) Duncan Murdoch

Re: [R] Query on linear mixed model

On Tue, May 18, 2010 at 4:52 AM, Vijayan Padmanabhan v.padmanab...@itc.in wrote: snip This does not answer your statistical question, but I did include some ideas to simplify your script. ##For ALL Product Comparison across All Time Points. options(contrasts=c('contr.treatment','contr.poly'))

Re: [R] looking for .. dec if vector if element x

Try this also: pmax(test - 1, 1) On Tue, May 18, 2010 at 8:34 AM, Knut Krueger r...@krueger-family.de wrote: Hi to all, I am just looking for more efficient ways ;-) is there a better way instead a loop to decrease x if greater y test - c(1,3,5,7,9) decrease if greater 1 to test2 -

Re: [R] lattice::panel.levelplot.raster too picky with unequal spacing

On 18 May 2010 15:30, Deepayan Sarkar deepayan.sar...@r-project.org wrote: Maybe a better test would be isTRUE(all.equal(diff(range(diff(ux))), 0)) I'll try that out for the next release. Sounds good (and works for me), thanks. baptiste __

Re: [R] Variable variables using R ... e.g., looping over data frames with a numeric separator

On 05/17/2010 03:51 PM, Monte Shaffer wrote: for(i in 1:L-1) { dataStr = gsub(' ','',paste(fData.,i)); dataVar = eval(dataStr); ## GOAL is to grab data frame 'fData.1' and do stuff with it, then in next loop grab data frame 'fData.2' and do stuff with it } As Dan Davison said, the more

[R] scaling with relative units in plots or retrieving axes limits in plots

Dears, a way to define x and y positions in plots in relative numbers (e.g in fractions between 0 and 1 referring to relative positions inside the plot region) would really help me. One example I would need this to would be to add text via text() to a plot always at a defined spot, e.g the

[R] scaling with relative units in plots or retrieving axes limits in plots

Dears, a way to define x and y positions in plots in relative numbers (e.g in fractions between 0 and 1 referring to relative positions inside the plot region) would really help me. One example I would need this to would be to add text via text() to a plot always at a defined spot, e.g the

Re: [R] scaling with relative units in plots or retrieving axes limits in plots

par(usr) #gives the extreme plotting coordinates ?par HTH, Josh On Tue, May 18, 2010 at 7:32 AM, Jannis bt_jan...@yahoo.de wrote: Dears, a way to define x and y positions in plots in relative numbers (e.g in fractions  between 0 and 1 referring to relative positions inside the plot

Re: [R] scaling with relative units in plots or retrieving axes limits in plots

Hi, For 2., I don't know if it's possible to retrieve the axis limits, but you can surely specify them in your call to plot (with arguments xlim and ylim). That's a cheap solution and others probably have better ones. Ivan Le 5/18/2010 16:23, Jannis a écrit : Dears, a way to define x and

[R] scaling with relative units in plots or retrieving axes limits in plots

Dears, a way to define x and y positions in plots in relative numbers (e.g in fractions between 0 and 1 referring to relative positions inside the plot region) would really help me. One example I would need this to would be to add text via text() to a plot always at a defined spot, e.g the

Re: [R] scaling with relative units in plots or retrieving axes limits in plots

Thanks for the replies! If anybody encounters a similar problem, the function that now does what I wanted is attached below. Best Jannis trnsf.coords = function(array_x,array_y) # This function transfers relative coordinates between 0 and 1 for two arrays with x # and y values into the

Re: [R] help color coding map in R

Thank you this was helpful Chris Anderson Data Analyst Medical Affairs wk: 925-677-4870 cell: 707-315-8486 Fax:925-677-4670 -Original Message- From: foolish.andr...@gmail.com [mailto:foolish.andr...@gmail.com] On Behalf Of Felix Andrews Sent: Monday, May 17, 2010 5:55 PM To: Anderson,

Re: [R] sample

Sorry, I made two mistakes. The first was matching the female with the male. The second was 2 variables should be selected randomly every time. Followed is a revised copy: ## Import data. moms - read.delim(females.txt, sep = , stringsAsFactors = FALSE, header = TRUE) dads -

[R] Counting Frequencies in Data Frame

Hi, I am sure there is an easy way to do it, but I can't find it. I have a data frame that has 15 columns and 7000 rows. The only values inside the data.frame are aa, ab, bb as you can see an example bellow. 1 2 3 1 aa ab ab 2 ab ab ab 3 aa aa aa 4 bb bb bb What I would like to do, is to

Re: [R] looking for .. dec if vector if element x

Henrique Dallazuanna schrieb: Try this also: pmax(test - 1, 1) O test - c(1,3,5,7,9,11,12,13,14) test test - pmax(test - 1, 1) test This works for 1 what about if I would dec 11: to 14 to close the gap between 9 and 10 ? I did not find the answer with the help file Thank you Knut

[R] C function call in R

dear all, we am trying to improve the performance of my R code, with the implentation of some function with custom C code. we found difficult to import and export/import data structure such us matrices or data.frame into the external C functions. we already tried the solution from Writing R

[R] Fatal error that doesn't let me start R

Hi, all. I have R installed in my computer. I guess I did something in my previous session, and now every time I start R, I find the following message: Fatal error: unable to restore saved data in .RData I uninstalled R and installed it again and I'm still getting this message. Can anyone

Re: [R] C function call in R

Hi, On Tue, May 18, 2010 at 10:50 AM, John Lande john.land...@gmail.com wrote: dear all, we am trying to improve the performance of my R code, with the implentation of some function with custom C code. we found difficult to import and export/import data structure such us matrices or

[R] Repeating Name for Rows

Hello, I have a large data frame (47:2186), where i want to label every 12th row. This command works, Day -rep(97:278, each = 12) However i need 97 to only labeled 11 rows and then from 98:278 can be labeled every 12 times. Thanks! __

Re: [R] Fatal error that doesn't let me start R

On 18/05/2010 10:51 AM, gbre...@ssc.wisc.edu wrote: Hi, all. I have R installed in my computer. I guess I did something in my previous session, and now every time I start R, I find the following message: Fatal error: unable to restore saved data in .RData I uninstalled R and installed it

Re: [R] Fatal error that doesn't let me start R

Hi, On Tue, May 18, 2010 at 10:51 AM, gbre...@ssc.wisc.edu wrote: Hi, all. I have R installed in my computer.  I guess I did something in my previous session, and now every time I start R, I find the following message: Fatal error: unable to restore saved data in .RData I uninstalled R

Re: [R] Counting Frequencies in Data Frame

Hi, Others will have fancier solutions, but is the way I would do it: dat - read.table(textConnection(1 2 3 1 aa ab ab 2 ab ab ab 3 aa aa aa 4 bb bb bb), header=TRUE) closeAllConnections() countAB - function(x) { aa - length(which(x == aa)) ab - length(which(x == ab)) bb - length(which(x

[R] (no subject)

Hello I have a data array with soil variables (caperf), in which the variable clay is factor (as I see entering str(caperf)) . I need to do a regression model, so I need to have arcilla (=clay) as a numeric variable.  For that I have entered as.numeric(as.character(arcilla)) and even entering

[R] avoiding reinstall already installed library

Dear R-experts, I am installing new libraries using install.packages(ggplot2,dependencies=T). But I perceive that many dependencies are already installed. As I am using a low-band internet, how can avoid reinstall installed libraries? cheers milton [[alternative HTML version deleted]]

Re: [R] Counting Frequencies in Data Frame

-Original Message- From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org] On Behalf Of M.Ribeiro Sent: Tuesday, May 18, 2010 7:13 AM To: r-help@r-project.org Subject: [R] Counting Frequencies in Data Frame Hi, I am sure there is an easy way to do it, but I

[R] Re: Change class factor to numeric

sorry I had a mistake sending my question without a subject. I do resend again. Please excuse me. Hello I have a data array with soil variables (caperf), in which the variable clay is factor (as I see entering str(caperf)) . I need to do a regression model, so I need to have arcilla

Re: [R] avoiding reinstall already installed library

Hi, On Tue, May 18, 2010 at 11:36 AM, milton ruser milton.ru...@gmail.com wrote: Dear R-experts, I am installing new libraries using install.packages(ggplot2,dependencies=T). But I perceive that many dependencies are already installed. As I am using a low-band internet, how can avoid

Re: [R] Re: Change class factor to numeric

Try arcilla-as.numeric(as.character(clay)) best milton On Tue, May 18, 2010 at 12:36 PM, Arantzazu Blanco Bernardeau aramu...@hotmail.com wrote: sorry I had a mistake sending my question without a subject. I do resend again. Please excuse me. Hello I have a data array with soil

Re: [R] avoiding reinstall already installed *package*

{ I've modified the subject; I can't stand it hitting square into my face ... } mr == milton ruser milton.ru...@gmail.com on Tue, 18 May 2010 12:36:23 -0300 writes: mr Dear R-experts, mr I am installing new libraries using mr install.packages(ggplot2,dependencies=T). mr

Re: [R] BRugs under Linux?

On 17.05.2010 18:46, Kevin E. Thorpe wrote: Hello. In this post: http://finzi.psych.upenn.edu/Rhelp10/2010-March/233815.html Uwe Ligges suggests using BRugs rather than R2WinBUGS under windows. He also notes that it is not in the main CRAN repository, but it is in extras which is a default

Re: [R] Re: Change class factor to numeric

Hi, I think that providing the output from str(data array or whatever you have) would help. Because, for now, we don't have much idea of what you really have. Moreover, some sample data is always welcomed (using the function dput for example) Ivan Le 5/18/2010 17:36, Arantzazu Blanco

Re: [R] survey package: weights used in svycoxph()

On Mon, 17 May 2010, Vinh Nguyen wrote: Dear R-help, Let me know if I should email r-devel instead of this list. This message is addressed to Professor Lumley or anyone familiar with the survey package. Does svycoxph() implement the method outlined in Binder 1992 as referenced in the help

Re: [R] (no subject)

Hi, Sorry, I'm not really getting what going on here ... perhaps having more domain knowledge would help me make better sense of our question. In particular: On Tue, May 18, 2010 at 11:35 AM, Arantzazu Blanco Bernardeau aramu...@hotmail.com wrote: Hello I have a data array with soil

Re: [R] Re: Change class factor to numeric

Hello so, here you have the output of the data frame. The data frame comes from a csv file. I could take Gr_2 instead of arcilla, because it is the same value... but curiously, it is a factor as well. str(caperf) 'data.frame':    556 obs. of  38 variables:  $ Hoja    : int  818 818 818

Re: [R] avoiding reinstall already installed *package*

On 18-May-10 15:49:37, Martin Maechler wrote: { I've modified the subject; I can't stand it hitting square into my face ... } mr == milton ruser milton.ru...@gmail.com on Tue, 18 May 2010 12:36:23 -0300 writes: mr Dear R-experts, mr I am installing new libraries using mr

Re: [R] Repeating Name for Rows

If I perceive the issue: Day - c(rep(97, each = 11), rep(98:278, each = 12)) - A R learner. -- View this message in context: http://r.789695.n4.nabble.com/Repeating-Name-for-Rows-tp2221457p2221492.html Sent from the R help mailing list archive at Nabble.com.

Re: [R] (no subject)

Hello Well, the problem is, that arcilla is the percentage of clay in the soil sample. So, for linear model, I need to work with that number or value. Now, R thinks that arcilla (arcilla means clay in spanish), is a factor, and gives me the value as a factor, so the output of the linear model

Re: [R] timing a function

Dear all, Just one last question. There seems to be no problem in writing z = system.time(y - f(x)) or z - system.time(y - f(x)) Then z contains the named vector of the elapsed times, and y the value of the function f(x). Am I right ? Thank you very much, Gustave 2010/5/17 Alexander

Re: [R] avoiding reinstall already installed *package*

Hi Martin, thanks for your reply, and very thanks for your kind tips about package and library So, I was trying to understand *why* we load packages using library(). I suggest that developers killl the problem on its root, deleting library function :-) Good to know already installed packages will

Re: [R] Re: Change class factor to numeric

Hi again, If you used the function read.table() to read from a csv file into a data.frame, it is weird that numeric data are converted into factors. I would check in the original data that you don't have a typo somewhere. I don't know all the possibilities, but a special character can

Re: [R] (no subject)

One last thing: before you take my advice on how to recode your nominal/categorical clay variable for your regression model, take some time to see how other people talk about this and do some searching on phrases like regression model with nominal variables (that's just the one I used). You'll

[R] unsigned 4 byte number

Does anybody know how to read in a unsigned 4 byte number from a binary file? According to the help for readBin, the signed argument only applies to size=1 or 2. But if you declare any size larger it assumes it is signed? Thanks -- View this message in context:

Re: [R] (no subject)

Arantzazu, Your problem is that the data were probably imported from Excel where you had at least one cell containing #N/A. You need to replace those cases in your dataframe with NA. Then you should be able to do as.numeric(as.character(arcilla)). -Peter Ehlers On 2010-05-18 10:07, Arantzazu

Re: [R] Using the zero-inflated binomial in experimental designs

Ivan Allaman ivanalaman at yahoo.com.br writes: I'm trying to use the inflated binomial distribution of zeros (since 75% of the values are zeros) in a randomized block experiment with four quantitative treatments (0, 0.5, 1, 1.5), but I'm finding it difficult, since the examples available

[R] Function that is giving me a headache- any help appreciated (automatic read )

note: whole function is below- I am sure I am doing something silly. when I use it like USGS(input=precipitation) it is choking on the precip.1 - subset(DF, precipitation!=NA) b - ddply(precip.1$precipitation, .(precip.1$gauge_name), cumsum) DF.precip - precip.1 DF.precip$precipitation -

[R] FW: Re: Change class factor to numeric

Hello everybody the problem has been solved. It was my mistake not to be ensured that any N/A had dissapeared. I do apologize for the inconveniences caused ;) friendly greeting sfrom Spain! Arantzazu Blanco Bernardeau Dpto de Química Agrícola, Geología y Edafología Universidad de

[R] Getting dates in an SPSS file in right format.

Dear all, I am trying to read an SPSS file into a data frame in R using method read.spss(), sample - read.spss(file.name,to.data.frame=TRUE) But dates in the data.frame 'sample' are coming as integers and not in the actual date format given in the SPSS file. Appreciate if anyone can help

Re: [R] avoiding reinstall already installed *package*

On 2010-05-18 10:05, (Ted Harding) wrote: On 18-May-10 15:49:37, Martin Maechler wrote: { I've modified the subject; I can't stand it hitting square into my face ... } mr == milton rusermilton.ru...@gmail.com on Tue, 18 May 2010 12:36:23 -0300 writes: mr Dear R-experts, mr

Re: [R] survey package: weights used in svycoxph()

On Tue, May 18, 2010 at 8:50 AM, Thomas Lumley tlum...@u.washington.edu wrote:   I don't believe so since svycoxph() calls coxph() of the survival package and weights are applied once in the estimating equation.  If the weights are implemented in the ratio, could you point me to where in

Re: [R] Getting dates in an SPSS file in right format.

On 5/18/2010 12:38 PM, Praveen Surendran wrote: Dear all, I am trying to read an SPSS file into a data frame in R using method read.spss(), sample - read.spss(file.name,to.data.frame=TRUE) But dates in the data.frame 'sample' are coming as integers and not in the actual date

Re: [R] Function that is giving me a headache- any help appreciated (automatic read )

I don't think you can do this precipitation!=NA) have a look at ?is.na --- On Tue, 5/18/10, stephen sefick ssef...@gmail.com wrote: From: stephen sefick ssef...@gmail.com Subject: [R] Function that is giving me a headache- any help appreciated (automatic read ) To: r-help@r-project.org

Re: [R] avoiding reinstall already installed *package*

On 18-May-10 16:42:40, Peter Ehlers wrote: On 2010-05-18 10:05, (Ted Harding) wrote: On 18-May-10 15:49:37, Martin Maechler wrote: { I've modified the subject; I can't stand it hitting square into my face ... } mr == milton rusermilton.ru...@gmail.com on Tue, 18 May 2010 12:36:23

[R] Sweave counter

Dear R users, I am using the Sweave package and I am doing some MCMC. I have a loop function for my MCMC. Every 100 iterations, I want the number of iterations already done to appear on my screen (but not on the final document). Is that possible ? Usually I can rely on if (i%%100==0) cat(i,

Re: [R] proportion of treatment effect by a surrogate (fitting multivariate survival model)

On Mon, May 17, 2010 at 7:42 PM, Vinh Nguyen vqngu...@uci.edu wrote: Dear R-help, I would like to compute the variance for the proportion of treatment effect by a surrogate in a survival model (Lin, Fleming, and De Gruttola 1997 in Statistics in Medicine).  The paper mentioned that the

Re: [R] unsigned 4 byte number

-Original Message- From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org] On Behalf Of tetonedge Sent: Tuesday, May 18, 2010 9:14 AM To: r-help@r-project.org Subject: [R] unsigned 4 byte number Does anybody know how to read in a unsigned 4 byte number from a

Re: [R] Sweave counter

Hi Jimmy, You can use message() instead of cat(): if (i%%100==0) message(i) Best, Ista On Tuesday 18 May 2010 1:10:51 pm Jimmy Söderly wrote: Dear R users, I am using the Sweave package and I am doing some MCMC. I have a loop function for my MCMC. Every 100 iterations, I want the number of

Re: [R] Sweave counter

Under UNIX I usually write something like system(paste(echo ',...,', sep=sep)) where I replace the '...' with whatever I need to show. This would need some adjustments for non-scalars, though. Benno Am 18.Mai.2010 um 19:10 schrieb Jimmy Söderly: Dear R users, I am using the

[R] Glasso question

The glassolist function says that if rholist is set to NULL, then 10 values in a (hopefully reasonable) range are used. I'm trying to figure out how these 10 values are chosen, I can't find this in any documentation. Thanks for your help! -- View this message in context:

Re: [R] (no subject)

On May 18, 2010, at 12:07 PM, Arantzazu Blanco Bernardeau wrote: Hello Well, the problem is, that arcilla is the percentage of clay in the soil sample. So, for linear model, I need to work with that number or value. Now, R thinks that arcilla (arcilla means clay in spanish), is a

[R] Res: Using the zero-inflated binomial in experimental designs

Hi Ben! First I thank you for your attention. Unfortunately, the ANOVA does not work with vglm. In another email, Rafael warned me that actually a lot of zeros does not necessarily imply a distribution of zeros binomail inflated. So how could I test if my variable is or not a binomial zero

Re: [R] C function call in R

John Lande john.land...@gmail.com writes: dear all, we am trying to improve the performance of my R code, with the implentation of some function with custom C code. we found difficult to import and export/import data structure such us matrices or data.frame into the external C functions.

[R] Maximization of quadratic forms

Dear R Help, I am trying to fit a nonlinear model for a mean function $\mu(Data_i, \beta)$ for a fixed covariance matrix where $\beta$ and $\mu$ are low- dimensional. More specifically, for fixed variance-covariance matrices $\Sigma_{z=0}$ and $\Sigma_{z=1}$ (according to a binary covariate

Re: [R] Getting dates in an SPSS file in right format.

On 05/18/2010 11:52 AM, Chuck Cleland wrote: On 5/18/2010 12:38 PM, Praveen Surendran wrote: Dear all, I am trying to read an SPSS file into a data frame in R using method read.spss(), sample- read.spss(file.name,to.data.frame=TRUE) But dates in the data.frame 'sample' are coming as

Re: [R] survey package: weights used in svycoxph()

On Tue, 18 May 2010, Vinh Nguyen wrote: Binder's estimating equations are the usual way of applying weights to a Cox model, so nothing special is done apart from calling coxph(). To quote the author of the survival package, Terry Therneau, Other formulae change in the obvious way, eg, the

[R] Regarding the 'R' Load Command

Hi, I'm new to 'R' and need some help on the Load command. Any responses will be highly appreciated. Thanks in advance! As per manuals, the Load command expects a binary file input that is saved using a save command. However it is required that we need to call the 'R' program from Java

Re: [R] GUI commands to call for a protein from protein data bank

I tried doing that and this is what I go: dlg - aDialog(items=list( ProtienCode=stringItem() ), OK_handler=function(.) { # . is reference to dlg object values - .$to_R() f - function(ProtienCode) pdb - read.pdb(.$get_ProteinCode()) #cat(ProteinCode is,ProtienCode,\n) do.call(f, values) } )

Re: [R] looking for .. dec if vector if element x

Thank you Adrian, its working fine. Knut __ R-help@r-project.org 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, self-contained,

  1   2   >