Re: [R] Looking up the directory a file is located in

On Thu, Nov 18, 2010 at 12:00 AM, Cliff Clive wrote: > > Basically I'm just looking for a command that can look up the name of the > directory of the script that is running.  If I move or copy the script to > another directory, it should be able to read the name of the new directory > without me h

Re: [R] stacking consecutive columns

On Nov 17, 2010, at 16:37 , Graves, Gregory wrote: > Follows is the exact solution to this: > > v <- NULL > > #note that decreasing is FALSE so preceding year preceeds > > for(i in 2:46) { > kk <- melt(yearmonth[, c(1, i, i+1)], id.vars="month", variable_name="year") > v[[i-1]] <- kk[order(k

[R] filling large matrix( or dist object) with values from small matrix.

The script below provides a trivial example of filling a large matrix with values from a small one with corresponding row and column details. this is fine for a small example but too slow for many large examples ( my small matrices are 65x65 and my large matrices up to 900x900, and I have many

[R] [R-pkgs] The *makesweave* Package for Using Make with Sweave Efficiently

Hi R Community, I've just completed the initial version of (Linux-based) R package, called *makesweave*, for using Make to build Sweave documents efficiently. The idea is that R is started once (per shell) as a background process, then Make builds each Sweave source file, using the same R instanc

[R] [R-pkgs] lubridate v2.2 available on cran

Version 2.2 of the lubridate package is now available. lubridate makes it easier to work with date-time data. For example, it provides: * simple functions to extract and modify components of a date-time, such as years, months, days, hours, minutes, and seconds: year(), month(), day(), ... * quick

Re: [R] How to catch warnings

Dear Michael, I would like to thank you for your reply. This line made things easier. One more question what If I want to halt or pause the program when a warning happens? Right now I get only a message printed but it would be nicer if the execution is paused so to try to print more values. Bes

Re: [R] How to catch warnings

Sorry Alex, I don't quite follow what you want. What do you mean by "it would be nicer if the execution is paused so to try to print more values" ? Michael On 18 November 2010 21:11, Alaios wrote: > > Dear Michael, > I would like to thank you for your reply. This line made things easier. One >

[R] Why some times R-halts?

Hello everyone. I have written a small function that returns a 3*n array where n is typicall between 5-55. I issued in R shweights<-shadowing_weights(x,xr) and I got the result in less that one sec. then I tried to see what is inside the shweights array and thus I issued the edit(shweights) which

Re: [R] plot vs print ??

Hello everybody. My question arised from the output of lattice's histogram. But might be extended to any other object that could be printed. I think I've understood your answers, print calls the plot function when the object to be printed is a trellis plot object. I guess I can always use plot

Re: [R] rgl, multiple graphics ??

Hello I just get one plot How can I concoct several rgl plots? -- View this message in context: http://r.789695.n4.nabble.com/rgl-multiple-graphics-tp3047006p3048434.html Sent from the R help mailing list archive at Nabble.com. __ R-help@r-project.

[R] why help ignores some values

Please take a look at the values below . [31,] 30.000 30.000 [32,] 31.000 31.000 [33,] 32.000 32.000 [34,] 33.000 33.000 [35,] 34.000 34.000 [36,] 35.0

Re: [R] why help ignores some values

Possibly because R has no idea what the 35.872-th element of a matrix is? However you are coming up with these elements, you are running into a floating-point arithmetic problem (FAQ 7.31). Using round() might help. Sarah On Thu, Nov 18, 2010 at 6:19 AM, Alaios wrote: > Please take

Re: [R] How to catch warnings

Like this ? f <- function() { # on first warning will print to console and return results so far out <- c(1,2,3,4) extras <- list(good=c(2,3,4,5), bad=c(9,9,9), ugly=c(9,9,9,9,9)) ok <- TRUE for (ex in extras) { tryCatch(out <- cbind(out, ex), warning=function(w) { ca

Re: [R] How to catch warnings

-BEGIN PGP SIGNED MESSAGE- Hash: SHA1 On 11/18/2010 11:11 AM, Alaios wrote: > Dear Michael, > I would like to thank you for your reply. This line made things easier. One > more question what If I want to halt or pause the program when a warning > happens? Right now I get only a message p

[R] classification question

Hi. Working with a data set like: "age", "demographic data (n fields)", "interests(n fields)" has performed X actions for event Y. I want to ask how likely is it for another person with his/her age, demographic data and interests to perform actions for that event. My query set might be part

Re: [R] where are my pspline knots?

Hello: I hope that someone more knowledgeable will confirm or correct what I'm about to say. I don't have time now to check this by studying the "pspline" code. From reading the pspline{survival} help page, I believe it uses standard B-splines, "penalising the integrated secon

Re: [R] why help ignores some values

Just to update this f[35.99,36.22] works fine  and returns the data from f[35,36] it does this rounding the R itself. the problem finally was a value f[0.99,3] that was returning numeric(0) which I understand as a matrix doesnot have  any f[0,0] Regards Alex --- On Thu, 11/18/10, Sarah Gosle

Re: [R] How to catch warnings

Thanks a lot. Unfortunately browser() works pretty bad for me. As I am using Rkward this gives me a prompt while at the same time R console is blocked which make things not very easy. Regards Alex --- On Thu, 11/18/10, Rainer M Krug wrote: From: Rainer M Krug Subject: Re: [R] How to catch wa

Re: [R] rgl, multiple graphics ??

On 18/11/2010 5:55 AM, skan wrote: Hello I just get one plot How can I concoct several rgl plots? It's basically up to you to handle everything. demo(stereo) shows you one way. Duncan Murdoch __ R-help@r-project.org mailing list https://stat.

Re: [R] Updata Rdata File

On 18/11/2010 1:49 AM, Jason Kwok wrote: How do I add data to a .rdata file? In my case, I have a time series that needs to get updated every day. Load it, update the variables, save them. Duncan Murdoch __ R-help@r-project.org mailing list https:/

Re: [R] How to catch warnings

-BEGIN PGP SIGNED MESSAGE- Hash: SHA1 On 11/18/2010 12:54 PM, Alaios wrote: > Thanks a lot. Unfortunately browser() works pretty bad for me. As I am > using Rkward this gives me a prompt while at the same time R console is > blocked which make things not very easy. First option: you know

[R] Odp: Why some times R-halts?

Hi r-help-boun...@r-project.org napsal dne 18.11.2010 11:29:28: > Alaios > Odeslal: r-help-boun...@r-project.org > > > Hello everyone. > I have written a small function that returns a 3*n array where n is typicall > between 5-55. > I issued in R > shweights<-shadowing_weights(x,xr) > and I

Re: [R] Give me all operator

works nice : lapply( Mystruct, "[", "xy")  # should return a list with "xy" elements Thanks --- On Wed, 11/17/10, David Winsemius wrote: From: David Winsemius Subject: Re: [R] Give me all operator To: "Alaios" Cc: r-help@r-project.org, ivan.calan...@uni-hamburg.de Date: Wednesday, November 17

Re: [R] randomForest parameters for image classification

1. Memory issue: You may want to try to increase nodesize (e.g., to 5, 11, or even 21) and see if that degrades performance. If not, you should be able to grow more trees with the larger nodesize. Another option is to use the sampsize argument to have randomForest() do the random subsampling for

[R] Help me with R plotting

Hello everyone. I want some help with plots. I have some robots in an area. Every robot is placed on x,y coordinates and every robot has a power consumption of some watts. I would like to show where are the robots by showing in a map dots (where every dots is the x,y coordinate). Below that dot

Re: [R] Help me with R plotting

?plot Will give a good solution #Example: set.seed(5) xx <- data.frame(x = runif(50), y = runif(50), z = sample(c(1:3), 50, T)) with(xx, plot(x, y, col = z, pch = 19)) Now you can start going into how to play with the colors, and how to add ?legend to the image... Tal Contact

[R] sweep by levels of a factor

Hi, I'd appreciate help with this. I have a data matrix with one column, called f in the example below, a factor. I'd like to subtract the means from each of other columns for each level of the factor. That is, in the example, to go from the first matrix below to the second. I know SWEEP will tak

Re: [R] format secondary axis for dates

Thanks for your help, David! Somehow I just could not find this information in the documentation! Jannis David Winsemius schrieb: On Nov 16, 2010, at 2:32 PM, Jannis wrote: Dear List, this may be a Newbi question and may have been asked several times, but i am too stupid to find the po

[R] Mixed multinomial logit model (mlogit script)

Dear all, I am trying to run a mixed multinomial logit model in R since my response variable has 4 non-ordinal categories. I am using the package mlogit that estimates the parameters by maximum likelihood methods. First of all, I prepared my data using the mlogit.data command. In the mlogit com

Re: [R] where are my pspline knots?

The pspline function uses P-splines (Eilers and Marx, Statistical Science, 1981), which are a spline basis using a regular set of knots. Looking at the code for pspline, which isn't so hard, let dx = (max(x) - min(x))/ ntermwhere nterm is round(2.5 * desired degrees of freedom) The

Re: [R] New Sampling question

This can of course be done, but before I make any attempt to do it I have to ask: why do you want this? On Wed, Nov 17, 2010 at 7:08 PM, wangwallace wrote: > > I have another question about drawing samples from a data frame. This might > sound really tricky. Let me use a data frame I have posted

Re: [R] plot does not work

Dear jim, just to update this. My problem is now fixed. Regards Alex --- On Mon, 10/25/10, jim holtman wrote: From: jim holtman Subject: Re: [R] plot does not work To: "Alaios" Cc: "Rhelp" Date: Monday, October 25, 2010, 12:59 PM try 'graphics.off()' to close any device that might be open a

[R] PS

Note that there is NOT an intercept term in my pspline basis. One of the features of psplines is that \sum beta_i f_i(x) (where f_i are the spline basis functions) is linear if and only if the coefficients beta are a linear sequence. This makes it easy to decompose the fit into linear and non-li

Re: [R] sweep by levels of a factor

Hi: d$r <- with(d, ave(x, f, FUN = function(u) u - mean(u))) > d f x r 1 1 2 1 2 1 0 -1 3 2 0 -2 4 2 4 2 5 2 2 0 HTH, Dennis On Thu, Nov 18, 2010 at 6:02 AM, Lancaster, Anthony < anthony_lancas...@brown.edu> wrote: > Hi, > I'd appreciate help with this. I have a data matrix with one colu

[R] problems subsetting

Dear all, I have searched the forums for an answer - and there is plenty of questions along the same line - but none of the paproaches shown worked to my problem: I have a data frame that I get from a csv: summarystats<-as.data.frame(read.csv(file=f_summary)); where I have the columns Datase

Re: [R] problems subsetting

Hi, I got a bit lost with your explanation for your second problem. A reproducible example would DEFINITELY help to understand what you have and what you're trying to get. For your first problem subset1 <- summarystats[summarystats$Class == 1 & summarystats$Type == 1 & summarystats$Category

[R] Plotting number of patients at risk below survival curve

Dear list, does anyone know of a R-package that has implemented the increasingly popular inclusion of the number of patients at risk below Kaplan-Meier curves like in http://bloodjournal.hematologylibrary.org/content/vol116/issue19/images/large/zh89991058760001.jpeg any hint (or negative answer)

Re: [R] Help me with R plotting

Thanks a lot for your help I tried first this test<-(myvalues)-mean(myvalues) so to put the values close to the 0...30 scale... then I passed this as an argument plot(x,y,col=test) which plots the places where the robots are with a small color. I think this is ok for now. What I want next is a

Re: [R] problems subsetting

Hi Gerrit, indeed, that works. Excellent tip! For reference, I did this: subset1<-subset(summarystats,(Type==1)&(Class==1)&(Category==1)) I am still not totally sure when one uses "&" amd when "&&" - I was under the impression that && stands for logical AND Thanks a lot. Martin On 11

[R] Exporting Dataframe to Excel spreadsheet

Hi, I have a  data frame contains 30,000 rows and 105 columns.  How can i extract each of the columns as separate  Excel spreadsheet of the same excel file.  Regards, Pankaj Barah Department of Biology, Norwegian University of Science & Technology (NTNU) Realfagbygget, N-7491 Trondheim

Re: [R] Updata Rdata File

Thanks for the response Duncan. I'm sorry I wasn't clear in my question. I need to add an additional value every day to the data and not update it. I know how to load and save but I dont' know how to add data. load(file="C:\\datafile.rdata") save(data,file="C:\\datafile.rdata") Thanks. Jason

Re: [R] Help me with R plotting

Hi, Try this, robots <- data.frame(id=letters[1:20], x=rnorm(20), y=rnorm(20), consumption=runif(20, 10, 100)) library(ggplot2) ggplot(robots) + geom_point(aes(x, y, colour=cut(consumption, c(0, 30, 50, 100 + geom_text(aes

Re: [R] Exporting Dataframe to Excel spreadsheet

On Thu, Nov 18, 2010 at 10:26 AM, pankaj borah wrote: > I have a  data frame contains 30,000 rows and 105 columns.  How can i > extract each of the columns as separate  Excel spreadsheet of the same > excel file. See: http://rwiki.sciviews.org/doku.php?id=tips:data-io:ms_windows -- Statistics

Re: [R] Help me with R plotting

(Assuming you aren't going into ggplot2) You can have a look here, it might offer a nice solution: http://rgm2.lab.nig.ac.jp/RGM2/R_man-2.9.0/library/shape/man/colorlegend.html Also, you can use ?cut to change your coloring to specific chunks (although having a full gradient is probably nicer)

Re: [R] ncdf4 for Windows/R-2.12.0

Hi David, This link looks like it points to the resources you need: http://www.unidata.ucar.edu/software/netcdf/docs/faq.html#windows_netcdf4 Perry -- View this message in context: http://r.789695.n4.nabble.com/ncdf4-for-Windows-R-2-12-0-tp3047807p3048935.html Sent from the R help mailing li

Re: [R] problems subsetting

Hi, On Thu, Nov 18, 2010 at 10:25 AM, Martin Tomko wrote: > Hi Gerrit, > indeed, that works. Excellent tip! > > For reference, I did this: > > subset1<-subset(summarystats,(Type==1)&(Class==1)&(Category==1)) > > I am still not totally sure when one uses "&" amd when "&&"  - I was under > the impr

Re: [R] problems subsetting

Thanks Steve, that explains it... Unfortunately, I did nto get that from my R docs... I am still searching for a solution for matching the entries in my matrix by matching the rownames to the entires in a subset I got using Gerrit's method. Any idea? Thanks Martin On 11/18/2010 4:35 PM, Ste

Re: [R] problems subsetting

On 2010-11-18 07:25, Martin Tomko wrote: Hi Gerrit, indeed, that works. Excellent tip! For reference, I did this: subset1<-subset(summarystats,(Type==1)&(Class==1)&(Category==1)) I am still not totally sure when one uses "&" amd when"&&" - I was under the impression that&& stands for logical

Re: [R] problems subsetting

On Nov 18, 2010, at 10:25 AM, Martin Tomko wrote: Hi Gerrit, indeed, that works. Excellent tip! For reference, I did this: subset1<-subset(summarystats,(Type==1)&(Class==1)&(Category==1)) I am still not totally sure when one uses "&" amd when "&&" - I was under the impression that && stan

Re: [R] Looking up the directory a file is located in

So it sounds like the best we can do in R is to keep track of the script in a sort of master file that runs the script, and set the working directory in the master. Is that accurate? In Python any time you run a script, there is a built-in "__file__" variable that can tell you the file name of t

Re: [R] New Sampling question

Dear Ista Zahn-2, If you can give me some advice, I really appreciate it. I have been working on it for days. it seems hard for some novice of R like me to write flexible functions myself. This is for my dissertation. CSE and WSE are two scales of the same construct. The sampling strategy I wan

Re: [R] Help me with R plotting

Thanks a lot :) nice one --- On Thu, 11/18/10, baptiste Auguié wrote: From: baptiste Auguié Subject: Re: [R] Help me with R plotting To: "Alaios" Cc: "Rhelp" Date: Thursday, November 18, 2010, 3:34 PM Hi, Try this, robots <- data.frame(id=letters[1:20],                      x=rnorm(20),

Re: [R] problems subsetting

On Nov 18, 2010, at 10:42 AM, Martin Tomko wrote: Thanks Steve, that explains it... Unfortunately, I did nto get that from my R docs... I am still searching for a solution for matching the entries in my matrix by matching the rownames to the entires in a subset I got using Gerrit's metho

Re: [R] Looking up the directory a file is located in

On Thu, Nov 18, 2010 at 10:26 AM, Cliff Clive wrote: > In Python any time you run a script, there is a built-in "__file__" > variable > that can tell you the file name of the script itself. It would be nice to > have a feature like this in R. > R has that feature. source("c:/myfullpath/myfile.r

Re: [R] Looking up the directory a file is located in

On Thu, Nov 18, 2010 at 3:26 PM, Cliff Clive wrote: > > So it sounds like the best we can do in R is to keep track of the script in a > sort of master file that runs the script, and set the working directory in > the master.  Is that accurate? Errr no. Maybe. What? Your script is somewhere. So

Re: [R] Looking up the directory a file is located in

On Wed, Nov 17, 2010 at 4:08 PM, Cliff Clive wrote: > > Hello everyone, > > This should be an easy question, I think. > > I'd like to write a command in a program to set the working directory to > whatever directory the file is currently stored in.  Suppose I have a file > called "myRscript.r", an

[R] how to find near neighbors?

I am looking for an efficient way to find near neighbors... More specifically... I have two sets of points: A & B and I want to find points in set B which are closer to set A than some cutoff (or n-closest) I will appreciate very much any pointers... Ryszard

Re: [R] aspect ratio 1 and blank space

Look at the squishplot function in the TeachingDemos package, it may do what you want, or if not quite then you could perhaps tweak the code to include the values you want and create the correct aspect ratio. -- Gregory (Greg) L. Snow Ph.D. Statistical Data Center Intermountain Healthcare greg.

[R] conditional mean between two data frames with different levels

Hi guys, I have two data frames: one referred to 2008 and one to 2009. Their structure is identical except for the different data in them. I need to create a vector alfa of the same length of the dataframe 2009 and fill each element with the mean of 2008$var1 conditional to the subgroup indicated

[R] genralized linear regression - function glm - number of

Hello, Performing a linear regression through the function glm ("yi ~ X$V1 + X$V2 + X$V3 + X$V4 + X$V5 + X$V6 + X$V7 + X$V8 + X$V9 + X$V10"), I then edit the information about the coefficients: print(coefficients(summary(fit))) I note that the number of coefficients (7) is lower than the num

Re: [R] New Sampling question

Hi Wallace, Have you tried playing with sample()? Note that you can apply this function both to whole dataframes, as well as specific items within a vector. If you play with applying the function to different ways of indexing your sample data, you will likely arrive at your solution. for example:

Re: [R] New Sampling question

Also, I need some function at the end which would enable me to draw 1000 such random samples. thanks! :) -- View this message in context: http://r.789695.n4.nabble.com/New-Sampling-question-tp3047885p3048958.html Sent from the R help mailing list archive at Nabble.com. _

Re: [R] Odp: Sampling problem

I tried. yours works too. thanks a bounch!! Man -- View this message in context: http://r.789695.n4.nabble.com/Sampling-problem-tp3043804p3049013.html Sent from the R help mailing list archive at Nabble.com. __ R-help@r-project.org mailing list https:

Re: [R] how to find near neighbors?

Have a look at nncross() from the spatstat package. Best regards, Thierry ir. Thierry Onkelinx Instituut voor natuur- en bosonderzoek team Biometrie & Kwaliteitszorg Gaverstraat 4 9500 Geraardsbergen Belgium Research

Re: [R] Updata Rdata File

?rbind "Jason Kwok" wrote: >Thanks for the response Duncan. > >I'm sorry I wasn't clear in my question. I need to add an additional >value >every day to the data and not update it. > >I know how to load and save but I dont' know how to add data. >load(file="C:\\datafile.rdata") >save(data,file=

Re: [R] Updata Rdata File

Hi Jason, You can add as many objects as you want to save(). For instance: save(data1, data2, data3, data4, file = "whatever"). Once you have loaded the data, just add whatever data you want, and then save it (or edit your original object and resave it). Here is another example: > ls() charact

Re: [R] genralized linear regression - function glm - number of

On Nov 18, 2010, at 11:00 AM, Christine SINOQUET wrote: Hello, Performing a linear regression through the function glm ("yi ~ X$V1 + X$V2 + X$V3 + X$V4 + X$V5 + X$V6 + X$V7 + X$V8 + X$V9 + X$V10"), I then edit the information about the coefficients: print(coefficients(summary(fit))) I n

Re: [R] problems subsetting

Hi, On Thu, Nov 18, 2010 at 10:42 AM, Martin Tomko wrote: > Thanks Steve, that explains it... Unfortunately, I did nto get that from my > R docs... > > I am still searching for a solution for matching the entries in my matrix by > matching the rownames to the entires in a subset I got using Gerri

Re: [R] New Sampling question

You could try writing a loop a<-data.frame(c(1:10),c(21:30)) M<-10 #number of iterations- scale up to 1000 once you get your sampling function working res<-NULL #place to store your results for i in (1:M) { ares<-sample(a[,2],1) res<-c(res, ares) } res It's up to you how to

Re: [R] conditional mean between two data frames with different levels

Hi Alberto, It would help if you could provide a small example. I might break the problem down into three parts: 1) create a vector that has the final subgroupings you want 2) find the conditional means by subgroup 3) replicate the means as needed. My first guess would be start with: "==" or

Re: [R] stacking consecutive columns

On Wed, Nov 17, 2010 at 10:37 AM, Graves, Gregory wrote: > Follows is the exact solution to this: > > v <- NULL > > #note that decreasing is FALSE so preceding year preceeds > > for(i in 2:46) { >  kk <- melt(yearmonth[, c(1, i, i+1)], id.vars="month", variable_name="year") >  v[[i-1]] <- kk[order

Re: [R] how exactly does 'identify' work?

I think that your problem comes from a misunderstanding. The general rule is that you give the plot command 2 vectors, x and y (though you can give it the vectors separately, or together in a list or matrix). If you give plot only a single vector then it will use this as the y vector and use t

Re: [R] problems subsetting

Hi Martin, I think Steve meant dput(df), as I already told you Ivan Le 11/18/2010 17:40, Steve Lianoglou a écrit : Hi, On Thu, Nov 18, 2010 at 10:42 AM, Martin Tomko wrote: Thanks Steve, that explains it... Unfortunately, I did nto get that from my R docs... I am still searching for a sol

Re: [R] kalman filter in sspir

To learn why sspir does not have a filter function, you need to ask the package maintainer, Claus Dethlefsen . My belief is that he, Soren Lundbye-Christensen and Anette Luther Christensen found other outlets for their time since they completed the package and the companion paper, "C. D

[R] dmultinomial

Hello All, I'm trying to run a maximum likelihood analysis using dmultinomial (i'm avoiding dmultinom as I'd like to run it with vectors for the ML stuff). However, I'm having a hard time getting even the simplest example running. Any help would be greatly appreciated. > library(mc2d) > dmultin

Re: [R] problems subsetting

On Nov 18, 2010, at 10:42 AM, David Winsemius wrote: On Nov 18, 2010, at 10:25 AM, Martin Tomko wrote: Hi Gerrit, indeed, that works. Excellent tip! For reference, I did this: subset1<-subset(summarystats,(Type==1)&(Class==1)&(Category==1)) I am still not totally sure when one uses "&" am

[R] lme Random Effects and Covariates

1. I'm attempting to test for Random Effects. I've grouped the data on subject (grid) but want to use lme to build the model without subject as a RE then add it and do anova between the 2 models. This is the result I get and it appears it's adding Random Effects. tmp.dat4 <- groupedData(Trials

[R] Logistic regression with factorial effect

Hello, I’d like to evaluate the temporal effect on the relationship between a continuous variable (e.g. size) and the probability of mate success. Initially I was trying to do a logistic regression model incorporating the temporal effect, but I don’t know if that is the best option. I simulated s

[R] predict() an rpart() model: how to ignore missing levels in a factor

I am using an algorigm to split my data set into two random sections repeatedly and constuct a model using rpart() on one, test on the other and average out the results. One of my variables is a factor(crop) where each crop type has a code. Some crop types occur infrequently or singly. when the d

[R] Sample covariance matrix in R

Hello everyone. I would like to find the sample covariance matrix using R. So far I read on the wikipedia what a sample_covariance is http://en.wikipedia.org/wiki/Sample_covariance according to wikipedia one vector is enough to calculate the sample covariance matrix. In R I tried cov(myvector) a

[R] Help with lmer, nested data and repeated measures

Dear all, I'm discovering the somehow confusing (for me) world of linear mixed models after having been advised it could be the best option to analyze my dataset. After several days of reading, I'm not sure that what I ended up with makes some sense and I'd greatly appreciate any help and expla

Re: [R] Help with lmer, nested data and repeated measures

You would be better off posting to R-sig-mixed-models or R-sig-ecology -- Bert On Thu, Nov 18, 2010 at 9:58 AM, Simon Garnier wrote: > Dear all, > > I'm discovering the somehow confusing (for me) world of linear mixed models > after having been advised it could be the best option to analyze my

Re: [R] Logistic regression with factorial effect

You would be better off posting to R-sig-mixed-models or R-sig-ecology -- Bert On Thu, Nov 18, 2010 at 9:32 AM, Billy.Requena wrote: > > Hello, > > I’d like to evaluate the temporal effect on the relationship between a > continuous variable (e.g. size) and the probability of mate success. > Init

Re: [R] sweep by levels of a factor

Thank you very much, Tony L. On 18 November 2010 14:02, Lancaster, Anthony wrote: > Hi, > I'd appreciate help with this. I have a data matrix with one column, called > f in the example below, a factor. I'd like to subtract the means from each > of > other columns for each level of the factor.

Re: [R] conditional mean between two data frames with different levels

Thank you very much Josh, I guess you`re right. So this is an example: data frame 1 has 2 columns and 10 rows. The first column is "score" a variable indicating the number of goals scored by a football team score<-c(1,2,0,2,1,1,3,2,1,0), column 2 contains the football "teams " where teams<-c(a,b,c

Re: [R] New Sampling question

Hi, Mike, thank you very much!!:) the following two functions are really helpful, which I didn't even know. Actually, I searched the forum for something like this, but failed. Now I am still trying to make up my own functions. :) sample(a[,2],2) #randomly draw two numbers from a column sample(

Re: [R] Looking up the directory a file is located in

Thanks, Gabor! So far I like this one best: https://stat.ethz.ch/pipermail/r-help/2005-November/082347.html So if my script is called "myRscript.r", I can do the following: this.file = parent.frame(2)$ofile this.dir = gsub("/myRscript.r", "", this.file) setwd(this.dir) This will set the work

[R] RowSums Question

I have a question on RowSums. Lets say i have a timeSeries table A B C 1/1/90 NA 1 1 1/2/90 NA 1 1 1/3/90 NA 1 1 1/4/90 NA 1 1 1/5/901 1 1 1/6/901 1 1 if i use RowSums, i will get 1/5/903 1/

Re: [R] sweep by levels of a factor

Try this also: x$x - ave(x$x, x$f) On Thu, Nov 18, 2010 at 12:02 PM, Lancaster, Anthony < anthony_lancas...@brown.edu> wrote: > Hi, > I'd appreciate help with this. I have a data matrix with one column, called > f in the example below, a factor. I'd like to subtract the means from each > of >

Re: [R] Sample covariance matrix in R

Alex: ?cov > -Original Message- > From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org] On > Behalf Of Alaios > Sent: Thursday, November 18, 2010 12:54 PM > To: Rhelp > Subject: [R] Sample covariance matrix in R > > Hello everyone. > I would like to find the sample cov

Re: [R] RowSums Question

Try this: rowSums(tsObj, na.rm = TRUE) On Thu, Nov 18, 2010 at 3:58 PM, cameron wrote: > > > I have a question on RowSums. > > Lets say i have a timeSeries table >A B C > 1/1/90 NA 1 1 > 1/2/90 NA 1 1 > 1/3/90 NA 1 1 > 1/4/90 NA 1

Re: [R] Sample covariance matrix in R

The help and its examples are very comprehensive here. The usage you cite shows exactly what you need to do From: Alaios [mailto:ala...@yahoo.com] Sent: Thursday, November 18, 2010 1:30 PM To: Doran, Harold Subject: RE: [R] Sample covariance matrix in R Checked that Usage covr(x, y = NULL, na.r

Re: [R] Looking up the directory a file is located in

On 18/11/2010 11:56 AM, Cliff Clive wrote: Thanks, Gabor! So far I like this one best: https://stat.ethz.ch/pipermail/r-help/2005-November/082347.html So if my script is called "myRscript.r", I can do the following: this.file = parent.frame(2)$ofile this.dir = gsub("/myRscript.r", "", this.fi

Re: [R] predict() an rpart() model: how to ignore missing levels in a factor

I don't think that, considering the mechanism behind recursive partitioning, that there is any way for you to ignore the crop factor if it is not in the original test set. What decision should be made if, for instance, the next split in a decision tree were on crops and output was 5 for apples,

Re: [R] how exactly does 'identify' work?

Hi, I think the problem is 1 - when a linear model is fitted, ploting the qqnorm( test.lm$ res ) we dont 'know' what values are actually being used on the y-axis; and how do we refer to the ‘Index’ on the x-axis?? therefore, i dont know how to refer to the x and y coordinates in the identi

[R] generalized linear regression - function glm - dismissed predictors - more information about simulated data

Hello, a) Thanks a lot for the answer relative to the dismissed coefficients. However, the result below has been obtained in the R console and there was no warning (please, see the full display below (***) if you wish). "yi ~ X$V1 + X$V2 + X$V3 + X$V4 + X$V5 + X$V6 + X$V7 + X$V8 + X$V9 + X$V

Re: [R] odfWeave/XML Windows issue

Dave_F battelle.org> writes: > > > I am getting the following error when using odfWeave > > Error in xmlEventParse(infile, handlers = handlers, trim = FALSE, state = > state) : > File content_1.xml does not exist > > > Thanks, > Dave Dear Dave! I was getting the same message while tryi

Re: [R] how exactly does 'identify' work?

On 18/11/2010 1:50 PM, casperyc wrote: Hi, I think the problem is 1 - when a linear model is fitted, ploting the qqnorm( test.lm$ res ) we dont 'know' what values are actually being used on the y-axis; and how do we refer to the ‘Index’ on the x-axis?? therefore, i dont know how to refer

Re: [R] how exactly does 'identify' work?

yes, i tried to modify the "2L" part in plot.lm ### if (show[2L]) { ylim <- range(rs, na.rm = TRUE) ylim[2L] <- ylim[2L] + diff(ylim) * 0.075 qq <- qqnorm(rs, main = main, ylab = ylab23, ylim = ylim, ...) if (qqline)

[R] R Courses***Fundamentals & Advanced Programming ***December 2010 - in San Francisco, New York City and Washington DC by XLSolutions Corporation

Check out our Fundamentals and Advanced Programming R courses for December 2010 in Washington DC, San Francisco and New York City http://www.xlsolutions-corp.com/Rcourses *** R Fundamentals and Programming Techniques New York, Dec 13-14, 2010 Washington DC, Dec 16-17, 2010 San

  1   2   >