similar to: how to specify starting values in varIdent() of lme()

Hi All I was wondering how to get the same results with gls and lme. In my lme, the design matrix for the random effects is (should be) a identity matrix and therefore G should add up with R to produce the R matrix that gls would report (V=ZGZ'+R). Added complexity is that I have 3 levels, so I have R, G and say H (V=WHW'+ZGZ'+R). The lme is giving me the correct results, I am

Do I miss here something? dtaa = read.table("http://www.ats.ucla.edu/stat/mplus/examples/ma_snijders/mlbook1.dat", sep=",") head(dtaa) # shows the data as it should be save(dtaa,"dtaa",file="c:/dtaa") d = load("c:/dtaa") head(d) # all data is lost, it only shows [1] "dtaa" "dtaa" Thanks for your hint on this.

I'm trying to use func_odbc to write to a MS SQL db. Here's my func_odbc conf: [OPTIN] dsn=MSSQL-Optin write=INSERT into OptIn (orgID) values (${VAL1}) Dial Plan exten => +18665551212,n,Set(ODBC_OPTIN()=dave) When I do an odbc show, it shows that I am connected to the db. If I use isql, I can write to the db, however, when I use func_odbc, a record will not write. I'm using

Just idly watching * in console mode and saw that someone from 50.49.54.102 tried to register with my *. whois gives:- OrgName: Internet Assigned Numbers Authority OrgID: IANA Address: 4676 Admiralty Way, Suite 330 City: Marina del Rey StateProv: CA PostalCode: 90292-6695 Country: US NetRange: 50.0.0.0 - 50.255.255.255 CIDR: 50.0.0.0/8 NetName: RESERVED-50

Hi Reading the help for ?unz I was wondering if I can read data into R from within an zipfile that is on some website, like maybe: dtaa = read.table(unz("http://www.ats.ucla.edu/stat/examples/alsm/alsm.zip","Ch01pr19.dat")) Thanks for letting me know if you came acros such a thing before. Toby

This morning someone tried to make sip call through my Asterisk. My server just drop these calls and record them in CDR with IP address: 2012-11-28 06:30:51 SIP/216... 1000 "1000" <1000> Hangup 999011972592249388 ANSWERED 00:01 Hacker: 168.63.67.239 2. 2012-11-28 06:30:49 SIP/216... 1000 "1000" <1000> Hangup 88011972592249388 ANSWERED 00:01 Hacker:

(See DAAG book, p. 173, ex. 3) I'm a new user of R, and I'm following the DAAG text. I want to create a subset of the races2000 data frame, but get errors because of a mismatch of values in some columns: > library(DAAG) > attach(races2000) > hills2000 <- races2000[races2000$type == 'hill'] Error in as.matrix.data.frame(x) : dim<- : dims [product 770] do not

I've been getting a lot of timeouts on non-critical invite transactions. I turned on sip debug. They were the result of SIP invites like this: Retransmitting #10 (NAT) to 185.107.94.10:13057: SIP/2.0 401 Unauthorized Via: SIP/2.0/UDP 215.45.145.211:5060;branch=z9hG4bK-524287-1---zg4cfkl50hpwpv4p;received=185.107.94.10;rport=13057 From:

I have a new device that takes measurements anywhere from every second, to every 15 minutes (depending on changes). The matrix has a date, time and Y column (Y is the measurement). For three days it is 25,000 rows. How do I average the measurements by every 30 minutes so my matrix is 48 rows per day? I have been working on this and cannot figure out a simple method. Any ideas? Thank you. ----- In

I just encountered a weird issue on my servers - Tinc was using a constant 10-50% CPU on several servers, and these servers were also receiving a constant ~3 Mb/s of data over the Tinc interface, which is usually otherwise pretty quiet. Example: https://d.sb/2019/04/firefox_11-15.54.22.png Grafana dashboard: https://dash.d.sb/dashboard/snapshot/6nWZqagpgxzxUrybDZkNbF6JSflLlKmO?orgId=1 This seems

Dear R-Users, I am currently working with LME to analyse repeated measures data. I encounter a problem when including both a random effect and a correlation structure with different grouping levels into the LME model. The error message is: Error in lme.formula(diameter ~ flowers*timef + competition*timef + population*timef, : Incompatible formulas for groups in "random" and

Hi I have a huge text file and .dat file from which I want to read data. I do not need all the columns in the files. I want to extract only some columns from the .txt file or the .dat file, because reading the entire file is becoming very difficult due to memory issues. Is it possible to extract a few columns from .txt or .dat file while reading the data in R? Thanks Upasna -- "The past

Can anyone help explain to me why the two codes below have different result? I thought I can use log(time)~. to replace log(time)~dist+climb+timef.I am using CVlm from DAAG package. I think nihills is preloaded with the package. Thanks in advance. > CVlm(df=nihills, form.lm=formula(log(time)~.),plotit="Observed",m=2)Analysis of Variance Table Response: log(time) Df Sum Sq

Have you ever thought about what makes a good trader? I have! Apart from skills and persistence, it?s information. Knowledge is everything, and after you read this to the end, you?ll be more likely to achieve! Get C T X E . P K First Thing Today This Is Going To Explode! C T X E . P K - C a n t e x E n e r g y C o r p Global score: Excellent| Mean Recommendation: Srtong Byu| Investor

I issue the following: > d <- density(rnorm(1000)) > d and get: Call: density.default(x = rnorm(1000)) Data: rnorm(1000) (1000 obs.); Bandwidth 'bw' = 0.2235 x y Min. :-3.5157 Min. :2.416e-05 1st Qu.:-1.6892 1st Qu.:1.129e-02 Median : 0.1373 Median :7.267e-02 Mean : 0.1373 Mean :1.367e-01 3rd Qu.: 1.9639

Hi All, I'm using the following code to produce some stacked bar graphs. *setwd("C:\\Users\\Tinus\\Documents\\NMMU\\R\\Seamounts")* *SChla <- read.csv("SM_Chla_data.csv")* * * *#Extract mean values from data file* * * *Coral <- SChla[185:223,] #Reduce SChla to Coral only* *coral <- with(Coral , aggregate(cbind(Pico, Nano, Micro), list(Depth),FUN=mean))*

Howdy, I'm working on a macro that authenticates the calling extension against a list of allowed extensions but it looks like the Expression I'm attempting to send of pipe separated extensions is showing up as additional arguments to my macro. I expected to have 4 arguments to the below macro, Instead it looks like I'm actually getting 6. I'm open to suggestions to other ways

We have been noticing flurries of sshd reject messages in which some system out there in the hinterlands hits us with a flood of ssh login attempts. An example: Apr 6 05:41:51 dc sshd[88763]: Did not receive identification string from 67.19.58.170 Apr 6 05:49:42 dc sshd[12389]: input_userauth_request: illegal user anonymous Apr 6 05:49:42 dc sshd[12389]: Failed password for illegal user

You've omitted a comma. races2000 is a data frame, which for purposes of extracting rows behaves like a 2-dimenional object. The following works fine: hills2000 <- races2000[races2000$type == 'hill', ] Additionally, you might like to ponder > type <- races2000[names(races2000)=="type"] > type[1:4] Error in "[.data.frame"(type, 1:4) :

On 12/30/2017 08:18 PM, Dovid Bender wrote: > Script kiddies trying to find vulnerable systems that they can make > calls on. Lock down the box with iptables and use fail2ban to block > them. The via is probably bogus unless a box at the DoD was comprimised. > > > > On Sat, Dec 30, 2017 at 6:49 PM, sean darcy <seandarcy2 at gmail.com > <mailto:seandarcy2 at