Re: [R] [R-sig-Geo] raster::levels() not working in packaged function.

2020-10-31 Thread Marcelino de la Cruz Rot
r-sig-...@r-project.org https://stat.ethz.ch/mailman/listinfo/r-sig-geo . -- Marcelino de la Cruz Rot Depto. de Biología y Geología Física y Química Inorgánica Universidad Rey Juan Carlos Móstoles España __ R-help@r-project.org mailing list -- To UNSUBSC

Re: [R] Re : Privacy rights of an old user of this list

2013-09-17 Thread Adan Leobardo Martinez Cruz
Dear all, I will express my opinion without knowing the details of the posts John would like to be removed. In the current state, people posting on this and other servers have no clear way to go when trying to remove their posts. It is a likely event that the number of people attempting the rem

[R] SEM and count data

2013-02-25 Thread Rick Cruz
Hi all- New to R, and I've managed to learn a bit about Lavaan package for SEM, and pscl package for analyzing count data. I would like to find a way to combine these analyses as I am hoping to run a model with a multi-trait, multi-reporter factor (2 reporters by 4 traits apiece, with one o

Re: [R] Expected number of events, Andersen-Gill model fit via coxph in package survival

2012-10-20 Thread Omar De la Cruz C.
standard errors for zz provided by predict.coxph, turns out to be too small. Thank you for any suggestions. Omar. On Thu, Oct 11, 2012 at 1:26 PM, Omar De la Cruz C. wrote: > Thank you, Dr. Therneau, that was very helpful. > > Best regards, > > Omar. > > > On Mon, Oct 8

Re: [R] Expected number of events, Andersen-Gill model fit via coxph in package survival

2012-10-11 Thread Omar De la Cruz C.
Thank you, Dr. Therneau, that was very helpful. Best regards, Omar. On Mon, Oct 8, 2012 at 9:58 AM, Terry Therneau wrote: > >> I am interested in producing the expected number of events, in a >> recurring events setting. I am using the Andersen-Gill model, as fit >> by the function "coxph" in

[R] Expected number of events, Andersen-Gill model fit via coxph in package survival

2012-10-05 Thread Omar De la Cruz C.
Hello, I am interested in producing the expected number of events, in a recurring events setting. I am using the Andersen-Gill model, as fit by the function "coxph" in the package "survival." I need to produce expected numbers of events for a cohort, cumulatively, at several fixed times. My ultim

[R] Problems with plm's pgmm: changing results

2012-01-16 Thread Rolan Paul Veron Cruz
Hello, I'm a newbie with R, and this is my first time to post here.  Apologies in advanced if I should be posting my question on another list. I'm working on a problem on Arellano and Bond's GMM estimator, and I've chosen to use plm/pgmm.  My issue is, I'm getting different results from pgmm()

Re: [R] help with hclust

2011-09-15 Thread Marcelino de la Cruz
3 4 3 Pennsylvania Utah 3 3 HTH, MArcelino > >I would appreciate your help! > >Best regards, >Madeleine > Marcelino de la Cruz Rot D

Re: [R] "useR! 2011" T-shirt competition

2011-05-27 Thread Jose-Marcio Martins da Cruz
David Firth wrote: Dear R users: It is planned that participants in the R User Conference, useR! 2011 (August 16-18, http://R-project.org/useR-2011) will each receive a Better : http://www.R-project.org/useR-2011 __ R-help@r-project.org mailin

Re: [R] How to extract information from the following dataset?

2011-05-12 Thread Jose-Marcio Martins da Cruz
Xin Zhang wrote: Hi all, I have never worked with this kind of data before, so Please help me out with it. I have the following data set, in a csv file, looks like the following: Jan 27, 2010 16:01:24,000 125 - - - Jan 27, 2010 16:06:24,000 125 - - - Jan 27, 2010 16:11:24,000 176 - - - Jan 2

Re: [R] use of ROCR package (ROC curve / AUC value) in a specific case versus integral calculation

2011-03-13 Thread Jose-Marcio Martins da Cruz
them, you can sort them by their scores and enumerate the number of couples out of order you can get. Hope this help. José-Marcio -- --- Jose Marcio MARTINS DA CRUZ http://j-chkmail.ensmp.fr Ecole des Mines de Paris

Re: [R] calculate phase/amplitude of fourier transform function in R

2011-02-13 Thread Jose-Marcio Martins da Cruz
sammyny wrote: I did a fourier transform on a function in time domain to get the following functions in frequency domain (in latex): $Y_1[\omega] = \frac{1}{1-\phi_1 e^{-jw}}$ $Y_2[\omega] = \frac{1}{1-(\phi_1 + \phi_2)e^{-jw} +\phi_1\phi_2e^{-2jw}}$ How do I find the spectrum of this functio

Re: [R] Time unit in ts() and arima() functions

2011-02-12 Thread Jose-Marcio Martins da Cruz
Thanks for the hint. Helped a lot Gabor Grothendieck wrote: On Sat, Feb 12, 2011 at 7:48 AM, Jose-Marcio Martins da Cruz Functions which work with ts typically assume that a full cycle is represented by 1 unit so if a full cycle is a week then a week must be one unit and a day must be 1

[R] Time unit in ts() and arima() functions

2011-02-12 Thread Jose-Marcio Martins da Cruz
This question is surely trivial, sorry. I'm afraid I'm misunterpreting the information I got with the documentation, and I'm a little bit confused. I'm just an engineer with some little skills in statistics. Well, I have a time series - 600 days long - with some weekly periodicity inside. So

Re: [R] predict and arima

2011-02-09 Thread Jose-Marcio Martins da Cruz
D=1 and select d using ndiffs and the KPSS test # for stationarity. # ... if no model can be found - return NULL -- --- Jose Marcio MARTINS DA CRUZ http://j-chkmail.ensmp.fr Ecole des Mines de Paris 60, bd Saint Mi

Re: [R] predict and arima

2011-02-08 Thread Jose-Marcio Martins da Cruz
I think ther's a bug here : bdp wrote: Some code I have been playing with to do this follows ... get.best.arima<- function(x.ts, minord=c(0,0,0,0,0,0), maxord=c(2,1,1,2,1,1)) { # function based on 'Introductory Time Series with R' best.aic<- 1e8 # a big number n<- leng

[R] Cross validation for Ordinary Kriging

2011-01-06 Thread pearl may dela cruz
ear ALL, The last part of my thesis analysis is the cross validation. Right now I am having difficulty using the cross validation of gstat. Below are my commands with the tsport_ace as the variable: nfold <- 3 part <- sample(1:nfold, 69, replace = TRUE) sel <- (part != 1) m.model <- x2[sel, ] m

[R] Prediction error for Ordinary Kriging

2011-01-05 Thread pearl may dela cruz
ce) ~ 1, x2, cutoff=100, width=9) te<- fit.variogram(v,vgm(0.0437, "Exp", 26, 0)) y <- krige(tsport_ace~1, x2, x2_SG, model = ve.fit) spplot(y, 1, col.regions = bpy.colors(100), sp.layout = list("sp.lines",as(rsa2, "SpatialLines"),no.clip = TRUE)) I'm

[R] Convert tileplot's trellis to shapefile/ Add values to voronoi polygon

2010-12-16 Thread pearl may dela cruz
In any way, can you help me on either converting tileplot's trellis to shapefile/image or (solving the problem of voronoi polygons and adding another column to represent the polygons in tripack). Thanks a lot in advance, Pearl dela Cruz ___

[R] Unable to load an object

2010-06-11 Thread Luis Ridao Cruz
R-help, I seem not to get an object saved neither with ".RData" extension nor output via dput. Whenever I try to import the above object in another worspace I just get nothing. > geoFeatures <- load("geoFeatures.RData") > geoFeatures [1] "geoFeatures" The "geoFeatures.RData" workspace contain

[R] Sort a data frame according to romans

2009-12-08 Thread Luis Ridao Cruz
R-help, I have a data frame: > mydata strata nh NhWhfh 1 I 10 26 0.048 0.385 2 II 32 84 0.154 0.381 3 III 16 42 0.077 0.381 4 IV 4 11 0.020 0.364 5 V 10 26 0.048 0.385 7 VII 64 168 0.309 0.381 8VIII 49 129 0.237 0.380 9 IX 22 58 0.107 0.37

[R] could not find function lapply<-

2009-12-08 Thread Luis Ridao Cruz
R-help, I have a list whose elements are data frames. I want to change the colnames attribute in each element of this list but an error message comes up: > lapply(LD_strataNew,function(x) dimnames(x)[[2]][-1]) <- as.roman(1:9)[-6] Error in lapply(LD_strataNew, function(x) dimnames(x)[[2]][-1]) <

[R] Optim(...) estimate of stDev far too low

2009-09-30 Thread Luis Ridao Cruz
R-help, I'm just trying to find the ML (maximum likelihood) estimates of the mean and standard deviation of a set of observations: >xx=c(2.5,3.5,4,6,6.5,7.5) fn<-function(params,x=xx) { media<-params[1] st <-params[2] pdf=-sum(dnorm(log(xx),log(media),st,TRUE)) return(pdf) } optim(c(mu,stdev

[R] Maximum likelihood estimation of parameters make no biological sense

2009-09-24 Thread Luis Ridao Cruz
R-help, I'm trying to estimate some parameters using the Maximum Likehood method. The model describes fish growth using a sigmoidal-type of curve: fn_w <- function(params) { Winf <- params[1] k <- params[2] t0 <- params[3] b <- params[4]

[R] Error: C stack usage is too close to the limit

2009-09-14 Thread Luis Ridao Cruz
R-help, I 'm trying to optimize a model to data using log-likelihoods but I encounter the following error message: > l= c(49.4, 57.7,64.8,70.9,78.7,86.6,88.3,91.6,99,115) > t=3:12 > fn <- function(params, l=l, t=t) { Linf <- params[1] k <- params[2]

[R] cannot upgrade to R 2.9.1 in Ubuntu

2009-08-07 Thread Luis Ridao Cruz
R-help, At the moment the R version installed on my machine is 2.8.1. (Ubuntu 9.04) I wish to upgrade to R 2.9.1. I did: $ sudo apt-get upgrade ..but R is not upgraded although the "sources.list" file is updated with: deb http://cran.ii.uib.no/bin/linux/ubuntu jaunty/ When I run from the term

[R] RODBC package to connect to Oracle database Linux

2009-08-04 Thread Luis Ridao Cruz
t on this I would need also help as there is no example on Oracle drivers (just MySql, PostgreSQLt) Thanks in advance ---- Luis Ridao Cruz Faroe Marine Research Institute Nóatún 1, P.O. Box 3051 FO-110 Tórshavn Faroe Islands Tel : (+298

[R] Oracle ODBC driver for Linux

2009-08-04 Thread Luis Ridao Cruz
R-help, I get the following error message when trying to connect to an Oracle database through R (2.8.1) under Linux (Ubuntu 9.04). > channel<-odbcConnect("magnus",uid="luisr",pwd="juanayzakarias") Warning messages: 1: In odbcDriverConnect(st, ...) : [RODBC] ERROR: state IM002, code 0, message

[R] Problem installing RODBC package under UBUNTU

2009-06-24 Thread Luis Ridao Cruz
R-help, I'm now on the process of installing RODBC on UBUNTU but an error occured. I copy-paste the output: > install.packages("RODBC") Warning in install.packages("RODBC") : argument 'lib' is missing: using '/home/luisridaocruz/R/x86_64-pc-linux-gnu-library/2.8' trying URL 'http://cran.ii.uib

[R] rgdal package -- Error in installing package

2009-06-24 Thread Luis Ridao Cruz
R-help, I'm trying to install the rgdal package under Ubuntu but I get the following warning: > install.packages("rgdal") Warning in install.packages("rgdal", configure.args = "--with-gdal-modules=/usr/local/lib") : argument 'lib' is missing: using '/home/luisridaocruz/R/x86_64-pc-linux-gnu-

[R] R under Ubuntu

2009-05-26 Thread Luis Ridao Cruz
R-help, I have installed R under Ubuntu and I'm very new to a Linux distribution. To open an empty R session I just type R on the Terminal aplication. But how can I open a saved workspace? At present I just start R and then load ("my_workspace") but it must be possible to do it all at oncer

Re: [R] [R-sig-Geo] Comparing spatial distributions - permutation test implementation

2009-05-21 Thread Marcelino de la Cruz
wo populations. Syrjala SE. Ecology. 1996;77(1):75­80. http://dl.getdropbox.com/u/1047321/Syrjala1996.pdf [2] https://stat.ethz.ch/pipermail/r-sig-geo/2008-February/ thread.html#3137 ___ R-sig-Geo mailing list r-sig-...@stat.math.ethz.ch https://stat.ethz.ch/mailman/

Re: [R] RES: set a legend with an "image" graph

2009-03-17 Thread Luis Ridao Cruz
R-help, I'm sorry that I said " I don't have idea ". What I mean t is that I just don't know how to set a legend when "image" is called. Thanks in advance __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do

[R] set a legend with an "image" graph

2009-03-17 Thread Luis Ridao Cruz
R-help, I have an image plot and I wish to set a legend to it but I don't have an idea how to do it. Below it's the code line. image(yyy, col = rev(heat.colors(10)), axes = T) Thanks in advance. __ R-help@r-project.org mailing list https://stat.ethz.c

Re: [R] cluster analysis: mean values for each variable and cluster

2009-02-20 Thread Marcelino de la Cruz
") dist cluster=hclust(dist,method="ward") cluster plot(cluster) rect.hclust(cluster, k=4, border="red") I extracted 4 clusters from the data. My question is: is it possible to produce a summary of every mean values for each variable of each of the 4 clusters? Thanks a lot

Re: [R] multidimensional scaling with long form data

2009-02-18 Thread Marcelino de la Cruz
existing functions. Thanks! Marcelino de la Cruz Rot Departamento de Biología Vegetal E.U.T.I. Agrícola Universidad Politécnica de Madrid 28040-Madrid Tel.: 91 336 54 35 Fax: 91 336 56 56 marcelino.delac...@upm.es _ [[alte

[R] distance betwenn axis and axis annotation

2009-02-18 Thread Luis Ridao Cruz
R-help, How can I specify the distance between the axis and axis annotation? Thanks in advance __ 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.htm

[R] Filter a big matrix

2009-02-11 Thread cruz
ight be syntax errors, I just want to clearly describe my problem Since paste("X$a>=5") gives type of character; whereas the type of X$a>=5 should be logical. How can I do this? All thoughts are greatly appreciated. Many Thanks, cruz __ R

Re: [R] error in my previous message

2009-02-10 Thread cruz
try this: temp2 <- temp2[order(temp2[,1]),] # PROBLEM IS HERE dim(temp2) <- c(5,1) print(temp2) On Wed, Feb 11, 2009 at 2:23 PM, wrote: > i'm sorry. i had an error in my previous code because i left out a letter in > the rownames. > while fixing that, i also found a solution. so i'm sorry fo

Re: [R] installing R on Ubuntu

2009-02-08 Thread cruz
x27;t even start a R session in Emacs when I have read the ESS manual several times, but it is getting clearer to me now as I have been writing R codes by Emacs for a couple of weeks... It is exactly the same as we code in Windows. Thanks, cruz PS: The R gurus may find this thread annoying cause it i

Re: [R] R on Mobile Devices (Android)

2009-02-06 Thread cruz
On Fri, Feb 6, 2009 at 8:12 PM, Harsh wrote: > Hell R-list, > > At the cost of sounding far-fetched and almost incredulous, I would > like to know if any R user is remotely considering the use of R on > Mobile devices, and Android in particular. > Check out FreeRunner from Openmoko (http://wiki.o

Re: [R] annotating a filled contours plot with a grid of points

2009-02-06 Thread Marcelino de la Cruz
A > Dpto. Sistemas y Recursos Forestales > Ctra. La Coruña, Km. 7.5 > E-28040 Madrid (Spain ) > Tel.: +34 91 347 1461 > e-mail: dmartin inia.es /// dmartinbenito gmail.com > > <http://dmartinbenito.googlepages.com>http://dmartinbenito.googlepages.com >

Re: [R] Text Mining

2009-02-05 Thread cruz
There is an interesting article "An Introduction to Text Mining in R" by Ingo Feinerer on R News Volume 8/2, October 2008 (http://www.r-project.org/doc/Rnews/Rnews_2008-2.pdf) Check it out On Fri, Feb 6, 2009 at 9:16 AM, spiketide wrote: > > hi everyone... > > i am a newbie to text mining.

[R] non-numeric argument to binary operator

2009-02-04 Thread cruz
(paste("x", 2, sep = "")) : non-numeric argument to binary operator > All comments are appreciated. Thanks, cruz __ 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, reproducible code.

[R] Alternate looping

2009-01-06 Thread Luis Ridao Cruz
R-help, I'm using the "for" control flow to graph plots continuously: par(mfrow=c(3,5), mar=c(0.6,.2,1.2,.2),yaxt="n",xaxt="n") for(j in 1:11) { for(i in 1:15) { species <- spAldur[spAldur$ar == 1993+i & spAldur$aldur == j,] plot(spec...@coords[,1], spec...@coords[,2], xlim=c(-10,-3.5),

Re: [R] Newbie question

2009-01-03 Thread cruz
abntest[abntest$abntesr > 90,] On Sun, Jan 4, 2009 at 11:55 AM, wrote: > Hi: > > I'm loading in students test scores with: > >> abntest <- read.table("scores.txt") > > if I type: > >>abntest > > I get ALL the values. I want to be able to filter it by various things such > as: > > > if( abntes

[R] example of using loess()

2009-01-03 Thread cruz
involved. Thanks a lot, cruz __ 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, reproducible code.

[R] Fitting weighted polynomial regression model

2008-12-29 Thread cruz
7; and 'y' lengths differ In addition: Warning messages: 1: pseudoinverse used at 0.52416 0.029047 2: neighborhood radius 3.8766 3: reciprocal condition number 4.3293e-017 4: There are other near singularities as well. 19.538 > Any idea woul

[R] R simulations

2008-12-04 Thread João Cruz
I need help with these R simulations ___ The generation of pseudo-random numbers with uniform distribution [0.1] can be done through a method called Congruencial Committee, which consists the use of a recursive relationship of the type: xi +1 = (a.xi + c) mod

[R] CRAN mirror for Taipeh vs Taipei

2008-11-22 Thread cruz
Hi, I am residing in Taipei, Taiwan. My windows has English R (2.7.2) installed. I noticed that the CRAN mirror for Taipei is Taiwan (Taipeh), Is Taipeh a Deutsch or English name? Thanks, cruz __ R-help@r-project.org mailing list https://stat.ethz.ch

Re: [R] help for code in jump diffusion

2008-11-22 Thread cruz
i tried to run your code, this is how/why you got NaN: > mu<-0.2 > sig<-0.2 > S0<-100 > j<-0.2 > dt<-1/252 > int<-0.1 > i<-0 > is.nan function (x) .Primitive("is.nan") > k<-rnorm(1,0,1) > k [1] 0.3214954 > theta<-ifelse((k<(int*dt)),1,0) > theta [1] 0 > m<-rnorm(1) > m [1] -0.4525731 > ga

Re: [R] problem with urca package

2008-11-11 Thread cruz
is this what you expect? > ur.df(PIBTUN$V1,type='none',lags=1) ### # Augmented Dickey-Fuller Test Unit Root / Cointegration Test # ### The value of the test statistic is: 6.252

[R] estimated variance for parametric fit

2008-11-09 Thread cruz
Hi, What formula is appropriate for calculating the estimate of variance from the parametric fit? I have a linear regression model i.e. lm(a~b) Thanks, cruz __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do

Re: [R] Re move specific rows from a matrix

2008-11-07 Thread cruz
use the trick from this thread: https://stat.ethz.ch/pipermail/r-help/2008-November/178985.html your problem shall be fixed On Fri, Nov 7, 2008 at 8:21 PM, mentor_ <[EMAIL PROTECTED]> wrote: > > Hi, > > imagine the following two matrix: > Matrix A: > a b c > 1 2 3 > 4 5 6 > 7 8 9 > > Matrix B:

Re: [R] grouping data together

2008-11-07 Thread cruz
hi friend, this is from your previous posts on Kruskal-Wallis test:) i came up with this one: A5 <- read.table('kew.dat' ,header=TRUE) plot(factor(A5$Month, levels=month.abb), A5$Rain) is that what you want? On Sat, Nov 8, 2008 at 7:03 AM, Swanton0822 <[EMAIL PROTECTED]> wrote: > > Hi. > i hav

Re: [R] missing value where TRUE/FALSE needed

2008-11-07 Thread cruz
is this what you want? > vector1 [1] 65 1 34 100 42 20 79 43 89 10 > vector2 [1] 34 65 47 91 48 32 23 74 92 86 > > for (i in 1:10) { + for (j in 1:10) { + if (vector1[i] == vector2[j]) + show(c(i,j)) + } + } [1] 1 2 [1] 3 1 On Sat, Nov 8, 2008

[R] chi square table

2008-11-06 Thread cruz
them give me 20.090 Thanks, cruz __ 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, reproducible code.

Re: [R] How to avoid "$ operator is invalid for atomic vectors"

2008-11-06 Thread cruz
> > Does that answer your question? > Thanks:) I received one from Erin: x <- NULL __ 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 provi

Re: [R] How to avoid "$ operator is invalid for atomic vectors"

2008-11-06 Thread cruz
n of this object should not be fixed because it will grow... so in MATLAB, i.e. we define x = [ ], what about in R? Thanks, cruz __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://w

[R] How to avoid "$ operator is invalid for atomic vectors"

2008-11-06 Thread cruz
Hi, I am writing this in a wrong way, can someone please correct me? > A <- matrix() > length(A) <- 6 > dim(A) <- c(3,2) > colnames(A) <- c("X","Y") > A X Y [1,] NA NA [2,] NA NA [3,] NA NA > A$X Error in A$X : $ oper

Re: [R] Data manipulation question

2008-11-06 Thread cruz
On Thu, Nov 6, 2008 at 4:23 PM, Peter Jepsen <[EMAIL PROTECTED]> wrote: > > Here is an example: > > id <- c(rep("a",4),rep("b",2), rep("c",5), rep("d",1)) > start <- c(c(0,6,17,20),c(0,1),c(0,5,10,11,50),c(0)) > stop <- c(c(6,12,20,30),c(1,10),c(3,10,11,30,55),c(6)) > data <- as.data.frame(cbind(id

Re: [R] Reshape a matrix

2008-11-05 Thread cruz
On Thu, Nov 6, 2008 at 8:53 AM, dinesh kumar <[EMAIL PROTECTED]> wrote: > Dear R users, > > I have a matrix like > > A X1 > B Y2 > C Z3 > > I want to reshape this matrix into this format > >X Y Z > A 1 > B 2 > C 3 > is this what you

Re: [R] FW: need some help

2008-11-05 Thread cruz
> > How can I assign df$counts values to dd$counts for the corresponding > dd$comp==df$comp > ?factor the $comp then assign accordingly __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide

Re: [R] use abline() for regression model in the plot

2008-11-04 Thread cruz
On Wed, Nov 5, 2008 at 6:47 AM, Swanton0822 <[EMAIL PROTECTED]> wrote: > > after that i have been ask to do the following: > 'The regression equation for EXPEND depending on ECAB is > E(EXPEND)=119.0+1.73xECAB when there is outliners and > E(EXPEND)=97.8+1.94xECAB when the outliners are omitted fro

[R] qr() and Gram-Schmidt

2008-11-03 Thread cruz
r) > Q [,1][,2] [,3] [1,] -0.7071068 -0.70710680 [2,] 0.000 0.000 1 [3,] -0.7071068 0.70710680 Thanks, cruz __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE

[R] Reconstruct the Q matrix from the QR object

2008-10-29 Thread cruz
Hi, Is it possible to construct a Q from qr() that some of the rows could be specified (to be fixed values)? Thanks, cruz __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R

[R] How to label a "vector"

2008-10-01 Thread cruz
Hi, I have an object, Z: > typeof(Z) [1] "list" It looks like: [,1] [,2] [,3] [1,] 111 [2,] 111 How can it label the Rows and Columns with Text? so that it looks like: Y1 Y2 Y3 X1 111 X2 111

[R] savePlot() does not save plot with format set

2008-08-26 Thread Luis Ridao Cruz
R-help, Whenever I try to save a plot with "savePlot" the file is not stored in my hard disk with the selected format. Several formats are set and none of them works. I just get the file name with missing extension and it can't be open with programs like Paint and Microsoft Photo Editor Th only on

[R] "remove.packages" gives error message

2008-04-02 Thread Luis Ridao Cruz
R-help, I'm trying to remove some packages with the command "remove.packages" but R 2.4.1 will come out with this: > remove.packages(pkgs="FLAssess", lib = "C:\\Program Files\\R\\R-2.4.1\\library\\FLAssess") Error in .find.package(pkgs, lib) : there is no package called 'FLAssess' Thanks in adv

[R] equivalent of excel "sumif"

2008-03-04 Thread Luis Ridao Cruz
R-help, I have two data frames df1 and df2: >head(df1) lgdcm 0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 2 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 3 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 4 0 0 0 0 0

[R] data manipulation

2008-02-29 Thread Luis Ridao Cruz
R-help, I have a data frame in which I compute maximum by rows and I wish to find out the column name at which this maximun occurs. Does anyone know how to do it? Thanks in advance __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/list

[R] shaded area graph and extra plot

2008-02-14 Thread Luis Ridao Cruz
R-help, I'm using the code below to plot a shaded area graph. At the same time I want to plot a second series on the y-axis (from par(new=T) on) but as the two series have different x-axis range (first 1994:2007 and second 1996:2007) the corresponding x's do not match. How can this be sorted ou

[R] shaded area graph and extra plot

2008-02-12 Thread Luis Ridao Cruz
R-help, I'm using the code below to plot a shaded area graph. At the same time I want to plot a second series on the y-axis (from par(new=T) on) but as the two series have different x-axis range (first 1994:2007 and second 1996:2007) the corresponding x's do not match. How can this be sorted ou

[R] sort a data frame according to roman characters

2008-02-12 Thread Luis Ridao Cruz
R-help, I have a data frame with one column containing roman numbers The data are not sorted as : III III IV VVI VII VIII IX XXI XII XIII XIV XV Using data[order(data$Roman),] does not do the job. How can this be done? Thanks in advance. __

[R] function code

2008-01-25 Thread Luis Ridao Cruz
R-help, Sorry for this question (I guess it has been addressed before but I could not find it in the archives) but how can I see a function code when the following comes up: > svymean function (x, design, na.rm = FALSE, ...) { .svycheck(design) UseMethod("svymean", design) } Thanks i

[R] reduce the code used by "for"

2007-12-04 Thread Luis Ridao Cruz
R-help, I have a 3-way array: > dim(bugvinP) [1] 13 14 3 The array looks something like this (object trimmed for readability) , , slag = 1 ar 1994 1995 1996 1997 1998 1 NA 0.000 0.000 0.000 0.000 2 0.036 0.059 0.027 0.000 0.000 3 0.276 0.475 0.491 0.510 0.559 10 1.0

[R] Fail to open R 2.6.1 - additional information

2007-11-28 Thread Luis Ridao Cruz
This is additional information to my previous email (pasted below) ## I opened an R session ( R 2.6.0 on Windows) and got the following pop-up window: ## "This application has failed to start because tk84.dll was not found. Re-installing the application may fix the problem" ## I have now upgr

[R] Fail to open R 2.6.1

2007-11-28 Thread Luis Ridao Cruz
R-help, I opened an R session ( R 2.6.0 on Windows) and got the following pop-up window: "This application has failed to start because tk84.dll was not found. Re-installing the application may fix the problem" I have now upgraded to R 2.6.1 but the same fatal error comes up. Can someone let me

[R] print " character

2007-11-20 Thread Luis Ridao Cruz
R-help, Sorry if this question has been discussed/posted before but I can't just find it myself. How can I print the comment character (") ? Thanks in advance __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do

[R] Unable to save plot as ".pdf" , "ps" , "eps"

2007-11-15 Thread Luis Ridao Cruz
R-help, Whenever I try to save a plot with extension ".pdf" , "ps" or "eps" I get the following error/warning message: plot(rnorm(10)) savePlot("test",type="pdf") Error in savePlot("test", type = "pdf") : Invalid font type In addition: Warning messages: 1: In savePlot("test", type = "pdf") : f

[R] manipulating 3-way array

2007-10-18 Thread Luis Ridao Cruz
R-Help, I have a 3 dimensional array with dimensions: > dim(test) [1] 241 104 12 I copy-paste an example of the object > test[1:5,1:5,1:2] , , ar = 1996 slagnr puntar12 34 5 I1 60.30303 36.50610 0 0.0 0 I10 58.19512 35.48209 0 75.0 0 I11 0.