similar to: Time series plot

Hi All, I created a two variable lm() model slm<-lm(y[1:3000,8]~y[1:3000,12]+y[1:3000,15]) I made two predictions predict(slm,newdata=y[201:3200,]) predict(slm,newdata=y[601:3600,]) there is no error message for either of these. the results are identical, and identical to slm$fitted as well. if this is not the right way to apply the model coefficients to a new set of inputs, what is

Dear List, I'm trying to use R under Windows on a huge database in MySQL via ODBC (technical reasons for this...). Now I want to read tables with some 160.000.000 entries into R. I would be lucky if anyone out there has some good hints what to consider concerning memory management. I'm not sure about the best methods reading such huge files into R. for the moment I spilt the whole

Dear Group, Having a character vector like this one: [1] "03:38:55" "07:42:38" "08:04:27" "08:17:13" "08:41:14" "08:46:58" [7] "08:47:11" "08:53:51" "08:57:51" "08:58:56" I try to do a histogram over times of a day. All I want to know, if my solution is proper or if there is another

Dear Group I'm still trying to bring many data into R (see older postings). After solving some troubles with the database I do most of the work in MySQL. But still I could be nice to work on some data using R. Therefore I can use a dedicated Server with Gentoo Linux as OS hosting only R. This Server is a nice machine with two CPU and 4GB RAM which should do the job: Dual Intel XEON 3.06 GHz

Dear Group, if anyone has experience with the RMySQL Package maybe this behaviour is know: Reading data from a table into R the fields with datatype bigint(20) unsigned are transformed in some way: e.g. the query "select * from orders where userid = 14929859848712890325" selects the correct case but in R the userid is changed to 14929859848712890368. What happened here? This

Dear List, I have two data.frame of the following form: A: n V1 V2 1 12 0 2 10 8 3 3 8 4 8 4 6 7 3 7 12 0 8 1 0 9 18 0 10 1 0 13 2 0 B: n V1 V2 1 0 2 2 0 3 3 1 9 4 12 8 5 2 9 6 2 9 8 2 0 10 4 1 11 7 1 12 0 1 Now I want to merge those frame to one data.frame with summing up the columns V1 and V2 but not the column n. So the result

In my (limited) experience R is more powerful concerning data manipulation. An example: I have a vector holding a user id. Some user ids can appear more than once. Doing SELECT COUNT(DISTINCT userid) on MySQL will take approx. 15 min. Doing length(unique(userid)) will take (almost) no time... So I think the other way round will serve best: Do everything in R and avoid using SQL on the database...

I found that read.table() does not work well on lines that is longer than 236 bytes. Is this know problem? Is it under fixing? Regards, Jincai Jiang (Office) 212-761-3984 -------------------------------------------------------- NOTICE: If received in error, please destroy and notify sender. Sender does not waive confidentiality or privilege, and use is prohibited. [[alternative HTML

Hi all, I have a question regarding connection of two Asterisk servers to our PBX. Each Asterisk server has one PCI E1 card, and they are in failover mode with Linux HA. On our PBX we have only one E1 card towards Asterisk servers. My question is how to connect these two Asterisks to one E1 card on PBX, and when primary Asterisk server fails not to have to manually pull out E1 cable from primary

Hi all, I'm testing Linksys SPA922 phone and I have strange issue. when call is finished on the phone I see "CallEnded" and normal silence for cca. 5 seconds and then I get fast busy for cca. 20 sec. So, this isn't automatic hangup as on other phones I have tried (Cisco 7940, grandstream, XLite,... ) and I have to manually hangup handset to finish a call. Is this normal behavior

Hi there, I am wondering if it is possible to do an xyplot with two y-axes. I'd like to print two parameters in a time series but they both have different scales. Which parameter in xyplot can I add to achieve this result? Thanks a lot for any help. Antje [[alternative HTML version deleted]]

Hi all, I have a question regarding using T38 for fax sending and here is my scenario: fax -> SIP ATA (T38 enabled) -> Asterisk #1 -> IAX TRUNK -> Asterisk #2 -> SIP ATA (T38 enabled) -> fax My question is, how can I know if I'm really using T38? is T38 information coming to the other side (because of SIP to IAX conversion) or just plain g711a data? I'm using Linksys

Hello all, I am working with a 1000x1000 matrix, and I would like to return a 1000x1000 matrix that tells me which value in the matrix is greater than a theshold value (1 or 0 indicator). i have tried mat2<-as.matrix(as.numeric(mat1>0.25)) but that returns a 1:100000 matrix. I have also tried for loops, but they are grossly inefficient. THanks for all your help in advance. Lanre

Hola amigos, esta es mi primera duda, espero que no sea demasiado fácil. Tengo unos datos de dos variables y quiero mostrar recta de regresión y valor de correlación serie0 <- c(0.651, 0.712, 0.614, 0.645, 0.559, 0.647, 0.642, 0.534, 0.616, 0.621, 0.623) serie1 <- c(0.572, 0.641, 0.565, 0.596, 0.518, 0.604, 0.602, 0.501, 0.58, 0.589, 0.596) data <- cbind(serie0, serie1) colnames(data)

Dear group, First I provide you with an example, I found in the newsgroup. Then I'd like to explain my problem to you by means of the output. enviro <- data.frame(Year = rep(2001:2002, each = 365), Day = rep(1:365, 2), Precip = pmax(0, rnorm(365 * 2)), Temp = 2 + 0.2 * rnorm(365 * 2)) xyplot(Precip + Temp ~ Day | Year,

Hi all, I have a question concerning xyplot. My data is a data.frame looking like that: In the first column I have numbers from 0 to 23 (hours of a day), the second column contains the name of a weekday (Day as factor) and the third column contains the number I am interested in. So as an example, the first five rows look like that: Hour Day Freq 1 0 Mo 23 2 1 Mo

Hi there, I have a question, which I thought is very easy to solve, but somehow I can't find a solution. Probably someone could help me quickly? Here it is: I have two matrices: a [,1] [,2] [,3] [1,] 1 4 9 [2,] 2 6 10 [3,] 3 6 11 [4,] 4 8 12 b [,1] [,2] [1,] 1 4 [2,] 2 5 [3,] 3 6 Now I want to find

Dear all, Does any one knows why should I get the following error message, when trying to do a simple data.frame?? DataF<-data.frame(Subject,BiomR,Spp,Capas,Litter,Herbs,LitterD,MaxCanH,DDifS p,DSSp,Slope, CanDens,NearestSp) Erro em data.frame(Subject, BiomR, Spp, Capas, Litter, Herbs, LitterD, : arguments imply differing number of rows: 202, 0 The data I am using

Hello, I need to perform a mixed-model (with nesting) MANCOVA, using Type III sums of squares. I know how to perform each of these types of tests individually, but I am not sure if performing a mixed-model MANCOVA is possible. Please let me know. Erika <>< <>< <>< <>< <>< <>< <>< Erika Crispo, PhD candidate

Hi: I am using a subset of the below dataset to predict PRED_SUIT for the whole dataset but I am having trouble with 'newdata'. The model was created with 153 records and want to predict for 208 records. wolf2 <- structure(list(gridcell = c(367L, 444L, 533L, 587L, 598L, 609L, 620L, 629L, 641L, 651L, 662L, 674L, 684L, 695L, 738L, 748L, 804L, 805L, 872L, 919L, 929L, 938L, 950L, 958L,