[R] Help with assigning values to a row or column

Hello, I am attempting to write a function that assigns a value to a row or column within the matrix that it runs on. I am, however, having trouble accessing the row or column within this dynamic variable. I have looked through the archives (nothing there) and read the rules for posting to this

Re: [R] Package plm heterogenous slopes

Giovanni, Thank you for your reply. pvcm is indeed what I was looking for. I have a follow-up questio: plm documentation says that pvcm with method 'within' is equivalent of estimating a separate model for each individual. Let's assume that I would like to estimate a panel data model with an

Re: [R] Help with assigning values to a row or column

Hello, With a little broader searching, I did find a solution by Fernando Saldanha (http://tolstoy.newcastle.edu.au/R/help/05/04/3559.html): eval(parse(text = paste('testMatrix', '[1,] - 4', sep = ''))) Thanks anyway!! Best, Dave On Mon, Feb 8, 2010 at 12:19 AM, Dave Deriso

Re: [R] Out-of-sample prediction with VAR

Hello Peter, by judging from your code snippet: | ts_Y - ts(log_residuals[1:104]); # detrended sales data | ts_XGG - ts(salesmodeldata$gtrends_global[1:104]); | ts_XGL - ts(salesmodeldata$gtrends_local[1:104]); | training_matrix - data.frame(ts_Y, ts_XGG, ts_XGL); | |

Re: [R] metafor package: effect sizes are not fully independent

Hello Gang, Mike has already given you some excellent advice and references. I just want to add some information about the metafor package. To be precise, the meta-analytic mixed-effects model is given by: y_i = X_i * beta_i + u_i + e_i where X_i and beta_i have the usual interpretation, u_i

[R] Error on start R in server

Hello all, (Thank for your reply) I have a web-application in Apache Tomcat, when i start R in this application, I used packe RSJava Code ROmegahatInterpreter interp; String [] rargs = {--no-save}; REvaluator e; interp = new ROmegahatInterpreter (

Re: [R] Help with apply()

On 8/02/2010 2:08 PM, David Winsemius wrote: On Feb 7, 2010, at 8:26 PM, Nathan S. Watson-Haigh wrote: I have a 2 column data.frame: d[1:5,] a b 180015 C 280016 B 380023 C 480062 B 580069 B I want to apply a function across each

[R] using setMethod or setGeneric to change S4 accessor symbol from @ to $

I created some S4 objects that are essentially data frame objects. The S4 object definitions were necessary to verify data integrity and force a standardized data format. I am, however, finding myself redefining all the typical generic functions so that I can still manipulate my S4 objects as if

Re: [R] Help with apply()

On 8/02/2010 4:33 PM, Jim Lemon wrote: On 02/08/2010 12:26 PM, Nathan S. Watson-Haigh wrote: I have a 2 column data.frame: d[1:5,] a b 1 80015 C 2 80016 B 3 80023 C 4 80062 B 5 80069 B I want to apply a function across each row: for(i in 1:nrow(d)) { + myFun(con, d[i,]$a,

[R] [R-pkgs] Update of interval package

Hi all, I have uploaded an update of the interval package, which calculates nonparametric maximum likelihood estimates (NPMLE) of distributions and exact or asymptotic weighted logrank-type tests for interval censored data. An important part of the update is that now the default initial

[R] [R-pkgs] Package np update (0.30-6) adds nonparametric entropy test functionality...

Dear R users, Version 0.30-6 of the np package has been uploaded to CRAN. See http://cran.r-project.org/package=np Note that the cubature package is now required in addition to the boot package. The recent updates in 0.30-4 through 0.30-6 provides additional functionality in the form of five

[R] R installation require proxy setting

Could somebody please assist with the problem I have now. At our university certain changes have been made to users regarding accessing the internet. I have re-installed the (current) version of R, using custom settings (for Windows XP), choosing amongst other things Internet2 somewhere. Now,

Re: [R] Memory Problem

Hi, Can I use macro variables in R. If we can use macro variables in R, where i can get that programs or macro in R books. -- View this message in context: http://n4.nabble.com/Memory-Problem-tp1459740p1472700.html Sent from the R help mailing list archive at Nabble.com.

Re: [R] x-axis plot problem

Thanks guys, it works now properly when I change axas=i to axat=n, but there is a space between the box from the left side (y-axis) and the beginning plotted curve. please could you tell me how avoid that. any suggestions thanks, -- View this message in context:

[R] R: Package plm heterogenous slopes

Dear Otto, great. Re your follow-up question, it is not directly feasible in 'plm'; although I like the partitioned regression solution you cite, which should work. As an alternative suggestion, I would check the packages 'nlme' and 'lme4', which allow very flexible specifications. E.g., you

[R] Hi problem with read.csv

Hi all, I am getting problem with read.csv() command the problem is I have one file which is saved in ANSI while I am trying to open with R-Project its opend and all the cloumn names are comming proper like for the command names(kk) FILEID UNIQUEIDLATESTFLG ACTIVEFLG

[R] the hat ^ in regular expression

Dear UseRs, I'm trying to find variable names (string after the mydata$) in a expression. For example, myexpr - expression( ( mydata$variable1 / mydata$variable2 ) ^ 2 - 1 + 3 * 4 ) I would like to get variable1 and variable2. The following few lines split the original character string into

Re: [R] the hat ^ in regular expression

christophe dutang wrote: Dear UseRs, I'm trying to find variable names (string after the mydata$) in a expression. For example, myexpr - expression( ( mydata$variable1 / mydata$variable2 ) ^ 2 - 1 + 3 * 4 ) I would like to get variable1 and variable2. The following few lines split the

[R] Counting by rows based on multiple criteria

Dear all, I have a data frame of 6 columns and ~6 rows which I hope to perform the following calculation on. For each row, I wish to determine whether there are a greater number of positive or negative numbers. Then, if there are more positive numbers in the row, count how many occur -

Re: [R] Counting by rows based on multiple criteria

Will this do it for you: x - matrix(sample(c(-1,0,1), 6000, TRUE), ncol=6) # compute the occurances gt0 - rowSums(x 0) lt0 - rowSums(x 0) eq0 - rowSums(x == 0) x - cbind(x, result=ifelse(eq0 == 6, -, + ifelse(gt0 lt0, gt0, + ifelse(gt0 == lt0, 0, -lt0 head(x,20)

Re: [R] contour persp

On 07.02.2010 22:46, Andrew Wang wrote: I have this data set that both x y are ordered vectors of length 600 700 respectively; z is a 600 by 700 matrix whose entry z[i,j] is either a missing value (indicated by 'NaN') or a real number between 0 and 1. The contour function contour(x,y,z)

[R] About scales in graphics

Hi, I'm building two graphics just changing the scale to show the graphic illusion. The first graph would be 'correct'. atend = c(410,430,450,408,408,405) names(atend) = c('Janeiro', 'Fevereiro', 'Março', 'Abril', 'Maio', 'Junho') barplot(atend, las=1, xlab='Meses', ylab='Número de

[R] Physically open Excel file from R

Hello, everyone I wonder if it is possible to PHYSICALLY open an Excel file from R. The reason I ask is, I produce regularly an Excel file in R, and then I want to make it look good, so I have a VBA routine in another Excel file that works on the regular Excel file. This formatting file executes

Re: [R] the hat ^ in regular expression

Try this: library(gsubfn) myexpr - expression( ( mydata$variable1 / mydata$variable2 ) ^ 2 - 1 + 3 * 4) strapply(as.character(myexpr), mydata\\$(\\w+))[[1]] [1] variable1 variable2 See http://gsubfn.googlecode.com for more info on strapply. Another approach is setdiff(all.vars(myexpr),

Re: [R] Physically open Excel file from R

On 08.02.2010 13:48, Sergey Goriatchev wrote: Hello, everyone I wonder if it is possible to PHYSICALLY open an Excel file from R. The reason I ask is, I produce regularly an Excel file in R, and then I want to make it look good, so I have a VBA routine in another Excel file that works on the

Re: [R] Physically open Excel file from R

Dear Uwe, Worked as a charm! :-) Thank you kindly. Now, after I've run the VBA with shell.exec(), and closed the formater and the formatted file, I still have Excel shell (=Excel with no files open) open. Is it possible to close it from R as well? Best, Sergey 2010/2/8 Uwe Ligges

Re: [R] About scales in graphics

Silvano wrote: Hi, I'm building two graphics just changing the scale to show the graphic illusion. The first graph would be 'correct'. atend = c(410,430,450,408,408,405) names(atend) = c('Janeiro', 'Fevereiro', 'Março', 'Abril', 'Maio', 'Junho') barplot(atend, las=1, xlab='Meses',

Re: [R] Counting by rows based on multiple criteria

Spot on, thanks very much indeed. Steve _ Got a cool Hotmail story? Tell us now __ R-help@r-project.org mailing list

[R] PriorControl MCLUST

Hi I have used MCLUST to fit a mixture Gaussian model to Data and now I want to do the same for a subgroup of data with the mean and variance of previous step. I mean, I want to use MCLUST only with known means and variances and determine only probability. I have seen Prior Control But I do not

Re: [R] Physically open Excel file from R

Try shell('\file.xls\') where 'file.xls' is the excel filename. The escaped quotes (\) are not strictly necessary if the filename contains no spaces, but they are essential if it does. Sergey Goriatchev serg...@gmail.com 08/02/2010 12:48:20 Hello, everyone I wonder if it is possible to

[R] evolution of Nelder-Mead process

Dear list,   I am looking for an R-only implementation of a Nelder-Mead process that can find local maxima of a spatially distributed variable, e.g. height, on a spatial grid, and outputs the coordinates of the new point during each evaluation. I have found two previous threads about this

Re: [R] Memory Problem

What exactly is your definition of macro? What to you want to do? What is the problem that you are trying to solve? Why to you think macros will help? Typically R does not have macros; I assume that idea is a holdover from SAS. On Mon, Feb 8, 2010 at 4:30 AM, Meenakshi

[R] Fast way to determine number of lines in a file

Hi all, Is there a fast way to determine the number of lines in a file? I'm looking for something like count.lines analogous to count.fields. Hadley -- http://had.co.nz/ __ R-help@r-project.org mailing list

[R] objects masked from packages

dear all, I have a problem with a masked object in a package we created here. we make a package for a workflow of internal analysis of microarray data. to create the package we used: install.packages(pkgs=affyAnalysis, repos=NULL) R CMD INSTALL affyAnalysis Erzeuge Verzeichnisse ... Erzeuge

Re: [R] Fast way to determine number of lines in a file

Hi, parser::nlines does it in C. Romain On 02/08/2010 03:16 PM, Hadley Wickham wrote: Hi all, Is there a fast way to determine the number of lines in a file? I'm looking for something like count.lines analogous to count.fields. Hadley -- Romain Francois Professional R Enthusiast +33(0)

Re: [R] Fast way to determine number of lines in a file

Hadley Wickham hadley at rice.edu writes: Hi all, Is there a fast way to determine the number of lines in a file? I'm looking for something like count.lines analogous to count.fields. Hadley How about something like length(readLines(fname)) Ken

Re: [R] Hi problem with read.csv

dvkirankumar wrote: I am getting problem with read.csv() ... but if I am converting those files to UTF-8 formate and after loading to R-project Object and again I checked for coloumn names then its giving like ... ï..FILEID Add fileEncoding=UTF-8 Dieter -- View this

Re: [R] Fast way to determine number of lines in a file

If you are willing to use an external program parse the result of: system(wc -l small.dat) 10 small.dat On Windows there is a wc.exe program in the Rtools distribution. On Mon, Feb 8, 2010 at 9:16 AM, Hadley Wickham had...@rice.edu wrote: Hi all, Is there a fast way to determine the number

Re: [R] evolution of Nelder-Mead process

Juliane Struve wrote: Dear list, I am looking for an R-only implementation of a Nelder-Mead process that can find local maxima of a spatially distributed variable, e.g. height, on a spatial grid, and outputs the coordinates of the new point during each evaluation. I have found two previous

Re: [R] Fast way to determine number of lines in a file

parser::nlines does it in C. Looks promising, but I need something that uses connections because I'm working with big bzipped files. Hadley -- http://had.co.nz/ __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help

Re: [R] Fast way to determine number of lines in a file

Hi Ken, How about something like length(readLines(fname)) I'm trying to avoid the overhead of reading the file in twice. (I'm trying to preallocate a data structure for a chunked read) Hadley -- http://had.co.nz/ __ R-help@r-project.org mailing

Re: [R] Fast way to determine number of lines in a file

On 02/08/2010 04:16 PM, Hadley Wickham wrote: parser::nlines does it in C. Looks promising, but I need something that uses connections because I'm working with big bzipped files. Hadley Ah... the lack of c-level api for connections again ;-) -- Romain Francois Professional R Enthusiast

Re: [R] Fast way to determine number of lines in a file

I was looking for a fast line counter as well a while ago and ended up writing a small function in R: countLines() in the R.utils package At least at the time, it was faster than readLines() [for unknown reasons]. It is also more memory efficient. It supports connections. I don't think it

[R] optimized R-selection and R-replacement inside a matrix

Hello, I need to modify some huge arrays (2000 individuals x 50 000 variables). To format the data, I think I should benefit from optimized R-selection and R-replacement inside a matrix and prohibite a naive use of loops. Thank you in advance for providing information about the following

Re: [R] evolution of Nelder-Mead process

On 08-Feb-10 14:50:59, Paul Hiemstra wrote: Juliane Struve wrote: Dear list, I am looking for an R-only implementation of a Nelder-Mead process that can find local maxima of a spatially distributed variable, e.g. height, on a spatial grid, and outputs the coordinates of the new point

Re: [R] [BioC] objects masked from packages

Hi Assa, First off, please don't cross-post. This is really an R-help question, as it is about the R language rather than a specific BioC package. Assa Yeroslaviz wrote: dear all, I have a problem with a masked object in a package we created here. we make a package for a workflow of

[R] confidence interval for negatively skewed, leptokurtic sample

Hello, I´ve got a statistical problem that I hope you can help me with. It doesn´t have to do directly with R, so if there´s another forum which would suit better, please tell me! Now here´s the problem: I want to derive confidence intervals for a variable X, which is - given the descriptive

Re: [R] Data views (Re: (Another) Bates fortune?)

On Sun, Feb 7, 2010 at 2:40 PM, Emmanuel Charpentier charp...@bacbuc.dyndns.org wrote: Note : this post has been motivated more by the hierarchical data subject than the aside joke of Douglas Bates, but might be of interest to its respondents. Le vendredi 05 février 2010 à 21:56 +0100, Peter

Re: [R] using a variable name stored in another variable?

On Sun, 7 Feb 2010, Sharpie wrote: Chris Seidel wrote: Hello, I'm trying to figure out how to create a data object, and then save it with a user-defined name that is input as a command line argument. I know how to create the object and assign it the new name, however, I can't figure out

Re: [R] x-axis plot problem

Many many thanks.it is working now very well. again thanks a lot. -- View this message in context: http://n4.nabble.com/x-axis-plot-problem-tp1472286p1473006.html Sent from the R help mailing list archive at Nabble.com. __

Re: [R] Memory Problem

jim holtman jholt...@gmail.com 08/02/2010 14:09:52 Typically R does not have macros; I know exactly why Jim Holtman said that; R doesn't have a separate 'macro' construct with separate 'macro variables'. But it is perhaps a bit misleading to say that R doesn't have macros without saying a bit

[R] data frames; matching/merging

Hi all, I'm feeling a little guilty to ask this question, since I've written a solution using a rather clunky for loop that gets the job done. But I'm convinced there must be a faster (and probably more elegant) way to accomplish what I'm looking to do (perhaps using the merge function?). I

Re: [R] data frames; matching/merging

On Mon, Feb 8, 2010 at 11:39 AM, Jonathan jonsle...@gmail.com wrote: Hi all,    I'm feeling a little guilty to ask this question, since I've written a solution using a rather clunky for loop that gets the job done.  But I'm convinced there must be a faster (and probably more elegant) way to

Re: [R] data frames; matching/merging

Hi! I'm definitely not an expert in R (and it's my first reply!), but if I understand right, I think the aggregate function might do what you're looking for. Try ?aggregate to get more info. You might find what you need! HTH Ivan Le 2/8/2010 17:39, Jonathan a écrit : Hi all, I'm

Re: [R] data frames; matching/merging

x - read.table(textConnection(V1 V2 + 1a3 + 2a2 + 3b9 + 4c4 + 5a7 + 6b11), header=TRUE) closeAllConnections() # close; matrix with rownames - easy enough to change into a dataframe if you want cbind(tapply(x$V2,

Re: [R] data frames; matching/merging

On Feb 8, 2010, at 11:39 AM, Jonathan wrote: Hi all, I'm feeling a little guilty to ask this question, since I've written a solution using a rather clunky for loop that gets the job done. But I'm convinced there must be a faster (and probably more elegant) way to accomplish what I'm

Re: [R] data frames; matching/merging

You could try aggregate: If we call your data frame df: aggregate(df[2], by=df[1], FUN=min) will get you what you asked for (if not necessarily what you need ;-) ) Switching the columns around is easy enough if you need to; proceeding stepwise: df.new-aggregate(df[2], by=df[1], FUN=min)

[R] Color intervals in image function (image.plot)

Hi, The script below is my current coding in order to produce a contour plot of temperature across altitude and time. In my case, time,altitude and temperature are represented by x, y and z variables. ## Brazilan.Pallete -

[R] How to reverse the axis direction in log plot?

The following command gives me a plot where the axes are in log scale but 'x' increases in the right direction. plot(1:10,1:10, log='xy') I want to change the plot such that it is still in log scale but 'x' decreases rather increases in the right direction. I'm wondering if there is an option to

Re: [R] data frames; matching/merging

Here are 3 solutions assuming DF contains the data frame: # 1. aggregate aggregate(DF[2], DF[1], min) V1 V2 1 a 2 2 b 9 3 c 4 # 2. aggregate.formula - requires R 2.11.x aggregate(V2 ~ V1, DF, min) V1 V2 1 a 2 2 b 9 3 c 4 # 3. SQL using sqldf library(sqldf) sqldf(select

Re: [R] How to reverse the axis direction in log plot?

reverse the xlim; plot(1:10, 1:10, log='xy', xlim=rev(range(1:10))) On Mon, Feb 8, 2010 at 12:07 PM, bluesky...@gmail.com wrote: The following command gives me a plot where the axes are in log scale but 'x' increases in the right direction. plot(1:10,1:10, log='xy') I want to change the

Re: [R] simulat from dirichet process

There are a couple of packages with rdirichlet functions, including gtools and bayesm (and probably others). If these do not do what you want, give us some more detail of what you are trying to do. -- Gregory (Greg) L. Snow Ph.D. Statistical Data Center Intermountain Healthcare

Re: [R] How to reverse the axis direction in log plot?

plot(1:10,1:10,log='xy',xlim=c(10,1)) - Try http://prettygraph.com Pretty Graph , the easiest way to make R-powered graphs on the web. -- View this message in context: http://n4.nabble.com/How-to-reverse-the-axis-direction-in-log-plot-tp1473161p1473176.html Sent from the R help mailing

Re: [R] How to reverse the axis direction in log plot?

On Mon, Feb 8, 2010 at 11:13 AM, Hrishi Mittal hrishimit...@gmail.com wrote: plot(1:10,1:10,log='xy',xlim=c(10,1)) What if I want to use the automatically computed limits? I could use the following command. But it is not a very clean user interface, as the user has to compute the min and the

Re: [R] How to reverse the axis direction in log plot?

On Feb 8, 2010, at 12:30 PM, bluesky...@gmail.com wrote: On Mon, Feb 8, 2010 at 11:13 AM, Hrishi Mittal hrishimit...@gmail.com wrote: plot(1:10,1:10,log='xy',xlim=c(10,1)) What if I want to use the automatically computed limits? I could use the following command. But it is not a very

[R] R wiki link ?

Is this a transient problem, or has the link to the R wiki on the R home page (www.r-project.org) to http://wiki.r-project.org/ been corrupted? I can find http://rwiki.sciviews.org that works. JN __ R-help@r-project.org mailing list

[R] Color intervals in image.plot function

Hi, The script below is my current coding in order to produce a contour plot of temperature across altitude and time. In my case, time,altitude and temperature are represented by x, y and z variables. ## Brazilan.Pallete -

Re: [R] Color intervals in image function (image.plot)

On Feb 8, 2010, at 12:06 PM, FMH wrote: Hi, The script below is my current coding in order to produce a contour plot of temperature across altitude and time. In my case, time,altitude and temperature are represented by x, y and z variables.

[R] Follow-up Question: data frames; matching/merging

Wow.. thanks for the deluge of responses! Aggregate seems like the way to go here. But, suppose that instead of integers in column V2, I actually have dates (and instead of keeping the minimum integer, I want to keep the earliest date): df =

Re: [R] R installation require proxy setting

Look at ?download.file and scroll down to the section titled Setting Proxies. -- Gregory (Greg) L. Snow Ph.D. Statistical Data Center Intermountain Healthcare greg.s...@imail.org 801.408.8111 -Original Message- From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-

Re: [R] R wiki link ?

Prof. John C Nash wrote: Is this a transient problem, or has the link to the R wiki on the R home page (www.r-project.org) to http://wiki.r-project.org/ been corrupted? I can find http://rwiki.sciviews.org that works. Yes, the problem is known. I have to fix it. Best, Philippe Grosjean JN

Re: [R] Follow-up Question: data frames; matching/merging

It works if you use character class rather than factors for your date strings. Note last arg to data.frame: DF - data.frame(V1=c(1,1,2,3,1,2), + V2=c('2002-03-13','1989-03-10','1988-01-20','1997-05-15','1996-11-18','2000-01-12'), + stringsAsFactors = FALSE) aggregate(DF[2], DF[1], min) V1

Re: [R] metafor package: effect sizes are not fully independent

Thanks for the confirmation and pointer, Mike! Dr. Viechtbauer, I'm looking forward to the new functionality of specifying covariance structures in metafor(). Thanks both again for the great help, Gang On Sun, Feb 7, 2010 at 8:59 PM, Mike Cheung mikewlche...@gmail.com wrote: Dear Gang, It

[R] Selecting single TAZ based on area proportion

Good day all, I am having an issue coercing my data. Below i have some data on taxlots and an associated TAZ(transportation analsysi zone) that each property is within. The main issue is that some properties cross TAZ boundaries so i need to make a decision as to which TAZ the property

[R] Problem with R on USB-drive

Hi, I installed R on USB-drive, but when I run Rgui.exe from bin folder, I get this error: --- R version 2.10.1 (2009-12-14) Copyright (C) 2009 The R Foundation for Statistical Computing ISBN 3-900051-07-0 R is free software and comes with ABSOLUTELY NO WARRANTY. You are

Re: [R] using setMethod or setGeneric to change S4 accessor symbol from @ to $

On 02/07/2010 08:31 PM, Markus Weisner wrote: I created some S4 objects that are essentially data frame objects. The S4 object definitions were necessary to verify data integrity and force a standardized data format. I am, however, finding myself redefining all the typical generic functions

[R] Poisson and neg. bin. regression with random effects

Hi there, I have relative abundance data for 13 mammal species that I collected at various sites that ranged in road density. I'm trying to determine the effect of road density on animal abundance across body sizes. For most species, I have data that was collected in one year but

[R] constraining functions for optim()

Hi list! I am optimizing one function with two data sets where the unknown values (Vcmax and gi) are the same in each dataset (see script below). The script works, but I would like to add other constraints to this function. In this function, the optimization is carried out like: ### code and

Re: [R] Error on start R in server

On 02/08/2010 12:51 AM, Jiiindo wrote: Hello all, (Thank for your reply) I have a web-application in Apache Tomcat, when i start R in this application, I used packe RSJava Code ROmegahatInterpreter interp; String [] rargs = {--no-save}; REvaluator e; interp = new

[R] Wilcoxon signed-ranks test using package coin ?

Given the following data, and hypothesized median M.0 I've found a method to implement the Wilcoxon signed-rank test. Data: (with one zero difference and tied ranks) x - c(136, 103, 91, 122, 96, 145, 140, 138, 126, 120, 99, 125, 91,142, 119, 137) M.0 - 119 library(exactRankTests) Package

Re: [R] using a variable name stored in another variable?

Solved! Duncan Murdoch got it right: I think you want save(list=myobjectname, file= ...) assuming that the object has already been created with that name. If it hasn't, you'll need two steps: assign( myobjectname, value) save(list=myobjectname, file=...) [...] This works great. If

[R] Hypercube in R

Dear all, Does anybody have an idea or suggestion how to construct (plot) 4-dimensional hypercube in R. Thanks in advance for any pointers. Regards, Andrej __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do

Re: [R] constraining functions for optim()

It makes no sense to set your constraints as the solution that you want. I also don’t like your non-smooth objective function: abs(y1 + y2). It is also wrong. It should be abs(y1) + abs(y2), but even that is not preferable since it is non-smooth. What was wrong with the approach that I already

Re: [R] Contributed packages

Satish Vadlamani wrote: Folks: If you wanted to find out about what are the contributed packages and classify them, how would you go about it? For someone new like me, I would like to know what the possibilities are. When I click on install packages on my Windows version of R, it gives me

Re: [R] specifying colors in a heatmap/image -like plot

David and Jim, thanks for your help. Your advice was exactly what I needed. I tinkled with the col argument before but I was trying to assign colors to numbers with a command like: col(1-green,2-blue) Usually I need to see an example before I can implement something correctly. Thanks

Re: [R] constraining functions for optim()

Hi Ravi, It makes no sense to set your constraints as the solution that you want. The thing is that y.1 and y.2 are too far away from the measured values. These values (photosynthesis assimilation rates) should not be as far away from each other. With the data I sent where you proposed a very

Re: [R] using setMethod or setGeneric to change S4 accessor symbol from @ to $

Worked like a charm!! Thank you so much. I just plugged the following into my code ... setMethod($, CADresponses, function(x, name) slot(x, name)) ... and it worked perfect. If you don't mind, I have a quick follow up question, using your example setClass(A, representation(a=numeric,

Re: [R] using setMethod or setGeneric to change S4 accessor symbol from @ to $

On 02/08/2010 01:22 PM, Markus Weisner wrote: Worked like a charm!! Thank you so much. I just plugged the following into my code ... setMethod($, CADresponses, function(x, name) slot(x, name)) ... and it worked perfect. If you don't mind, I have a quick follow up question, using your

[R] Mixed logit models with a random coefficient

Hi All, Sorry to bother you. I'm trying to estimate a set of discrete choice data in R with mixed logit models where one coefficient is random and normally distributed. I've searched on the R help archive and don't see much information very specific to what I'm doing, so I write the code

[R] Invitation to connect on LinkedIn

LinkedIn Scott Hyde requested to add you as a connection on LinkedIn: -- Arnaud, I'd like to add you to my professional network on LinkedIn. - Scott Accept invitation from Scott Hyde

[R] interp function in package akima

hello! does anyone know how to get values out of an interp object (from akima package)? ive constructed an elevation map and would like to numerically see the calculated z values of data points, which were not empirically sampled (basically: id like to see the altitude at locations x,y which

[R] mle.stepwise

Hello, I have a basic question on mle.stepwise, which seems to return strange answers, even on a basic example. I posted this last May but never got an answer. I resend it here hoping it will this time. # BEGIN require(wle) x1 = runif(100) x2 = runif(100) y = 3 + 2*x1 + rnorm(100)

[R] What is the equivalent function in R to xyplot in S?

Page 140 of MASS uses the function xyplot. But I don't find it in R. Is there a package that I should load to use xyplot. Or there is a function with a different name in R that does the same thing as xyplot in S. xyplot(Gas ~ Temp | Insul, whiteside, panel = function(x, y, ...) { panel.xyplot(x,

Re: [R] What is the equivalent function in R to xyplot in S?

Tena koe library(lattice) ?xyplot HTH . Peter Alspach -Original Message- From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org] On Behalf Of bluesky...@gmail.com Sent: Tuesday, 9 February 2010 11:31 a.m. To: r-h...@stat.math.ethz.ch Subject: [R] What is

Re: [R] What is the equivalent function in R to xyplot in S?

Try help.search(xyplot). If nothing comes up, try RSiteSearch(xyplot). -Ista On Mon, Feb 8, 2010 at 10:31 PM, bluesky...@gmail.com wrote: Page 140 of MASS uses the function xyplot. But I don't find it in R. Is there a package that I should load to use xyplot. Or there is a function with a

Re: [R] interp function in package akima

On Feb 8, 2010, at 5:06 PM, Sacha Viquerat wrote: hello! does anyone know how to get values out of an interp object (from akima package)? ive constructed an elevation map and would like to numerically see the calculated z values of data points, which were not empirically sampled

Re: [R] mle.stepwise

eariasca eariasca at math.ucsd.edu writes: Hello, I have a basic question on mle.stepwise, which seems to return strange answers, even on a basic example. I posted this last May but never got an answer. I resend it here hoping it will this time. # BEGIN [snip] # END I am

[R] Dividing one column of form xx-yy into two columns, xx and yy

I have a data set where one column consists of two numerical factors, separated by a -. So my data looks something like this: 43-156 43-43 1267-18 . . . There are additional columns consisting of single factors as well, so reading the csv file (where the data is stored) with the sep=- addition

[R] object 'xxx' not found

The following line of code seems fairly straight forward, yet it is kicking back an error message: for (i in 1:length(mean.natveg.frac)){month.observed[i]=as.numeric(names(mean.natveg.frac[i]))%%12} Error message: Error in month.observed[i] = as numeric(names(mean.natveg.frac[i]))%%12 : object

Re: [R] using setMethod or setGeneric to change S4 accessor symbol from @ to $

On 02/08/2010 02:54 PM, Markus Weisner wrote: Thanks. Used getGeneric([) to figure out the general format for the setMethod, but am having some problem with how to set up the actual function: getGeneric([) standardGeneric for [ defined from package base function (x, i, j, ..., drop =

  1   2   >